仿真高阶时域系统.doc

现代工程控制理论 ——高阶系统的时域分析 第 PAGE10页 /共 NUMPAGES14页 现代工程控制理论 实验报告 一、实验目的 了解高阶系统的传递函数，对高阶系统进行近似简化 分析高阶系统的时域，探究影响高阶系统动态性能的因素 二、实验原理 高阶系统的闭环传递函数的一般形势可表示为： 当输入为阶跃函数即 左半平面一对非常接近的零点和极点（偶极子）可以相消 假设系统零点与极点相距很近，即很小则有： 左半平面距离虚轴很远的极点可以忽略 实验内容 对传递函数为的高阶系统运用 matlab进行仿真进行以下操作，观察响应曲线的变化： 1、减少极点，观察系统的响应曲线 2、增加偶极子，观察系统的响应曲线 3、增加极点，观察系统的响应曲线 4、增加零点，观察系统的响应曲线 实验方案 运用matlab对传递涵数为的系统 进行仿真并在分别传递函数去掉极点-2和-8后的进行仿真，比较仿真输出的响应曲线 仿真程序如下： den1=conv([1 1 1],conv([0.5 1],[0.125 1])); den1=conv([1 1 1],conv([0.5 1],[0.125 1])); num1=[1.05]; den2=conv([1 1 1],[0.5 1]); den3=conv([1 1 1],[0.125 1]); y1=step(num1,den1,t); %原来系统传涵 y2=step(num1,den2,t); %减少极点-8 y3=step(num1,den3,t); %减少极点-2 [a1,b1,c1,d1,e1,f1,g1]=value(y1,t(2)-t(1)); [a2,b2,c2,d2,e2,f2,g2]=value(y2,t(2)-t(1)); [a3,b3,c3,d3,e3,f3,g3]=value(y3,t(2)-t(1)); figure (1); plot(t,y1,k-,t,y2,r--,t,y3,linewidth,2); s1=[原来的传涵， g1]; s2=[去掉极点-8， g2]; s3=[去掉极点-2， g3]; legend(s1,s2,s3,4); xlabel(t(s)); ylabel(y); title(去掉极点的比较); hold on; 给系统附加附加极点 给系统附加一个极点-0.25，传递函数变为 给系统附加一个极点-16，传递函数变为 仿真程序如下： den1=conv([1 1 1],conv([0.5 1],[0.125 1])); den1=conv([1 1 1],conv([0.5 1],[0.125 1])); den4=conv(den1,[4 1]); den5=conv(den1,[0.0625 1]); num1=[1.05]; t=0:0.01:15; den4=conv([1 1 1],conv(den1,[2 1])); y1=step(num1,den1,t); %原来系统传涵 y4=step(num1,den4,t); %增加极点-0.25 y10=step(num1,den5,t); %增加极点-16 [a1,b1,c1,d1,e1,f1,g1]=value(y1,t(2)-t(1)); [a4,b4,c4,d4,e4,f4,g4]=value(y4,t(2)-t(1)); [a10,b10,c10,d10,e10,f10,g10]=value(y10,t(2)-t(1)); figure (7); plot(t,y1,k-,t,y4,r--,t,y10,linewidth,2); s1=[原来的传涵， g1]; s4=[附加极点-0.25， g4]; s10=[附加极点-16， g10]; legend(s1,s4,s10,4); xlabel(t(s)); ylabel(y); title(附加极点的比较); hold on; 给系统附加附加零点 （1）给系统附加一个零点-1，传递函数变为 （2）给系统附加一个零点-0.5，传递函数变为 仿真程序如下： den1=conv([1 1 1],conv([0.5 1],[0.125 1])); den1=conv([1 1 1],conv([0.5 1],[0.125 1])); num1=[1.05]; num2=[1.05 1.05]; num4=[