数字信号处理实验

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

>>grid,xlabel('NORMALIZED FREQUENC'),ylabel('幅度')

得到的结果为：

25

20

15

幅度

1050-0.5

-0.4-0.3

-0.2-0.100.10.2

NORMALIZED FREQUENC

0.30.40.5

图2-1-c L=21

结果分析：（1）由上三图可知，当L为奇数时，零点间距是固定的。 （2）而其峰值即为L的大小。

d.根据以上分析，asinc函数零点的间距为 2 /L，峰值为L。

2.、asinc函数的M文件

编写一个MATLAB的函数如asinc（W，L），直接从（3.6）式计算在频率格上的asinc（W，L）。该函数应有两个输入：长度L和频率W的向量。函数必须检查被零除的情形，如w=0时。

直接计算混叠sinc函数（3.14）式得到脉冲信号的DTFT。绘出幅度，保存该图以便与dtft得到的结果比较。 编写的asinc函数如下： function a=asinc(w,L) if (w==0) a=L;

else a=sin(1/2*w*L)/sin(1/2*w); end

利用混叠sinc函数式（3.14）得到DTFT，M文件： function R=a(w,L)

% w:数字域频率基本频率在[-pi,pi]之间 % L：矩形脉冲的长度

b=sin(1/2*w*L)./sin(1/2*w); % 混叠sinc函数 c=exp(-j*w*(L-1)/2);