控制系统工具箱.ppt

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

1、第6章 控制系统工具箱,在MATLAB的Control System Toolbox(控制系统工具箱)中提供了许多仿真函数与模块，用于控制系统的仿真和分析。本章着重介绍控制系统的模型、时域分析方法和频域分析方法等内容。,对线性定常系统，式中s的系数均为常数，且a1不等于零，这时系统在MATLAB中可以方便地由分子和分母多项式系数构成的两个向量唯一地确定，这两个向量分别用num和den表示。num=b1,b2,bm,bm+1den=a1,a2,an,an+1,连续系统的传递函数模型（tf对象）,单输入单输出（SISO)LTI系统的传递函数：,注意：它们都是按s的降幂进行排列的！,MATLAB中，

2、用函数对象tf()来建立系统的传递函数模型。函数的调用格式为：,G=tf(num,den),例：系统传递函数为：,num=12,24,0,20;den=2 4 6 2 2;G=tf(num,den)Transfer function:12 s3+24 s2+20-2 s4+4 s3+6 s2+2 s+2,零极点模型实际上是传递函数模型的另一种表现形式，其原理是分别对原系统传递函数的分子、分母进行分解因式处理，以获得系统的零点和极点的表示形式。,在MATLAB中零极点增益模型用z,p,K矢量组表示：z=z1,z2,zm；p=p1,p2,.,pn；KGain=k；,零极点增益模型（zpk对象）,D

3、(s)为特征多项式，D(s)=0为特征方程，K为系统增益，zi为零点，pj为极点,零极点决定了系统的性能（系统的阶次和多项式表达时的系数）。,MATLAB中，用函数对象zpk()来建立系统的零极点增益模型。函数的调用格式为：,G1=zpk(z,p,k),例：系统的零极点模型为：,则该模型可以由以下语句输入到MATLAB的工作空间中：,z=-1.9294;-0.0353+0.9287i;-0.0353-0.9287i;p=-0.9567+1.2272i;-0.9567-1.2272i;0.0433+0.6412i;0.0433-0.6412i;kGain=6;G1=zpk(z,p,kGain)Z

4、ero/pole/gain:6(s+1.929)(s2+0.0706s+0.8637)-(s2-0.0866s+0.413)(s2+1.913s+2.421),可见，系统的零点和极点分别为分子、分母多项式的根，增益是传递函数分子和分母最高项系数的比值。,LTI系统可以用一组一阶微分方程来描述，其矩阵形式即为现代控制理论中常用的状态空间表示法。状态方程与输出方程的组合称为状态空间表达式，又称为动态方程，经典控制理论用传递函数将输入输出关系表达出来，而现代控制理论则用状态方程和输出方程来表达输入输出关系，揭示了系统内部状态对系统性能的影响。,在MATLAB中，系统状态空间用（A,B,C,D)矩阵组

5、表示。,x为状态向量；u为输入向量；y为输出向量A为系统矩阵；B为输入矩阵；C为输出矩阵；D为前馈矩阵；,状态空间模型（ss对象）*,然后使用函数对象ss()来建立系统的状态空间模型。函数的调用格式为：,G3=ss(A,B,C,D),MATLAB中只需将各个系数矩阵按常规矩阵的方式输入到工作空间中即可:,例：两输入两输出系统：,A=1 6 9 10;3 12 6 8;4 7 9 11;5 12 13 14;B=4 6;2 4;2 2;1 0;C=0 0 2 1;8 0 2 2;D=zeros(2,2);G3=ss(A,B,C,D),c=x1 x2 x3 x4 y1 0 0 2 1 y2 8 0

6、 2 2 d=u1 u2 y1 0 0 y2 0 0 Continuous-time model.,a=x1 x2 x3 x4 x1 1 6 9 10 x2 3 12 6 8 x3 4 7 9 11 x4 5 12 13 14 b=u1 u2 x1 4 6 x2 2 4 x3 2 2 x4 1 0,6.1.4 线性离散时间系统的数学模型,在MATLAB中，离散系统的数学模型和连续系统的数学模型表达函数相同，只是需要输入系统的采样周期T,具体格式如下。状态空间模型：sys=ss(A,B,C,D，Ts)零极点模型：sys=zpk(z,p,k,Ts)传递函数模型：sys=tf(num,den，Ts)

