Выполнение пункта задания 7 на Matlab

>> S=fft(s);

>> MODS=(2/N)*abs(S); MODS(1)=(1/N)*abs(S(1));

>> k=0:(N-1);

>> subplot(2,1,2), stem(k,MODS, 'MarkerSize',3), grid, xlabel('k'), title('fft N=1024')

Для определения значений дискретных нормированных частот k, соответствующих гармоникам сигнала (рис. 3.4.2), можно воспользоваться внешней функцией fft_e1.

При использовании функции fft_e1 значение порога ε₁ следует выбирать из условия:

При обработке реальных сигналов значение ε₁, как правило, задается и зависит от конкретной задачи. Для обоснованного выбора ε₁ в данном примере вычисляется числитель и знаменатель – MAX_NOISE и MAXS – и находится их отношение:

Рис. 3.4.2. 1024-точечная последовательность – смесь вещественной периодической последовательности с нормальным белым шумом и модуль ее спектральной плотности.