Подгонка кривой ряда Фурье n-го порядка в Python - PullRequest
1 голос
/ 20 мая 2011

Я искал способ кодирования фрагмента в Python, который вычислял бы для любого n-го порядка подбора кривой ряда Фурье.Чтобы рассчитать определенный порядок подбора кривой ряда Фурье, скажем, 3-й порядок довольно прост, однако сделать это, когда порядок n является переменным, все еще не работаетВозможно, кто-то сделал это, но мой поиск пока не может его найти.Интересно, кто-нибудь мог бы помочь?Благодарю.

1 Ответ

0 голосов
/ 20 мая 2011

Ну, формула

n-th cos_coeff = (2/T)*integral(-T/2,T/2, f(t)*cos(n*t*2*pi/2)dt)
n-th sin coeff = (2/T)*integral(-T/2,T/2, f(t)*sin(n*t*2*pi/2)dt)

Проверьте scypi и scipy.integrate для получения подробной информации об интеграции.

Здесь должно быть

cos_coeff(f, T, N) = (2/T)*quad(lambda t: f(t)*cos(N*t*2*math.pi/2),-T/2,T/2)

(но не проверено)

Я не знаком с Дискретным преобразованием Фурье, но, возможно, вы также можете рассчитать указанный коэффициент по нему. Проверьте http://docs.scipy.org/doc/scipy/reference/tutorial/fftpack.html

...