7、其中Ts为采样周期。,连续系统状态空间转换为离散系统状态空间形式采用函数c2d，其调用格式为sysd=c2d(sys,Ts)或Ad,Bd,Cd,Dd=c2dm(A,B,C,D,Ts,method)%将带选项的连续系统%转换到离散系统Ad,Bd,Cd,Dd=c2dt(A,B,C,T,lambda)%带有输入纯延迟的连续形式转换到离散形式,在一些场合下需要用到某种模型，而在另外一些场合下可能需要另外的模型，这就需要进行模型的转换。,模型转换函数包括：tf2zp：传递函数模型转换为零极点增益模型；tf2ss：传递函数模型转换为状态空间模型；zp2tf：零极点增益模型转换为传递函数模型；zp2ss：零

8、极点增益模型转换为状态空间模型；,6.1.5 模型的转换,函数tf2zp()可以用来求传递函数的零极点和增益。,ss2tf：状态空间模型转换为传递函数模型ss2zp：状态空间模型转换为零极点增益模型,例：已知某系统的传递函数零：,求：同一系统所对应的极点增益模型。,num=1,11,30,0;den=1,9,45,87,50;z,p,k=tf2zp(num,den),z=0-6-5,p=-3.0000+4.0000i-2.0000-1.0000,k=1,即，结果表达式为：,已知一个单输入三输出系统的传递函数模型为：num=0 0-2;0-1-5;1 2 0;den=1 6 11 6;A,B,C

9、,D=tf2ss(num,den)A=-6-11-6 B=1 C=0 0-2 D=0 1 0 0 0 0-1-5 0 0 1 0 0 1 2 0 0,系统的零极点增益模型：,注意：零极点的输入可以写成行向量，也可以写成列向量。,z=-3;p=-1,-2,-5;k=6;num,den=zp2tf(z,p,k)num=0 0 6 18 den=1 8 17 10a,b,c,d=zp2ss(z,p,k)a=-1.0000 0 0 b=1 2.0000-7.0000-3.1623 1 0 3.1623 0 0 c=0 0 1.8974 d=0,控制系统常用到并联系统，这时就要对系统函数进行分解，使其表

10、现为一些基本控制单元（子传递函数）的和的形式。同时部分分式展开也可以用来对系统求解。,6.1.6 部分分式展开,MATLAB中，Residue()函数可以实现传递函数模型与部分分式模型的互换。,函数r,p,k=residue(num,den)对两个多项式的比进行部分分式展开，即把传函分解为微分单元的形式。b,a=residue(r,p,k)可以将部分分式转化为多项式比p(s)/q(s)。,其中，向量num和den是按s的降幂排列的多项式系数向量（即传递函数的分子和分母）。R、P分别为各个子传递函数的增益和极点，K为部分分式展开后的余项（常数项）。,部分分式展开：num=2,0,9,1;den=

11、1,1,4,4;r,p,k=residue(num,den),p=0.0000+2.0000i-1.0000,k=2,r=0.0000+0.2500i-2.0000,结果表达式：,其解可以表示为：,例：对传递函数模型的串联：num,den=series(num1,den1,num2,den2)将串联连接的传递函数进行相乘。,可以使用series()函数直接对串联环节进行化简：,G=series(sys1,sys2),格式：,例：对状态空间模型的串联：a,b,c,d=series(a1,b1,c1,d1,a2,b2,c2,d2)串联连接两个状态空间系统。,6.1.7 模型的连接,例：对传递函数模

12、型的并联：num,den=parallel(num1,den1,num2,den2),可以使用parallel()函数直接对并联环节进行化简：,G=parallel(G1,G2),格式：,例：并联连接两个状态空间：a,b,c,d=parallel(a1,b1,c1,d1,a2,b2,c2,d2),更简单地：,G=G1G2,MATLAB中使用Gfeedback(G1,G2,sign)直接对反馈系统进行化简，sign用来指示系统2输出到系统1输入的连接符号(反馈极性)，sign省略时，默认为负，即sign=-1，sign=1表示正反馈。,1、传递函数的形式表示的反馈系统：num,den=feedb

