自动控制原理实验五 利用MATLAB绘制系统根轨迹
实验五 利用MATLAB绘制系统根轨迹
实验五利用MATLAB绘制系统根轨迹一、实验目的（1）熟练掌握使用MA TLAB绘制控制系统零极点图和根轨迹图的方法；（2）熟练使用根轨迹设计工具SISO；（2）学会分析控制系统根轨迹的一般规律；（3）利用根轨迹图进行系统性能分析；（4）研究闭环零、极点对系统性能的影响。
二、实验原理及内容1、根轨迹与稳定性当系统开环增益从变化时，若根轨迹不会越过虚轴进入s右半平面，那么系统对所有的K值都是稳定的；若根轨迹越过虚轴进入s右半平面，那么根轨迹与虚轴交点处的K值，就是临界开环增益。
应用根轨迹法，可以迅速确定系统在某一开环增益或某一参数下的闭环零、极点位置，从而得到相应的闭环传递函数。
2、根轨迹与系统性能的定性分析1）稳定性。
如果闭环极点全部位于s左半平面，则系统一定是稳定的，即稳定性只与闭环极点的位置有关，而与闭环零点位置无关。
2）运动形式。
如果闭环系统无零点，且闭环极点为实数极点，则时间响应一定是单调的；如果闭环极点均为复数极点，则时间响应一般是振荡的。
3）超调量。
超调量主要取决于闭环复数主导极点的衰减率，并与其它闭环零、极点接近坐标原点的程度有关。
4）调节时间。
调节时间主要取决于最靠近虚轴的闭环复数极点的实部绝对值；如果实数极点距虚轴最近，并且它附近没有实数零点，则调节时间主要取决于该实数极点的模值。
5）实数零、极点影响。
零点减小闭环系统的阻尼，从而使系统的峰值时间提前，超调量增大；极点增大闭环系统的阻尼，使系统的峰值时间滞后，超调量减小。
而且这种影响将其接近坐标原点的程度而加强。
【自我实践5-1】在实验内容（2）中控制系统的根轨迹上分区段取点，构造闭环系统传递函数，分别绘制其对应系统的阶跃响应曲线，并比较分析。
1：阻尼比=0.00196，k=5.942：阻尼比=0.246,k=2.163：阻尼比=0.669 k=0.6694：阻尼比=1.0,k=0,3855：阻尼比=1.2（无此阻尼，取-0.2）,k=24.5将数据填入实验数据记录表中。
自动控制原理 matlab实验报告
自动控制原理实验（二）一、实验名称：基于MATLAB的控制系统频域及根轨迹分析二、实验目的：（1）、了解频率特性的测试原理及方法；（2）、理解如何用MATLAB对根轨迹和频率特性进行仿真和分析；（3）、掌握控制系统的根轨迹和频率特性两大分析和设计方法。
三、实验要求：（1）、观察给定传递函数的根轨迹图和频率特性曲线；（2）、分析同一传递函数形式，当K值不同时，系统闭环极点和单位阶跃响应的变化情况；（3）、K值的大小对系统的稳定性和稳态误差的影响；（4）、分析增加系统开环零点或极点对系统的根轨迹和性能的影响。
四、实验内容及步骤（1）、实验指导书：实验四（1）、“rlocus”命令来计算及绘制根轨迹。
会出根轨迹后，可以交互地使用“rlocfind”命令来确定点击鼠标所选择的根轨迹上任意点所对应的K值，K值所对应的所有闭环极点值也可以使用形如“[K, PCL] = rlocfind(G1)”命令来显示。
（2）、波特图：bode(G1, omga)另外，bode图还可以通过下列指令得出相位和裕角：[mag,phase,w] = bode(sys)（3）、奈奎斯特图：nuquist(G, omega)（2）课本：例4-1、4-2、4-7五实验报告要求（1）、实验指导书：实验四思考题请绘制下述传递函数的bode图和nyquist图。
1. 根据实验所测数据分别作出相应的幅频和相频特性曲线；2. 将思考题的解题过程(含源程序)写在实验报告中。
幅频特性曲线相频特性曲线Gs = zpk([10], [-5; -16; 9], 200)subplot(1, 2, 1)bode(Gs)gridsubplot(1, 2, 2)nyquist(Gs)grid（2）课本：例4-1、4-2、4-7图像结果：程序：Gs = zpk([-1], [0; -2; -3],1) rlocus(Gs)图像结果：程序：Gs = zpk([-2], [-1-j; -1+j],1) rlocus(Gs)程序：K=[0.5 1 2]for i=1:1:3num=[1,1,0,0]; den=[1,1,K(i)]; sys=tf(num,den); rlocus(sys); hold ongrid onend图像结果：目标：改变增益K和转折频率依次调节源程序：k1=[4.44,10,20];num=[1,2];den=conv([1,1],[1,2,4]);%一阶转折频率 1/T（wn1=2,wn2=1）二阶转折频率 wn3=wn'=2，伊布西塔=1/2 num1=[1,1];den1=conv([1,2],[1,2,4]);%一阶转折频率 1/T（wn1=1,wn2=2）二阶转折频率 wn3=wn'=2，伊布西塔=1/2 t=[0:0.1:7]; %for i=1:3g0=tf(k1(i)*num,den);g=feedback(g0,1);[y,x]=step(g,t);c(:,i)=y;g1=tf(k1(i)*num1,den1);g(1)=feedback(g1,1);[y1,x]=step(g(1),t);c1(:,i)=y1;endplot(t,c(:,1),'-',t,c(:,2),'-',t,c(:,3),'-',t,c1(:,1),'-',t,c1(:,2), '-',t,c1(:,3),'-');gridxlabel('Time/sec'),ylabel('out')结果分析：在本题中（1）改变k值：k值越大，超调量越大，调节时间越长，峰值时间越短，稳态误差越小（2）改变转折频率：超调量，调节时间，峰值时间，稳态误差同样有相应的变化。
自动控制原理MATLAB仿真实验报告
实验一 MATLAB 及仿真实验（控制系统的时域分析）一、实验目的学习利用MATLAB 进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性； 二、预习要点1、 系统的典型响应有哪些2、 如何判断系统稳定性3、 系统的动态性能指标有哪些 三、实验方法（一） 四种典型响应1、 阶跃响应：阶跃响应常用格式：1、)(sys step ；其中sys 可以为连续系统，也可为离散系统。
2、),(Tn sys step ；表示时间范围0---Tn 。
3、),(T sys step ；表示时间范围向量T 指定。
4、),(T sys step Y =；可详细了解某段时间的输入、输出情况。
2、 脉冲响应：脉冲函数在数学上的精确定义：0,0)(1)(0〉==⎰∞t x f dx x f其拉氏变换为：)()()()(1)(s G s f s G s Y s f ===所以脉冲响应即为传函的反拉氏变换。
脉冲响应函数常用格式： ① )(sys impulse ； ②);,();,(T sys impulse Tn sys impulse③ ),(T sys impulse Y =（二） 分析系统稳定性 有以下三种方法：1、 利用pzmap 绘制连续系统的零极点图；2、 利用tf2zp 求出系统零极点；3、 利用roots 求分母多项式的根来确定系统的极点 （三） 系统的动态特性分析Matlab 提供了求取连续系统的单位阶跃响应函数step 、单位脉冲响应函数impulse 、零输入响应函数initial 以及任意输入下的仿真函数lsim.四、实验内容 (一) 稳定性1． 系统传函为()27243645232345234+++++++++=s s s s s s s s s s G ，试判断其稳定性2． 用Matlab 求出253722)(2342++++++=s s s s s s s G 的极点。
%Matlab 计算程序num=[3 2 5 4 6];den=[1 3 4 2 7 2];G=tf(num,den);pzmap(G);p=roots(den)运行结果： p =+ - + -P ole-Zero MapReal AxisI m a g i n a r y A x i s-2-1.5-1-0.500.5-1.5-1-0.50.511.5图1-1 零极点分布图由计算结果可知，该系统的2个极点具有正实部，故系统不稳定。
利用MATLAB绘制系统根轨迹
利用MATLAB绘制系统根轨迹————————————————————————————————作者：————————————————————————————————日期：第4章 利用MATLAB 绘制系统根轨迹一、 利用MATLAB 绘制系统根轨迹相关知识假设闭环系统中的开环传递函数可以表示为：)()())(()())(()(021********s KG p s p s p s z s z s z s K den numK a s a s a s b b s b s K s G n m nn n n m m m m k =+⋅⋅⋅+++⋅⋅⋅++==++⋅⋅⋅++++⋅⋅⋅++=---- 则闭环特征方程为： 01=+dennumK特征方程的根随参数K 的变化而变化，即为闭环根轨迹。
控制系统工具箱中提供了rlocus()函数，可以用来绘制给定系统的根轨迹，它的调用格式有以下几种：rlocus(num ，den) rlocus(num ，den ，K) 或者 rlocus(G) rlocus(G ，K)以上给定命令可以在屏幕上画出根轨迹图，其中G 为开环系统G 0(s)的对象模型,K 为用户自己选择的增益向量。
如果用户不给出K 向量，则该命令函数会自动选择K 向量。
如果在函数调用中需要返回参数，则调用格式将引入左端变量。
如[R ，K]=rlocus(G)此时屏幕上不显示图形，而生成变量R 和K 。
R 为根轨迹各分支线上的点构成的复数矩阵，K 向量的每一个元素对应于R 矩阵中的一行。
若需要画出根轨迹，则需要采用以下命令：plot(R ，¹¹)plot()函数里引号内的部分用于选择所绘制曲线的类型，详细内容见表1。
控制系统工具箱中还有一个rlocfind()函数，该函数允许用户求取根轨迹上指定点处的开环增益值，并将该增益下所有的闭环极点显示出来。
这个函数的调用格式为：[K ，P]=rlocfind(G)这个函数运行后，图形窗口中会出现要求用户使用鼠标定位的提示，用户可以用鼠标左键点击所关心的根轨迹上的点。
实验 五 用MATLAB绘制系统根轨迹(2)
实验数据分析时可围绕要达成的实验目标，借助公式、模型或必要软件，结合作图、计算等方法对实验数据进行数据处理后，得出实验结果。实验数据分析过程应条理清晰，层次清楚。
rlocus(num,den) %直接在s平面上绘制系统的根轨迹图，[k,r]=rlocfind(num,den)%在作好的根轨迹图上，确定被选的闭环极点位置的增益值k和此时的闭环极点r(向量)的值。
在作出根轨迹图后，再执行该命令，命令窗口会出现提示语，“Select a point in the graphics windows”,此时将鼠标移至根轨迹图并选定位置，单击左键确定，出现“＋”标记，在MATLAB窗口上即得到该点的根轨迹开环增益K值和对应的所有闭环根r(列向量)。
rlocus(num,den),
[k,r]=rlocfind(num,den)
num=[1 3];
den=[1 2 0];
sgrid
rlocus(num,den);
[k,r]=rlocfind(num,den)
5、实验现象、实验数据记录：
6、实验现象、实验数据的分析：
借助适当原理、公式、模型、利用必要的软件或其它分析方法、手段，对实验过程中出现的现象、所获取的实验数据进行具体、充分的分析。本部分的撰写应完整呈现出分析过程。
num=[1 5 5];
den=[1 3 4 2 0]
G=tf(num,den);
G1=zpk(G);
z=G1.z;
p=G1.p;
Z=z{:};
MATLAB的根轨迹分析
基于MATLAB的根轨迹分析一.实验目的:1. 学习利用MATLAB的语言绘制控制系统根轨迹的方法。
2. 学习利用根轨迹分析系统的稳定性及动态特性。
二.实验内容：1. 应用MATLAB语句画出控制系统的根轨迹。
2. 求出系统稳定时，增益k的范围。
3. 分析系统开环零点和极点对系统稳定性的影响。
三.实验步骤1•给定某系统的开环传递函数G(s)H(s)=k/s(s*s+4s+16)用MATLAB与语言绘出该系统的根轨迹。
程序如下：num=[1];den=[1,4,16,0];G=tf(nu m,de n)G1=zpk(G)Z=tzero(G)P=pole(G)pzmap( nu m,de n);title('pole-zero Map')rlocus( nu m,de n)根轨迹如图-2-6-8-8Root Locus 结论：由上图可知增益k的取值范围：0<k<642.将系统的开环传递函数改为：G(s)H(s)=k/s(s*s+4s+5)绘出该系统根轨迹图，观察增加了开环零点后根轨迹图的变化情况。
程序如下：num=[1,1];den=[1,4,5,0];G=tf(nu m,de n)G1=zpk(G)Z=tzero(G)P=pole(G)pzmap( nu m,de n);title('pole-zero Map')rlocus( nu m,de n)根轨迹如图Root LocusrRrxA-2-2.5 -2 -1.5 -1-0.5Real Axis结论：增加了开环零点后根轨迹的变化tax A■『r r L ■iJ 1「「⑴改变了根轨迹在实轴上的分布。
(2) 改变了根轨迹渐近线的条数、倾角及截距口(3) 若增加的廿坏零点和某个极点車合或距离很近，构成开坏偶极子，则两者和互抵消。
因此，可加入•个零点來抵消有损于系统性能的极点耳(4) 根轨迹曲线将向左倔移，有利于改善系统的动态件能，而且，所加的零点越靠近虚轴，彫响越大。
根轨迹分析的MATLAB实现
函数命令调用格式：
[k,poles]=rlocfind(sys) [k,poles]=rlocfind(sys,p)
函数命令使用说明：
rlocfind（）函数命令可计算出与根轨迹上极点（位 置为poles）相对应的根轨迹增益（k）。rlocfind（）函 数既适用于连续时间系统，也适用于离散时间系统。
pzmap(a,b,c,d)函数可以在复平面内绘制用状态空间 模型描述系统的零极点图，对于MIMO系统，可绘制所有 输入到输出间的传递零点。
pzmap(sys)函数可以在复平面里绘制以传递函数模 型sys表示开环系统的零极点。传递函数模型sys即 G(s)=num(s)/den(s)。
pzmap(p,z)函数可在复平面里绘制零极点图，其中 行矢量p为极点位置，列矢量z为零点位置。这个函数命令 用于直接绘制给定的零极点图。
[k,poles]=rlocfind(sys,p)函数可对指定根计算对应 的增益与根矢量p。
[k,poles]=rlocfind(sys)函数输入参量sys可以是由函 数tf（）、zpk（）、ss（）中任何一个建立的LTI对象模 型。函数命令执行后，可在根轨迹图形窗口中显示十字形 光标，当用户选择根轨迹上一点时，其相应的增益由k记 录，与增益相关的所有极点记录在poles中。
2．求系统根轨迹的函数rlocus()
函数命令调用格式：
[r,k]=rlocus(a,b,c,d) [r,k]=rlocus(sys) [r,k]=rlocus(a,b,c,d,k) [r,k]=rlocus(num,den,k)
函数命令使用说明：
rlocus()函数命令用来绘制SISO系统的根轨迹图。 rlocus()函数既可适用于连续时间系统，也适用于离散时 间系统。
自动控制原理MATLAB仿真实验报告
实验一 典型环节的MATLAB 仿真 一、实验目的1．熟悉MATLAB 桌面和命令窗口，初步了解SIMULINK 功能模块的使用方法。
2．通过观察典型环节在单位阶跃信号作用下的动态特性，加深对各典型环节响应曲线的理解。
3．定性了解各参数变化对典型环节动态特性的影响。
二、实验内容① 比例环节1)(1=s G 和2)(1=s G ；Simulink 图形实现：示波器显示结果：② 惯性环节11)(1+=s s G 和15.01)(2+=s s GSimulink 图形实现：示波器显示结果：③ 积分环节s s G 1)(1Simulink 图形实现：示波器显示结果：④ 微分环节s s G )(1Simulink 图形实现：波器显示结果：⑤ 比例+微分环节（PD ）2)(1+=s s G 和1)(2+=s s G1）、G1（s ）=s+2Simulink 图形实现：示波器显示结果：2）、G2（s）=s+1 Simulink图形实现：示波器显示结果：⑥ 比例+积分环节（PI ）s s G 11)(1+=和s s G 211)(2+=1）、G1（1）=1+1/sSimulink 图形实现：示波器显示结果：2）G2（s）=1+1/2s Simulink图形实现：示波器显示结果：三、心得体会通过这次实验我学到了很多，对课本内容加深了理解，熟悉MATLAB桌面和命令窗口，初步了解SIMULINK功能模块的使用方法，加深对各典型环节响应曲线的理解，这为对课程的学习打下了一定基础。
实验二线性系统时域响应分析一、实验目的1．熟练掌握step( )函数和impulse( )函数的使用方法，研究线性系统在单位阶跃、单位脉冲及单位斜坡函数作用下的响应。
2．通过响应曲线观测特征参量ζ和nω对二阶系统性能的影响。
3．熟练掌握系统的稳定性的判断方法。
二、实验内容1．观察函数step( )的调用格式，假设系统的传递函数模型为243237()4641s s G s s s s s ++=++++绘制出系统的阶跃响应曲线？2．对典型二阶系统222()2n n n G s s s ωζωω=++1）分别绘出2(/)n rad s ω=，ζ分别取0,0.25,0.5,1.0和2.0时的单位阶跃响应曲线，分析参数ζ对系统的影响，并计算ζ=0.25时的时域性能指标,,,,p r p s ss t t t e σ。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验五利用MATLAB绘制系统根轨迹
一、实验目的
（1）熟练掌握使用MA TLAB绘制控制系统零极点图和根轨迹图的方法；
（2）熟练使用根轨迹设计工具SISO；
（2）学会分析控制系统根轨迹的一般规律；
（3）利用根轨迹图进行系统性能分析；
（4）研究闭环零、极点对系统性能的影响。
二、实验原理及内容
1、根轨迹与稳定性
当系统开环增益从变化时，若根轨迹不会越过虚轴进入s右半平面，那么系统对所有的K值都是稳定的；若根轨迹越过虚轴进入s右半平面，那么根轨迹与虚轴交点处的K值，就是临界开环增益。
应用根轨迹法，可以迅速确定系统在某一开环增益或某一参数下的闭环零、极点位置，从而得到相应的闭环传递函数。
2、根轨迹与系统性能的定性分析
1）稳定性。
如果闭环极点全部位于s左半平面，则系统一定是稳定的，即稳定性只与闭环极点的位置有关，而与闭环零点位置无关。
2）运动形式。
如果闭环系统无零点，且闭环极点为实数极点，则时间响应一定是单调的；如果闭环极点均为复数极点，则时间响应一般是振荡的。
3）超调量。
超调量主要取决于闭环复数主导极点的衰减率，并与其它闭环零、极点接近坐标原点的程度有关。
4）调节时间。
调节时间主要取决于最靠近虚轴的闭环复数极点的实部绝对值；如果实数极点距虚轴最近，并且它附近没有实数零点，则调节时间主要取决于该实数极点的模值。
5）实数零、极点影响。
零点减小闭环系统的阻尼，从而使系统的峰值时间提前，超调量增大；极点增大闭环系统的阻尼，使系统的峰值时间滞后，超调量减小。
而且这种影响将其接近坐标原点的程度而加强。
【自我实践5-1】
在实验内容（2）中控制系统的根轨迹上分区段取点，构造闭环系统传递函数，分别绘制其对应系统的阶跃响应曲线，并比较分析。
1：阻尼比=0.00196，k=5.94
2：阻尼比=0.246,k=2.16
3：阻尼比=0.669 k=0.669
4：阻尼比=1.0,k=0,385
5：阻尼比=1.2（无此阻尼，取-0.2）,k=24.5
3、基于SISO 设计工具的系统根轨迹设计
用根轨迹法进行系统校正过程中，分析补偿增益和附加实数（或复数）零极点之间匹配的规律。
在MA TLAB 命令窗口中输入命令【rltool 】，然后回车，或者输入【rltool （sys ）】函数命令，就可打开系统根轨迹的图形界面。
【综合实践】绘制180°根轨迹。
请绘制： （1）12()(1)(1)g
K G s T s T s =
++；（2）12(1)()(1)(1)
g K as G s T s T s +=
++；
（3）12()(1)(1)(1)
g
K G s T s T s ps =
+++的根轨迹，其中T 1=0.2，T 2=2，a=0.1，p=1，分析
附加零点、极点对根轨迹的影响；固定T 值，分别改变a 和p 的值看附加零、极点位置的变化对根轨迹形状的影响。
将结果填入下表。
根轨迹图
根轨迹的绘制及参量分析。
分别绘制（1）12(1)()(1)(1)
g K as G s T s T s +=
++；（2）12()(1)(1)(1)
g
K G s T s T s ps =
+++
的0°根轨迹，比较其与180°根轨迹不同。
其中T 1=0.2，T 2=2，a=0.1，p=1。
（1）0°根轨迹： 180°根轨迹：
（2）0°根轨迹： 180°根轨迹：
对于上面的180°和0°根轨迹，求系统临界稳定时的K g 值，求K g =5时系统极点的位置（在根轨迹上的小红块上点击鼠标右键显示极点坐标值，该小红块可以用鼠标拖动）；分析此时系统的阶跃动态响应和Bode 图。
可取其它值，再求出对应的闭环零、极点响应曲线）
g
伯德图的绘制：
伯德图
三、思考题 1)
附加开环零点总对系统的稳定性是否有利？ 不一定，引入负的零点会使系统更加稳定，但是引入一个负的零点可能会破坏系统的稳定性
2)
附加开环极点总对系统的稳定性是否不利？ 不一定，引入负的极点会使系统更加稳定，但是引入一个负的极点可能会破坏系统的稳定性
3)
对实际系统，如何通过附加零点和极点来改善系统性能？举例说明。
对传递函数12(1)()(1)(1)
g K as G s T s T s +=
++；通过增加一阶惯性环节来增加系统的极点，
通过增加一阶微分环节来增加系统零点，即可增加系统的稳定性。
4)
如何绘制系统等效根轨迹？
绘制根轨迹主要通过引入开环传递函数G 等效，在G 等效将参数变量置于常规根轨迹所对应的开环传递函数G 中K 的位置上，然后按照常规根轨迹的作图法进行作图。
四、实验能力要求
（1）熟练掌握使用MA TLAB 绘制控制系统零极点图和根轨迹图的方法。
（2）通过根轨迹图能够确定有用的参数，比如：分离点坐标及相应参数、临界开环增益点及相应参数、任意一点对应的开环增益K 值和闭环极点。
（3）利用根轨迹图进行系统性能分析，以阻尼比为依据分区段比较不同闭环极点对应系统性能的变化。
能够确定系统稳定的开环增益范围。
（4）了解闭环零、极点对系统性能的影响。
（5）学会使用根轨迹设计工具校验系统的动态性能和稳态性能。