北理工自动控制实验报告

1、北京理工大爭北京理工大爭Seeing砒 Tevfmology本科实验报告实验名称：控制理论基础实验课程名称：实验时间：任课教师：实验地点：实验教师：实验类型：原理验证综合设计自主创新学生姓名：学号班级组号：学院：同组搭档：专业：成绩：实验 1 控制系统的模型建立一、实验目的1. 掌握利用 MATLAB 建立控制系统模型的方法。2. 掌握系统的各种模型表述及相互之间的转换关系。3. 学习和掌握系统模型连接的等效变换。二、实验原理1. 系统模型的 MATLAB 描述 系统的模型描述了系统的输入、输出变量以及内部各变量之间的关系，表征一个系统的 模型有很多种，如微分方程、传递函数模型、状态空间模型等

2、。这里主要介绍系统传递函数(TF)模型、零极点增益(ZPK)模型和状态空间(SS)模型的MATLAB描述方法。传递函数(TF)模型 传递函数是描述线性定常系统输入-输出关系的一种最常用的数学模型，其表达式一般为b sm + bsml + b si + b厂(u) = mmi10-anSn + ani Sni + + 1$1 + %在 MATLAB 中，直接使用分子分母多项式的行向量表示系统，即 num = bm, bm-1, b1, bden = an, an-1, a1, a0调用廿函数可以建立传递函数TF对象模型，调用格式如下：Gtf = tf(num,den) Tfdata函数可以从TF

