数字信号处理实验二离散时间福利叶变换.docxVIP

实验目的：理解数值计算在离散时间傅里叶变换中的作用。掌握MATLAB求解DTFT的方法。实验内容：经由正、负离散时间傅立叶变换表达式是信号分析的一个关键部分。当LTI系统用于滤波的时候，作为冲激响应离散时间傅立叶的频率响应，提供了LTI系统间接的描述。离散时间傅立叶变换X是w的周期复值函数，周期总是2π，并且基周期通常选在区间[-π，π]上。对离散时间傅立叶变换DTFT来说有两个问题：1. DTFT的定义对无限长信号是有效的。2. DTFT是连续变量w的函数。在MATLAB中任何信号向量必须是有限长度的，所以DTFT无法用MATLAB直接计算。当能从变换定义式推导出解析式并计算它时，可以用MATLAB计算。第二个问题是频率抽样问题。MATLAB擅长在有线网络点上计算DTFT，通常选足够多的频率以使绘出的图平滑，逼近真实的DTFT。1、脉冲信号的DTFT设矩形脉冲r[n]由下式定义使用dtft函数计算12点脉冲信号的DTFT。会出区间π≤ω＜π上对ω的DTFT。把实部和虚部分开绘出，但是注意这些图不是很有用。另绘出DTFT的幅度。选择频率样本的数量是脉冲长度的5到10倍。以使绘出的图看上去平滑。用不同数量的频率样本做实验。绘图时，注意要正确的标注频率坐标轴的变量ω。2、指数信号对于信号x[n]=(0.9)nu[n],使用freqz函数计算其DTFT X(ejw)。对ω在区间π≤ω＜π上绘出幅度与相位特性。这需要从freqz返回的[X,W]向量的移位。解释为什么幅度特性是ω的偶函数，而相位特性是ω的奇函数。 [HH,WW]=freqz(b,a,N,’whole’)与dtft类似，freqz有两个输出：变换数值(HH)和频率个点(WW)。第四个输出参数是可选择的，但如将其设定为’whole’，则输出量WW制定频率格点的范围是从ω=0到ω=2π。如果省略第四个参数，频率格点由0≤ω＜π区间上等间距的N个点组成。3、复指数信号如果去a=z0=rejθ为复数，则可应用相同的变换。这一情形之所以重要是因为由此可更深入了解这个复数的幅度(r)和相位(θ)是如何影响DTFT的。再一次取z0=0.95ej3π/11，计算其DTFT，并对ω绘出幅度。注意作为ω的函数，幅度响应的尖峰位于何处。把尖峰位置与z0的极坐标形式联系起来。如果的角度变为θ=3π/5，粗略绘出预期的DTFT幅度，并通过对freqz的计算进行绘图来验证。实验结果1 脉冲信号DTFTfunction [H,W]=dtft(h,N)N=fix(N);L=length(h); h=h(:);if(NL)error(DTFT:#data samples cannot exceed # freq samlples)endW=(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)); xn=ones(1,12);stem(xn); [X,W]=dtft(xn,100);subplot(211),plot(W/2/pi,real(X));grid,title(MAGNITUDE RESPONSE)xlabel(NORMALIZED FREQUENCY),ylabel(Re(H(w))) subplot(212),plot(W/2/pi,imag(X));grid,title(MAGNITUDE RESPONSE)xlabel(NORMALIZED FREQUENCY),ylabel(Im(H(w))) subplot(211),plot(W/2/pi,abs(X)); grid,title(MAGNITUDE RESPONSE) xlabel(NORMALIZED FREQUENCY),ylabel(|H(w)|) subplot(212),plot(W/2/pi,180/pi*angle(X)); grid ,xlabel(NORMALIZED FREQUENCY),ylabel(DEGREES), title(PHASE RESPONSE)2、指数信号b=[1]; a=[1,-0.9]; N=1024; [H,W]=freqz(b,a,N,whole); W=(2*pi/N)*[-(N-1)/2:(N-1)/2]; H=fftshift(H); subplot(211),plot(W/pi,abs(H)); grid,title(Magnitude), xlabel(Normalized Frequency), ylabel(|H(w)|); subplot(212),plot(W/pi,180/pi*angle(H)); grid,tit