13、ack(num1,den1,num2,den2,sign)2、两个状态函数表示的系统按反馈方式连接：a,b,c,d=feedback(a1,b1,c1,d1,a2,b2,c2,d2)a,b,c,d=feedback(a1,b1,c1,d1,a2,b2,c2,d2,sign)系统1的所有输出连接到系统2的输入，系统2的所有输出连接到系统1的输入总系统的输入/输出数等同于系统1。a,b,c,d=feedback(a1,b1,c1,d1,a2,b2,c2,d2,inp1,out1)部分反馈连接，将系统1的指定输出out1连接到系统2的输入，系统2的输出连接到系统1的指定输入inp1，以此构成 闭环系

14、统。,MATLAB中使用Gccloop(G,sign)直接对单位反馈系统进行化简，sign用来指示系统的反馈极性，sign省略时，默认为负，即sign=-1。,1、由传递函数表示的开环系统构成闭环系统：numc,denc=cloop(num,den,sign)当sign=1时采用正反馈；当sign=-1时采用负反馈；sign缺省时，默认为负反馈。2、闭环系统的状态空间模型ac,bc,cc,dc=cloop(a,b,c,d,sign)通过将所有的输出反馈到输入，从而产生闭环系统的状态空间模型。ac,bc,cc,dc=cloop(a,b,c,d,outputs,inputs)表示将指定的输出out

15、puts反馈到指定的输入inputs，以此构成闭环系统。,应用举例：1）系统1为：系统2为：求按串联、并联、正反馈、负反馈连接时的系统状态方程及系统1按单位负反馈连接时的状态方程。,2）系统1、系统2方程如下所示。,求部分并联后的状态空间，要求u11与u22连接，u13与u23连接，y11与y21连接。,ctrb和obsv函数可以求出状态空间系统的可控性和可观性矩阵。格式：co=ctrb(a,b)ob=obsv(a,c)对于nn矩阵a，nm矩阵b和pn矩阵cctrb(a,b)可以得到nnm的可控性矩阵co=b ab a2b an-1bobsv(a,c)可以得到nmn的可观性矩阵ob=c ca

16、ca2 can-1当co的秩为n时，系统可控；当ob的秩为n时，系统可观。,6.1.8 模型的属性,例6.12 已知某系统的闭环传递函数为：,要求判断系统的稳定性及系统是否为最小相位系统。,7 s4-32 s3-684 s2+2689 s+10365-s4-5 s3-86 s2+345 s+1236,%exp4_1.mclear;clc;close all;%系统描述n=7,-32,-684,2689,10365;d=1,-5,-86,345,1236;%求系统的零极点z,p,k=tf2zp(n,d)%检验零点的实部；求取零点实部大于零的个数ii=find(real(z)0)n1=length