3、对象模型中提取分子分母多项式，调用格式如下：num,den = tfdata(Gtf) 返回 cell 类型的分子分母多项式系数num,den = tfdata(Gtf,v) 返回向量形式的分子分母多项式系数零极点增益(ZPK)模型传递函数因式分解后可以写成g = * Zi)(szj(s Pi)(S 卩 2)(S 霭)式中，Zi,Z2 -Zm称为传递函数的零点，卩，卩2卩n称为传递函的极点，k为传递 系数(系统增益)。在MATLAB中，直接用z,p,k矢量组表示系统，其中z, p，k分别表示系统的零极点及 其增益，即：Z= Zi，Z2Zm；P=Pi，P2卩；k=k；调用zpk函数可以创建ZPK

4、对象模型，调用格式如下：G = zpk(z, p, k)同样，MATLAB提供了 zpkdata命令用来提取系统的零极点及其增益，调用格式如下：z, p, k = zpkdata(Gzpk)返回cell类型的零极点及增益z, p, k = zpkdata(Gzpk,v)返回向量形式的零极点及增益函数 pzmap 可用于求取系统的零极点或绘制系统得零极点图，调用格式如下： pzmap(G)在复平面内绘出系统模型的零极点图。p,z = pzmap(G)返回的系统零极点，不作图。状态空间(SS )模型 由状态变量描述的系统模型称为状态空间模型，由状态方程和输出方程组成：rx = Ax + buy =

5、 Ax + Du其中：x为n维状态向量；u为r维输入向量；y为m维输出向量；A为nXn方 阵，称为系统矩阵；B为nXr矩阵，称为输入矩阵或控制矩阵；C为mXn矩阵，称为 输出矩阵；D为mXr矩阵，称为直接传输矩阵。在MATLAB中，直接用矩阵组A,B,C,D表示系统，调用SS函数可以创建ZPK对象模型， 调用格式如下：GSS = SS(A,B,C,D)同样， MATLAB 提供了 SSdata 命令用来提取系统的 A、B、C、D 矩阵，调用格式如下： A,B,C,D = ssdata (Gss)。它返回系统模型的A、B、C、D矩阵。三种模型之间的转换 上述三种模型之间可以互相转换， MATLA

6、B 实现方法如下TF 模型f ZPK 模型：zpk(SYS)或廿2zp(num,den)TF 模型SS 模型：ss(SYS)或 tf2ss(num,den)ZPK 模型TF 模型：tf(SYS)或 zp2tf(z,p,k)ZPK 模型SS 模型：ss(SYS)或 zp2ss(z,p,k)SS 模型TF 模型：tf(SYS)或 ss2tf(A,B,C,D)SS 模型ZPK 模型：zpk(SYS)或 ss2zp(A,B,C,D)系统模型的连接 在实际应用中，整个控制系统是由多个单一的模型组合而成，基本的组合方式有串联连 接、并联连接和反馈连接。图 1-2 分别为串联连接、并联连接和反馈连接的结构框

7、图和等 效总传递函数。在 MATLAB 中可以直接使用“*”运算符实现串联连接，使用“”运算符实现并联连 接。反馈系统传递函数求解可以通过命令feedback实现，调用格式如下：T = feedback(G,H)T = feedback(G,H,sign)其中，G为前向传递函数，H为反馈传递函数；当sign = +1时，GH为正反馈系统传 递函数；当sign = -1时，GH为负反馈系统传递函数；默认值是负反馈系统。三、实验内容1.已知控制系统的传递函数如下G (s) 2s2+18s+40s3+5s2+8s+6试用 MATLAB 建立系统的传递函数模型、零极点增益模型及系统的状态空间方程模型，

8、并绘制系统零极点图。 代码： num=2 18 40; den=1 5 8 6; Gtf=tf(num,den); Gzpk=zpk(Gtf); pzmap(Gzpk);G-ss = a G-ss = a - Klz2b =xlsi TOC o 1-5 h z -110-1-12.51500-3uly21.778x31.414Klx3yl3.3741. 1250d =ulyl 0Coni: inuoust iute state-space model.零极点图：Real Axis (seconds )sdnoces(sixAyranigam2.已知控制系统的状态空间方程如下0 1 Real A

9、xis (seconds )sdnoces(sixAyranigam2.已知控制系统的状态空间方程如下0 1 00 & 1 a o o-1 -2 _3试用 MATLAB 建立系统的传递函数模型、零极点增益模型及系统的状态空间方程模型 并绘制系统零极点图。代码：A=0 1 0 0;0 0 1 0;0 0 0 1;-1 -2 -3 -4;B=0;0;0;1;C=10 2 0 0;D=0;Gss=ss(A,B,C,D);Gtf=tf(Gss);Gzpk=zpk(Gss);pzmap(Gzpk);grid on;实验结果：Gss =a =k2a701a0k200107300a1-1-2-3-4b =1

10、11X0k20 x30k41c =k2k3护102Q0d =111y0Gont innous-timE st at bspace model GtfGtf =2 s + 10s4 + 4 s3 + 3 s2 + 2 s + 1Continuciiis-tzLuie! transfer functzLon. GzpltGzpk =2 (s+5?Cs+3.234) Cs+0.6724)佔2 + 0. 0936s + 0. 45995Continucus-tim.e sero/pole/gairi Jiiodel.零极点图：Role-Zero1sdnoces(sixAyranigamI8 6 4 2

11、 co00000.9840.9930.9640.935J_0.9994r-3-0.870.740.452 4 60000.999_0.9930.9350.870.740.45r 0.984 rr0.964-0.8-5-4.5-4-3.5-3-2.5-2-1.5-1-0.50Real Axis (seconds-1)3 已知三个系统的传递函数分别为5 -1-45 Sj +2乳s + 3）（47）（5十D（$十叫）（十g）试用 MATLAB 求上述三个系统串联后的总传递函数。 代码：num1=2 6 5;den1=1 4 5 2;G1=tf(num1,den1);num2=1 4 1;den2=1

12、 9 8 0;G2=tf(num2,den2);num3=conv(conv(1 3,1 7),0 5);den3=conv(conv(1 1,1 4),1 6);G3=tf(num3,den3);G=G1*G2*G3; 实验结果：22十6 s十5s3 + 4+ 5 s + 2CantinuouE-t:LnLE transfer function. G2G2 =s2 + 4 s + 1s3 + 9 s2 + 8 aCont muous-11110 t ransfer function.) 03G3 =E 32 + 50 s + 105s3 斗 IL + 34 s + 24C outinuous

13、-t lute t ransfer funct ion. GG =106 + 1705 + 1065 sn4 + 3150 s3 + 45SO2 + 2930 s + 525sA9 十 24 sA8 十 226 s7 十s5 十 2905 s5 十 4516 s4 十 4044 s5 十 19対点2 十 384 sCoutinuoiistiiKB transf&r function.4.已知如图E2-1所示的系统框图试用 MATLAB 求该系统的闭环传递函数。 代码：num1=1; den1=1 1;G1=tf(num1,den1); num2=1;den2=0.5 1;G2=tf(num2,

14、den2); num3=3;den3=1 0; G3=tf(num3,den3);G=(G1+G2)*G3; T=feedback(G,G2);U=G/(1+G2*G);实验结果： Tr =2. 25 32 + 7. 5 s + 50. 25 4 + 1.25 s3 + 2 52 + 5. 5 s + 6Coritinuous-time transfer functioii. UU =1. 125 s5 十 7. 125 s i 十 1弓事3 十 lfi.5 e2 + fi s0. 1257 + s6 + 3. 125+ 7 s4 + 13. 25 s3 + 14.5 s2 + 6 sCant

15、inuous-time transfer function.5已知如图 E2-2 所示的系统框图试用 MATLAB 求该系统的闭环传递函数。 代码： num1=10;den1=1 1;G1=tf(num1,den1);num2=2; den2=conv(1 0,1 1);G2=tf(num2,den2);num3=1 3; den3=1 2;H1=tf(num3,den3);num4=5 0; den4=1 6 8;H2=tf(num4,den4); T1=feedback(G2,H1,1);G=G1*T1; T=feedback(G,H2);实验结果： T20+ 160 吕2 + 400 s

16、 + 3206 + 1Q5 + 35 吕乜 + 443 + 32 s2 + 116 a -慕Continiious一time transfer functiari.四实验收获与心得经过上学期的MATLAB实验，对MATLAB软件已经有一定的了解了。这次比较简单的几 个实验让我学到了怎么将传递函数通过MATLAB来表达，怎么画传递函数的零极点，怎么设 置软件比较容易看。倒腾了半天，最大的收获是知道了 conv函数是怎么用的，以及有三个 因式的该怎样书写。实验 2控制系统的暂态特性分析一、实验目的学习和掌握利用 MATLAB 进行系统时域响应求解和仿真的方法。 考察二阶系统的时间响应，研究二阶系统

17、参数对系统暂态特性的影响。二、实验原理系统的暂态性能指标控制系统的暂态性能指标常以一组时域量值的形式给出，这些指标通常由系统的单位 阶跃响应定义出来，这些指标分别为：延迟时间 td ：响应曲线首次到达稳态值的 50%所需的时间。上升时间 tr ：响应曲线从稳态值的 10%上升到 90%所需要的时间长，对于欠阻尼 系统，通常指响应曲线首次到达稳态值所需的时间。峰值时间 tp ：响应曲线第一次到达最大值的时间。调整时间 ts ：响应曲线开始进入并保持在允许的误差(2%或5%)范围内所 需要的时间。超调量。：响应曲线的最大值和稳态值之差，通常用百分比表示其中y(t)为响应曲线。在MATLAB中求取单

18、位阶跃响应的函数为step，其使用方法如下 step(sys)在默认的时间范围内绘出系统响应的时域波形 step(sys,T)绘出系统在0- T范围内响应的时域波形step(sys,ts:tp:te)绘出系统ts - te范围内，以tp为时间间隔取样的响应波形。y,t=step()该调用格式不绘出响应波形，而是返回响应的数值向量及其对应的时间向量。系统的暂态性能指标可以根据上述定义，在响应曲线上用鼠标读取关键点或通过搜索曲 线对应的数值向量中关键点来确定。2.LTIViewer 工具在MATLAB中提供了线性是不变系统仿真的工具LT Viewer，可以方便地观察系统的响应 曲线和性能指标。在命

19、令窗口中键入litview即可启动LTIViewer。这里简要介绍LTI Viewer工 具。【File】菜单Import选项：可以从Workspace或文件中导入系统模型。Export选项：将当前窗口中的对象模型保存到Workspace或文件中。Toolbox preferences选项：属性设置功能，可以设置控制系统中得各种属性值。Page Setup选项：页面设置功能，可以对打印输出和显示页面进行设置。【Edit】菜单Plot Configuration 选项：对显示窗口及显示内容进行配置。Line Style 选项：线型设置功能，可以对输出响应曲线的线型进行设置。Viewer Pref

20、erences 选项：对当前窗口的坐标、颜色、字体、响应曲线的特性参数等属 性进行设置。3 ）右键菜单 在运行界面上点击鼠标右键，将会弹出一个弹出式菜单，菜单上个选项的功能分别为：Plot Types：选择绘制的系统曲线类型，可选的类型有单位阶跃响应、单位冲击响应、波特 图、奈奎斯特图、零极点图等。System：选择需要仿真的系统。Characteristic：系统的性能指标选项。Grid ：显示和关闭网格。Normalize :正常显示模式。Full View :满界面显示模式。Properties :性能编辑器选项，可以设置画面的标题、坐标标志、坐标范围、线型、颜 色、性能指标等。三、实验

21、内容1.已知单位负反馈系统前向通道的传递函数为试用MATLAB绘制系统的单位阶跃响应曲线。代码：num1=80; den1=1 2 0;G1=tf(num1,den1);T=feedback(G1,1);step(T);实验结果：Sep Resp onsg1.81.6 - .1.121.6 - .1.1234Time (seconds)5672. 已知二阶系统2. 已知二阶系统(1)Z=6.0,sn=5，试用MATLAB绘制系统单位阶跃响应曲线，并求取系统的暂态性 能指标。代码：Wn=5;E=0.6; num1=WnA2;den1=1 2*E*Wn num1;G=tf(num1,den1);s

22、tep(G);edutilpmA1.41.21Step ResponseSystem: GPeak amplitude: 1.09edutilpmA1.41.21Step ResponseSystem: GPeak amplitude: 1.09Overshoot (%): 9.48At time (seconds): 0.783System: GRise time (seconds): 0.371System: GSettling time (seconds): 1.19System: GFinal value: 1000.20.40.60.811.21.41.6Time (seconds)

23、由图中可以看到：上升时间tr:响应曲线上升到90%所需要的时间为0.371s。峰值时间tp:响应曲线第一次到达最大值的时间为0.787s。调整时间ts:响应曲线开始进入并保持在允许的误差范围内所需要的时间为1.19s。超调量。：响应曲线的最大值和稳态值只差的百分比为。=9%(2)n=1,Z从0变化到2,求此系统的单位阶跃响应。代码：Wri=l ：for i=fl;5E=0.a#i； numl=WnA2;denl=L 2*EWn num.1:G=tfdeni):si:ep(Gj 30);hold on;-end.Legend C E二CT ,J E=0. aJ / E=0. 8J / E= 1.

24、2J E=1. 6J/ E二2)实验结果：0 5 10 15 20 25 30Time (seconds)(3)Z=5.0 ,sn从0 5 10 15 20 25 30Time (seconds)(3)Z=5.0 ,sn从0变化到1 ( n工0)，求此系统的单位阶跃响应。 代码：E=0.5：for i=0;5Step Resp onse2864286420000edutilpmA02+i;riml=WrL2;denl=l 2*E*Wn nujnl;G=tf tnmnlj deni):st ep tGj 30):hold on;?ndlegend C Wn=0J/ Wn=0. 2J/Wn=a.

