Подгонка кривой ряда Фурье n-го порядка в Python

Question

Подгонка кривой ряда Фурье n-го порядка в Python

1

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

umeshu 20 май 2011, в 12:12

Источник

0

Что именно вы хотите? Вы хотите рассчитать коэффициенты Фурье произвольной функции? (ну, периодическая произвольная функция)
Evpok 20 май 2011, в 09:50
0

Хай евпок. Да, это правильно.
umeshu 20 май 2011, в 10:11
0

Например, у меня есть известные периодические данные, x или произвольная периодическая функция, и я хочу вычислить коэффициенты. Я бы изобразил функцию, cosCoeffs, sinCoeffs = fourier (x, T, N), где T - период, а N - произвольный порядок. Благодарю.
umeshu 20 май 2011, в 10:22

Показать ещё 1 комментарий

Теги:

python

fft

curve-fitting

1 ответ

Ещё вопросы

Что именно вы хотите? Вы хотите рассчитать коэффициенты Фурье произвольной функции? (ну, периодическая произвольная функция)
Например, у меня есть известные периодические данные, x или произвольная периодическая функция, и я хочу вычислить коэффициенты. Я бы изобразил функцию, cosCoeffs, sinCoeffs = fourier (x, T, N), где T - период, а N - произвольный порядок. Благодарю.

Evpok · Answer 1 · 2011-05-20T08-05-00.000Z

Ну, формула

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