北理工自动控制实验报告

北理工自动控制实验报告本科实验报告实验名称： 控制理论基础实验信息与电子学院信息与电子学院课程名称：实验时间：任课教师：实验地点：实验教师：实验类型：□原理验证■综合设计□自主创新学生姓名：学号/班级：组 号：学院：同组搭档：专业：成 绩：实验1控制系统的模型建立一、 实验目的掌握利用MATLAB建立控制系统模型的方法。掌握系统的各种模型表述及相互之间的转换关系。学习和掌握系统模型连接的等效变换。二、 实验原理系统模型的MATLAB描述系统的模型描述了系统的输入、输出变量以及内部各变量之间的关系，表征一个系统的模型有很多种，如微分方程、传递函数模型、状态空间模型等。这里主要介绍系统传递函数（TF）模型、零极点增益（ZPK模型和状态空间（SS）模型的MATLAE描述方法。1）传递函数（TF）模型传递函数是描述线性定常系统输入-输出关系的一种最常用的数学模型，其表达式一般为?????+???1???-1+?+????+??????=—??????+???-1???-1+?+????+??在MATLAB中，直接使用分子分母多项式的行向量表示系统，即num=[bm,bm-1,…b1,b]den=[an,an-1,…a1,a0]调用tf函数可以建立传递函数TF对象模型，调用格式如下：Gtf=tf(num,den)Tfdata函数可以从TF对象模型中提取分子分母多项式，调用格式如下：[num,den]=tfdata(Gtf)返回cell类型的分子分母多项式系数返回向量形[num,den]=tfdata(Gtf,'v')返回向量形式的分子分母多项式系数零极点增益(ZPK模型传递函数因式分解后可以写成一“？??？？?)(？？??)・・・(?2???)????=(????)(????)・・・(??■???)式中，？?,？？…??？称为传递函数的零点，??,??…?金？称为传递函的极点，k为传递系数（系统增益）。在MATLAB中，直接用[z,p,k]矢量组表示系统，其中z,p,k分别表示系统的零极点及其增益，即：z=[??,??…???]；p=[??,??…???；k=[k];调用zpk函数可以创建ZPK对象模型，调用格式如下：??=??????????）同样，MATLAB提供了zpkdata命令用来提取系统的零极点及其增益，调用格式如下：[??????]=??????????????（返?回？？？@?）类型的零极点及增益[??????]=??????????????（??????返?回向量形式的零极点及增益函数pzmap可用于求取系统的零极点或绘制系统得零极点图，调用格式如下：??????????在复平面内绘出系统模型的零极点图。[????]=??????????（?回）的系统零极点，不作图。状态空间(SS模型由状态变量描述的系统模型称为状态空间模型，由状态方程和输出方程组成：???????????

{??=????????■■■■■■■■■■其中：X为n维状态向量；U为r维输入向量；y为m维输出向量；A为nXn方阵，称为系统矩阵；B为nXr矩阵，称为输入矩阵或控制矩阵；C为mXn矩阵，称为输出矩阵；D为mXr矩阵，称为直接传输矩阵。在MATLA沖，直接用矩阵组[A,B,C,D]表示系统，调用ss函数可以创建ZPK对象模型，调用格式如下：Gss=ss(A,B,C,D)同样，MATLAB提供了ssdata命令用来提取系统的AB、CD矩阵，调用格式如下：[A,B,C,D]=ssdata(Gss)。它返回系统模型的A、B、CD矩阵。三种模型之间的转换上述三种模型之间可以互相转换， MATLAB实现方法如下TF模型fZPK模型：zpk(SYS)或tf2zp(num,den)TF模型fSS模型：ss(SYS)或tf2ss(num,den)ZPK模型fTF模型：tf(SYS)或zp2tf(z,p,k)ZPK模型fSS模型：ss(SYS)或zp2ss(z,p,k)SS模型fTF模型：tf(SYS)或ss2tf(A,B,C,D)SS模型fZPK模型：zpk(SYS)或ss2zp(A,B,C,D)系统模型的连接在实际应用中，整个控制系统是由多个单一的模型组合而成，基本的组合方式有串联连接、并联连接和反馈连接。图1-2分别为串联连接、并联连接和反馈连接的结构框图和等效总传递函数。在MATLAB中可以直接使用“*”运算符实现串联连接，使用“+”运算符实现并联连接。反馈系统传递函数求解可以通过命令 feedback实现，调用格式如下：T=feedback(G,H)T=feedback(G,H,sign)其中，G为前向传递函数，H为反馈传递函数；当sign=+1时，GH为正反馈系统传递函数；当sign=-1时，GH为负反馈系统传递函数；默认值是负反馈系统。三、实验内容1.已知控制系统的传递函数如下????=2??+18??+40??+5?*+8??+6试用MATLAB!立系统的传递函数模型、零极点增益模型及系统的状态空间方程模型，并绘制系统零极点图。代码：num=[21840];den=[1586];Gtf=tf(num,den);Gzpk=zpk(Gtf);pzmap(Gzpk);

gridon;Gss=ss(Gzpk);实验结果：G-ss-a=xlx2x3xl-110玄2-1-12.515Q0卡b=ul龙10x2玄31.414c-xlx2x3yi3.374n250d=ulyi0Coniinuous-tim.*siat?-spacsitodel.零极点图：Pole-ZeroMap■51'5c5

ra0JrSanocesctpXAvra^aa0.9450.890.810.680.50.30.9760.994-0.9940.9760.9450.890.81Pole-ZeroMap■51'5c5

ra0JrSanocesctpXAvra^aa0.9450.890.810.680.50.30.9760.994-0.9940.9760.9450.890.810.680.5-1.5-5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5-1RealAxis(seconds )0.3E-1 -0.52.已知控制系统的状态空间方程如下'0100_0_00L00X=工+00010-l—2-3-+Iy=200卜试用MATLAB建立系统的传递函数模型、零极点增益模型及系统的状态空间方程模型，并绘制系统零极点图。代码：A=[0100;0010;0001;-1-2-3-4];B=[0；0；0；1]；C=[10200];D=0;Gss=ss(A,B,C,D);Gtf=tf(Gss);Gzpk=zpk(Gss);pzmap(Gzpk);gridon;实验结果：Gas二注-Kl 就xl010t)x20019x30001^4-1 ~2-3-4b-ulxl0x20x30xiIxlx2 就yl10 2 0 0d=Lilyl0»Gtf&tf=2s+104s^3+33^2 25+ICoatinuous-tinuttansf^rfwutipih»GzjiGzpk=+0l+0l4599)Cs+3.234)Cs+0.6721)(s2+0.0936sCentinuous-tim.ezero/pole/^ainm.odel-零极点图：<snocesccEXAVTa^aa卩'0.984‘0.964<snocesccEXAVTa^aa卩'0.984‘0.9640.9350.870.74 0.450.9930.999\\4321Pole-ZeroMap0.99925-0.9930.9930.984 0.964 0.935 0.87 0.74 0.45-5-3.5 -3 -2.5 -2-5-3.5 -3 -2.5 -2-1.5 -1 -0.5 0RealAxis(seconds-1)3已知三个系统的传递函数分别为G©Af(j+1Xj+4Xs+6)试用MATLAB求上述三个系统串联后的总传递函数。代码：num仁[265];den1=[1452];G仁tf(num1,den1);num2=[141];den2=[1980];G2=tf(num2,den2);num3=conv(conv([13],[17]),[05]);den3=conv(conv([11],[14]),[16]);G3=tf(num3,den3);G=G1*G2*G3;实验结果：G1二2s^24-6s+5s'3十4e*2+5e+2匚outuiuou£-time■tiansffrrfunction.»C2C2三sR2+4s+1J3+9s^2+S5Carrtiznuous—七亠血亡tran3£er■fimv七ion.>>G303=5s*2+50s+1OBs'3411s'2+34s+24Con±inuoiis-timetransfertuncticn.>>GC—10$6+170$5+1065f4+3150s3+-孤QT2+2980$+5253^9十24sAB十226s^7十10S4s"o十2905s^5■+451Gs\+4044s^3十19J6事2十J34sCoutinuous'timetransferfunction.4.已知如图E2-1所示的系统框图试用MATLAB求该系统的闭环传递函数。代码：num1=1;den仁[11];G仁tf(num1,den1);num2=1;den2=[0.51];G2=tf(num2,den2);num3=3;den3=[10];G3=tf(num3,den3);G=(G1+G2)*G3;T=feedback(G,G2);U=G/(1+G2*G);实验结果：»TI=2-25 +7-5=+6Q.25 +I•旣s'3+2s'2+5,55+6Cantinuous-tinetransferfunctiorii»U□=1.L25s*5r7.125e'4+16.5£'3+16.5e'2+6£0,.135s"7■+s'0+3.125F5+7了4+13.255^3+14.55*2+6s匚orrtini_ioti3-timetrarts-ferfunction.5•已知如图E2-2所示的系统框图试用MATLAB求该系统的闭环传递函数。代码：num仁10;den仁[11];G仁tf(num1,den1);num2=2;den2=conv([10],[11]);G2=tf(num2,den2);num3=[13];den3=[12];H仁tf(num3,den3);num4=[50];den4=[168];H2=tf(num4,den4);T仁feedback(G2,H1,1);G=G1*T1;T=feedback(G,H2);实验结果：»TT=20s'3+160 +400s+320s'G■*10s5+35亠44 +■82s^24116s-13Centinuous-timetransferfundion,四.实验收获与心得经过上学期的MATLABS验，对MATLA软件已经有一定的了解了。这次比较简单的几个实验让我学到了怎么将传递函数通过 MATLAB来表达，怎么画传递函数的零极点，怎么设置软件比较容易看。倒腾了半天，最大的收获是知道了conv函数是怎么用的，以及有三个因式的该怎样书写。实验2控制系统的暂态特性分析、实验目的学习和掌握利用MATLAB进行系统时域响应求解和仿真的方法。考察二阶系统的时间响应，研究二阶系统参数对系统暂态特性的影响。二、实验原理系统的暂态性能指标控制系统的暂态性能指标常以一组时域量值的形式给出，这些指标通常由系统的单位阶跃响应定义出来，这些指标分别为：（1）延迟时间td:响应曲线首次到达稳态值的50%所需的时间。（2）上升时间tr:响应曲线从稳态值的10%上升到90%所需要的时间长，对于欠阻尼系统，通常指响应曲线首次到达稳态值所需的时间峰值时间tp:响应曲线第一次到达最大值的时间。调整时间ts:响应曲线开始进入并保持在允许的误差(土2%或±5%范围内所需要的时间。(5)超调量6：响应曲线的最大值和稳态值之差，通常用百分比表示其中y(t)为响应曲线。在MATLAB中求取单位阶跃响应的函数为step，其使用方法如下step(sys)在默认的时间范围内绘出系统响应的时域波形step(sys,T)绘出系统在0-T范围内响应的时域波形step(sys,ts:tp:te) 绘出系统ts-te范围内，以tp为时间间隔取样的响应波形。[y,t]=step(…)该调用格式不绘出响应波形，而是返回响应的数值向量及其对应的时间向量。系统的暂态性能指标可以根据上述定义，在响应曲线上用鼠标读取关键点或通过搜索曲线对应的数值向量中关键点来确定。LTIViewer工具在MATLAB^提供了线性是不变系统仿真的工具LTViewer，可以方便地观察系统的响应曲线和性能指标。在命令窗口中键入litview即可启动LTIViewer。这里简要介绍LTIViewer工具。1） 【File】菜单Import选项：可以从Workspace或文件中导入系统模型。Export选项：将当前窗口中的对象模型保存到Workspace或文件中。Toolboxpreferences选项：属性设置功能，可以设置控制系统中得各种属性值。PageSetup选项：页面设置功能，可以对打印输出和显示页面进行设置。2） 【Edit】菜单PlotConfiguration 选项：对显示窗口及显示内容进行配置。LineStyle选项：线型设置功能，可以对输出响应曲线的线型进行设置。ViewerPreferences选项：对当前窗口的坐标、颜色、字体、响应曲线的特性参数等属性进行设置。3）右键菜单在运行界面上点击鼠标右键，将会弹出一个弹出式菜单，菜单上个选项的功能分别为：PlotTypes:选择绘制的系统曲线类型，可选的类型有单位阶跃响应、单位冲击响应、波特图、奈奎斯特图、零极点图等。System:选择需要仿真的系统。Characteristic :系统的性能指标选项。Grid:显示和关闭网格。Normalize:正常显示模式。Full View:满界面显示模式。Properties:性能编辑器选项，可以设置画面的标题、坐标标志、坐标范围、线型、颜色、性能指标等。三、实验内容1.已知单位负反馈系统前向通道的传递函数为G(沪试用MATLAB^制系统的单位阶跃响应曲线。代码：num仁80;den仁[120];G仁tf(num1,den1);T=feedback(G1,1);step(T);实验结果：StepResponse1.8 r r ( ■ , t n1.6--IEi1.4- -1.2 -I1：f\f\e1 upm0.8ii0.60.40.2r0L I I E J E E :0 1 2 3 4 5 6 7Time(seconds)2.已知二阶系统Z=6.0,3n=5,试用MATLAB^制系统单位阶跃响应曲线，并求取系统的暂态性能指代码：Wn=5;E=0.6;num仁WnT;den1=[12*E*Wnnum1];G=tf(num1,den1);step(G);oyniDi11.9 fl■Risetime(seconds):0.371 System:GSettlingtime(seconds):1.19System:GPeakamplitude:1.09Overshoot(%):9.48System:GSettlingtime(seconds):1.19System:GPeakamplitude:1.09Overshoot(%):9.48Attime(seconds):0.783System:GFinalvalue:1660.40.2-0L Ir L I J L [ L [ :0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6Time(seconds)由图中可以看到：上升时间tr:响应曲线上升到90%所需要的时间为0.371s。峰值时间tp:响应曲线第一次到达最大值的时间为0.787s调整时间ts:响应曲线开始进入并保持在允许的误差范围内所需要的时间为 1.19s。超调量6：响应曲线的最大值和稳态值只差的百分比为6=9%3n=1 ,Z从0变化到2,求此系统的单位阶跃响应。畑1;forE=C.4*i,nunl=Wn'1:denl=[l2*E*WnntiiiLl]:(numljdeni);step^30):holdon;endlegend(E=0J/E=Q.4J/E=On8/E=1+2J/E=1,6J/E=2J):实验结果:StepResponse2[ L斗 1jF1.E=01.8I\.1\\1\\E=0.41 J\E=0.81.6GJ1 1 / .E=1.2,1 I i1 1 \ E=1.61.4- j i 1 1•E=21.21 1j 11-；/N i Ii i| \ i \1P 111/ 1 I0.8r ? 1 I 1 I-0.6jl/,/. :rJj//Z-l: i11Lr■0.4-!?./ Ii1 I |f-0.2拙 1J ! ij \i h \j1fI1■IJ0/ r rJ「 I,fc1'0 5 10 15 20 25 30Time(seconds)Z=5. 0 ,3n从0变化到(3n半0)，求此系统的单位阶跃响应。代码:E=0.5：fori=0;5Wn=0.2*i:nujn :denl=[l2*E*Wnnum];G=tfCnrrmljdeni):step(G1,30}:heldon;endlegendCVofQjj5 『J*n=CL4’JWn=CL3Wn=Q»『」'恤=1‘)：实验结果:

StepResponse1.41.210.80.6StepResponse1.41.210.80.60.40.200 5 10 15 20 25 30Time(seconds)edm观察上述实验结果，分析这两个特征参数对系统暂态特性的影响。由上面的两个图可以看到，Z参数主要影响系统的稳定性，超调量，上升时间。Z越大，上升时间，峰值时间变长，而超调量变小。3n参数主要影响时间，3n越大，上升时间，峰值时间越小，而对超调量等其他因素近似没有影响。四、实验心得通过这个实验，对二阶系统的3n，Z对系统的影响有更深的理解，也知道怎么系统更稳定，达到稳定的时间更短，更快，更好，对以后的实际电路有深远的影响实验3根轨迹分析一、实验目的学习和掌握利用 MATLAB绘制根轨迹图的方法。学习和掌握利用系统根轨迹图分析系统的性能。二、实验原理根轨迹分析的 MATLAB现根轨迹是指系统某一参数变化时，闭环特征根在s平面上运动的轨迹。在MATLAB^，提供了用于根轨迹分析的专门函数。rlocus函数该函数的使用方法如下：rlocus(sys) 绘制单输入单输出LTI系统的根轨迹图。rlocus(sys,k) 使用用户指定的根轨迹增益k来绘制系统的根轨迹图。[r,k]=rlocus(sys)返回根轨迹增益值和闭环极点值，不绘制根轨迹图。rlocfind 函数该函数的使用方法如下：[k,poles]=rlocfind(sys) 计算鼠标选取点处的根轨迹增益值和闭环极点值，可在图形窗口根轨迹图中显示出十字光标，当用户选择其中一点时，相应的增益值和极点值记录在 k和poles中。[k,poles]=rlocfind(sys,p) 计算最靠近给定闭环极点p处的根轨迹增益。sgrid函数该函数的使用方法如下：Sgrid可在连续系统根轨迹或零极点图上绘制出栅格线，栅格线由等阻尼系数和等自然频率线构成。sgrid('new')先清除当前的图形，然后绘制出栅格线，并将坐标轴属性设置成holdon。sgrid(z,Wn) 指定阻尼系数 z和自然频率Wnsgrid(z,Wn,'new')指定阻尼系数z和自然频率Wn在绘制栅格线之前清除当前的图形并将坐标轴属性设置成holdon.Rltool工具MATLAB提供了一个根轨迹设计工具Rltool， 在命令窗口输入 rltool 命令即可启动该工具，也可输入rltool (sys)命令打开带系统模型sys的根轨迹设计器，运行界面如图3-1所示。Rltool 工具既可以分析系统根轨迹也能对系统进行设计， 具体使用方法请参考MATLAB^助或查阅相关资料三、实验内容1.已知系统开环传递函数为(1)使用MATLAB绘制系统的根轨迹图。代码：hujuf[15]:d*n^coiw(eon*j([1l]j[13])?[112])：&=tf「nun,den>:zpk1=zpk■G);rlocus1zpk1j实验结果:RootLocus-12-10-8 -6 -4RealAxis(seconds-1)-202o2ooo20RootLocus-12-10-8 -6 -4RealAxis(seconds-1)-202o2ooo20-30丿SGnocescRTXAvyrnrig卩-30-14(2)求根轨迹的两条分支离开实轴时的K值，并确定该K值对应的所有闭环极点。代码及结果:»[kjpoLes]=rloefind.1zpkl)a inthegraphicswiiwiowselecled_poirrt=-2.1213-0-000713.3812-u.7575-2-1213+0.00071-2.1213-0,C007i

;3O-30-14-12-10-8 -6 -4 -2-1RealAxis(seconds )020oooJr-sanocesctpXAWUFanLyam;3O-30-14-12-10-8 -6 -4 -2-1RealAxis(seconds )020oooJr-sanocesctpXAWUFanLyamo-2分析可得:k=3.3812;poles=-11.7575,-2.1213-0.0007i,-2.1213+0.0007i。(3)以区间[-40, -5]之间的值替代12- =s处的极点，重新绘制根轨迹图,观察其对根轨迹图的影响。代码：fori=5:5:40num=[15];den=conv(conv([11],[13]),[1i]);G=tf(num,den);zpk1=zpk(G);rlocus(zpk1);holdon;end

40-40-45 -40 -35 -30 -25 -20 -15 -10 -51RealAxis(seconds )03020oo-w如sanocesforxAWUFa^^aam40-40-45 -40 -35 -30 -25 -20 -15 -10 -51RealAxis(seconds )03020oo-w如sanocesforxAWUFa^^aamo-3从图中可以看到，轨迹的分支数并没有增加或减少，渐进线也并没变化，极值越大，轨迹分离之后趋近于无穷大的时间越快，分离点没有变化。2.已知系统开环传递函数为r(5+s)使用MATLAB绘制系统的根轨迹图代码::nuui»=[1S]:d*n=conveonv([12],[10])r[1S32])；fnum,den'1:zpkl=zplt(G):