25、4J/Wn=0. 6J/Wn=0. 8J/Wn=LJ ):实验结果：Time (seconds)421861 0 04200Time (seconds)421861 0 042000edutilpmA（4）观察上述实验结果，分析这两个特征参数对系统暂态特性的影响。由上面的两个图可以看到，Z参数主要影响系统的稳定性，超调量，上升时间。Z越大， 上升时间，峰值时间变长，而超调量变小。 n参数主要影响时间，3n越大，上升时间，峰值时间越小，而对超调量等其他因素 近似没有影响。四、实验心得通过这个实验，对二阶系统的sn,Z对系统的影响有更深的理解，也知道怎么系统更 稳定，达到稳定的时间更短，更快，更好

26、，对以后的实际电路有深远的影响。实验 3根轨迹分析一、实验目的学习和掌握利用MATLAB绘制根轨迹图的方法。学习和掌握利用系统根轨迹图分析系统的性能。二、实验原理根轨迹分析的MATLAB实现根轨迹是指系统某一参数变化时，闭环特征根在s平面上运动的轨迹。在MATLAB中, 提供了用于根轨迹分析的专门函数。1 ）rlocus 函数该函数的使用方法如下：rlocus(sys)绘制单输入单输出LTI系统的根轨迹图。rlocus(sys,k)使用用户指定的根轨迹增益k来绘制系统的根轨迹图。r,k= rlocus(sys)返回根轨迹增益值和闭环极点值，不绘制根轨迹图。rlocfind 函数该函数的使用方法