17、(ii);%检验极点的实部；求取极点实部大于零的个数jj=find(real(p)0)n2=length(jj);%判断系统是否稳定,if(n20)disp(the system is unstable)disp(the unstable pole are:)disp(p(jj)%显示不稳定极点p(jj)else disp(the system is stable)end%判断系统是否为最小相位系统if(n10)disp(the system is a nonminimal phase one)else disp(the syetem is a minimal phase one)end%绘制

18、零极点图pzmap(p,z),源程序,ii=find(条件式)：用来求取满足条件的向量的下标向量，以列向量表示。,例如：条件式为real(p0)，其含义就是找出极点向量p中满足实部的值大于0的所有元素下标，并将结果返回到ii向量中去。这样如果找到了实部大于0的极点，则会将该极点的序号返回到ii下。如果最终的结果里ii的元素个数大于0，则认为找到了不稳定极点，因而给出系统不稳定的提示，反之得出系统稳定的结论。,pzmap(p,z)：根据系统已知的零极点p和z绘制出系统的零极点分布图,说明：,z=-8.9995-2.5260 8.0485+0.5487i 8.0485-0.5487ip=-8.26

19、56-2.4242 7.8449+0.3756i 7.8449-0.3756ik=7,ii=3 4jj=3 4the system is unstablethe unstable pole are:7.8449+0.3756i 7.8449-0.3756ithe system is a nonminimal phase one,运行窗口,例6.13 系统闭环模型如下所示：,判断系统的稳定性，以及系统是否为最小相位系统。,clearclcclose all%系统描述num=3 16 41 28;den=1 14 110 528 1494 2117 112;%求系统的零极点z,p,k=tf2zp(

20、num,den)%检验零点的实部；求取零点实部大于零的个数ii=find(real(z)0)n1=length(ii);%检验极点的实部；求取极点实部大于零的个数jj=find(real(p)0)n2=length(jj);%判断系统是否稳定,if(n20)disp(the system is unstable)disp(the unstable pole are:)disp(p(jj)else disp(the system is stable)end%判断系统是否为最小相位系统if(n10)disp(the system is a nonminimal phase one)else dis

21、p(the system is a minimal phase one)end%绘制零极点图pzmap(p,z),源程序,z=-2.1667+2.1538i-2.1667-2.1538i-1.0000 p=-1.9474+5.0282i-1.9474-5.0282i-4.2998-2.8752+2.8324i-2.8752-2.8324i-0.0550,k=3ii=jj=the system is stablethe system is a minimal phase one,运行窗口,1、step()函数的用法,y=step(num,den,t)：其中num和den分别为系统传递函数描述中的

22、分子和分母多项式系数，t为选定的仿真时间向量，可以由t=0:step:end产生。该函数返回值y为系统在仿真时刻各个输出所组成的矩阵。,y,x,t=step(A,B,C,D,iu)：其中A,B,C,D为系统的状态空间描述矩阵，iu用来指明输入变量的序号。x为系统返回的状态轨迹。,y,x,t=step(num,den)：此时时间向量t由系统模型的特性自动生成，对传递函数模型，状态变量x返回为空矩阵。,在MATLAB的控制系统工具箱中提供了求取单位阶跃响应的函数step()和冲激响应的函数impulse()。,6.2 控制系统的时域分析,如果对具体的响应值不感兴趣，而只想绘制系统的阶跃响应曲线，可

23、调用以下的格式，自动输出响应曲线：step(num,den)；step(num,den,t)；step(A,B,C,D,iu,t)；step(A,B,C,D,iu)；,线性系统的稳态值可以通过函数dcgain()来求取，其调用格式为：dc=dcgain(num,den)或dc=dcgain(a,b,c,d),clc；clear；close all%开环传递函数描述num=20;den=1 8 36 40 0;%求闭环传递函数numc,denc=cloop(num,den);%绘制闭环系统的阶跃响应曲线t=0:0.1:10;y=step(numc,denc,t);y1,x,t1=step(num

24、c,denc);%对于传递函数调用，状态变量x返回为空矩阵plot(t,y,r:,t1,y1)title(the step responce)xlabel(time-sec)%求稳态值disp(系统稳态值dc为：)dc=dcgain(numc,denc),源程序,运行窗口,系统稳态值dc为：dc=1,2、impulse()函数的用法,求取脉冲激励响应的调用方法与step()函数基本一致。y=impulse(num,den,t)；y,x,t=impulse(num,den)；y,x,t=impulse(A,B,C,D,iu,t)impulse(num,den)；impulse(num,den,t

25、)impulse(A,B,C,D,iu)；impulse(A,B,C,D,iu,t),clc;clear;%系统模型建立num=3;den=conv(1 1+3i,1 1-3i);%求系统的单位阶跃响应y,x,t=step(num,den);%求响应的稳态值finalvalue=dcgain(num,den)%求响应的峰值及对应的下标yss,n=max(y);%计算超调量及峰值时间percentovershoot=100*(yss-finalvalue)/finalvaluetimetopeak=t(n),例:已知系统G(s)=3/(s+1+3i)(s+1-3i)计算系统在阶跃输入时的峰值时间

26、、超调量、调整时间等系统瞬态性能指标（稳态误差允许正负2）,%计算上升时间n=1;while y(n)0.98*finalvalue)endsettlingtime=t(k),运行结果,finalvalue=0.3000percentovershoot=35.0914timetopeak=1.0491risetime=0.4417settlingtime=3.5337,6.2.2 控制系统的根轨迹分析方法,pzmap：绘制线性系统的零极点图;rlocus：求系统根轨迹。rlocfind：计算给定一组根的根轨迹增益。sgrid：在连续系统根轨迹图和零极点图中绘制出阻尼系数和自然频率栅格。,根轨迹

27、图绘制：rlocus(),rlocus(num,den)或rlocus(a,b,c,d)：根据SISO开环系统的传递函数模型和状态空间描述模型，直接在屏幕上绘制出系统的根轨迹图。开环增益的值从零到无穷大变化。rlocus(a,b,c,d,k)或rlocus(num,den,k)：通过指定开环增益k的变化范围来绘制系统的根轨迹图。r=rlocus(num,den,k)或者r,k=rlocus(num,den)：不在屏幕上直接绘出系统的根轨迹图，而根据开环增益变化矢量k，返回闭环系统特征方程1k*num(s)/den(s)=0的根r，它有length(k)行，length(den)-1列，每行对应

28、某个k值时的所有闭环极点。或者同时返回k与r。对于参数根轨迹，可以通过传递函数的等效变换而进行绘制。若给出传递函数描述系统的分子项num为负，则利用rlocus函数绘制的是系统的零度根轨迹（正反馈系统或非最小相位系统）。,rlocfind()函数,rlocfind()：找出给定的一组根（闭环极点）对应的根轨迹增益。,k,p=rlocfind(a,b,c,d)或者k,p=rlocfind(num,den)它要求在屏幕上先已经绘制好有关的根轨迹图。然后，此命令将产生一个光标以用来选择希望的闭环极点。命令执行结果：k为对应选择点处根轨迹开环增益；p为对应K的系统闭环特征根。不带输出参数项k,p时，同

29、样可以执行，只是此时只将k的值返回到缺省变量ans中。,源程序,%根轨迹图的绘制clcclearclose all%已知系统开环传递函数模型num=1;den=conv(0.01 1 0,0.02 1);rlocus(num,den)k1,p=rlocfind(num,den)k2,p=rlocfind(num,den)title(root locus),Select a point in the graphics windowselected_point=-21.1302-0.8630ik1=9.6354p=1.0e+002*-1.0774-0.2113+0.0086i-0.2113-0.0

30、086iSelect a point in the graphics windowselected_point=0.1023+68.1800ik2=139.7356p=1.0e+002*-1.4809-0.0095+0.6868i-0.0095-0.6868i,说明：重根即两条根轨迹曲线汇合处，纯虚根即根轨迹与虚轴的交点。,sgrid()函数（等阻尼比等振动频率栅格线）,sgrid：在现存的屏幕根轨迹或零极点图上绘制出自然振荡频率wn、阻尼比矢量z对应的格线。sgrid(new)：是先清屏，再画格线。sgrid(z,wn)：则绘制由用户指定的阻尼比矢量z、自然振荡频率wn的格线。,由控制理论知

31、，离虚轴近的稳定极点(主导极点）对整个系统的响应贡献最大，且在系统处于临界阻尼状态时，相角为45o，即根轨迹图上，应当选择相角为45o的射线附近的K值。通过rlocfind，配合sgrid()函数可以找出合适的增益值。,%系统开环传递函数为G(s)=k/s(s+1)(s+2)%试寻找一个合适的k值使得闭环系统具%有较理想的阶跃响应。%开环系统描述clc;clear;close allnum=1;den=conv(1 0,conv(1 1,1 2);z=0.1:0.2:1;wn=1:6;sgrid(z,wn);text(-0.3,2.4,z=0.1)text(-0.8,2.4,z=0.3)tex

32、t(-1.2,2.1,z=0.5)text(-1.8,1.8,z=0.7)text(-2.2,0.9,z=0.9),例6.21,%通过sgrid指令可以绘出指定阻尼比z和%自然振荡频率wn的栅格线hold onrlocus(num,den)axis(-4 1-4 4)k,p=rlocfind(num,den)%由控制理论知，离虚轴近的稳定极点对整个系统的响应贡献大%通过rlocfind，配合前面所画的z及wn%栅格线%从而可以找出能产生主导极点阻尼比%z=0.707的合适增益numc,denc=cloop(k,den);figure(2)step(numc,denc),运行窗口,Select

33、a point in the graphics windowselected_point=-0.4988+0.4348ik=0.6862p=-2.2454-0.3773+0.4040i-0.3773-0.4040i,6.3 控制系统的频域分析,MATLAB提供了函数bode()来绘制系统的波特图，其用法如下：,bode(num,den)：可绘制出以连续时间多项式传递函数表示的系统的波特图。bode(a,b,c,d,iu,w)或bode(num,den,w)：可利用指定的角频率矢量绘制出系统的波特图。由于横坐标按对数分度，因此w必须由logsapce（对数生成语句）生成。,6.3.1 函数bod

34、e,bode(a,b,c,d,iu)：可得到从系统第iu个输入到所有输出的波特图。当带输出变量mag,pha,w或mag,pha引用函数时，可得到系统波特图相应的幅值mag、相角pha及角频率点w矢量或只是返回幅值与相角。相角以度为单位，幅值可转换为分贝单位：magdb=20log10(mag),bode(a,b,c,d)：自动绘制出系统的一组Bode图，它们是针对连续状态空间系统a,b,c,d的每个输入的Bode图。其中频率范围由函数自动选取，而且在响应快速变化的位置会自动采用更多取样点。,MATLAB提供了函数nyquist()来绘制系统的极坐标图，其用法如下：,nyquist(num,d

35、en)：可绘制出以连续时间多项式传递函数表示的系统的极坐标图。nyquist(num,den,)或nyquist(a,b,c,d,iu,)：可利用指定的角频率矢量绘制出系统的极坐标图。,6.3.2 奈奎斯特图（幅相频率特性图）,nyquist(a,b,c,d,iu)：可得到从系统第iu个输入到所有输出的极坐标图。当不带返回参数时，直接在屏幕上绘制出系统的极坐标图（图上用箭头表示的变化方向，负无穷到正无穷）。当带输出变量re,im,引用函数时，可得到系统频率特性函数的实部re和虚部im及角频率点w矢量（为正的部分）。可以用plot(re,im)绘制出对应w从负无穷到零变化的部分。,nyquist

36、(a,b,c,d)：绘制出系统的一组Nyquist曲线，每条曲线相应于连续状态空间系统a,b,c,d的输入/输出组合对。其中频率范围由函数自动选取，而且在响应快速变化的位置会自动采用更多取样点。,（1）源程序,%clearclose allk=26;z=;p=-6 1;num,den=zp2tf(z,p,k);figure(1)subplot(211)nyquist(num,den)subplot(212)pzmap(p,z)figure(2)numc,denc=cloop(num,den);step(numc,denc),逆时针绕行圈数为1，不稳的极点数，所以系统稳定。,从闭环系统的阶跃响应

37、看系统的却是稳定的。,6.3.3 其他常用频域分析函数,MATLAB除了提供前面介绍的基本频域分析函数外，还提供了大量在工程实际中广泛应用的库函数，由这些函数可以求得系统的各种频率响应曲线和 特征值。如：,margin：求幅值裕度和相角裕度及对应的转折频率freqs：模拟滤波器特性nichols：求连续系统的尼科尔斯频率响应曲线（即对数幅相曲线）ngrid：尼科尔斯方格图,margin()函数,margin函数可以从频率响应数据中计算出幅值裕度、相角裕度以及对应的频率。幅值裕度和相角裕度是针对开环SISO系统而言，它指示出系统闭环时的相对稳定性。当不带输出变量引用时，margin(sys)可在

38、当前图形窗口中绘制出带有裕量及相应频率显示的Bode图，其中幅值裕度以分贝为单位。sys为系统模型描述,幅值裕度是在相角为-180度处使开环增益为1的增益量，如在-180度相频处的开环增益为g，则幅值裕度为1/g；若用分贝值表示幅值裕度，则等于：-20*log10(g)。类似地，相角裕度是当开环增益为1.0时，相应的相角与180度角的和。,margin(mag,phase,w)：由bode指令得到的幅值mag（不是以dB为单位）、相角phase及角频率w矢量绘制出带有裕量及相应频率显示的bode图。margin(num,den)：可计算出连续系统传递函数表示的幅值裕度和相角裕度并绘制相应波特图

39、。类似，margin(a,b,c,d)可以计算出连续状态空间系统表示的幅值裕度和相角裕度并绘制相应波特图。gm,pm,wcg,wcp=margin(mag,phase,w)：由幅值mag（不是以dB为单位）、相角phase及角频率w矢量计算出系统幅值裕度和相角裕度及相应的相角交界频率wcg、截止频率wcp，而不直接绘出Bode图曲线。,源程序,例6.24,%margin绘制出幅值裕度和相角裕度的波特图。clearclcclose all%某系统的开环传递函数为：G(s)=k/s(s+1)(0.2s+1)%求k分别为2和20时的幅值裕度与相角裕度num1=2;num2=20;den=conv(1 0,conv(1 1,0.2 1);figure(1)margin(num1,den)figure(2)margin(num2,den);,