北京邮电大学-现代信号处理-作业.docxVIP

2012111203班 AaronHwang 学号：2012140619 1.2 设为有限长序列。 （1）计算16点DFT，并画出幅度谱序列。 解：程序代码如下 n=0:15; x=5*cos(0.25*pi*n); figure(1); stem(n,x); xlabel(n);ylabel(x(n)); X=fft(x); X=abs(X); figure(2); stem(n,X); xlabel(k);ylabel(X(K)); 所得图像如下 （2）在给序列后面补16个零后，计算32点DFT，并画出DFT幅度谱序列。 解：程序代码如下 n=0:31; n1=0:15; x1=5*cos(0.25*pi*n1); x=[x1 zeros(1,16)]; figure(1); stem(n,x);xlabel(n);ylabel(x(n)); X=fft(x); X=abs(X); figure(2);stem(n,X); stem(n,Xk);xlabel(k);ylabel(X(k)); 所得图像如下： （3）把DFT的点数扩大为64，然后重复（2） 解：程序代码如下 n=0:63; n1=0:15; x1=5*cos(0.25*pi*n1); x=[x1 zeros(1,48)]; figure(1); stem(n,x);xlabel(n);ylabel(x(n)); X=fft(x); X=abs(X); figure(2);stem(n,X); stem(n,X);xlabel(k);ylabel(X(k)); 所得图像如下： （4）依据DTFT与DFT之间的关系，解释补零操作对DFT的影响。 解：通过补零操作和DFT点数的增加可以更加准确的得出原信号的DFS图像。 1.5（1）设为单零点系统的传递函数，求其幅度响应的解析式。令，画出幅度响应的草图。 解：由已知，可知响应为 程序代码如下 r0=0.95; w0=0.25*pi; w=linspace(0,2*pi,100); H=1-r0*exp(-j*w0)*exp(-j*w); plot(w,abs(H)); title(H（jw）); 所得图像如下： （2）设为单极点系统的传递函数，求其幅度响应的解析式。令，画出幅度响应的草图。 解：由已知， 可知其幅度响应的解析式为 程序代码如下 r0=0.95;w0=0.25*pi; w=linspace(0,2*pi,100); H=1-r0*exp(-j*w0)*exp(-j*w);Hjw=1./H; plot(w,abs(Hjw));title(H（jw）); 所得图像如下 （3）在程序1_5_1中修改的值，然后观察零极点位置的变化，以及幅度响应曲线的变化情况，观察有何规律。 附：程序1_5_1 用MATLAB计算、显示数字滤波器的零极点和幅度响应。 r1=0.95,w1=0.25*pi; r2=0.9,w2=0.5*pi; K=2; b=K*[1 -2*cos(w2) r2*r2]; a=[1 -2*r1*cos(w1) r1*r1]; zplane(b,a); pause; w=linspace(0,pi,512); H=freqz(b,a,w); plot(w/pi,20*log10(abs(h))); xlabel(Normalized frequency(\omega/pi)); ylabel(Mabnitude(dB)); 解：原程序所得图像如下 当时，其他条件不变，可得图像 而当时，其他条件不变，可得图像 而当时，其他条件不变，可得图像 与此类似，通过对其他参数的修改试验，可得以下结论： 当 时，极点在单位圆上； 当 时，极点在圆外，随着 的增大，越来越远离单位圆； 当 时，极点在圆内，随着 的减小，越来越靠近原点； 当 时，零点在单位圆上； 当 时，零点在圆外，随着 的增大，越来越远离单位圆； 当 时，零点在圆内，随着 的减小，越来越靠近原点； 当 时 增大，极点越向（0,1）靠拢； 当 时 增大，零点越向（1,0）靠拢； 幅度响应曲线的极点频率处幅度值出现峰值，而零点频率处幅度值出现谷值 1.8 设全极点滤波器具有下列格型反射系数： 画出该滤波器的直接型结构图和格型结构图。 解：有程序 k1=0.2;k2=0.3;k3=0.5;k4=0.7; K=[k1,k2,k3,k4]; [b,a]=latc2tf(K,allpole) 得到 b = 1 0 0 0 0 a = 1.0000 0.