实习报告(MATLAB).docx

本科生实验报告实习类型 认识实习 题 目 系统的零极点及频率响应特性 学院名称 信息科学与技术学院 专业名称 信息工程 学生姓名 欧建全 学生学号 201513010118 指导教师 贾勇 实习地点 5715 实习成绩 二〇一六年十二月 二〇一六年十二月实验目的1．掌握系统函数零极点的定义。2．熟悉零极点与频率响应的关系。3．掌握极点与系统稳定性的关系。4．状态方程与系统函数的关系。5．在 MATLAB 中实现系统函数与状态方程间的转换。实验原理1. 对(序列)信号x(n)进行ZT, 得X(z), 从而得到它的零极点分布. 2. 对(离散)系统, 求出它的系统函数H(z) , 也可得到它的零极点分布. 3. 按教材(3.6.13)式, 信号或系统的幅度特性由零点至单位圆周上的矢量长度和极点至单位圆周上的矢量长度之比. 4. 极点影响频率特性的峰值, 零点影响频率特性的谷值. 零极逾靠近单位圆, 这些特征越明显. 如有极点, 则频率特性曲线在4πω处出现峰值. 5. 本实验借助于计算机分析信号或系统的频率响应, 目的是掌握用极、零点分布的几何分析法分析频率响应, 实验时需代入相应的X(z) 或H(z) 中, 再在0~2π中等间隔的取点. 如100等分:w=[0:2*pi/100:2*pi], 再用plot等函数作出图形.在MATLAB中，可以用函数[z，p，K]=tf2zp（num，den）求得有理分式形式的系统转移函数的零、极点，用函数zplane（z，p）绘出零、极点分布图；也可以用函数zplane（num，den）直接绘出有理分式形式的系统转移函数的零、极点分布图。?另外，在MATLAB中，可以用函数?[r，p，k]=residuez（num，den）完成部分分式展开计算；可以用函数sos=zp2sos（z，p，K）完成将高阶系统分解为2阶系统的串联。计算零、极点可以用roots函数，若参数为传统函数H(s)的分子多项式系数b，则得到零点；若为分母多项式系数a，则得到极点。MATLAB还提供了zero(sys)和pole(sys)函数直接计算零、极点，其中sys表示系统传递函数。另外，【p,z】=pzmap(sys)函数也具有计算极点p和零点z的功能。不带返回值的pzmap(sys)则绘制出系统的零、极点分布图。?零极点和传递函数的多项式系数一样，可以作为LTI系统的描述方法。MATLAB提供了（b,a）=zp2tf(z,p,k)和【z,p,k】=tf2zp(b,a)两个函数用于在上述两种描述方法之间进行转换，其中k为用零、极点表示传函数时的系统增益。?MATLAB提供了freqs（b，a）函数用来绘制系统的频率响应，包括幅度响应和相位响应，其中b和a分别对应传递函数的分子和分母多项式系数。如果将调用方式改为H=freqs（b，a，omega），则不绘图输出，而是计算抽样点omega处的频响并传递到H中。MATLAB 语言提供了系统函数，零极点和状态方程之间的相互转换语句，也提供了得到系统频率特性的语句： tf2zp：从系统函数的一般形式求出其零点和极点。zp2tf：从零极点求出系统函数的一般式。ss2zp：从状态方程式求系统的零极点。zp2ss：从零极点求系统的状态方程。freqs：由 H(s)的一般形式求其幅频特性和相频特性。实验内容已知系统函数H(s)=s2-0.5s+2/s2+0.4s+1求其零极点图。MATLAB 程序如下：num=[1 －0.5 2]； den=[1 0.4 1]；[z，p]=tf2zp（num，den）；Pzplane (z, p) 运行结果如下：图 1已知系统和状态方程和输出方程X’=(1 0 1 -3)X+(1 0)f Y=[-0.25 1]X求其系统的零极点。 MATLAB 程序如下：A = [1，0；1，-3]； B = [1，0]； C = [-0.25，1]； D = 0;[z，p] =ss2zp (A,B,C,D) zplane (z，p) 已知系统函数H(s)=0.2s2-0.3s+1/s2+0.4s+1求其零极点图。MATLAB 程序如下：num = [0.20.31]； den = [10.41]； w =logspace (-1，1)； freqs(num，den，w)运行结果如下：图 2第四章 心得与体会 ? 就对MATLAB相关的命令操作而言，通过这次实验的亲身操作和实践，学习掌握了许多原本不知道的或者不太熟悉的命令。掌握系统函数零极点的定义，悉零极点与频率响应的关系。掌握极点与系统稳定性的关系，态方程与系统函数的关系。在 MATLAB 中实现系统函数与状态方程间的转换。实习报告6