27、如下：k,poles=rlocfind(sys)计算鼠标选取点处的根轨迹增益值和闭环极点值，可在图形窗口 根轨迹图中显示出十字光标，当用户选择其中一点时，相应的增益值和极点值记录在k和 poles 中。k,poles=rlocfind(sys,p)计算最靠近给定闭环极点p处的根轨迹增益。sgrid函数该函数的使用方法如下：Sgrid可在连续系统根轨迹或零极点图上绘制出栅格线，栅格线由等阻尼系数和等自然 频率线构成。sgrid( new先清除当前的图形，然后绘制出栅格线，并将坐标轴属性设置成hold on。 sgrid(z,Wn)指定阻尼系数z和自然频率Wn。sgrid(z,Wn, new )指

28、定阻尼系数z和自然频率Wn，在绘制栅格线之前清除当前的图 形并将坐标轴属性设置成hold on.Rltool 工具MATLAB提供了一个根轨迹设计工具Rltool，在命令窗口输入rltool命令即可启动该工 具，也可输入rltool (sys)命令打开带系统模型sys的根轨迹设计器，运行界面如图3 - 1所 示。Rltool工具既可以分析系统根轨迹也能对系统进行设计，具体使用方法请参考MATLAB 帮助或查阅相关资料。三、实验内容1已知系统开环传递函数为使用MATLAB绘制系统的根轨迹图。 代码：k=l ;num=1 5:den=conv (conv 1 1L 1 3)? 1 12):niuR

