哈工程数字信号处理2013版 实验二 离散时间傅里叶变换.docVIP

实验二 离散时间傅里叶变换 实验目的 离散时间傅里叶变换应该如何建立函数和应用 学习一些新函数如fftshift（），ceil(),fix()以及asinc()是如何使用的 实验原理 1、经由正、逆离散时间傅里叶变换表达的信号傅里叶表示式是信号分析的一个关键部分。 X()=（3.9） （3.10） 类似地，当LTI系统用于滤波时，作为冲击响应离散时间傅里叶变换的频率响应，提供了LTI系统简介的描述。离散时间傅里叶变换X()是的周期复值函数，周期总是2π，并且基周期通常选在区间[-π，π）上。对离散时间傅里叶变换DTFT来说有两个问题： DTFT的定义对无限长信号是有效的。 DTFT是连续变量的函数。 在MATLAB中，任何信号（向量）必须是有限长度的，仅此就是第一点成为问题。因此，不可能使用MATLAB计算无限长信号的DTFT。有一个值得注意的例外情形，当能从变换定义式推导出解析式并只是计算它时，可以使用MATLAB计算无限长信号的DTFT。 2、对于频率抽样问题。MATLAB擅长在有限网格点上计算DTFT。通常选择足够多的频率以使绘出的图平滑，逼近真实的DTFT。对计算有利的最好选择是在（-π，π）区间上一组均匀地隔开的频率，或者对共轭对称变换选择[0,π]区间。采用上述抽样办法，DTFT式变成X()= DTFT的周期性意味着在-π≤0区间上的数值是那些对kN/2的数值。因为上市是在有限数量的频率点=2πk/N处计算，并在有限范围内求和，因此它是可计算的。由于信号长度必须是有限的（0≤nL），这个求和式不适用于x[n]=u[n]的情形。 在对DTFT进行抽样时，并不要求N=L，尽管通常经由DFT进行计算。在正确应用FFT计算N点DFT前，需要对x[n]进行时间混叠。 计算DTFT需要两个函数，MATLAB的freqz函数计算无限长信号，dtft（h，H）函数计算有限长信号的DTFT。 实验内容及结果分析 dtft函数 function[H,W]=dtft(h,N) N=fix(N); L=length(h); h=h(:); if(NL) error(DTFT:#data samples cannot exceed#freq samples) end W=(2*pi/N)*[0:N-1]; mid=ceil(N/2)+1; W(mid:N)=W(mid:N)-2*pi; W=fftshift(W); H=fftshift(fft(h,N)); end asinc函数 function [Y]=asinc(w,L) if(w==[0]) Y=L; else Y=sin(0.5*L.*w)./sin(0.5.*w) end 1.a首先假设r(n)是一个30点的u(n),对u（n）求dtft得到如下 L=30； n=0:L； y=heaviside(n)； [X,W]=dtft(y,128)； subplot(311),plot(W/2/pi,real(X)); xlabel(W),ylabel(real(H1)); grid,title(实部) subplot(312),plot(W/2/pi,imag(X)); xlabel(W),ylabel(img(H1)) grid,title(虚部) subplot(313),plot(W/2/pi,abs(X)); grid,title(幅度) xlabel(W),ylabel(|H1(w)|); 再对作图 n=-pi:0.01*pi:pi;L=30; r=asinc(n,L).*exp(-j*n*(L-1)/2); subplot(211);plot(n,abs(r));grid ylabel(|r|);gtext(w) n=0:29;r=ones(30,1); [X,W]=dtft(r,100); subplot(212);plot(W,abs(X)); ylabel(|x|);gtext(w);grid 结果分析：得到两种图像一致，验证了r(n)的dtft与内插函数的图像一致。 1.b nn=0:11; u=ones(1,12); [X,W]=dtft(u,60); subplot(331),plot(W/2/pi,real(X)); xlabel(W),ylabel(real(H1)); grid,title(实部) subplot(332),plot(W/2/pi,imag(X)); xlabel(W),ylabel(img(H1)) grid,title(虚部) subplot(333),plot(W/2/pi,abs(X)); grid,title(幅度) xlabel(W),ylabel(|H1(w)|); nn=0:11; u=ones(1,12