matlab的典型应用,MATLAB8.X程序设计及典型应用第八章 MATLAB的典型应用.ppt

【例8-33】图8.46为一典型的反馈控制系统结构，其中 .试分别求系统的开环和闭环单位脉冲响应。 编写文件名为exm8_33的脚本文件： clear G=tf(4,[1 2 3 4]);Gc=tf([1 -3],[1,3]);H=tf(1,[.01,1]); G_o=G*Gc; G_c=feedback(G_o,H); subplot(1,2,1),impulse(G_o); subplot(1,2,2),impulse(G_c); 执行结果如图所示： (a)开环系统 (b)闭环系统 由图可知，开环系统最终稳定于 ， 而闭环系统却不收敛于0；由此可知，控制器和闭环系统结构并不总能够改进系统的控制效果。 ?step(sys) ?step(sys,t) ?[y,t]=step(sys) ?[y,t,x]=step(sys) ?y=step(sys,t) ?step(A,B,C,D) ?step(A,B,C,D,iu) ?step(A,B,C,D,iu,t) ?step(num,den) ?step(num,den,t) ?step(sys1,sys2,...,sysn) ?step(sys1,sys2,...,sysn,t) ?step(sys1,'PlotStyle1',sys2,'PlotStyle2',…,sysn,'PlotStylen') ?用于绘制线性系统的单位阶跃响应曲线，各参数意义与指令impulse()相同。 一、线性系统的时域分析 【例8-34】已知某系统的传递函数为 其中 试绘制其单位阶跃响应曲线。 编写文件名为exm8_34的脚本文件： clear z=[];p=[-.2+.6j -.2-.6j -3];k=2; G=zpk(z,p,k); step(G) title('例8-34：单位阶跃响应曲线实例') grid on 执行结果如图所示： 【例8-35】典型的二阶系统传递函数为 试分析不同参数下的系统单位阶跃响应。 分析：该系统中的参数为系统的固有频率 ， 为系统的阻尼比。 (一)若系统的固有频率为 ，改变阻尼比的值分别为 编写文件名为exm8_35_1的脚本文件： clear,wn=1; xi=[0 0.4 0.8 1,2]; t=0:0.1:12; hold on for i=1:length(xi) Gc=tf(wn^2,[1,2*xi(i)*wn,wn^2]); step(Gc,t) end hold off title('不同阻尼比下系统阶跃响应曲线图');grid 在指令窗中执行该文件结果如图所示： (2)将阻尼比的值固定，取 ，改变系统的固有频率。 编写文件名为exm8_35_2的脚本文件： clear,wn=[0.1:0.2:1,2]; xi=0.55; t=0:0.1:12; hold on for i=1:length(wn) Gc=tf(wn(i)^2,[1,2*xi*wn(i),wn(i)^2]); step(Gc,t) end hold off title('不同固有频率系统阶跃响应曲线图'); grid 在指令窗中执行该文件结果如图所示： 由两图可知，系统的响应特性很大程度上取决于系统的阻尼比和固有频率。固有频率越高，系统响应越快；阻尼比直接影响系统超调量和振荡次数。 【例8-35】典型的二阶系统传递函数为 试分析不同参数下的系统单位阶跃响应。 ?initial(sys,x0) ?initial(sys,x0,t) ?[y,t,x]=initial(sys,x0) ?[y,t,x]=initial(sys,x0,t) ?initial(sys1,sys2,...,sysn,x0) ?initial(sys1,sys2,...,sysn,x0,t) ?initial(sys1,'PlotStyle1',sys2,'PlotStyle2',…,sysn,'PlotStylen',x0) ?计算LTI对象sys的零输入响应，调用格式、参数含义与指令impulse()相似。输入变量x0为系统的初始状态列向量。 注意：该指令只适用于系统的状态空间模型。 一、线性系统的时域分析 【例8-36】已知控制系统的状态空间表达式为 当初始状态为 时，绘制系统的零输入响应曲线。 编写文件名为exm8_36的脚本文件： clear,A=[-0.5573 -.7814;.7814 0]; B=[1 0]';C=[