实验三-零极点分布对系统频率响应的影响(数字信号实验-)

实验名称专业、年级学号姓名零极点分布对系统频率响应的影响以下内容由实验指导教师填写〔实验内容请以批注的形式批阅〕实验工程完成情况实验工程成绩指导教师时间年月日备注：〔1〕、按照要求独立完成实验内容。〔2〕、实验结束后，把电子版实验报告按要求格式改名〔例：09号_张三_实验七.doc〕后，实验室统一刻盘留档。实验三零极点分布对系统频率响应的影响一、实验目的学习用分析零极点分布的几何方法分析研究信号和系统频率响应。二、实验原理如果知道信号的Z变换以及系统的系统函数H(z)，可以得到它们的零极点分布，由零极点分布可以很方便地对它们的频率响应进行定性分析。信号的幅度特性由零点矢量长度之积除以极点矢量的长度之积，当频率ω从0变化到2π时，观察零点矢量长度和极点矢量长度的变化，重点观察那些矢量长度较短的情况。另外,由分析知道,极点主要影响频率响应的峰值，极点愈靠近单位圆，峰值愈锋利；零点主要影响频率特性的谷值，零点愈靠近单位圆，谷值愈深，如果零点在单位圆上，那么频率特性为零。根据这些规律可以定性画出频率响应的幅度特性。峰值频率和谷值频率可以近似用响应的极点和零点的相角表示，例如极点π/4，峰值频率近似为π/4，极点愈靠近单位圆，估计法结果愈准确。本实验借助计算机分析信号和系统的频率响应，目的是掌握用极、零点分布的几何分析法分析频率响应，实验时需要将z=ejω代入信号的Z变换和系统函数中，再在0~2π之间，等间隔选择假设干点，并计算它的频率响应。三、实验内容〔包括代码与产生的图形〕要求：不仅打印幅度特性曲线，而且要有系统频率特性的文字分析。1.假设系统用下面差分方程描述：y(n)=x(n)+ay(n-1) 假设a=0.7,0.8,0.9，分别在三种情况下分析系统的频率特性，并打印幅度特性曲线。代码：A=[1,-a];subplot(3,1,3);zplane(B,A);xlabel('ʵ²¿Re');ylabel('Ð鲿Im');title('y(n)=x(n)-ay(n-1)´«Ê亯ÊýÁã¡¢¼«µã·Ö²¼');gridon[H,w]=freqz(B,A,'whole');subplot(3,1,2);plot(w/pi,abs(H),'linewidth',2);gridon;xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|');title('·ùƵÏìÓ¦ÌØÐÔ');axis([0,2,0,6]);subplot(3,1,1);plot(w/pi,angle(H),'linewidth',2);gridon;axis([-0.1,2.1,-3,3]);xlabel('\omega/\pi');ylabel('\phi(\omega)');title('ÏàƵÏìÓ¦ÌØÐÔ');图像：代码：A=[1,-a];subplot(3,1,3);zplane(B,A);xlabel('ʵ²¿Re');ylabel('Ð鲿Im');title('y(n)=x(n)-ay(n-1)´«Ê亯ÊýÁã¡¢¼«µã·Ö²¼');gridon[H,w]=freqz(B,A,'whole');subplot(3,1,2);plot(w/pi,abs(H),'linewidth',2);gridon;xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|');title('·ùƵÏìÓ¦ÌØÐÔ');axis([0,2,0,6]);subplot(3,1,1);plot(w/pi,angle(H),'linewidth',2);gridon;axis([-0.1,2.1,-3,3]);xlabel('\omega/\pi');ylabel('\phi(\omega)');title('ÏàƵÏìÓ¦ÌØÐÔ');图像：代码：B=1;a=0.9A=[1,-a];subplot(3,1,3);zplane(B,A);xlabel('ʵ²¿Re');ylabel('Ð鲿Im');title('y(n)=x(n)-ay(n-1)´«Ê亯ÊýÁã¡¢¼«µã·Ö²¼');gridon[H,w]=freqz(B,A,'whole');subplot(3,1,2);plot(w/pi,abs(H),'linewidth',2);gridon;xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|');title('·ùƵÏìÓ¦ÌØÐÔ');axis([0,2,0,6]);subplot(3,1,1);plot(w/pi,angle(H),'linewidth',2);gridon;axis([-0.1,2.1,-3,3]);xlabel('\omega/\pi');ylabel('\phi(\omega)');title('ÏàƵÏìÓ¦ÌØÐÔ');图像：分析：由y(n)=x(n)+ay(n-1)可知:H[z]=B[z]/A[z]=1/(1-az^(-1))系统极点z=a，零点z=0，当B点从w=0逆时针旋转时，在w=0点，由于极点向量长度最短，形成波峰,并且当a越大，极点越接近单位圆，峰值愈高愈锋利；在w=pi点形成波谷；z=0处零点不影响幅频响应。2.假设系统用下面差分方程描述：y(n)=x(n)+ax(n-1)假设a=0.7,0.8,0.9，分别在三种情况下分析系统的频率特性，并打印幅度特性曲线。a=代码：B=[1,a];subplot(3,1,3);zplane(B,A);xlabel('ʵ²¿Re');ylabel('Ð鲿Im');title('y(n)=x(n)-ax(n-1)´«Ê亯ÊýÁã¡¢¼«µã·Ö²¼');gridon[H,w]=freqz(B,A,'whole');subplot(3,1,2);plot(w/pi,abs(H),'linewidth',2);gridon;xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|');title('·ùƵÏìÓ¦ÌØÐÔ');axis([0,2,0,6]);subplot(3,1,1);plot(w/pi,angle(H),'linewidth',2);gridon;axis([-0.1,2.1,-3,3]);xlabel('\omega/\pi');ylabel('\phi(\omega)');title('ÏàƵÏìÓ¦ÌØÐÔ');图像：代码：B=[1,a];subplot(3,1,3);zplane(B,A);xlabel('ʵ²¿Re');ylabel('Ð鲿Im');title('y(n)=x(n)-ax(n-1)´«Ê亯ÊýÁã¡¢¼«µã·Ö²¼');gridon[H,w]=freqz(B,A,'whole');subplot(3,1,2);plot(w/pi,abs(H),'linewidth',2);gridon;xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|');title('·ùƵÏìÓ¦ÌØÐÔ');axis([0,2,0,6]);subplot(3,1,1);plot(w/pi,angle(H),'linewidth',2);gridon;axis([-0.1,2.1,-3,3]);xlabel('\omega/\pi');ylabel('\phi(\omega)');title('ÏàƵÏìÓ¦ÌØÐÔ');图像：a=0.9代码：B=[1,a];subplot(3,1,3);zplane(B,A);xlabel('ʵ²¿Re');ylabel('Ð鲿Im');title('y(n)=x(n)-ax(n-1)´«Ê亯ÊýÁã¡¢¼«µã·Ö²¼');gridon[H,w]=freqz(B,A,'whole');subplot(3,1,2);plot(w/pi,abs(H),'linewidth',2);gridon;xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|');title('·ùƵÏìÓ¦ÌØÐÔ');axis([0,2,0,6]);subplot(3,1,1);plot(w/pi,angle(H),'linewidth',2);gridon;axis([-0.1,2.1,-3,3]);xlabel('\omega/\pi');ylabel('\phi(\omega)');title('ÏàƵÏìÓ¦ÌØÐÔ');图像：分析：由y(n)=x(n)+ax(n-1)可知:H[z]=B[z]/A[z]=(1-az^(-1))/1系统极点z=0，零点z=a，当B点从w=0逆时针旋转时，在w=0点，由于零点向量长度最长，形成波峰：在w=pi点形成波谷；z=a处极点不影响相频响应。 3.假设系统函数用下式描述：y(ny(ny(n-2)+x(n)+x(n-1)试分析它的频率特性，要求打印其幅度特性曲线，并求出峰值频率和谷值频率。代码：A=[1,-1.273,0.81];B=[1,1];subplot(3,1,3);zplane(B,A);xlabel('ʵ²¿Re');ylabel('Ð鲿Im');title('y(n)=1.273y(n-1)-0.81y(n-2)+x(n)+x(n-1)´«Ê亯ÊýÁã¡¢¼«µã·Ö²¼');gridon[H,w]=freqz(B,A,'whole');subplot(3,1,2);plot(w/pi,abs(H),'linewidth',2);gridon;xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|');title('·ùƵÏìÓ¦ÌØÐÔ');axis([0,2,0,20]);subplot(3,1,1);plot(w/pi,angle(H),'linewidth',2);gridon;axis([-0.1,2.1,-3,3]);xlabel('\omega/\pi');ylabel('\phi(\omega)');title('ÏàƵÏìÓ¦ÌØÐÔ');图像：分析：由y(ny(ny(n-2)+x(n)+x(n-1)可知:z^(-2))系统极点z1=0.79+j0.62*1.62^(-2)，z2=0.79-j0.62*1.62^(-2)零点z1=-1，z2=0当B点从w=0逆时针旋转时，当旋转到接近极点z1=0.79+j0.62*1.62^(-2)是极点向量长度最短，幅度特性出现峰值。当转到w=pi点形成波谷；z=a处零点不影响幅频响应。当旋转到接近极点z2=0.79-j0.62*1.62^(-2)是极点向量长度再次最短，幅度特性再次出现峰值。四、总结1