демонстрирует расчет
Листинг 14.1 демонстрирует расчет Фурье-спектра по N=128 точкам, причем предполагается, что интервал дискретизации данных yi равен ?. В середине листинга применяется встроенная функция fft, а его оставшаяся часть предназначена для корректного пересчета соответствующих значений частот ?i (они вычисляются в последней строке листинга). Обратите внимание, что результаты расчета представляются в виде модуля Фурье-спектра (Рисунок 14.2), поскольку сам спектр является комплексным. Очень полезно сравнить полученные амплитуды и местоположение пиков спектра (Рисунок 14.3) с определением синусоид в листинге 14.1.