29、j den,.1 ;zpkl = spk (G):rlocus (spkl);实验结果:Root Locus3020100-10-20-10-6-4-202Real Axis (seconds-1)(2)求根轨迹的两条分支离开实轴时的K值，并确定该K值对应的所有闭环极点。代码及结果：-14 -12-8sdnoces(sixAyranigam-30 kj poles = rlncfind.(zpkl)Select a point in the graphics windowselected_point =-2.1213 -2.1213 -0.000713.3812poles -11.7575-2

30、.1213 十 0.00071-2.1213 - 0.0007130-20-14 -12-10 -8 -6 -4Real Axis (seconds-1)-2Root locus200200011)1sdnoces(sixAyraniga-30分析可得：k=3.3812; poles=-11.7575,-2.1213-0.0007i,-2.1213+0.0007i。重新绘制根轨迹图观察其(3)以区间-40,-5 之间的值替代12 - = s处的极点, 对根轨迹图的影响。重新绘制根轨迹图观察其代码：for i=5:5:40num=1 5;den=conv(conv(1 1,1 3),1 i);G

31、=tf(num,den);zpk1=zpk(G);rlocus(zpk1);hold on;end实验结果：-45 -40 -35 -30 -25 -20 -15 -10 Real Axis (seconds-1)-5 0 50 0 0 0 01sdnoces(sixAyranigamI从图中可以看到，轨迹的分支数并没有增加或减少，渐进线也并没变化，极值越大，轨 迹分离之后趋近于无穷大的时间越快，分离点没有变化。2.已知系统开环传递函数为(1)使用MATLAB绘制系统的根轨迹图。代码：k=l :nun= 1 8:den=conv (ronv( 1 2., 1 0 ) 1 S 32)G=tf (

32、nunij den):zpk1 = zpk (G):rlocus(zpk1);-25实验结果：-25 一一- /u 1 z-25 -20 -15-10 -5 0Real Axis (seconds-1)5 10 155 0 5 0 5 0 5 01 1 - 1 1 2 )1sdnoces(sixAyranigamI计算两条分支进入右半平面和两条分支复数极点出发在实轴相交处的K值。 代码及结果： lijpoles = rlocfindfzpk0Select a point in the graphics windowsalected_point =-0.0237 十 3.4938i54.7660

33、pole吕=-5.0253 + 3.354Si-5.0253 - 3.3548i0.0253 十 3.4641i0.0263 - 3.4641iL11 L1/-i” riir25L11 L1/-i” riir252015105Root locus0-5-10-15-20-25 -20 -15-10 -5 0Real Axis (seconds-1)s-dnoces(sixAyranigam5 10 15-25代码及结果： kjpcles =rlorfirid(zpkl)Select a point in the graphics windowselect ed_point =-10. 1659

34、 - 0.077612.0693e+03poles =5. 3414 +11.417515.2414 -11.4175i-10.2414 + 0.02571-10.2414 - 0.0257i252015/ 252015/ -1+- /L-十划- -rri15 10 15Root Locus0 5 0 5 0 5 01 - 1 1 2 )1s-dnoces(sixAyranigamI-25-25 -20 -15 -10 -5 0 Real Axis (seconds-1)以区间-20,-1 之间的值替代零点的位置，重新绘制根轨迹图，观察其对根轨 迹图的影响。代码：for i=0.25:5z=-

35、2i:p= 0 -2 -4+4j 44j；2pil=2pkm k):rlocus(zpkl);h.oId on;1- end结果：-30-25-20-15-10-5051015Real Axis (seconds-1)Root Locus00001 1 2 -30-25-20-15-10-5051015Real Axis (seconds-1)Root Locus00001 1 2 )1sdnoces(sixAyranigamI-30分析：零点绝对值越大，左半支汇合点越靠近左边，绝对值越大。3.已知单位负反馈系统的开环传递函数为(1)使用 MATLAB 绘制系统的根轨迹图。 代码：nuni=l

36、;den= 1 4 5 0;G=tfnumj den);zpkl=zpk(G):rlocus Czpk1):结果;. / rirri-Root Locus3210-1-2-3-6 -5-4 -3 -2 -1 Real Axis (seconds-1)01sdnoces(sixAyranigam-4(2)分析使系统稳定的K值范围和使系统无超调的K值范围，并通过观察系统的单位 阶跃响应曲线加以验证。-4代码：for k=0:d:20;num.=i:den.= 1 4 5 0:G-=t f (nunij den.):E=feedback CGj 1): step (g, 30): hold on;-

