MATLAB7控制系统分析.ppt

《MATLAB7控制系统分析.ppt》由会员分享，可在线阅读，更多相关《MATLAB7控制系统分析.ppt（43页珍藏版）》请在三一办公上搜索。

1、第八章 MATLAB7控制系统分析,6.1 系统的时域分析,6.1.1 阶跃响应分析6.1.2 冲激响应分析6.1.3 任意输入的时域响应分析,例6.1 典型二阶系统如下所示：绘制出当 时系统的单位阶跃响应。解：w=2:2:12;kosai=0.5;figure(1);hold onfor wn=w num=wn2;den=1,2*kosai*wn,wn2;step(num,den);endhold off，grid on，title(单位阶跃响应);xlabel(时间);ylabel(振幅);,例6.2 现有一连续二阶系统如下所示：绘制出系统的单位阶跃响应。解：A=-0.7524-0.726

2、8;0.7268 0;B=1-1;0 2;C=2.8776 8.9463;D=0;sys=ss(A,B,C,D)step(sys);grid on;title(单位阶跃响应);xlabel(时间);ylabel(振幅);,例6.4 已知某二阶系统如下所示：试求该二阶系统的阶跃响应。解：num=2-3.4 1.5;den=1-1.6 0.8;sys(num,den,0.1)step(sys)dstep(num,den);grid on;title(离散系统阶跃响应);label(时间);ylabel(振幅)；,例6.6 某多输入系统如下所示：试求该二阶系统的单位冲击响应。解：A=-2.8-1.5

3、377 0 1.4;1.5377 0 0 0;1 2.6601-4.8-4.2136;0 0 4.2136 0;B=4 1;2 0;2 1;0 0;C=1 1 0 3;0 2 0 1;D=0-2;-2 0;impulse(A,B,C,D);grid on;title(单位冲击响应);xlabel(时间);ylabel(振幅);,例6.7 某离散二阶系统如下所示：试求该系统的单位冲击响应。解：A=-1.2-0.75;0.75 1.075;B=1;1;C=4 3.7505;D=1;sys=ss(A,B,C,D,0.1);impulse(sys)dimpulse(A,B,C,D);grid on;t

4、itle(离散冲击响应);xlabel(时间);ylabel(振幅);,例6.9 已知某系统如下所示：以t=0.6为取样周期，将系统转换成离散系统，并求出离散系统的单位阶跃响应、冲击响应和零输入响应(x0=1 1 1 1T)。,解：A1=-1.5-0.8 0 0;0.8 0 0 0;0.3 0.4-4.0-1.25;0 0-1.25 0;B1=1 0 1 0;C1=1 2 1 2;D1=0;t=0.6;A,B,C,D=c2dm(A1,B1,C1,D1,t,tusin);subplot(2,2,1)dstep(A,B,C,D);grid on;title(离散阶跃响应);xlabel(时间);y

5、label(振幅);subplot(2,2,2);dimpulse(A,B,C,D);grid on;title(离散冲击响应);xlabel(时间);ylabel(振幅);subplot(2,2,3);x0=1 1 1 1;dinitial(A,B,C,D,x0);grid on;title(零输入响应);xlabel(时间);ylabel(振幅);axis(0 6-0.5 2.5);subplot(2,2,4);z,p,k=ss2zp(A,B,C,D,1),zplane(z,p);grid on;title(离散零极点图);xlabel(实部);ylabel(虚部);,例6.10 对于二阶

6、系统如下所示：试求出周期为4秒的方波输出响应。解：u,t=gensig(square,4,10,0.1);H=tf(2 5 1,1 2 3);tf(0 1-1,1 1 5);lsim(H,u,t);grid on;title(周期为4秒的方波输出响应);xlabel(时间);ylabel(振幅);,例6.11 有二阶系统如下所示：试求出系统对100点随机噪声的响应曲线。解：num=2-6.8 3.6;den=3-4.3 1.75;u=rand(100,1);dlsim(num,den,u);grid on;title(随机噪声响应);xlabel(时间);ylabel(振幅);,例6.12 设

