《北京理工大学自动控制原理实验报告(共36页).doc》由会员分享，可在线阅读，更多相关《北京理工大学自动控制原理实验报告(共36页).doc（36页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上 本科实验报告实验名称： 控制理论基础实验 课程名称：控制理论基础实验时间：任课教师：实验地点：实验教师：实验类型： 原理验证 综合设计 自主创新学生姓名：学号/班级：组 号：学 院：同组搭档：专 业：成 绩：实验1 控制系统的模型建立一、实验目的1、掌握利用MATLAB 建立控制系统模型的方法。2、掌握系统的各种模型表述及相互之间的转换关系。3、学习和掌握系统模型连接的等效变换。二、实验原理1、系统模型的 MATLAB描述系统的模型描述了系统的输入、输出变量以及内部各变量之间的关系，表征一个系统的模型有很多种，如微分方程、传递函数模型、状态空间模型等。这里主要介绍系
2、统传递函数（TF）模型、零极点增益（ZPK）模型和状态空间（SS）模型的MATLAB 描述方法。1）传递函数（TF）模型传递函数是描述线性定常系统输入-输出关系的一种最常用的数学模型，其表达式一般为 在MATLAB 中，直接使用分子分母多项式的行向量表示系统，即num = bm, bm-1, b1, b0den = an, an-1, a1, a0调用tf 函数可以建立传递函数TF 对象模型，调用格式如下：Gtf = tf(num,den)Tfdata 函数可以从TF 对象模型中提取分子分母多项式，调用格式如下：num,den = tfdata(Gtf) 返回cell 类型的分子分母多项式系数
3、num,den = tfdata(Gtf,v) 返回向量形式的分子分母多项式系数2）零极点增益（ZPK）模型传递函数因式分解后可以写成 式中, z1 , z2 , ,zm 称为传递函数的零点， p1,p2,pn称为传递函数的极点，k 为传递系数（系统增益）。在MATLAB 中，直接用z,p,k矢量组表示系统，其中z，p，k 分别表示系统的零极点及其增益，即：z=z1,z2,zm;p=p1,p2,pn;k=k;调用zpk 函数可以创建ZPK 对象模型，调用格式如下：Gzpk = zpk(z,p,k)同样，MATLAB 提供了zpkdata 命令用来提取系统的零极点及其增益，调用格式如下：z,p,
4、k = zpkdata(Gzpk) 返回cell 类型的零极点及增益z,p,k = zpkdata (Gzpk,v) 返回向量形式的零极点及增益函数pzmap 可用于求取系统的零极点或绘制系统得零极点图，调用格式如下：pzmap(G) 在复平面内绘出系统模型的零极点图。p,z = pzmap(G) 返回的系统零极点，不作图。3）状态空间（SS）模型由状态变量描述的系统模型称为状态空间模型，由状态方程和输出方程组成： 其中：x 为n 维状态向量；u 为r 维输入向量； y 为m 维输出向量； A 为nn 方阵，称为系统矩阵； B 为nr 矩阵，称为输入矩阵或控制矩阵；C 为mn 矩阵，称为输出矩
5、阵； D为mr 矩阵，称为直接传输矩阵。在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 矩阵4）三种模型之间的转换上述三种模型之间可以互相转换，MATLAB 实现方法如下TF 模型ZPK 模型：zpk(SYS)或tf2zp(num,den)TF 模型SS 模型：ss(SYS)或tf2ss(num,den)ZPK 模型TF 模型：tf
6、(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)2、系统模型的连接在实际应用中，整个控制系统是由多个单一的模型组合而成，基本的组合方式有串联连接、并联连接和反馈连接。下图分别为串联连接、并联连接和反馈连接的结构框图和等效总传递函数。在MATLAB 中可以直接使用“*”运算符实现串联连接，使用“”运算符实现并联连接。反馈系统传递函数求解可以通过命令feedback 实现，调用格式如下：T = feedback(
7、G,H)T = feedback(G,H,sign)其中，G 为前向传递函数，H 为反馈传递函数；当sign = +1 时，GH 为正反馈系统传递函数；当sign = -1 时，GH 为负反馈系统传递函数；默认值是负反馈系统。三、实验内容1、已知控制系统的传递函数如下试用MATLAB 建立系统的传递函数模型、零极点增益模型及系统的状态空间方程模型，并绘制系统零极点图。实验代码：num=2 18 40;den=1 5 8 6;Gtf=tf(num,den)Gzpk=zpk(Gtf)Gss=ss(Gtf)pzmap(Gzpk);grid on; 实验结果：（1）首先建立系统的传递函数模型描述，上述
8、程序的运行结果为：Gtf = 2 s2 + 18 s + 40 - s3 + 5 s2 + 8 s + 6（2）零极点增益模型为：Gzpk = 2 (s+5) (s+4) - (s+3) (s2 + 2s + 2)系统零极点图为： （3）状态空间方程模型：Gss = a = x1 x2 x3 x1 -5 -2 -1.5 x2 4 0 0 x3 0 1 0 b = u1 x1 4 x2 0 x3 0 c = x1 x2 x3 y1 0.5 1.125 2.5 d = u1 y1 02、已知控制系统的状态空间方程如下 试用MATLAB 建立系统的传递函数模型、零极点增益模型及系统的状态空间方程模型
9、，并绘制系统零极点图。实验代码：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=(Gss)pzmap(Gzpk);grid on;实验结果：（1）系统矩阵 a = x1 x2 x3 x4 x1 0 1 0 0 x2 0 0 1 0 x3 0 0 0 1 x4 -1 -2 -3 -4 b = u1 x1 0 x2 0 x3 0 x4 1 c = x1 x2 x3 x4 y1 10 2 0 0 d = u1 y1 02）再创建ZPK对象模型：Gzpk = a
10、= x1 x2 x3 x4 x1 0 1 0 0 x2 0 0 1 0 x3 0 0 0 1 x4 -1 -2 -3 -4 b = u1 x1 0 x2 0 x3 0 x4 1 c = x1 x2 x3 x4 y1 10 2 0 0 d = u1 y1 0（3）传递函数：Gtf = 2 s + 10 - s4 + 4 s3 + 3 s2 + 2 s + 1（4）零极点图：3、已知三个系统的传递函数分别为 试用MATLAB 求上述三个系统串联后的总传递函数。 实验代码：num1=2 6 5;den1=1 4 5 2;G1=tf(num1,den1)num2=1 4 1;den2=1 9 8 0;
11、G2=tf(num2,den2)z=-3 -7;p=-1 -4 -6;k=5;G3=zpk(z,p,k)G=G1*G2*G3实验结果：G1 = 2 s2 + 6 s + 5 - s3 + 4 s2 + 5 s + 2G2 = s2 + 4 s + 1 - s3 + 9 s2 + 8 sG3 = 5 (s+3) (s+7) - (s+1) (s+4) (s+6) G = 10 (s+3.732) (s+3) (s+7) (s+0.2679) (s2 + 3s + 2.5) - s (s+8) (s+6) (s+4) (s+2) (s+1)44、已知如图E2-1 所示的系统框图 试用MATLAB
12、求该系统的闭环传递函数。实验代码：num1=1;den1=1 1;G1=tf(num1,den1);num2=1;den2=0.5 1;G2=tf(num2,den2);num3=3;den3=1 0;G3=tf(num3,den3);H=G2;G=(G1+G2)*G3;Gtf=feedback(G,H,-1)实验结果：Gtf = 2.25 s2 + 7.5 s + 6 - 0.25 s4 + 1.25 s3 + 2 s2 + 5.5 s + 65、已知如图E2-2 所示的系统框图 试用MATLAB 求该系统的闭环传递函数。实验代码：num1=10;den1=1 1;G1=tf(num1,de
13、n1)num2=2;den2=1 1 0;G2=tf(num2,den2);num3=1 3;den3=1 2;H2=tf(num3,den3);num4=5 0;den4=1 6 8;H1=tf(num4,den4);G=G1*feedback(G2,H2,+1);Gtf=feedback(G,H1,-1) 实验结果：Gtf = 20 s3 + 160 s2 + 400 s + 320 - s6 + 10 s5 + 35 s4 + 44 s3 + 82 s2 + 116 s - 48四、体会和建议本次实验比较基础，学习如何创建传递函数模型，并得到对应的零极点模型和状态空间方程。零极点图之前在
14、理论课上画过，这次实验用MATLAB得以实现，非常直观和准确。在做框图问题时，要细心和耐心，注意连接顺序和反馈的部分。 实验2 控制系统的暂态特性分析 一、实验目的1、学习和掌握利用MATLAB 进行系统时域响应求解和仿真的方法。2、考察二阶系统的时间响应，研究二阶系统参数对系统暂态特性的影响。二、实验原理1、系统的暂态性能指标控制系统的暂态性能指标常以一组时域量值的形式给出，这些指标通常由系统的单位阶跃响应定义出来，这些指标分别为：（1）延迟时间 td ：响应曲线首次到达稳态值的50%所需的时间。（2）上升时间 tr ：响应曲线从稳态值的10%上升到90%所需要的时间长，对于欠阻尼系统，通常
15、指响应曲线首次到达稳态值所需的时间。（3）峰值时间tp ：响应曲线第一次到达最大值的时间。（4）调整时间ts ：响应曲线开始进入并保持在允许的误差（2%或5%）范围内所需要的时间。（5）超调量 ：响应曲线的最大值和稳态值之差，通常用百分比表示其中y(t )为响应曲线。在MATLAB 中求取单位阶跃响应的函数为step，其使用方法如下step(sys) 在默认的时间范围内绘出系统响应的时域波形step(sys,T) 绘出系统在0 T 范围内响应的时域波形step(sys,ts:tp:te) 绘出系统在ts te 范围内，以tp 为时间间隔取样的响应波形y,t = step() 该调用格式不绘出响
16、应波形，而是返回响应的数值向量及其对应的时间向量。系统的暂态性能指标可以根据上述定义，在响应曲线上用鼠标读取关键点或通过搜索曲线对应的数值向量中关键点来确定。2、LTI Viewer工具在MATLAB 中提供了线性是不变系统仿真的工具LTI Viewer，可以方便地观察系统的响应曲线和性能指标。在命令窗口中键入litview 即可启动LTI Viewer。这里简要介绍LTIViewer 工具的使用方法。1）【File】菜单Import 选项：可以从Workspace 或MAT 文件中导入系统模型。Export 选项：将当前窗口中的对象模型保存到Workspace 或文件中。Toolbox pr
17、eferences 选项：属性设置功能，可以设置控制系统中得各种属性值。Page Setup 选项：页面设置功能，可以对打印输出和显示页面进行设置。2）【Edit】菜单Plot Configuration 选项：对显示窗口及显示内容进行配置。Line Style 选项：线型设置功能，可以对输出响应曲线的线型进行设置。Viewer Preferences 选项：对当前窗口的坐标、颜色、字体、响应曲线的特性参数等属性进行设置。3）右键菜单在运行界面上点击鼠标右键，将会弹出一个弹出式菜单，菜单上个选项的功能分别为：Plot Types：选择绘制的系统曲线类型，可选的类型有单位阶跃响应、单位冲击响应、
18、波特图、奈奎斯特图、零极点图等。System：选择需要仿真的系统。Characteristic：系统的性能指标选项。Grid：显示和关闭网格。Normalize：正常显示模式。Full View：满界面显示模式。Properties：性能编辑器选项，可以设置画面的标题、坐标标志、坐标范围、线型、颜色、性能指标等。三、实验内容1、已知单位负反馈系统前向通道的传递函数为试用MATLAB 绘制系统的单位阶跃响应曲线。 实验代码：num=80;den=1 2 0;G=tf(num,den); /建立系统传递函数模型T=feedback(G,1); /负反馈连接step(T)实验结果 ：2、已知二阶系统
19、（1）= 0.6，=5，试用MATLAB 绘制系统单位阶跃响应曲线，并求取系统的暂态性能指标。（2）=1 ，从0 变化到2，求此系统的单位阶跃响应。（3）= 0.5， 从0 变化到 1（ 0），求此系统的单位阶跃响应。（4）观察上述实验结果，分析这两个特征参数对系统暂态特性的影响。（1）使用step函数可以绘制系统的单位阶跃响应曲线实验代码：num1=25;den1=1 6 25;G1=tf(num1,den1)step(G1) 实验结果：求系统的暂态性能指标 ：在workspace中输入ltiview，弹出LTI Viewer1）找到纵坐标为0.5的点，对应的时间为0.271s，即td=0.
20、271s。2）找到纵坐标为首次为1的点，对应的时间为0.556s，即tr=0.556s。3）找到纵坐标首次首次达到最大值的点（最大值为1.09），对应的时间为0.784s，即tp=0.784s。4）找到纵坐标为1.02的点，对应的时间为1.16s，即允许的误差取2%时，ts=1.16s。找到纵坐标为首次为1.05的点，对应的时间为1.03s，即ts=1.03s。5）超调量=（1.09-1）/1*100%=9%（2）取分别为0，0.5，1，1.5，2，分别求出单位阶跃响应实验代码：num2=1;den2=1 0 1;G2=tf(num2,den2); num3=1;den3=1 1 1;G3=t
21、f(num3,den3); num4=1;den4=1 2 1;G4=tf(num4,den4); num5=1;den5=1 3 1;G5=tf(num5,den5); num6=1;den6=1 4 1;G6=tf(num6,den6);ltiview然后加入G2G6，用不同颜色绘制出图像。放大，得到如下图形：（3）分别取=0.2，0.4，0.6，0.8，1实验代码：num7=0.04;den7=1 0.2 0.04;G7=tf(num7,den7); num8=0.16;den8=1 0.4 0.16;G8=tf(num8,den8); num9=0.36;den9=1 0.6 0.36
22、;G9=tf(num9,den9); num10=0.64;den10=1 0.8 0.64;G10=tf(num10,den10); num11=1;den11=1 1 1;G11=tf(num11,den11);系统单位阶跃响应的变化情况如图所示：（4）实验结果分析为阻尼比，称为无阻尼自振荡角频率。可以看出，保持不变，从0变化到2时，系统由无阻尼邓福振荡变成欠阻尼振荡，再到临界阻尼，最后达到过阻尼；而保持=0.5不变，从0变化到1时，系统处于欠阻尼状态，而越大，系统的上升时间、峰值时间和过渡时间越小，超调量不变，越快到达稳定状态。四、体会和建议这次实验是研究系统的暂态特性，主要分析了上升时
23、间、峰值时间、过渡时间、超调量等，并分析了当和n变化时，系统的单位阶跃响应会发生怎样的变化。这次实验还学习了使用LTIVIEWER，这个工具使用方便很有用处。 实验3 根轨迹分析一、实验目的1、学习和掌握利用MATLAB 绘制根轨迹图的方法。2、学习和掌握利用系统根轨迹图分析系统的性能。二、实验原理1、根轨迹分析的 MATLAB实现根轨迹是指系统某一参数变化时，闭环特征根在s 平面上运动的轨迹。在MATLAB 中，提供了用于根轨迹分析的专门函数。1）rlocus函数该函数的使用方法如下：rlocus(sys) 绘制单输入单输出LTI 系统的根轨迹图。rlocus(sys,k) 使用用户指定的根
24、轨迹增益k 来绘制系统的根轨迹图。r,k = rlocus(sys) 返回根轨迹增益值和闭环极点值，不绘制根轨迹图2）rlocfind函数该函数的使用方法如下：k,poles=rlocfind(sys) 计算鼠标选取点处的根轨迹增益值和闭环极点值，可在图形窗口根轨迹图中显示出十字光标，当用户选择其中一点时，相应的增益值和极点值记录在k和poles 中。k,poles=rlocfind(sys,p) 计算最靠近给定闭环极点p 处的根轨迹增益。3）sgrid函数该函数的使用方法如下：Sgrid 可在连续系统根轨迹或零极点图上绘制出栅格线，栅格线由等阻尼系数和等自然频率线构成。sgrid(new)
25、先清除当前的图形，然后绘制出栅格线，并将坐标轴属性设置成hold on。sgrid(z,Wn) 指定阻尼系数z 和自然频率Wn。sgrid(z,Wn,new) 指定阻尼系数z 和自然频率Wn，在绘制栅格线之前清除当前的图形并将坐标轴属性设置成hold on。2、Rltool工具MATLAB 提供了一个根轨迹设计工具Rltool，在命令窗口输入rltool 命令即可启动该工具，也可输入rltool（sys）命令打开带系统模型sys 的根轨迹设计器。Rltool 工具既可以分析系统根轨迹也能对系统进行设计。三、实验内容1、已知系统开环传递函数为（1）使用MATLAB 绘制系统的根轨迹图。（2）求根
26、轨迹的两条分支离开实轴时的K 值，并确定该K 值对应的所有闭环极点。（3）以区间-40，-5之间的值替代s = 12处的极点，重新绘制根轨迹图，观察其对根轨迹图的影响。根轨迹图可以用rlocus函数来实现，k,poles=rlocfind(sys)可以显示出相应的增益值和极点值。（1）实验代码：num=1,5;den=conv(1,1,1,3);den=conv(den,1,12);G=tf(num,den);rlocus(G) /绘制根轨迹图实验结果：（2）在（1）的基础上，输入k,poles=rlocfind(G)实验结果：selected_point = -2.1137 - 0.0932
27、ik =3.4099poles = -11.7554 + 0.0000i -2.1223 + 0.0926i -2.1223 - 0.0926i两条分支离开实轴时的K值为3.4099，该K值对应的所有闭环极点为-11.7554 + 0.0000i，-2.1223 + 0.0926i，-2.1223 - 0.0926i。（3）实验代码：ss=input(ss=);num=1,5;den=conv(1,1,1,3);den=conv(den,1,-ss);G=tf(num,den);rlocus(G)实验结果：a、ss=-402、已知系统开环传递函数为（1）使用MATLAB 绘制系统的根轨迹图。（
28、2）计算两条分支进入右半平面和两条分支复数极点出发在实轴相交处的K 值。（3）以区间-20，-1之间的值替代零点的位置，重新绘制根轨迹图，观察其对根轨迹图的影响。 （1）实验代码：num=1,8;den=conv(1,0,1,2);den=conv(den,1,8,32);G=tf(num,den);rlocus(G)实验结果：（2）在（1）的基础上输入k,poles=rlocfind(G)，并选定两条分支复数极点出发在实轴相交处的点。实验结果：k =2.0695e+03poles = 5.2416 +11.4178i 5.2416 -11.4178i -10.2608 + 0.0000i -
29、10.2225 + 0.0000i则这时的K值为2.0695e+03.再选定两条分支进入右半平面的点，结果为：k =54.0578poles = -5.0144 + 3.3568i -5.0144 - 3.3568i 0.0144 + 3.4462i 0.0144 - 3.4462i则此时的K值为54.0578.（3）实验代码：num=1,-ze;den=conv(1,0,1,2);den=conv(den,1,8,32);G=tf(num,den);rlocus(G) 实验结果：a、ze=-20b、ze=-1实验结果分析：可以看出，当零点越来越大时，汇合点越来越接近原点。3. 已知单位负反馈
30、系统的开环传递函数为（1）使用MATLAB 绘制系统的根轨迹图。（2）分析使系统稳定的K 值范围和使系统无超调的K 值范围，并通过观察系统的单位阶跃响应曲线加以验证。（1）实验代码：num=1;den=conv(1,0,1,4,5);G=tf(num,den);rlocus(G)实验结果：（2）在（1）的基础上，输入k,poles=rlocfind(G)，选择根轨迹与虚轴的交点，读出K值为：k= 19.7427.那么，使系统稳定的K值范围是k19.7427.下面用单位阶跃响应曲线进行验证：实验代码：K = input(K = );num = K;den = conv(1,0,1,4,5);G1
31、 = tf(num,den);G = feedback(G1,1,-1);step(G)实验结果： a、 K=20b、K=18使系统无超调的K值范围。实验代码：K = input(K = );num = K;den = conv(1,0,1,4,5);G1= tf(num,den);G= feedback(G1,1,-1);rltool(G)step(G)先输入K=1,得到根轨迹图，测出根轨迹分离点和会合点的K值，分别为1和0.852.因此，系统无超调的K值范围是0.852，1。然后用单位阶跃响应进行验证：实验结果：a、K=0.5b、 K=0.9c、K=2.8四、体会和建议这次实验主要研究了根
32、轨迹图，及能够影响根轨迹图的参数。MATLAB可以用不同的颜色清晰的画出不同的根轨迹，比手动画图要清晰的多。有了根轨迹图，分析系统的参数也更加容易了。 实验4 系统的频率特性分析一、实验目的1、学习和掌握利用MATLAB 绘制系统Nyquist 图和Bode 图的方法。2、学习和掌握利用系统的频率特性分析系统的性能。二、实验原理系统的频率特性是一种图解方法，分析运用系统的开环频率特性曲线，分析闭环系统的性能，如系统的稳态性能、暂态性能常用的频率特性曲线有Nyquist 图和Bode 图。在MATLAB中，提供了绘制Nyquist 图和Bode 图的专门函数。1、Nyquist图nyquist
33、函数可以用于计算或绘制连续时间LTI 系统的Nyquist 频率曲线，其使用方法如下：nyquist(sys) 绘制系统的Nyquist 曲线。nyquist(sys,w) 利用给定的频率向量w 来绘制系统的Nyquist 曲线。re,im=nyquist(sys,w) 返回Nyquist 曲线的实部re 和虚部im，不绘图。2、Bode图bode 函数可以用于计算或绘制连续时间LTI 系统的Bode 图，其使用方法如下：bode(sys) 绘制系统的Bode 图。bode(sys,w) 利用给定的频率向量w 来绘制系统的Bode 图。mag,phase=bode(sys,w) 返回Bode
34、图数据的幅度mag 和相位phase，不绘图。3、幅度和相位裕度计算margin 函数可以用于从频率响应数据中计算出幅度裕度、相位裕度及其对应的角频率，其使用方法如下：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、已知系统开环传递函数为绘制系统的Nyquist 图，并讨论其稳定性。实验代码：num
35、=1000;den=conv(1 3 2,1 5);Gtf=tf(num,den);nyquist(Gtf)实验结果：该系统不稳定。2、已知系统的开环传递函数为（1）绘制系统的零极点图，根据零极点分布判断系统的稳定性。（2）绘制系统Bode 图，求出幅度裕度和相位裕度，判断闭环系统的稳定性。bode 函数可以用于计算或绘制连续时间LTI 系统的Bode 图，margin 函数可以用于从频率响应数据中计算出幅度裕度、相位裕度及其对应的角频率。（1）实验代码：num=10*25/16 5/4 1;den=conv(1 0 0 ,10/3 1);den=conv(den, 0.2/3 1);den=
36、conv(den, 1/40 1);Gtf=tf(num,den);Gzpk=zpk(Gtf);pzmap(Gzpk);grid on; 实验结果：可以看出极点均位于左半平面，所以该系统稳定。（2）实验代码：num=10*25/16 5/4 1;den=conv(1 0 0 ,10/3 1);den=conv(den, 0.2/3 1);den=conv(den, 1/40 1);Gtf=tf(num,den);bode(Gtf);Gm,Pm,Wcg,Wcp=margin(Gtf)实验结果：Gm =0.1198Pm =60.6660Wcg =0.6614Wcp = 4.3878实验结果分析：在
37、开环幅频特性大于0dB的所有频段内，相频特性曲线对-180度线的正负穿越次数之差为0，可以判断出系统是稳定的。再使用工具读出对应关键点的坐标，可以求出相位裕度=60度，幅值裕度：-20lgK=-6.41，K=2.093、已知系统的开环传递函数为分别判断当开环放大系数K = 5和K = 20时闭环系统的稳定性，并求出幅度裕度和相位裕度。 （1）K=5实验代码：num=5;den=conv(1 0,1 1);den=conv(den,0.1 1);Gtf=tf(num,den);bode(Gtf);Gm,Pm,Wcg,Wcp=margin(Gtf);实验结果：Gm =2.2000Pm =13.57
38、09Wcg =3.1623Wcp =2.1020实验结果分析：在开环幅频特性大于0dB的所有频段内，相频特性曲线对-180度线没有穿越，可以判断出系统是稳定的。读出关键点坐标，求得相位裕度=15度，幅值裕度：-20lgK=5.63，K=0.52（2）K=20实验代码：num=20;den=conv(1 0,1 1);den=conv(den,0.1 1);Gtf=tf(num,den);bode(Gtf);Gm,Pm,Wcg,Wcp=margin(Gtf);实验结果：实验结果分析：在开环幅频特性大于0dB的所有频段内，相频特性曲线对-180度线的正负穿越次数之差为-1，可以判断出系统是不稳定的。四、体会和建议这次实验通过绘制频率特性图：乃奎斯特图和波特图，来分析系统的稳定性。这需要用到理论课上学到的知识，如何用着两种图来判断系统的稳定性。经过实验后，我对频率特性图有了更深的认识，对MATLAB的操作也更熟练了。专心-专注-专业