第8章 matlab在自动控制原理中的应用.ppt

《第8章 matlab在自动控制原理中的应用.ppt》由会员分享，可在线阅读，更多相关《第8章 matlab在自动控制原理中的应用.ppt（81页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第第8章章 Matlab在自动控在自动控制原理中的应用制原理中的应用n研究线性系统中的时不变系统u系统复杂，具有各种数学描述方法，由各个数学模型求出总的系统数学模型，研究系统特性的同时，注意研究系统中各个部件的状态。u控制系统多为闭环反馈系统，一般研究其开环系统如何反映和影响系统的闭环特性u由于控制系统的应用广泛性，在Matlab的控制工具箱有很多，如鲁棒控制工具箱，神经网络工具箱等线性控制系统的特殊之处8.1 控制工具箱中的LTI对象n线性系统可以采用四种不同方法来进行描述，每种方法又有几个参数矩阵，因此对系统进行调用和计算时都很不方便，基于面向对象的设计思想，建立专用的数据结构类型，将各种

2、模型封装成统一的LTI对象。在一个名字下包含了该系统的所有属性。n8.1.1 LTI 对象的类型和属性对象的类型和属性控制工具箱中一般有三个对象ss、tf、zpk共同属性p239 表8.1 Ts0，Ts1 Td 输入时延 特有属性p240 表8.28.1.2 LTI模型的建立n ndss(a,b,c,d)生成描述状态空间模型n nfilt(num,den)生成DSP形式的离散传递函数(z-1)n nss(a,b,c,d)生成状态空间模型n ntf(num,den,)生成传递函数模型n nzpk(z,p,k)生成零极增益模型n ns1=tf(3,4,5,1,3,5,7,9)n ns2=tf(3,

3、4,5,1,3,5,7,9,0.1,InputName,电流,OutputName,转速)n ns3=tf(4,5,1,5,7,9,InputDelay,0.1,InputName,u,OutputName,y)n ns4=filt(3,4,5,1,3,5,7,9,0.1)n ns5=ss(s1)不同行代表不同输出，不同列代表不同输入n nz=,-0.5;z=,-0.5;n np=0.3,0.1+2i,0.1-2i;p=0.3,0.1+2i,0.1-2i;n nk=2,3;k=2,3;n ns6=zpk(z,p,k,-1)s6=zpk(z,p,k,-1)双输入，单输出，离散系统双输入，单输出，

4、离散系统n nz=;-0.5;z=;-0.5;%单元阵列单元阵列n np=0.3;0.1+2i,0.1-2i;p=0.3;0.1+2i,0.1-2i;n nk=2;3;k=2;3;n ns7=zpk(z,p,k)s7=zpk(z,p,k)单输入，双输出，连续系统单输入，双输出，连续系统8.1.3对象属性的获取和修改对象属性的获取和修改n n对象属性的获取uuget uuset n n模型的参数转换和提取uu dssdatauu ssdatauu tfdatauu zpkdatan nget(s1)get(s1)n nget(s5)get(s5)n ns5.as5.an nset(s1,num,

5、0,1,2,3,4,den,2,4,6,8,10)set(s1,num,0,1,2,3,4,den,2,4,6,8,10)n nget(s1)get(s1)n n用单元阵列的访问方法提取单项属性和对它单独赋值用单元阵列的访问方法提取单项属性和对它单独赋值n ns1.nums1.numn ns1.num:s1.num:n ns1.num=0,5,4,3,2s1.num=0,5,4,3,2n n n ns6.ps6.pn ns6.p2s6.p2n ns6.p2=0.5;0.7s6.p2=0.5;0.72模型类型的参数转换和提取第六章中采用转换命令：第六章中采用转换命令：ss2tfss2tf，sst

6、zpsstzp，tf2zptf2zp，tf2sstf2ss，zp2tfzp2tf，zp2sszp2ss等。用这些命令时，输入变等。用这些命令时，输入变元中要键入系数矩阵，不太方便。在采用元中要键入系数矩阵，不太方便。在采用LTILTI模型以后，就不再用这些命令来进行模型变换模型以后，就不再用这些命令来进行模型变换了，而用能直接调用系统的了，而用能直接调用系统的LTILTI名称的命令来名称的命令来实现这些转换。这些命令就是实现这些转换。这些命令就是dssdatadssdata，ssdatassdata，tfdatatfdata和和zpkdatazpkdata，它们分别用来获得转换后的系，它们分别