7、一高阶系统开环的传递函数为：试求出当输入为幅值1的方波信号时系统的输出响应。解：num=1.064;den=2-3.685 1.791;u1=ones(1,50),-1*ones(1,50);u=u1,u1,u1;figure(1);dlsim(num,den,u);grid on;title(离散系统仿真);xlabel(时间);ylabel(振幅);,6.2 系统的根轨迹分析,6.2.1 函数指令方式6.2.2 给予根轨迹的设计工具,例6.13 某离散二阶系统如下：试绘制该系统的零极点图。解：n1=0.0001 0.0218 1.0436 9.3599;d1=0.0006 0.0268 0

8、.6365 6.2711;sys=tf(n1,d1);pzmap(sys);p,z=pzmap(sys)title(零极点图);xlabel(实部);ylabel(虚部);,例6.14 某离散二阶系统如下：试绘制该系统闭环的根轨迹图。解：n1=0.0001 0.0218 1.0436 9.3599;d1=0.0006 0.0268 0.6365 6.2711;sys=tf(n1,d1);rlocus(sys);title(根轨迹图);xlabel(实部);ylabel(虚部);,例6.15 已知一个单位反馈系统开环传递函数如下：试绘制系统闭环的根轨迹图，并在根轨迹图上任选一点，计算该点的增益k

