用MATLAB进行离散系统的Z域分析.docVIP

数字信号处理 课程实验报告 实验指导教师：*** 实验名称 应用MATLAB进行离散系统的Z域分析 专业、班级 - 姓 名 - 实验地点 - 实验日期 - 一、实验内容 给定一个系统函数，求出该系统零极点。画出零极点分布图，单位冲击响应图。 二、实验目的 1、理解零极点图，会通过零极点位置判断该系统是否稳定 2、分析系统的单位冲击响应与幅频响应。 三、实验原理 对于离散系统其系统函数可由差分方程变换，分子的零点为零点，分子的零点为极点，则系统函数的零极点可用MATLAB的多项式求根函数roots( )实现。 求系统函数零极点时，H(Z)是以Z降幂次序排列，系数向量一定要有多项式的最高幂次开始，一直到常数项，缺项要用0补齐。 四、涉及实验的相关情况介绍（包含使用软件或实验设备等情况） MATLAB7.0 此软件是美国MathWorks公司出品的商业数学软件。中文名为“矩阵实验室”，用于算法开发，数据可视化，数据分析以及数值计算的高级技术计算语言和交互式环境。 操作系统为Windows XP 五、实验记录（程序、相关的图形、相关数据记录及分析） 程序： a=[1 -0.5 -0.005 0.03]; b=[0 1 2 1]; figure(1) p=roots(a) q=roots(b) pa=abs(p) p=p; q=q; x=max(abs([p q 1])); x=x+0.1; y=x; clf hold on w=0:pi/300:2*pi; t=exp(i*w); plot(t) axis(square) plot([-x x],[0 0]) plot([0 0],[-y y]) text(0.1,x,jIm[z]) text(y,1/10,Re[z]) plot(real(p),imag(p),x) plot(real(q),imag(q),o) title(pole-zero diagram for discrete system) hold off figure(2) h=impz(b,a); stem(h);xlabel(k);title(μ￥????3??ìó|); [H,w]=freqz(b,a); figure(3) plot(w/pi,abs(H)); xlabel(?μ?ê\omega); title(?μ?ê?ìó|) 相关图形： 六、实验总结 通过matlab仿真输出的图，可以看出这个系统的系统函数存在且连续，极点在圆内，系统稳定。