7、用来获得转换后的系统状态空间、传递函数和零极增益参数。与统状态空间、传递函数和零极增益参数。与ssss，tf tf，zpkzpk命令的不同在于这些带命令的不同在于这些带datadata的命令仅的命令仅仅用来转换参数，但并不生成新的系统。仅用来转换参数，但并不生成新的系统。n nf1,g1=tfdata(s1)n nf11,g11n n n nz1,p1,k1,T1s=zpkdata(s1)n nz11,p11n n n na2,b2,c2,d2,Ts2,InputDelay=ssdata(s2)n n模型类型的检验cs1=class(s1)isa(s1,tf)isctisdtissison n

8、输出反馈n n反馈联接n n系统串联联接n n系统并联联接n n系统增广联接系统模型的连接例6.19 系统的串联、并联和反馈 系统的串联系统的串联由图所示，由图所示，Y YB B=WWB BU UB B=WWB BWWA AU UA A=WUWU故故W(s)=WW(s)=WA A(s)W(s)WB B(s)(s)多项式相乘由卷积函数多项式相乘由卷积函数convconv实现，其表示式为：实现，其表示式为：f=conv(fAf=conv(fA，fB)fB)，g=conv(gAg=conv(gA，gB)gB)系统的并联系统的并联 Y Y=WWA AU U+WWB BU U=(=(WWA A+WWB

9、B)U U=WUWU故故 W(s)=WA(s)+WB(sW(s)=WA(s)+WB(s)f=polyadd(conv(fA f=polyadd(conv(fA，gB)gB)，conv(fBconv(fB，gA)gA)g=conv(gA g=conv(gA，gB)gB)例6.19 系统的串并联和反馈 系统的反馈系统的反馈系统的连接方法如图系统的连接方法如图6.18-36.18-3。复合系统的传递函。复合系统的传递函数数故故MATLABMATLAB表达式为表达式为 f=conv(fAf=conv(fA，gB)gB)g=polyadd(conv(fA g=polyadd(conv(fA，fB)fB)

10、，conv(gAconv(gA，gB)gB)n n%用传递函数法写出这两个系统的描述参数：用传递函数法写出这两个系统的描述参数：n nfA=5,10;gA=5,2,1;fA=5,10;gA=5,2,1;n nfB=4;gB=1,1,0;fB=4;gB=1,1,0;n n%两环节串联后合成的传递函数两环节串联后合成的传递函数fh1,gh1fh1,gh1为为n nfh1=conv(fA,fB);fh1=conv(fA,fB);n ngh1=conv(gA,gB);gh1=conv(gA,gB);n ndisp(disp(串联后的传递函数串联后的传递函数)n nprintsys(fh1,gh1,s)

11、printsys(fh1,gh1,s)%两环节并联后合成的传递函数两环节并联后合成的传递函数fh2,gh2fh2,gh2为为n nfh2=polyadd(conv(fA,gB),conv(gA,fB);fh2=polyadd(conv(fA,gB),conv(gA,fB);n ngh2=conv(gA,gB);gh2=conv(gA,gB);n ndisp(disp(并联后的传递函数并联后的传递函数)n nprintsys(fh2,gh2,s)printsys(fh2,gh2,s)n n%将将B B环节放在负反馈支路上后合成的传递函数环节放在负反馈支路上后合成的传递函数fh3,gh3fh3,g

12、h3为为n nfh3=conv(fA,gB);fh3=conv(fA,gB);n ngh3=polyadd(conv(fA,fB),conv(gA,gB);gh3=polyadd(conv(fA,fB),conv(gA,gB);n ndisp(disp(将将B B环节放在反馈支路上后的传递函数环节放在反馈支路上后的传递函数)n nprintsys(fh3,gh3,s)printsys(fh3,gh3,s)输出反馈系统的结构图下图所示。可以使用函数cloop得到闭环系统的数学模型。格式1：numc,denc=cloop(num,den,sign)输入开环系统的传递函数，左变量返回系统的闭环参数，

13、numc为分子多项式系数向量，denc为分母多项式系数向量；右变量中的sign=1为正反馈，sign=-1为负反馈。符号向量sign的缺省值为-1。+-uy输出反馈格式2：Ac,Bc,Cc,Dc=cloop(A,B,C,D,sign)MIMO系统时，状态空间模型的全输出反馈。格式3：Ac,Bc,Cc,Dc=cloop(A,B,C,D,outputs,inputs)Ac,Bc,Cc,Dc=cloop(A,B,C,D,outputs,inputs)MIMO系统时，状态空间模型的选择输出反馈。向量outputs，inputs分别指明输出输入端口号数，正反馈时inputs取正值，负反馈时inputs取

14、负值。例：单位反馈系统结构图如图所示，求闭环系统的数学模型。+-uynum=10;den=1,2,0;nc,dc=cloop(num,den,-1);printsys(nc,dc)num/den=10-s2+2 s+10一般反馈系统的结构图如图所示。可以使用函数feedback得到闭环系统的数学模型。格式1：numc,denc=feedback(num1,den1,num2,den2,sinumc,denc=feedback(num1,den1,num2,den2,sign)gn)系统1和系统2均为传递函数模式时，左变量为返回的闭环系统参数；右变量中，num1 和den1为系统1参数向量，nu

15、m2和den2为系统2参数向量，符号变量sign=1为正反馈，+-uy反馈联接反馈联接符号向量符号向量sign=-1sign=-1为负反馈。符号向量为负反馈。符号向量signsign的缺省的缺省值为值为-1-1。格式格式2 2：Ac,Bc,Cc,Dc=feedback(A1,B1,C1,D1,A2,B2,CAc,Bc,Cc,Dc=feedback(A1,B1,C1,D1,A2,B2,C2,D2,sign)2,D2,sign)系统系统1 1与系统与系统2 2均为状态空间模式时的使用格式。均为状态空间模式时的使用格式。格式格式3 3：Ac,Bc,Cc,Dc=feedback(A1,B1,C1,D1

16、,A2,B2,CAc,Bc,Cc,Dc=feedback(A1,B1,C1,D1,A2,B2,C2,D2,outputs1,inputs1)2,D2,outputs1,inputs1)MIMOMIMO系统时，系统系统时，系统1 1与系统与系统2 2为状态空间模型的为状态空间模型的反馈联接。向量反馈联接。向量outputs1outputs1和和inputs1inputs1分别指明当系分别指明当系统统2 2作为反馈模块时，联接系统作为反馈模块时，联接系统1 1的输出输入端口的输出输入端口号数，正反馈时号数，正反馈时inputs1inputs1取正值，负反馈时取正值，负反馈时inputs1input

17、s1取负值。取负值。例：反馈控制系统结构图如下图所示，求例：反馈控制系统结构图如下图所示，求闭环系统的数学模型。闭环系统的数学模型。n1=10;d1=1,2,0;n1=10;d1=1,2,0;n2=0.2,1;d2=0.01,1;n2=0.2,1;d2=0.01,1;n,d=feedback(n1,d1,n2,d2,-1);n,d=feedback(n1,d1,n2,d2,-1);printsys(n,d)printsys(n,d)num/den=num/den=0.1 s+10 0.1 s+10 -0.01 s3+1.02 s2+4 s+10 0.01 s3+1.02 s2+4 s+10+-

18、uy系统串联联接的结构图如下图所示。在使用过程中可用函数series得到串联系统的数学模型。格式1：num,den=series(num1,den1,num2,den2)系统1和系统2均为多项式模型时，左变量为返回的闭环系统参数。右变量中num1和den1为系统1参数向量；右变量中num2和den2为系统2参数向量。uy系统串联联接格式2：A,B,C,D=series(A1,B1,C1,D1,A2,B2,C2,D2)A,B,C,D=series(A1,B1,C1,D1,A2,B2,C2,D2)系统1与系统2均为状态空间模式时的使用格式。格式3：A,B,C,D=series(A1,B1,C1,D

19、1,A2,B2,C2,D2,outputs1,inputs2)MIMO系统时，系统1与系统2之间的串联联接方式，向量outputs1指明系统1的输出号数，inputs2指明系统2的输入号数。例：已知系统串联校正结构图如下所示，求系统的传递函数。n1=0.5,1;d1=0.1,1;n1=0.5,1;d1=0.1,1;n2=10;d2=1,2,0;n2=10;d2=1,2,0;n0,d0=series(n1,d1,n2,d2);n0,d0=series(n1,d1,n2,d2);n3=1;d3=1;n3=1;d3=1;nc,dc=feedback(n0,d0,n3,d3)nc,dc=feedbac

20、k(n0,d0,n3,d3)printsys(nc,dc)printsys(nc,dc)num/den=num/den=5 s+10 5 s+10 -0.1 s3+1.2 s2+7 s+10 0.1 s3+1.2 s2+7 s+10uy+-系统并联联接的结构图如图所示，可以使用函数parallel得到闭环系统的数学模型。具体使用格式如下：格式1：num,den=parallel(num1,den1,num2,den2)num,den=parallel(num1,den1,num2,den2)系统1和系统2均为传递函数时，左变量为返回的闭环系统参数；右变量中，num1,den1为系统1参数向量，

21、num2,den2为系统2参数变量。uy系统并联联接格式2：A,B,C,D=parallel(A1,B1,C1,D1,A2,B2,C2,D2)系统1与系统2均为状态空间模式时的使用格式。格式3：A,B,C,D=parallel(A1,B1,C1,D1,A2,B2,C2,D2,inp1,inp2,out1,out2)MIMO系统时，系统1与系统2之间的并联联接方式，系统1的inp1和系统2的inp2连接，系统1的out1和系统2的out2连接。例：系统结构图如图所示，求并联系统的数学模型。Review:LTIReview:LTI模型的简单组合与复杂模型组合模型的简单组合与复杂模型组合(1)若假定

22、两环节均为单输入单输出的系统若假定两环节均为单输入单输出的系统SA和和SB。两个环节级联：两个环节级联：sysseries(SA，SB)两个环节并联：两个环节并联：sys=parallel(SA，SB)A环节前向，环节前向，B环节反馈：环节反馈：S=feedback(SA,SB)(2)当当在在多多输输入入多多输输出出系系统统中中，必必须须增增加加输输入入变变量量和和输输出出变变量量的编号：的编号：级联：级联：sys=series(SA，SB，outputA,inputB)并联：并联：sys=parallel(SA,SB,InputA,InputB,OutputA,OutputB)反馈：反馈：s

23、ys=feedback(SA,SB,feedout,feedin,sign)例例 计算图5.1所示的系统的传递函数。MATLAB源程序为：s1=tf(2,5,1,1,2,3)%系统s1的传递函数模型s2=zpk(-2,-10,5)%系统s2的零极点增益模型sys=feedback(s1,s2)%s1环 节 前 向，s2环 节 反 馈5(s+2)/(s+10)程序运行结果为：程序运行结果为：Transfer function:系统s1的传递函数模型2 s2+5 s+1-s2+2 s+3Zero/pole/gain:系统s2的零极点增益模型5(s+2)-(s+10)Zero/pole/gain:系

24、统s1、s2的反馈零极点增益模型0.18182(s+10)(s+2.281)(s+0.2192)-(s+3.419)(s2+1.763s+1.064)2LTI模型的复杂模型组合模型的复杂模型组合 对复杂系统的任意组合，在MATLAB中，则采用集成的软件包，让机器自动去完成复杂的组合，人们只要输入各环节的LTI模型和相应的联接矩阵与输入矩阵，指定输出变量，软件包会自动判别输入的模型表述方式，作出相应的运算并最后给出组合后系统的状态方程。在求解过程中，主要涉及append()函数和函数和connect()函数函数。通常，由以下五个步骤来完成：通常，由以下五个步骤来完成：对方框图中的各个环节进行编号

25、，建立它们的对象模型。利用append函数命令建立无连接的状态空间模型。sap=append(s1,s2,sm)按规定写出系统的互联接矩阵q互联矩阵q中的每一行由组合系统的一个输入编号和构成该输入的其它输出编号组成，其中该行的第一个元素为该输入的编号，接下来的元素则由构成该输入的其它子框的输出编号组成，如果为负反馈，则编号应取负号。选择组合系统中需保留的对外的输入和输出端的编号并列出。Inputs=i1,i2,outputs=j1,j2,用connect命令生成组合后的系统。实验 9 模型转换及连接n nP256 P256 例例8.1 SIMO8.1 SIMO系统模型转换系统模型转换n n注：

26、注：MISOMISO系统模型转换（双输入单输出）系统模型转换（双输入单输出）na=-0.5572-0.7814;0.7814 0;nb=1-1;0 2;nc=1.9691 6.4493;nd=0 0;nnum,den=ss2tf(a,b,c,d,1)nnum,den=ss2tf(a,b,c,d,1)ns1=ss(a,b,c,d)ns2=tf(s1)ns3=zpk(s1)n nP258 例8.2二、控制系统分析(课本8.2，8.3节)uu控制系统时域分析uu控制系统根轨迹分析uu控制系统频域分析uu控制系统稳定性分析n n控制系统时域响应分析阶跃响应函数功能：给定系统数学模型，求系统的单位阶跃响

27、应。格式1：step(num,den)or step(A,B,C,D)格式2：step(num,den,t)or step(A,B,C,D,t)格式3：y=step(num,den)or y,x=step(A,B,C,D)格式4：y,t=step(num,den)or y,x=step(A,B,C,D,t)其中，num和den对应TF模型的系数向量。A,B,C,D对应SS模型的系数矩阵。格式1：给定num,den（或者A,B,C,D），求系统的阶跃响应并作图，时间向量t的范围自动设定。格式2：时间向量t的范围由人工设定，等间隔。例：t=0:0.1:3。格式3：返回变量格式。返回输出变量y、状态

28、变量x，不作图。格式4：返回变量格式。包括或不包括时间向量t，不作图。例：系统的传递函数为：求阶跃响应并作系统性能分析。%绘制阶跃响应曲线num=4;den=1 1 4;step(num,den)%计算峰值y,t=step(num,den);max(y)ans=1.4419%求峰值时间tp=spline(y,t,max(y)tp=1.5708脉冲响应函数函数功能：给定系统数学模型，求系统的单位脉冲响应。impulse(num,den)or impulse(A,B,C,D)impulse(num,den,t)or impulse(A,B,C,D,t)y,x=impulse(num,den)or

29、y,x=impulse(A,B,C,D)y,x,t=impulse(num,den)or y,x=impulse(A,B,C,D,t)其中，num和den对应TF的系数向量。A,B,C,D对应SS模型的系数矩阵。格式1：给定num,den(A,B,C,D)，求系统的脉冲响应并作图，时间向量t的范围自动设定。格式2：时间向量t的范围由人工设定，等间隔。例：t=0:0.1:3。格式3：返回变量格式。返回输出变量y、状态变量x至MATLAB窗口，不作图。格式4：返回变量格式。包括或不包括时间向量t，不作图。例：系统的传递函数为：例：系统的传递函数为：求脉冲响应并作系统性能分析。求脉冲响应并作系统性能

30、分析。绘制脉冲响应num=4;den=1 1 4;impulse(num,den)计算误差面积y,x,t=impulse(num,den);trapz(t,y)ans=0.9927根轨迹n当开环系统中某个(或几个)参数从0到+连续变化 时,系统闭环特征方程的根(即闭环极点)在根平面(S平面)上连续移动而 形成的轨迹称为系统的根轨迹.根轨迹绘制的基本法则根轨迹绘制的基本法则n n法则法则法则法则1 1 根轨迹的起点和终点根轨迹的起点和终点根轨迹的起点和终点根轨迹的起点和终点:始于开环极点，始于开环极点，始于开环极点，始于开环极点，终于开环零点或无穷远。终于开环零点或无穷远。终于开环零点或无穷远。

31、终于开环零点或无穷远。n n法则法则法则法则2 2 根轨迹的分支数和对称性：根轨迹的分支数和对称性：根轨迹的分支数和对称性：根轨迹的分支数和对称性：max(m,n)max(m,n)n n法则法则法则法则3 3 实轴上的根轨迹实轴上的根轨迹实轴上的根轨迹实轴上的根轨迹n n法则法则法则法则4 4 根轨迹的渐近线根轨迹的渐近线根轨迹的渐近线根轨迹的渐近线n n法则法则法则法则5 5 根轨迹的分离点根轨迹的分离点根轨迹的分离点根轨迹的分离点n n法则法则法则法则6 6 起始角与终止角起始角与终止角起始角与终止角起始角与终止角n n法则法则法则法则7 7 根轨迹与虚轴的交点根轨迹与虚轴的交点根轨迹与虚

32、轴的交点根轨迹与虚轴的交点n n根轨迹作图与系统根轨迹分析函数功能：给定系统的开环模型，绘制根轨迹图或计算绘图变量。用法：rlocus(num,den)or rlocus(num,den,k)r=rlocus(num,den,k)orr,k=rlocus(num,den)rlocus(a,b,c,d)or rlocus(a,b,c,d,k)r=rlocus(a,b,c,d,k)or r,k=rlocus(a,b,c,d)+-uy格式1：根据传递函数模型绘制系统的根轨迹图。格式rlocus(num,den)中，k为机器自适应产生的从0的增益向量；格式rlocus(num,den,k)中，k为人工

33、给定的增益向量。格式2：返回变量格式，不作图。k为返回增益值，r为返回的闭环根矩阵，矩阵行数为length(k)，其列数为(length(den)-1)。格式3：绘制SISO系统的根轨迹图；格式rlocus(a,b,c,d)中，k为机器自适应产生的从0的增益向量；格式rlocus(a,b,c,d,k)中，k为人工给定的增益向量。格式4：状态空间模型时返回变量格式。例：控制系统为绘制系统的根轨迹图.。%绘制系统的根轨迹图num=1;den=1,2,10,0;rlocus(num,den)系统性能分析：在复平面上绘出闭环零、极点随某个可调参数变化的轨迹，根据零、极点的分布定性分析出来可能合适的零极

34、点位置，以便选定可调的参数，使闭环系统不仅稳定，而且满足动态性能指标例：已知系统的开环传递函数为试绘制该系统的根轨迹，并确定使闭环系统稳定的Kg范围。绘制根轨迹n=1 1;d1=1 0;d2=1-1;d3=1 4 16;d=conv(d1,conv(d2,d3);rlocus(n,d)k1,p=rlocfind(n,d)Select a point in the graphics windowselected_point=-0.0046+1.5439ik1=23.1109p1=-1.5058+2.7146i -1.5058-2.7146i 0.0058+1.5486i 0.0058-1.548

35、6ik2,p2=rlocfind(n,d)Select a point in the graphics windowselected_point=-0.0046+2.5965ik2=36.0909p2=0.0138+2.5855i 0.0138-2.5855i -1.5138+1.7627i -1.5138-1.7627i通过上述分析，可以得出系统稳定的Kg范围为23.1109Kg36.0909函数rlocfind()功能：用鼠标确定根轨迹上某一点的增益值和该点对应的n个闭环根。k,poles=rlocfind(num,den)or r,poles=rlocfind(a,b,c,d)k,pol

36、es=rlocfind(num,den,p)or k,poles=rlocfind(a,b,c,d,p)格式1：在已经作出的根轨迹图上，用鼠标选择闭环极点的位置后，返回对应闭环极点的根轨迹增益值k和对应的n个闭环根，并在图上标注十字。使用该命令前，需要用rlocus()指令作出根轨迹图。格式2：给定右变量p为根轨迹上某点坐标值，返回对应该点的根轨迹增益值k和对应的n个闭环根。函数pzmap()功能：给定系统数学模型，作出零极点位置图。pzmap(num,den)or pzmap(a,b,c,d)or pzmap(p,z)p,z=pzmap(num,den)or p,z=pzmap(a,b,c,

37、d)格式1：零极点绘图指令。零点标记为，极点标记为。格式2：返回零极点值，不作图。结论结论:(1)系统的频率特性系统的频率特性与传递函数和微分方程与传递函数和微分方程 一一对应一一对应,它从频率的角度描述系统的特性它从频率的角度描述系统的特性.(2)当系统或环节的输入信号是正弦信号时当系统或环节的输入信号是正弦信号时,其稳其稳 态输出仍为与输入同频率的正弦信号态输出仍为与输入同频率的正弦信号.(3)此同频率的正弦输出信号的幅值与输入正弦信此同频率的正弦输出信号的幅值与输入正弦信号的幅值之比等于幅频特性号的幅值之比等于幅频特性(4)稳态同频率的正弦输出信号的初相角与输入正稳态同频率的正弦输出信号

38、的初相角与输入正弦信号的初相角之差为相频特性弦信号的初相角之差为相频特性频率特性频率特性是指线性系统或环节在正弦信号作用下是指线性系统或环节在正弦信号作用下是指线性系统或环节在正弦信号作用下是指线性系统或环节在正弦信号作用下,稳态输出与稳态输出与稳态输出与稳态输出与输入之比对频率的关系特性输入之比对频率的关系特性输入之比对频率的关系特性输入之比对频率的关系特性.n n控制系统频域响应分析MATLAB提供了多个用于系统频域分析的函数。波特函数Bodebode(num,den)or bode(a,b,c,d,iu)bode(num,den,w)or bode(a,b,c,d,iu,w)m,p,w=

39、bode(num,den)or m,p,w=bode(a,b,c,d,iu)函数功能：对数频率特性作图函数，即波特图绘图。格式1：给定开环系统的数学模型作波特图，频率向量w自动给出。当开环模型以状态空间模型给出时须指定第几个输入。格式2：给定开环系统的数学模型作波特图，频率向量w由人工给出。w的单位为rad/s，可以由命令logspace得到对数等分的w值。格式3：返回变量式，不作图。其中，m为频率特性G(jw)的幅值变量，m=|G(jw)|；p为频率特性G(jw)的幅角向量，p=argG(jw)，单位为()；w为频率向量，单位为rad/s。例：控制系统的开环传递函数为作波特图，并确定谐振峰值

40、的大小Mr与谐振频率wr。n=10;d=1 2 10;bode(n,d);m,p,w=bode(n,d);mr=max(m)mr=1.6667wr=spline(m,w,mr)wr=2.8284奈奎斯特稳定判据n n根据开环频率特性在复平面上绘出幅相轨迹，根据开环的nyquist曲线，可判断闭环系统的稳定性。n n反馈控制系统稳定的充要条件是nyquist曲线按逆时针包围临界点(-1,j0)p圈，其中p为开环传递函数位于右半s平面的极点数，否则闭环系统不稳定。奈奎斯特函数nyquist()函数功能：奈奎斯特轨线作图命令，即极坐标图。nyquist(num,den)or nyquist(a,b,

41、c,d)nyquist(num,den,w)or nyquist(a,b,c,d,w)re,im,w=nyquist(num,den)or re,im,w=nyquist(a,b,c,d)格式1：给定开环系统的数学模型作nyquist图。频率w的范围自动给定。格式2：给定开环系统的数学模型作nyquist图。频率w的范围人工给定。格式3：返回极坐标图参数变量，不作图。其中re为复变函数G(jw)的实部变量，re=ReG(jw);im为复变函数G(jw)的虚部变量，im=ImG(jw);w为频率向量，单位为rad/s。例：二阶系统令n=1，分别作出=2,1,0.707,0.5时的奈奎斯特轨线。n

42、=1;d1=1 4 1;d2=1 2 1;d3=1 1.414 1;d4=1 1 1;nyquist(n,d1);hold onnyquist(n,d2);nyquist(n,d3);nyquist(n,d4);axis(equal);x=0:0.1:2*pi;plot(sin(x),cos(x),:)例：分别由例：分别由w w的自动变量与人工变量作下的自动变量与人工变量作下述系统的奈奎斯特轨线。述系统的奈奎斯特轨线。n=1;n=1;d=1 1 0;d=1 1 0;nyquist(n,d);nyquist(n,d);figure(2)figure(2)w=0.5:0.2:3;w=0.5:0.2

43、:3;nyquist(n,d,w);nyquist(n,d,w);axis(equal);axis(equal);有时为了使曲线变化趋势看起来更明显，常有时为了使曲线变化趋势看起来更明显，常采用人工变量。采用人工变量。稳定裕度及其分析方法n n奈奎斯特稳定判据能够根据开环频率特性来判别闭环系统的稳定性，同时给出了稳定的边界条件，但设计系统时，仅仅知道是否稳定还不够，还要了解稳定的程度，使设计留有一定的余地，在系统实际受到扰动或自身参数发生变化时，仍然能够保持稳定性，这就是研究稳定裕度的意义稳定裕度n n既然既然(-1,j0)(-1,j0)点在判断系统闭环稳定性方面有着点在判断系统闭环稳定性方面

44、有着特殊的意义，那么就可以利用开环频率特性离特殊的意义，那么就可以利用开环频率特性离(-1,j0)(-1,j0)点的距离（远近程度）表示系统的相对点的距离（远近程度）表示系统的相对稳定性稳定性稳定裕度稳定裕度n n当当 时，相位差与时，相位差与-180-180度差多少度差多少 n n当当 角角 时，幅值比与时，幅值比与1 1差多少差多少 稳定裕度函数margin函数功能：计算系统的稳定裕度，包括相位裕度Gm和幅值裕度Pm。margin(num,den)or margin(a,b,c,d)Gm,Pm,wg,wp=margin(num,den)or Gm,Pm,wg,wp=margin(a,b,c

45、,d)Gm,Pm,wg,wp=margin(m,p,w)格式1：给定开环系统的数学模型作波特图，并在图上标注幅值裕度Gm和对应的频率wg，相位裕度Pm和对应的频率wp。幅值裕度Gm定义为：单位为dB(分贝)。相位裕度Pm定义为单位为度()。格式2：返回变量格式，不作图。返回幅值裕度Gm和对应的频率wg，相位裕度Pm和对应的频率wp。格式3：给定频率特性的参数向量、幅值m、相位p和频率w，由插值法计算幅度裕度幅值裕度Gm和对应的频率wg，相位裕度Pm和对应的频率wp。例：控制系统结构图如图所示。求系统的稳定裕度，并分别采用格式2与格式3比较计算误差。n=2;d=1 3 2 0;margin(n,

46、d)gm1,pm1,wg1,wp1=margin(n,d);m,p,w=bode(n,d);gm2,pm2,wg2,wp2=margin(m,p,w);gm1,gm2;pm1,pm2;wg1,wg2;wp1,wp2+-uyn n控制系统的稳定性分析借助于MATLAB语言工具，控制系统的稳定性分析采用直接方法，成为异常方便，异常简单的问题。因此，控制理论中经常使用的几种间接分析方法，如劳斯判据、赫尔维斯判据和奈奎斯特判据等，这里从时域和频域两个方面介绍MATLAB语言的应用。uu时域稳定性分析uu频域稳定性分析时域稳定性分析在数学模型的基础上，可采用直接求根法确定系统的稳定性。在前述课程的学习过

47、程中，可以知道系统的稳定性取决于系统特征根的正负，这样要求取系统的根，在一般的求解过程中比较繁琐，而MATLAB恰恰提供了这方面的便利。系统的传递函数模型一般表示为系统的特征方程为函数pzmap(num,den)or pzmap(a,b,c,d)函数功能：绘出零极点在s平面上的位置，以图形可视方式显示系统的稳定性。例：控制系统结构图如图所示，试分别确定例：控制系统结构图如图所示，试分别确定k=2k=2，k=10k=10时系统的稳定性。时系统的稳定性。nnc=2;ndc=conv(1 0,conv(1 1,1 2)nn0=1;nd0=1;nn,d=feedback(nc,dc,n0,d0)nro

48、ots(d)ans=ans=-2.5214 -2.5214 -0.2393+0.8579i -0.2393+0.8579i -0.2393-0.8579i -0.2393-0.8579i+-uy因此k=2时，由于系统的闭环特征根全部具有负实部，系统是稳定的。同理可求k=10时闭环特征根ans=-3.3089 0.1545+1.7316i 0.1545-1.7316i因此，当k=10时，由于系统有一对共轭复数根的实部为正值，系统不稳定。利用函数pzmap(num,den)也可判断稳定性。nnc=2;ndc=conv(1 0,conv(1 1,1 2)nn0=1;nd0=1;nn,d=feedback(nc,dc,n0,d0)npzmap(n,d)nnc=2;ndc=conv(1 0,conv(1 1,1 2)nn0=1;nd0=1;nn,d=feedback(nc,dc,n0,d0)nmargin(n,d)K2(注意幅值交角频率注意幅值交角频率w_c和相位交角频率和相位交角频率w_g,w_c w_g稳定稳定)nnc=10;ndc=conv(1 0,conv(1 1,1 2)nn0=1;nd0=1;nn,d=feedback(nc,dc,n0,d0)nmargin(n,d)K10