数字信号处理实验报告3 DSP信号与系统实验报告 信号加窗及谱分析 电子科技大学 2018版

toc tic

X2 = dft(x2); toc

LINES = 2; ROWS = 1; ID = 0; STEP = 1; ID = ID+STEP; subplot(LINES,ROWS,ID) stem(X1)

ID = ID+STEP; subplot(LINES,ROWS,ID) stem(X2)

（3）利用MATLAB中提供的FFT函数，计算2个要求序列的DFT值，并调用MATLAB中提供的计算程序运行时间函数，对比与直接利用DFT定义式计算所花费的时间的程序。

直接调用理论上和上面的效果相同，只要把dft改成fft即可。但是时间上差距较大：

代码如下：

n = 0:255;

L = length(n);

Fs1 = 1;

x1 = cos(5*pi*n/16);

Fs2 = 25000;

x2 = square(2*pi*1000*n/Fs2,50); tic

X1 = fft(x1); toc tic

X2 = fft(x2); toc

LINES = 2; ROWS = 1; ID = 0; STEP = 1; ID = ID+STEP; subplot(LINES,ROWS,ID) stem(X1)

ID = ID+STEP; subplot(LINES,ROWS,ID) stem(X2)

（4）对要求的2个序列，分别采用不同的窗（Bartlett窗，Hann窗，Hamming窗和Blackman窗），利用FFT函数计算DFT值，对比其结果的程序。

直接调用理论上和上面的效果相同，只要把dft改成fft即可。

hold on grid on

n = 0:255; Fs = 1;

L = length(n);

x = cos(5*pi*n/16);

w1 = bartlett(L)'; w2 = hann(L)'; w3 = hamming(L)'; w4 = blackman(L)';

LINES = 5; ROWS = 2; ID = -1; STEP = 2; ID = ID+STEP; subplot(LINES,ROWS,ID) stem(x)

ID = ID+STEP; subplot(LINES,ROWS,ID) stem(x.*w1)

ID = ID+STEP; subplot(LINES,ROWS,ID) stem(x.*w2)

ID = ID+STEP; subplot(LINES,ROWS,ID) stem(x.*w3)

ID = ID+STEP; subplot(LINES,ROWS,ID) stem(x.*w4)

ID = 0;

ID = ID+STEP; subplot(LINES,ROWS,ID) fdtft(x,L,Fs)

ID = ID+STEP; subplot(LINES,ROWS,ID) fdtft(x.*w1,L,Fs)

ID = ID+STEP; subplot(LINES,ROWS,ID) fdtft(x.*w2,L,Fs)

ID = ID+STEP; subplot(LINES,ROWS,ID) fdtft(x.*w3,L,Fs)

ID = ID+STEP; subplot(LINES,ROWS,ID) fdtft(x.*w4,L,Fs)

（5）（拓展要求）分析利用上面的方法画出的信号频谱与理论计算出来的频谱