37、endlegend C k=(T J k=4J.J k=8J / k=12J / k= ir / k=Z(T):结果：Step Response21.510.552530k=0k=4k=8 k=12 k=16 k=20分析代码10 15 20 Time (seconds) 由上图可以看出，使得系统稳定的 K 值范围为： 0k20Step Response21.510.552530k=0k=4k=8 k=12 k=16 k=20分析代码for k=0:2:20：nuin= 0. 2; den=1 450;G=tf (rvuuij den); g=f eedback (G-j 1): step(E

38、j 30):hold on:endlegend(JI=0J; Jk=2J/k=4J/k=10J; Jk=6J/k=8J/k=12J; Jk=18 Jk=2OJ5 ;结果：Step Response1.4 |cc11.28600edutilpmA0.40.2k=08600edutilpmA0.40.2k=0一 k=2k=4k=10一 k=14k=6k=8一 k=12k=16k=18k=2000 5 10 15 20 25 30Time (seconds)分析：由上图可以看出，使得系统无超调的K值范围为：0k=2四、实验心得刚刚学的自控的根轨迹画法，由于用手画，开始觉得根轨迹其实不是那么的好用，但

39、是 通过这次试验，通过用MATLAB试验快速准确的画出根轨迹之后，发现可以很直观的看出使 系统具有各种各样的性质需要的K值范围，从而可以很快的找到符合要求的值。实验 4 系统的频率特性分析一、实验目的学习和掌握利用MATLAB绘制系统Nyquist图和Bode图的方法。学习和掌握利用系统的频率特性分析系统的性能。二、实验原理系统的频率特性是一种图解方法，分析运用系统的开环频率特性曲线，分析闭环系统的 性能，如系统的稳态性能、暂态性能常用的频率特性曲线有Nyquist图和Bode图。在MATLAB 中，提供了绘制 Nyquist 图和 Bode 图的专门函数。1Nyquist 图nyquist

40、函数可以用于计算或绘制连续时间 LTI 系统的 Nyquist 频率曲线，其使用方法 如 下： nyquist(sys) 绘制系统的 Nyquist 曲线。nyquist(sys,w)利用给定的频率向量w来绘制系统的Nyquist曲线。re,im=nyquist(sys,w)返回Nyquist曲线的实部re和虚部im,不绘图。Bode 图bode函数可以用于计算或绘制连续时间LTI系统的Bode图，其使用方法如下： bode(sys) 绘制系统的 Bode 图。bode(sys,w)利用给定的频率向量w来绘制系统的Bode图。mag,phase=bode(sys,w)返回Bode图数据的幅度mag和相位phase，不绘图。幅度和相位裕度计算margin 函数可以用于从频率响应数据中计算出幅度裕度、相位裕度及其对应的角频率， 其使用方法如下： margin(sys)margin(sys)margin(mag,phas