MATLAB7控制系统分析.ppt

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 on for wn=w num=wn2; den=1,2*kosai*wn,wn2; step(num,den); end hold off，grid on，title(单位阶跃响应); xlabel(时间);ylabel(振幅);,例6.2 现有一连续二阶系统如下所示： 绘制出系统的单位阶跃响应。 解： A

2、=-0.7524 -0.7268;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 某多输入系统

3、如下所示： 试求该二阶系统的单位冲击响应。 解： A=-2.8 -1.5377 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

4、,C,D,0.1);impulse(sys) dimpulse(A,B,C,D); grid on; title(离散冲击响应); 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); s

5、ubplot(2,2,1) dstep(A,B,C,D);grid on;title(离散阶跃响应); xlabel(时间);ylabel(振幅); 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

6、,D,1),zplane(z,p);grid on; title(离散零极点图);xlabel(实部);ylabel(虚部);,例6.10 对于二阶系统如下所示： 试求出周期为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

7、 1.75; u=rand(100,1); dlsim(num,den,u); grid on; title(随机噪声响应); xlabel(时间); ylabel(振幅);,例6.12 设一高阶系统开环的传递函数为： 试求出当输入为幅值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

8、.1 函数指令方式 6.2.2 给予根轨迹的设计工具,例6.13 某离散二阶系统如下： 试绘制该系统的零极点图。 解： n1=0.0001 0.0218 1.0436 9.3599; d1=0.0006 0.0268 0.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

9、; sys=tf(n1,d1); rlocus(sys); title(根轨迹图); xlabel(实部); ylabel(虚部);,例6.15 已知一个单位反馈系统开环传递函数如下： 试绘制系统闭环的根轨迹图，并在根轨迹图上任选一点，计算该点的增益k，以及所有的极点位置。 解： 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),可对指定根计算的增益和根

10、矢量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(虚部);,6.2.2 基于根轨迹的设计工具 打开根轨迹设计器的指令： rltool rltool(sys) 在MATLAB命令窗口中输入“rltool”后，直接打开系统根轨迹设计器；而rltool(sys)函数执行后则打开系统模型sys带根轨迹图的设计器。,例6.18 已知系统开环传递函数如下： 试用根轨迹设计器来查

11、看该系统在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.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); plo

12、t(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 on title(Nyquist曲线图);xlabel(实数轴);ylabel(虚数轴); figure(2);nyquist(num2,den,w);grid on title(Nyquist曲线图);xlabel(实数轴);ylabel(虚数轴);,例6.22 已知离散系统： 绘制出系统的

13、Nyquist曲线，判别闭环系统的稳定性，并绘制出闭环系统的单位冲激响应。 解： num=0.692;den=1 -1.758 0.375; z,p,k=tf2zp(num,den); p figure(1);dnyquist(num,den,0.1); title(离散Nyquist曲线图);xlabel(实数轴);ylabel(虚数轴); figure(2);num2,den2=cloop(num,den); dimpulse(num2,den2);title(离散冲击相应曲线); xlabel(时间);ylabel(幅值);,6.3.2 Bode图分析 例6.23 有典型二阶系统求取不同

14、值时的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); hold on; semilogx(w1,pha); end subplot(2,1,1);grid ontitle(Bode图); xlabel(频率 (rad/sec);ylabel(增益 (dB);

15、 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(-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);den

16、t=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 on title(Bode图);xlabel(频率 (rad/sec);ylabel(增益 (dB); subplot(2,1,2);grid on, xlabel(频率 (rad/sec);ylabel(相位 (deg);hold off;,例6.25 系统传递函数模型为：

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

18、chols 曲线图); 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 系统稳定性的概念 对线性定常系统稳定性的概念是这样定义的：若系统由于受到扰动作用而偏离了原来的平衡状态，当扰动去除后，如果能恢复到原来的平衡状态，则称该系统是

19、稳定的，否则该系统就是不稳定的。 与稳定性有关的概念还包括：控制系统的特征多项式与特征方程，系统稳定性判别的特征方程法，Hurwitz稳定判据。 结构不稳定系统：结构不稳定系统即不改变系统的结构仅仅调整其各个元部件的参数而无法保证稳定的系统。,6.4.2 系统稳定性分析法概述,基于判别特征方程根的方法（或代数稳定判据）：Hurwitz稳定判据，劳斯判据 根轨迹法 Nyquist稳定判据 基于Bode图的对数频率稳定判据,6.4.3 系统稳定性分析的MATLAB实现 例6.28 已知系统的开环传递函数为： 试对闭环系统判别其稳定性。 解： k=100;z=-2; p=0 -1 -20; n1,d

20、1=zp2tf(z,p,k); G=tf(n1,d1); P=n1+d1 roots(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=tf(num,den); rl

21、ocus(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=

22、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(Sy

23、stem is controlled); 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 controlled); 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); end P=s1(3,:);s1(3,:)*A;s1(3,:)*A*A; P1=inv(P),A1=P*A*P1,B1=P*B,例6.37 已知线性系统的动态方程如下，试确