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

e1=1e-7;

[MODm,m]=fft_e1(N,MOD,e1)

MODm =

0.7000 0.5000 0.5000 0.7000

m =

257 385 641 769

Пример 3. 2

Вычислить спектральную плотность 1024-точечной последовательности, представляющей собой смесь вещественной периодической последовательности x(n) с нормальным белым шумом. Построить графики последовательности и модуля ее спектральной плотности.

Сформируем 1024-точечную последовательность (вектор s) – сумму последовательности x(n) (3.4.1) (вектора х) и нормального белого шума (вектора noise).

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

>> N=1024; n=0:(N-1);

>> x=0.7*sin(0.5*pi.*n)+0.5*sin(0.75*pi.*n);

>> noise=randn(1,length(x));

>> s=x+noise;

>> subplot(2,1,1), stem(n,s,'MarkerSize',3), grid, xlabel('n'), title ('s=x+rand N=1024')

Вычислим модуль спектральной плотности на основе ДПФ с учетом постоянных множителей.