-25-25Avi r i -25-25Avi r i rii-20-20 -15-10 -5 0-1RealAxis(seconds )5 10 15巧205o5o-50JrsanocesforxAvyanLyam—5计算两条分支进入右半平面和两条分支复数极点出发在实轴相交处的 K值。代码及结果：»[kjpolfls]=rlocfindtipk1)SuLtct盘paintinthugraphicswiadQwsel«ded_paint=-0.0237+3.4938i54.7660poles=-5.0253+3.354Si-5.0253-3.354Si0.0253+3.4641i0.0363-3U641i

I 1 1,r 厂一+xri L+k \i r r t\\i r-20-25-20-15 -10 -5 0-1RealAxis(seconds )5 10 15巧205O5o-5o5)sanocescKKXAyapowm-25代码及结果:-25>>[It,,pflles] ind.'.rpkliSelectapaintinthegraphicswindows^l«cted_point二-ID.1650-0.077Bik=2,0Q93e+0Jpales-5,2414+11.41751氐3414-11.41751-10.2414十0.0257i-10.2414-0.02571JrSonocetpXAwvranlaaml-25-25+*-JrSonocetpXAwvranlaaml-25-25+*-/ T—■ Ji°1 [fT f-20 -15-10 -5 05 10 152015050505-20-1RealAxis(seconds)以区间［-20,-1］之间的值替代零点的位置，重新绘制根轨迹图，观察其对根轨迹图的影响。代码：for1=0.25:5p=10—2-H-lj=』亠曲]；zpkl^zpktz,p』):rlocus(zpkl':holdon;30sanoceSvCEXAycazjgL L L L I ! ——I L'■y-K.r r r f r、*一i r30sanoceSvCEXAycazjgL L L L I ! ——I L'■y-K.r r r f r、*一i rRootLocus-30 -25 -20 -15 -10 -505 10 151RealAxis(seconds)20ooo20--30分析：零点绝对值越大，左半支汇合点越靠近左边，绝对值越大。已知单位负反馈系统的开环传递函数为(1)使用MATLAB^制系统的根轨迹图。代码：nm=l；dcn=CI450】；G=tf(nun,dmj;zpkl=zpkfG):rlocus'.zpk1):结果;-3 -2 -1-5 -4-2-3-4-6)sdnocesctpxAvran—9aml-1RealAxis(seconds)(2)分析使系统稳定的K值范围和使系统无超调的K值范围，并通过观察系统的单位阶跃响应曲线加以验证。代码：fork=0：4:20：num*;den=[1450]；Gnf(nun,d&n);g=feedbackftj1)；step3Q>:holdon：end応endCk=0‘,*k=4-/k=8#>^=1?/k=16\Jk=20J):结果: k=0 k=0 k=4 k=8 k=12 k=16 k=201.5eup1m0.500 5 10 15 20 25 30Time(seconds)分析：由上图可以看出，使得系统稳定的 K值范围为：0<k<20代码："fork=0:2:20\nujnFOi2#k;den=[1450]:G二tffnum,den);(G-,1):s+ep(gj30);holdon:endlegend Jlc=2JJ1k=『J1^10*f'k=14’Jk=E"J Jk=15;''k=I6';'k=1g'Jk=2Cl‘)：StepResponse1.41.210.80.60.40.200 5 10 15 20 25 30Time(seconds)I出II分析：由上图可以看出，使得系统无超调的K值范围为：0<k<=2四、实验心得刚刚学的自控的根轨迹画法，由于用手画，开始觉得根轨迹其实不是那么的好用，但是通过这次试验，通过用MATLAB^验快速准确的画出根轨迹之后，发现可以很直观的看出使系统具有各种各样的性质需要的K值范围，从而可以很快的找到符合要求的值。实验4系统的频率特性分析一、 实验目的学习和掌握利用MATLAB绘制系统Nyquist图和Bode图的方法。学习和掌握利用系统的频率特性分析系统的性能。二、 实验原理系统的频率特性是一种图解方法，分析运用系统的开环频率特性曲线，分析闭环系统的 性能，如系统的稳态性能、暂态性能常用的频率特性曲线有Nyquist图和Bode图。在MATLAB中提供了绘制Nyquist图和Bode图的专门函数。Nyquist图nyquist函数可以用于计算或绘制连续时间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,phase,w)[Gm,Pm,Wcg,Wcp]=margin(sys)[Gm,Pm,Wcg,Wcp]=margin(mag,phase,w)

其中不带输出参数时，可绘制出标有幅度裕度和相位裕度值的Bode图，带输出参数时，返回幅度裕度Gm相位裕度Pm及其对应的角频率Wcg和Wcp、实验内容1.已知系统开环传递函数为1000绘制系统的Nyquist图'并讨论其稳定性。代码：num=1000;den=conv([132],[15]);G=tf(num,den);nyquist(G);8060NyquistDiagramoo42ooo24--oPXAVFanTQap8060NyquistDiagramoo42ooo24--oPXAVFanTQap00di806020o分析：系统没有开环正实部极点，但是上图在(-1，j0)左边有负穿越，所以系统不稳定。2.已知系统的开环传递函数为*)=1fio*)=1fio>(0.2V1 \5~——S+11?)13人40 )101(1)绘制系统的零极点图，根据零极点分布判断系统的稳定性。代码:num=conv([1.251.251],[010]);0 0],[10/3den=conv(conv(conv([20 0],[10/31]),[0.2/31]),[1/401]);G=tf(num,den);G仁zpk(G);pzmap(G1);margin(G);结果：

Pole-ZeroMap-0-0-