数字信号处理实验

R=b.*c; % 矩形脉冲的DTFT subplot(211)

plot(w/2/pi,abs(R)) % 频谱 grid

xlabel('NORMALIZED FREQUENCY'),ylabel('|H(w)|'),title('由数学式得出的幅频特性') subplot(212)

plot(w/2/pi,180/pi*angle(R)) %相频特性 grid

xlabel('NORMALIZED FREQUENCY'),ylabel('DEGREES'),title('由数学式得出的相频特性') 取L=5，在命令窗口输入： >> w=[-pi:pi/50:pi]; >> L=5; >>a(w,L); 得到的结果为：

由数学式得出的幅频特性

64

20-0.5

|H(w)|

-0.4-0.3

-0.2-0.100.10.2NORMALIZED FREQUENCY由数学式得出的相频特性

0.30.40.5

200

DEGREES

100

0-100-200-0.5

-0.4

-0.3

-0.2-0.100.10.2NORMALIZED FREQUENCY

0.3

0.4

0.5

图2-2 由数学式得到的结果

通过dtft得到DTFT： >> n=0:4;

>> r=ones(5,1);

>> [H,W]=dtft(r,128); >>subplot(211)

>>plot(W/2/pi,abs(H));

>>grid,title('幅频响应'),xlabel('NORMALIZED FREQUENCY'),ylabel('|H(w)|') >>subplot(212)

>>plot(W/2/pi,180/pi*angle(H));

>>grid,title('相频响应'),xlabel('NORMALIZED FREQUENC'),ylabel('DEGREES')

得到的结果为：