>> NOISE=fft(noise);
>> MODNOISE=(2/N)*abs(NOISE); MODNOISE(1)=(1/N)*abs(NOISE(1));
>> MAX_NOISE=max(MODNOISE);
>> MAXS=max(MODS)
>> [MAX_NOISE MAXS MAX_NOISE/MAXS]
ans =
0.1659 0.7514 0.2208
Таким образом, диапазон возможных значений ε1 равен:
0.2207< ε1<1.
Выберем ε1=0.5 и определим амплитуды и частоты гармоник, которые, относятся к сигналу:
Выполнение пункта задания 9 на Matlab
>> e1=0.5;
>> [MODm,m]=fft_e1(N,MODS,e1)
MODm =
0.7514 0.4879 0.4879 0.7514
m =
257 385 641 769
Получились значения k1 и k2, такие же, как в примере 3.1.