实验二 离散时间傅里叶变换

1、实验二 离散时间傅里叶变换1 实验原理1、经由正、逆离散时间傅里叶变换表达的信号傅里叶表示式是信号分析的一个关键部分。X()=（3.9）（3.10）类似地，当LTI系统用于滤波时，作为冲击响应离散时间傅里叶变换的频率响应，提供了LTI系统简介的描述。离散时间傅里叶变换X()是的周期复值函数，周期总是2，并且基周期通常选在区间-，）上。对离散时间傅里叶变换DTFT来说有两个问题：（1）DTFT的定义对无限长信号是有效的。（2）DTFT是连续变量的函数。在MATLAB中，任何信号（向量）必须是有限长度的，仅此就是第一点成为问题。因此，不可能使用MATLAB计算无限长信号的DTFT。有一个值得注意的

2、例外情形，当能从变换定义式推导出解析式并只是计算它时，可以使用MATLAB计算无限长信号的DTFT。2、对于频率抽样问题。MATLAB擅长在有限网格点上计算DTFT。通常选择足够多的频率以使绘出的图平滑，逼近真实的DTFT。对计算有利的最好选择是在（-，）区间上一组均匀地隔开的频率，或者对共轭对称变换选择0,区间。采用上述抽样办法，DTFT式变成X()= DTFT的周期性意味着在-<0区间上的数值是那些对k>N/2的数值。因为上市是在有限数量的频率点=2k/N处计算，并在有限范围内求和，因此它是可计算的。由于信号长度必须是有限的（0n<L），这个求和式不适用于xn=un的情形

3、。在对DTFT进行抽样时，并不要求N=L，尽管通常经由DFT进行计算。在正确应用FFT计算N点DFT前，需要对xn进行时间混叠。3、 计算DTFT需要两个函数，MATLAB的freqz函数计算无限长信号，dtft（h，H）函数计算有限长信号的DTFT。2 实验要求理解数值计算在离散时间傅里叶变换中的作用。3 实验内容1. 脉冲信号的DTFT(1)要求：设矩形脉冲rn= 1 0n<L 0 其他a.证明rn的DTFT可由（3.13）得出，记asinc（，L）（3.14）b.使用dtft函数计算12点脉冲信号的DTFT。绘出在区间-<上对的DTFT。把实部和虚部分开绘出。另绘出DTFT的

4、幅度。选择频率样本的数量是脉冲长度的5到10倍，以使绘出的图看上去平滑。用不同数量的频率样本做试验。c.注意asinc函数零点的位置是规律分布的。对奇数长脉冲，比如L=15的脉冲重复进行DTFT计算并绘出幅度，同样再次检验零点位置，注意峰值高度。d.对于asinc函数零点的间距与asinc函数的直流值，确定出通用规则。（2）程序M文件function H,W = dtft(h,N)N=fix(N);L=length(h);h=h(:);if(N<L) errorendW=(2*pi/N)*0:(N-1)'mid=ceil(N/2)+1;W(mid:N)=W(mid:N)-2*pi

5、;W=fftshift(W);H=fftshift(fft(h,N);%bnn=0:11;u=ones(1,12);X,W=dtft(u,72);subplot(221),plot(W,real(X);grid,title('REAL RESPONSE')xlabel('FREQUENCY W'),ylabel('REAL A')subplot(222),plot(W,imag(X);grid,title('IMAGE RESPONSE')xlabel('FREQUENCY W'),ylabel('IMAG

6、E A')subplot(223),plot(W,abs(X);grid,title('MAGNITUDE RESPONSE')xlabel('FREQUENCY W'),ylabel('|H(w)|')subplot(224),plot(W,angle(X);grid,title('PHASE RESPONSE')xlabel('FREQUENCY W'),ylabel('DEGREES')运行结果%cnn=0:14;u=ones(1,15);X,W=dtft(u,90);Y,W=dtft

7、(X,90);subplot(111),plot(W,abs(Y);grid,title('MAGNITUDE RESPONSE')xlabel('FREQUENCY W'),ylabel('|H(w)|')运行结果%d如图L=12时由R(ejw)=0得sin(wL/2)=0即wL/2=k*pi 则w=k*pi/36所以零点间距为pi/6直流值：12零点间距*直流值=(pi/6)*12=2*pi（3） 结果分析使用dtft函数可以快速准确的计算出脉冲信号的DTFT，频率样本的数量越大时，绘出的图形越平滑。2. asinc的M文件（1） 内容编写一

8、个MATLAB的函数如asinc（，L） ，直接从（3.14）式计算在频率格上的asinc(,L)。该函数应有两个输入：长度L和频率的向量。函数必须检查被零清除的情形，如=0时。直接计算混叠sinc函数（3.13）式得到的脉冲信号的DTFT。绘出幅度。保存该图以便将其与用dtft得到的结果进行比较。(2)程序%asinc函数function y=asinc(w,L) N=length(w); for i=1:N if w(i)=0 y(i)=L; else y(i)=sin(1/2*w(i)*L)/sin(1/2*w(i); end end L=12;N=84; W=(2*pi/N)*0:(N

9、-1)' W=W-pi;H=asinc(W,L); figure(3)plot(W,abs(H); grid,title('MAGNITUDE RESPONSE');xlabel('NORMALIZED FREQUENCY'),ylabel('|H(W)|')运行结果3. 无限长信号的DTFT通常不可能计算一个无限长信号的DTFT，但指数信号计算比较容易。当|a|<1, 有 利用freqz函数可以计算上式： 与dtft类似，freqz有两个输出：交换数值（HH）和频率格点（WW），第四个输入参数是可以选择的，但如将其设定为whole

10、,则输出向量WW指定频率格点的范围是从到。如果省略第四个参数，频率格点由区间上等间距的N个点组成。4.指数信号（1）内容：对于信号xn=un,使用freqz函数计算其DTFT。A. 对在区间-<上绘出幅度与相位特性。这需要从freqz返回的X,W向量的移位。解释为什么幅度特性是的偶函数，而相位特性是的奇函数。B. 推算一阶系统的幅度特性与相位特性的表示式。C. 直接以这些表达式来计算幅度特性与相位特性，并与用freqz函数计算出得结果相对比。（2） 程序%aN=1000;a=1,-0.9;b=1;X,W=freqz(b,a,N);W=-pi:0.1:pi;X=freqz(b,a,W);s

11、ubplot(211),plot(W,abs(X);grid,title('MAGNITUDE RESPONSE')xlabel('FREQUENCY W'),ylabel('|X(w)|')subplot(212),plot(W,angle(X);grid,title('PHASE RESPONSE')xlabel('FREQUENCY W'),ylabel('DEGREES')运行结果b假设一阶差分方程：yn=xn+ayn-1,其系统方程H()=其幅度,相位%cw=-pi:0.1:pi;x=1./(1-exp(-j.*w);subplot(211),plot(w,abs(x);grid,title('MAGNITUDE RESPONSE')xlabel('FREQUENCY W'),y