9、，以及所有的极点位置。解：n1=1;d1=conv(1 0,conv(0.5 1,4 1);s1=tf(n1,d1);rlocus(s1);k,poles=rlocfind(s1)title(根轨迹图);xlabel(实部);ylabel(虚部);,根据上面选择的极点，利用rlocfind(sys,p),可对指定根计算的增益和根矢量p。n1=1;d1=conv(1 0,conv(0.5 1,4 1);s1=tf(n1,d1);rlocus(s1);k,poles=rlocfind(s1,0.0009+0.7044i)sgrid;title(根轨迹图);xlabel(实部);ylabel(虚部)

10、;,6.2.2 基于根轨迹的设计工具打开根轨迹设计器的指令：rltoolrltool(sys)在MATLAB命令窗口中输入“rltool”后，直接打开系统根轨迹设计器；而rltool(sys)函数执行后则打开系统模型sys带根轨迹图的设计器。,例6.18 已知系统开环传递函数如下：试用根轨迹设计器来查看该系统在k=31时的闭环阶跃给定响应曲线、Bode图和单位冲击给定响应曲线。解：n1=1 1;d1=conv(conv(1 0,1-1),1 4 16);s1=tf(n1,d1);rltool(s1);,6.3 系统的频域分析,6.3.1 频域响应与Nyquist图6.3.2 Bode图分析6.

11、3.3 Nichols图,6.3 系统的频域分析6.3.1 频域响应与Nyquist图例6.19 系统的闭环传递函数如下：请画出系统的幅频特性。解：num=4;den=1 2 4;w=1:0.01:3;g=freqs(num,den,w);mag=abs(g);plot(w,mag);,例6.20 系统的传递函数如下：求当K分别取1700和6300时，系统的极坐标频率特性图。解：k1=1700;k2=6300;w=8:1:80;num1=k1;num2=k2;den=1 52 100 0;figure(1);nyquist(num1,den,w);grid ontitle(Nyquist曲线图

12、);xlabel(实数轴);ylabel(虚数轴);figure(2);nyquist(num2,den,w);grid ontitle(Nyquist曲线图);xlabel(实数轴);ylabel(虚数轴);,例6.22 已知离散系统：绘制出系统的Nyquist曲线，判别闭环系统的稳定性，并绘制出闭环系统的单位冲激响应。解：num=0.692;den=1-1.758 0.375;z,p,k=tf2zp(num,den);pfigure(1);dnyquist(num,den,0.1);title(离散Nyquist曲线图);xlabel(实数轴);ylabel(虚数轴);figure(2);

13、num2,den2=cloop(num,den);dimpulse(num2,den2);title(离散冲击相应曲线);xlabel(时间);ylabel(幅值);,6.3.2 Bode图分析例6.23 有典型二阶系统求取不同值时的Bode图。解：wn=7;kosai=0.1:0.3:2.0;w=logspace(-1,1,100);figure(1);num=wn2;for kos=kosai den=1 2*kos*wn wn2;mag,pha,w1=bode(num,den,w);subplot(2,1,1);hold on;semilogx(w1,mag);subplot(2,1,2

14、);hold on;semilogx(w1,pha);endsubplot(2,1,1);grid ontitle(Bode图);xlabel(频率(rad/sec);ylabel(增益(dB);subplot(2,1,2);grid on,xlabel(频率(rad/sec);ylabel(相位(deg);hold off;,例6.24 系统传递函数模型为：求出有理传递函数的频率响应，然后在同一张图上绘出以四阶Pade近似表示的系统频率响应。解：Pade近似可以表示e的指数次幂，MATLAB代码如下：num=1 1;den=conv(1 2,conv(1 2,1 2);w=logspace(

15、-1,2);t=0.5;m1,p1=bode(num,den,w);p1=p1-t*w*180/pi;n2,d2=pade(t,4);numt=conv(n2,num);dent=conv(den,d2);m2,p2=bode(numt,dent,w);subplot(2,1,1);semilogx(w,20*log(10*m1),w,20*log(10*m2),r-);subplot(2,1,2);semilogx(w,p1,w,p2,r-);subplot(2,1,1);grid ontitle(Bode图);xlabel(频率(rad/sec);ylabel(增益(dB);subplot

16、(2,1,2);grid on,xlabel(频率(rad/sec);ylabel(相位(deg);hold off;,例6.25 系统传递函数模型为：求该离散系统当Ts=0.1s时的Bode图。解：num=0.8 0.15;den=1-0.6 0.52;dbode(num,den,0.1);title(Bode图);ylabel(幅值);grid onxlabel(频率);ylabel(相位);hold off;,6.3.3 Nichols图例6.26 系统的开环传递函数为：试绘制出该系统的Nichols图。解：k1=16.7/0.125;z1=0;p1=-1.25-4-16;num,den

17、=zp2tf(z1,p1,k1);ngrid(new);nichols(num,den);title(Nichols 曲线图);xlabel(开环相位);ylabel(开环增益);,例6.27 有四阶系统：试绘制出该系统的Nichols图，Ts=0.05。解：num=1.6;den=1 1.2 1.24 0.92 0.30;ngrid(new);dnichols(num,den,0.05);title(离散Nichols 曲线图);xlabel(开环相位);ylabel(开环增益);,6.4 系统的稳定性分析,6.4.1 系统稳定性的概念 对线性定常系统稳定性的概念是这样定义的：若系统由于受到

18、扰动作用而偏离了原来的平衡状态，当扰动去除后，如果能恢复到原来的平衡状态，则称该系统是稳定的，否则该系统就是不稳定的。与稳定性有关的概念还包括：控制系统的特征多项式与特征方程，系统稳定性判别的特征方程法，Hurwitz稳定判据。结构不稳定系统：结构不稳定系统即不改变系统的结构仅仅调整其各个元部件的参数而无法保证稳定的系统。,6.4.2 系统稳定性分析法概述,基于判别特征方程根的方法（或代数稳定判据）：Hurwitz稳定判据，劳斯判据根轨迹法Nyquist稳定判据基于Bode图的对数频率稳定判据,6.4.3 系统稳定性分析的MATLAB实现例6.28 已知系统的开环传递函数为：试对闭环系统判别其

19、稳定性。解：k=100;z=-2;p=0-1-20;n1,d1=zp2tf(z,p,k);G=tf(n1,d1);P=n1+d1roots(P),例6.29 已知系统的闭环传递函数为：试对系统闭环判别其稳定性。解：n1=4 300;d1=0.002 0.442 8 300;s1=tf(n1,d1);G=zpk(s1);G.p1,例6.30 已知一个单位负反馈系统开环传递函数为G(s)，试在系统闭环的根轨迹图上选择一点，求出该点的增益k及其系统的闭环极点位置，并判断在该点系统闭环的稳定性。解：num=1 3;den=conv(conv(conv(1 0,1 5),1 6),1 2 2);sys=

20、tf(num,den);rlocus(sys);k,poles=rlocfind(sys)title(根轨迹图);xlabel(实轴);ylabel(虚轴);,例6.31 已知一个单位负反馈系统开环传递函数为G(s)，试计算当k在3337范围内时系统的闭环极点位置，并判断系统闭环的稳定性。解：num=1 3;den=conv(conv(conv(1 0,1 5),1 6),1 2 2);cpole=rlocus(num,den,33:1:37);range=33:1:37;range,cpole,例6.32 已知一个单位负反馈系统开环传递函数为：试用Bode图法判断系统闭环的稳定性。解：num

21、=0 0 0 3.6;den=2 8 5 0;s1=tf(num,den);Gm,Pm,Wcp,Wcg=margin(s1)margin(s1);title(Bode 图);ylabel(幅值);grid on;xlabel(频率)ylabel(相位);,6.5 系统的可控性与可观性分析,6.5.1 系统的可控性分析6.5.2 系统的可观测性分析,例6.33 已知三维状态方程如下：试确定系统的可控性。解：A=1 2 0;1 1 0;0 0 1;B=0 1;1 0;1 1;n=3;CAM=ctrb(A,B);rcam=rank(CAM);if rcam=n disp(System is cont

22、rolled);elseif rcamn disp(System is not controlled);end,例6.34 已知离散系统方程为：假设方程中有：采样周期T=0.1s，试确定离散系统的可控性。,解：f=0.8760 0 0;0.2546 0.6621-0.5701;0.1508 0.4221 1;g=0.2105;0.1033;0.1768;c=0 1 3.5;d=0;n=3;CAM=ctrb(f,g);rcam=rank(CAM);if rcam=n disp(System is controlled);elseif rcamn disp(System is not contro

23、lled);end,例6.35 现有一系统，其状态方程为：其中：试将该系统状态方程转换为可控标准型。解：A=-2 2-2;0-1 0;2-6 1;B=0;1;2;S=ctrb(A,B);if det(S)=0 s1=inv(S);endP=s1(3,:);s1(3,:)*A;s1(3,:)*A*A;P1=inv(P),A1=P*A*P1,B1=P*B,例6.37 已知线性系统的动态方程如下，试确定系统的可控性及可观测性。解：A=-2 3;3-2;B=1 1;1 1;C=2 1;1-2;D=0;n=2;CAM=ctrb(A,B);rcam=rank(CAM);if rcam=n disp(Sys

24、tem is controlled);else disp(Systen is no controlled);endob=obsv(A,C);roam=rank(ob);if roam=n disp(System is observable);elseif roam=n disp(Systen is no observable);end,例6.39 已知连续系统的传递函数模型为：试确定当 时系统的可控性。解：for alph=-1:1 alph num=1,alph;den=2 7 24 19;A,B,C,D=tf2ss(num,den);cam=ctrb(A,B);rcam=rank(cam)

25、,ob=obsv(A,C);roam=rank(ob);n=3;if rcam=n disp(System is controlled);else disp(Systen is no controlled);end if roam=n disp(System is observable);elseif roam=n disp(Systen is no observable);endend,例6.41 已知系统的动态方程为：试将系统动态方程转化为可观测标准型，并求出其变换矩阵。解：A=2-1;1 2;B=-2;1;C=2 3;V=obsv(A,C);rank(V);V1=inv(V);V10=V1;V11=V10(2,:);VA=V11*A;N=V11;VA;M=N,A1=inv(M)*A*MC1=C*M,