数字信号处理实验

W=fftshift(W);

H=fftshift(fft(h,N)); %move negative freq components

要计算12点脉冲信号的DTFT,则取L=12，并依此取不同的N值，N=60,72,84,96,108,120， 当N=60时，在命令窗口输入: >> r=ones(12,1); >> [H,W]=dtft(r,60); >>subplot(311) >>plot(W,real(H));

>>grid,xlabel('w'),ylabel('real(H)') >>subplot(312) >>plot(W,imag(H));

>>grid,xlabel('w'),ylabel('imag(H)') >>subplot(313) >>plot(W,abs(H));

>>grid,xlabel('w'),ylabel('幅度') 得到结果为：

1510

real(H)

50-5-4105

-3

-2

-1

0w

1

2

3

4

imag(H)

0-5-10-41510

-3

-2

-1

0w

1

2

3

4

幅度

50-4

-3-2-1

0w

1234

图2-1-b N=60

当N=72时，

>> [H,W]=dtft(r,72); >>subplot(311)