实验三-控制系统的根轨迹分析
实验五 基于MATLAB控制系统的根轨迹及其性能分析
实验五基于MATLAB控制系统的根轨迹及其性能分析一、实验目的1、熟练掌握使用MATLAB绘制控制系统零极点图和根轨迹图的方法。
2、学会分析控制系统根轨迹的一般规律。
3、利用根轨迹图进行系统性能分析。
4、研究闭环零、极点对系统性能的影响。
二、实验原理1、根轨迹与稳定性当系统开环增益从变化时,若根轨迹不会越过虚轴进入s右半平面,那么系统对所有的K 值都是稳定的;若根轨迹越过虚轴进入s右半平面,那么根轨迹与虚轴交点处的K值,就是临界开环增益。
应用根轨迹法,可以迅速确定系统在某一开环增益或某一参数下的闭环零、极点位置,从而得到相应的闭环传递函数。
2、根轨迹与系统性能的定性分析1)稳定性。
如果闭环极点全部位于s左半平面,则系统一定是稳定的,即稳定性只与闭环极点的位置有关,而与闭环零点位置无关。
2)运动形式。
如果闭环系统无零点,且闭环极点为实数极点,则时间响应一定是单调的;如果闭环极点均为复数极点,则时间响应一般是振荡的。
3)超调量。
超调量主要取决于闭环复数主导极点的衰减率,并与其它闭环零、极点接近坐标原点的程度有关。
4)调节时间。
调节时间主要取决于最靠近虚轴的闭环复数极点的实部绝对值;如果实数极点距虚轴最近,并且它附近没有实数零点,则调节时间主要取决于该实数极点的模值。
5)实数零、极点影响。
零点减小闭环系统的阻尼,从而使系统的峰值时间提前,超调量增大;极点增大闭环系统的阻尼,使系统的峰值时间滞后,超调量减小。
而且这种影响将其接近坐标原点的程度而加强。
三、实验内容1、绘制系统的零极点图直接在s 复平面上绘制系统对应的零极点位置,极点用“×”表示,零点用“○”表示。
例1、已知系统的开环传递函数,绘制系统的零极点图。
《图5》22s 5s 5G(s)H(s)s(s 1)(s 2s 2)++=+++2、绘制控制系统的根轨迹图并分析根轨迹的一般规律例2、若已知系统开环传递函数,绘制控制系统的根轨迹图,并分析根轨迹的一般规律。
Matlab实验三 绘制根轨迹
for
Principles of Automatic Control
实验三 绘制根轨迹 1: 绘制根轨迹 2: 参量分析
① 绘制根轨迹 rlocus(sys)
Gk
(s)
K *(s 1) (s 2)(s 3)
rlocus(num, den)
rlocus(sys,k)
②参量分析(根轨迹图上一顿乱点即可)
K1
(s a )nm
1
根轨迹渐进线的方程是新的根轨迹方程。
Байду номын сангаас
• 例: 绘制根轨迹及其渐近线
G(s)
K1
s(s 1)(s 2)
⑥讨论增加零点对根轨迹的影响 试试-2至-4之间的零点
G(s)
K1
s(s 1)(s 2)
r=rlocus(sys)
Gk
(s)
K *(s 1) (s 2)(s 3)
r=rlocus(num, den)
[r,k]=rlocus(sys)
[r,k]=rlocus(num, den)
③测量出根轨迹增益和对应闭环极点坐标,
在窗口显示
Gk
(s)
K *(s 1) (s 2)(s 3)
[k,poles] = rlocfind(sys)
[k,poles] = rlocfind(sys,p) P为已知的要研究的闭环极点。
④绘制零、极点以及在窗口显示零极点
pzmap(sys)
Gk
(s)
K *(s 1) (s 2)(s 3)
[p,z]=pzmap(sys) 求解零极点的好方法
⑤绘制根轨迹渐近线
一般人我不告诉他:
当根轨迹渐进线与实轴的交点已求出后, 可得到方程,这是根轨迹渐进线的方程。
《自动控制原理》实验指导书
《自动控制原理》实验指导书梅雪罗益民袁启昌许必熙南京工业大学自动化学院目录实验一典型环节的模拟研究--------------------------1 实验二典型系统时域响应和稳定性-------------------10 实验三应用MATLAB进行控制系统根轨迹分析----------15 实验四应用MATLAB进行控制系统频域分析------------17 实验五控制系统校正装置设计与仿真-----------------19 实验六线性系统校正-------------------------------22 实验七线性系统的频率响应分析---------------------26 附录:TDN—ACP自动控制原理教学实验箱简介----------31实验一 典型环节的模拟研究一. 实验目的1.熟悉并掌握TD-ACC +设备的使用方法及各典型环节模拟电路的构成方法。
2.熟悉各种典型环节的理想阶跃响应曲线和实际阶跃响应曲线。
对比差异、分析原因。
3.了解参数变化对典型环节动态特性的影响。
二.实验内容下面列出各典型环节的方框图、传递函数、模拟电路图、阶跃响应,实验前应熟悉了解。
1.比例环节 (P)A 方框图:如图1.1-1所示。
图1.1-1B 传递函数:K S Ui S Uo =)()( C 阶跃响应:)0()(≥=t Kt U O 其中 01/R R K =D 模拟电路图:如图1.1-2所示。
图1.1-2注意:图中运算放大器的正相输入端已经对地接了100K 的电阻,实验中不需要再接。
以后的实验中用到的运放也如此。
E 理想与实际阶跃响应对照曲线:① 取R0 = 200K ;R1 = 100K 。
② 取R0 = 200K ;R1 = 200K 。
2.积分环节(I)A .方框图:如右图1.1-3所示。
图1.1-3B .传递函数:TSS Ui S Uo 1)()(=C .阶跃响应: )0(1)(≥=t t Tt Uo 其中 C R T 0=D .模拟电路图:如图1.1-4所示。
自控实验报告实验三线性系统的根轨迹
实验三 线性系统的根轨迹一、实验目的1. 熟悉MATLAB 用于控制系统中的一些基本编程语句和格式。
2. 利用MATLAB 语句绘制系统的根轨迹。
3. 掌握用根轨迹分析系统性能的图解方法。
4. 掌握系统参数变化对特征根位置的影响。
二、实验报告1.根据内容要求,写出调试好的MATLAB 语言程序,及对应的结果。
2. 记录显示的根轨迹图形,根据实验结果分析根轨迹的绘制规则。
3. 根据实验结果分析闭环系统的性能,观察根轨迹上一些特殊点对应的K 值,确定闭环系统稳定的范围。
4.写出实验的心得与体会。
三、实验内容请绘制下面系统的根轨迹曲线同时得出在单位阶跃负反馈下使得闭环系统稳定的K 值的范围。
一、 )136)(22()(22++++=s s s s s Ks G1、程序代码:G=tf([1],[1,8,27,38,26]); rlocus (G); [k,r]=rlocfind(G)G_c=feedback(G,1); step(G_c)2、实验结果:-8-6-4-22468Root LocusReal AxisI m a g i n a r y A x i sselected_point = -8.8815 + 9.4658i k =1.8560e+04 r =-10.2089 + 8.3108i -10.2089 - 8.3108i 6.2089 + 8.2888i6.2089 - 8.2888iTime (seconds)A m p l i t u d eselected_point =-9.5640 - 7.6273i k =1.3262e+04 r =-9.5400 + 7.6518i -9.5400 - 7.6518i 5.5400 + 7.6258i5.5400 - 7.6258iTime (seconds)A m p l i t u d eTime (seconds)A m p l i t u d eselected_point =-0.0095 + 2.1118i k =73.9872 r =-3.9617 + 2.4724i -3.9617 - 2.4724i -0.0383 + 2.1409i -0.0383 -2.1409iTime (seconds)A m p l i t u d e3、结果分析:根轨迹与虚轴有交点,所以在K 从零到无穷变化时,系统的稳定性会发生变化。
自动控制原理 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)改变转折频率:超调量,调节时间,峰值时间,稳态误差同样有相应的变化。
根轨迹实验报告
根轨迹实验报告根轨迹实验报告引言:根轨迹是控制系统理论中的一个重要概念,它描述了系统在参数变化下的稳定性和响应特性。
本实验旨在通过实际操作和数据分析,深入理解根轨迹的原理和应用。
通过对比不同系统的根轨迹,可以更好地理解系统的稳定性和控制性能。
一、实验目的本实验的目的是通过实际操作和数据分析,加深对根轨迹的理解,掌握根轨迹的绘制方法和分析技巧。
同时,通过对比不同系统的根轨迹,分析系统参数对根轨迹的影响,进一步认识系统的稳定性和控制性能。
二、实验装置与方法实验所需的装置包括控制系统实验台、计算机和相应的控制软件。
实验过程中,首先将系统接入实验台,通过控制软件设置系统参数,然后进行数据采集和分析。
根据实验要求,可以改变系统参数、增加干扰等,观察根轨迹的变化。
三、实验结果与分析在实验过程中,我们分别绘制了不同系统的根轨迹,并进行了数据分析。
通过观察根轨迹的形状和位置,我们可以判断系统的稳定性和响应特性。
以一个简单的一阶系统为例,我们改变了系统的比例增益和时间常数,绘制了对应的根轨迹。
通过观察根轨迹的位置和形状,我们可以发现以下规律:当比例增益增大时,根轨迹向左移动,系统的稳定性增强;当时间常数增大时,根轨迹变得更加平缓,系统的响应速度变慢。
在另一个二阶系统的实验中,我们改变了系统的阻尼比和自然频率,绘制了对应的根轨迹。
通过观察根轨迹的形状和分布,我们可以得出以下结论:当阻尼比增大时,根轨迹变得更加收敛,系统的稳定性提高;当自然频率增大时,根轨迹变得更加散布,系统的响应速度增加。
通过对比不同系统的根轨迹,我们可以进一步分析系统的稳定性和控制性能。
例如,当两个系统的根轨迹重合或者相似,可以认为它们具有相似的稳定性和响应特性;而当根轨迹相交或者离散较大时,可能存在系统不稳定或者不良的控制性能。
四、实验总结通过本次实验,我们深入了解了根轨迹的原理和应用。
通过实际操作和数据分析,我们掌握了根轨迹的绘制方法和分析技巧。
线性系统的根轨迹-自动控制原理实验报告
自动控制原理实验报告实验题目:线性系统的根轨迹班级:学号:姓名:指导老师:实验时间:一、实验目的1. 熟悉MATLAB 用于控制系统中的一些基本编程语句和格式。
2. 利用MATLAB 语句绘制系统的根轨迹。
3. 掌握用根轨迹分析系统性能的图解方法。
4. 掌握系统参数变化对特征根位置的影响。
二、实验内容同时得出在单位阶跃负反馈下使得闭环系统稳定的K 值的范围。
2.1绘制下面系统的根轨迹曲线)136)(22()(22++++=s s s s s Ks G程序:G=tf([1],[1 8 27 38 26 0]); rlocus (G); %绘制系统的根轨迹[k,r]=rlocfind(G) %确定临界稳定时的增益值k 和对应的极点r G_c=feedback(G,1); %形成单位负反馈闭环系统 step(G_c) %绘制闭环系统的阶跃响应曲线-12-10-8-6-4-20246-10-8-6-4-20246810Root LocusReal AxisI m a g i n a r y A x i s0204060801001201400.10.20.30.40.50.60.70.80.91Step ResponseTime (sec)A m p l i t u d e得出在单位阶跃负反馈下使得闭环系统稳定的K 值的范围:K>28.74252.2绘制下面系统的根轨迹曲线)10)(10012)(1()12()(2+++++=s s s s s K s G 程序:G=tf([1 12],[1 23 242 1220 1000]); rlocus (G); %绘制系统的根轨迹[k,r]=rlocfind(G) %确定临界稳定时的增益值k 和对应的极点r G_c=feedback(G,1); %形成单位负反馈闭环系统 step(G_c) %绘制闭环系统的阶跃响应曲线-60-50-40-30-20-100102030-50-40-30-20-1001020304050Root LocusReal AxisI m a g i n a r y A x i s01234560.0020.0040.0060.0080.010.012Step ResponseTime (sec)A m p l i t u d e得出在单位阶跃负反馈下使得闭环系统稳定的K 值的范围: K>1.1202e+032.3绘制下面系统的根轨迹曲线)11.0012.0)(10714.0()105.0()(2++++=s s s s s K s G 程序:G=tf([5 100],[0.08568 1.914 17.14 100 0]); rlocus (G); %绘制系统的根轨迹[k,r]=rlocfind(G) %确定临界稳定时的增益值k 和对应的极点r G_c=feedback(G,1); %形成单位负反馈闭环系统step(G_c) %绘制闭环系统的阶跃响应曲线-60-50-40-30-20-10010203040-60-40-200204060Root LocusReal AxisI m a g i n a r y A x i s012345670.10.20.30.40.50.60.70.80.91Step ResponseTime (sec)A m p l i t u d e得出在单位阶跃负反馈下使得闭环系统稳定的K 值的范围:K> 7.8321根据实验结果分析根轨迹的绘制规则:⑴绘制根轨迹的相角条件与系统开环根轨迹增益 值的大小无关。
控制系统的根轨迹分析实验报告
一、实验目的1. 熟悉控制系统根轨迹的基本概念和绘制方法。
2. 掌握利用MATLAB软件绘制和分析控制系统根轨迹的方法。
3. 通过根轨迹分析,了解系统参数变化对系统性能的影响。
4. 培养实验操作能力和数据处理能力。
二、实验原理根轨迹是指当系统的某一参数(如开环增益K)从0变化到无穷大时,闭环系统的特征根在s平面上的变化轨迹。
通过分析根轨迹,可以了解系统在参数变化时的稳定性、瞬态响应和稳态误差等性能。
三、实验设备1. 计算机2. MATLAB软件3. 控制系统实验箱四、实验内容1. 绘制控制系统根轨迹(1)首先,根据实验要求,搭建控制系统的数学模型。
(2)利用MATLAB中的rlocus函数绘制系统的根轨迹。
(3)观察根轨迹的变化规律,分析系统在不同参数下的稳定性。
2. 分析系统性能(1)根据根轨迹,确定系统的稳定裕度,包括增益裕度和相位裕度。
(2)分析系统在不同参数下的瞬态响应,如上升时间、调整时间、超调量等。
(3)分析系统在不同参数下的稳态误差,如稳态误差和稳态误差系数。
3. 改变系统参数,观察根轨迹变化(1)改变系统的参数,如增益、时间常数等。
(2)重新绘制根轨迹,观察根轨迹的变化规律。
(3)分析系统参数变化对系统性能的影响。
五、实验结果与分析1. 绘制控制系统根轨迹(1)根据实验要求,搭建控制系统的数学模型,得到开环传递函数。
(2)利用MATLAB中的rlocus函数绘制系统的根轨迹。
(3)观察根轨迹的变化规律,分析系统在不同参数下的稳定性。
2. 分析系统性能(1)根据根轨迹,确定系统的稳定裕度,包括增益裕度和相位裕度。
(2)分析系统在不同参数下的瞬态响应,如上升时间、调整时间、超调量等。
(3)分析系统在不同参数下的稳态误差,如稳态误差和稳态误差系数。
3. 改变系统参数,观察根轨迹变化(1)改变系统的参数,如增益、时间常数等。
(2)重新绘制根轨迹,观察根轨迹的变化规律。
(3)分析系统参数变化对系统性能的影响。
根轨迹校正实验报告
根轨迹校正实验报告一、实验目的本实验旨在通过观察系统的根轨迹,对系统进行校正,以达到控制系统的稳定性、快速性和精确性要求。
二、实验原理1. 根轨迹根轨迹是指在极坐标系下,由系统特征方程的根在复平面内的运动轨迹。
2. 根轨迹的性质- 当系统的开环传递函数中,理论上根轨迹的起点是传递函数零点的位置。
- 根轨迹对称于实轴。
- 根轨迹总是从系统的零点出发,逐渐趋向于系统的极点。
3. 根轨迹设计的基本要求- 所有根轨迹应该位于左半平面。
- 根轨迹的密度越大,系统的稳定性越好。
- 根轨迹与虚轴的交点个数为系统开环传递函数的极点数与零点数之差。
- 根轨迹经过的区域越小,系统的快速性越好。
三、实验步骤本次实验使用了MATLAB软件进行根轨迹校正实验,具体步骤如下:1. 给定开环控制系统的传递函数,并画出其对应的零极点分布图。
通过观察零极点的位置,确定系统的初始根轨迹起点。
2. 使用MATLAB的rlocus函数,绘制出开环根轨迹。
通过该函数,我们可以根据系统传递函数的特点,得到根轨迹的形状。
3. 根据根轨迹的形状和性质,校正系统。
可以通过调整控制器的参数或改变系统的结构等方式,来使根轨迹满足系统的要求。
4. 经过多次调整和校正,得到符合要求的根轨迹。
通过观察根轨迹的形状和分布,判断系统是否稳定、快速和准确。
四、实验结果与分析经过根轨迹校正,我们得到了一条符合要求的根轨迹。
通过分析根轨迹的形状和性质,我们可以得出以下结论:1. 系统的稳定性由于根轨迹位于左半平面,且大部分根轨迹较为密集,因此系统的稳定性较好。
没有根轨迹位于右半平面,避免了系统的不稳定性。
2. 系统的快速性根轨迹的起点与旁边的极点较近,根轨迹与虚轴的交点附近也没有极点,因此根轨迹经过的区域较小。
这意味着系统的快速性较好,能够快速响应输入变化。
3. 系统的准确性根轨迹与实轴的交点个数与系统的极点数与零点数之差相符,说明系统的准确性较好。
这样的根轨迹设计使得系统能够准确响应输入信号,实现精确控制。
控制系统的稳定性分析实验报告
控制系统的稳定性分析实验报告一、实验目的1.了解控制系统的稳定性分析方法。
2.通过实验,掌握系统稳态误差、系统阻尼比、系统根轨迹等稳态分析方法。
3.掌握控制系统的稳定性分析实验步骤。
二、实验原理1.系统稳态误差分析系统稳态误差是指系统在达到稳态时,输出与输入之间的偏差。
对于稳态误差的分析,可以采用开环传递函数和闭环传递函数进行分析。
开环传递函数:G(s)闭环传递函数:G(s)/(1+G(s)H(s))其中,H(s)为系统的反馈环节,G(s)为系统的前向传递函数。
稳态误差可以分为静态误差和动态误差。
静态误差是指系统在达到稳态时,输出与输入之间的偏差;动态误差是指系统在达到稳态时,输出与输入之间的波动。
2.系统阻尼比分析系统阻尼比是指系统在达到稳态时,振荡的阻尼程度。
阻尼比越大,系统越稳定;阻尼比越小,系统越不稳定。
系统阻尼比的计算公式为:ζ=1/(2ξ)其中,ξ为系统的阻尼比,ζ为系统的阻尼比。
3.系统根轨迹分析系统根轨迹是指系统的极点随着控制参数变化而在复平面上的轨迹。
根轨迹分析可以用来判断系统的稳定性和性能。
系统的根轨迹可以通过以下步骤进行绘制:(1)确定系统的传递函数G(s)(2)将G(s)写成标准形式(3)计算系统的极点和零点(4)绘制系统的根轨迹三、实验步骤1.系统稳态误差分析实验(1)将系统的开环传递函数和闭环传递函数写出。
(2)通过实验,测量系统的静态误差和动态误差。
(3)根据静态误差和动态误差的测量结果,计算系统的稳态误差。
2.系统阻尼比分析实验(1)通过实验,测量系统的振荡频率和衰减周期。
(2)根据振荡频率和衰减周期的测量结果,计算系统的阻尼比。
3.系统根轨迹分析实验(1)将系统的传递函数写成标准形式。
(2)计算系统的极点和零点。
(3)绘制系统的根轨迹,并根据根轨迹的形状,判断系统的稳定性和性能。
四、实验结果分析通过实验,我们可以得到系统的稳态误差、阻尼比和根轨迹等数据。
根据这些数据,我们可以分析系统的稳定性和性能,并对系统进行优化。
自动控制原理第四章根轨迹法
根轨迹法可用于仿真和实验研究,通过模拟和实验 验证系统的性能和稳定性,为实际系统的设计和优 化提供依据。
根轨迹法的历史与发展
历史
根轨迹法最早由美国科学家威纳于1940年提出,经过多年的 发展与完善,已经成为自动控制领域中一种重要的分析和设 计方法。
发展
随着计算机技术和数值分析方法的不断发展,根轨迹法的应 用范围和精度得到了进一步拓展和提高。未来,根轨迹法有 望与其他控制理论和方法相结合,形成更加完善和高效的控 制系统分析和设计体系。
根轨迹的性能分析
根轨迹的增益敏感性和鲁棒性
通过分析根轨迹在不同增益下的变化情况,可以评估系统的性能和鲁棒性。
根轨迹与性能指标的关系
通过比较根轨迹与某些性能指标(如超调量、调节时间等),可以评估系统的 性能。
04
根轨迹法与其他控制方法的比较
根轨迹法与PID制根轨迹图,直观地分析系统的稳定性、响应速度和超调量等性
特点
根轨迹法具有直观、简便、易于掌握等优点,特别适合用于分析 开环系统的稳定性和性能。
根轨迹法的应用场景
控制系统设计
根轨迹法可用于控制系统设计,通过调整系统参数 ,优化系统的性能指标,如稳定性、快速性和准确 性等。
故障诊断与排除
根轨迹法可用于故障诊断与排除,通过观察系统根 轨迹的变化,判断系统是否出现故障,以及故障的 类型和程度。
在绘制根轨迹时,需要遵循一定 的规则,如根轨迹与虚轴的交点 、根轨迹的分离点和汇合点等。
03
根轨迹分析方法
根轨迹的形状分析
根轨迹的起点和终点
根轨迹的起点是开环极点的位置,而 终点是闭环极点的位置。通过分析起 点和终点的位置,可以判断根轨迹的 形状。
根轨迹的分支数
实验3 控制系统的根轨迹作图
自动控制原理实验报告课程名称 自动控制原理 成 绩 实验项目 控制系统的根轨迹作图 指导教师 齐立省 学生姓名 赵儒桐 学号 201100805035 班级专业 11电子信息工程 实验地点 综合楼226 实验日期 年 月 日一、实验目的1.利用计算机完成控制系统的根轨迹作图2.了解控制系统根轨迹图的一般规律3.利用根轨迹进行系统分析及校正二、实验步骤1.在Windows 界面上用鼠标双击matlab 图标,即可打开MATLAB 命令平台。
2.练习相关M 函数根轨迹作图函数:rlocus(sys)rlocus(sys,k)r=rlocus(sys)[r,k]=rlocus(sys)函数功能:绘制系统根轨迹图或者计算绘图变量。
图1-1 格式1:控制系统的结构图如图1-1所示。
输入变量sys 为LTI 模型对象,k 为机器自适应产生的从0→∞的增益向量, 绘制闭环系统的根轨迹图。
格式2:k 为人工给定的增益向量。
格式3:返回变量格式,不作图。
R 为返回的闭环根向量。
格式4:返回变量r 为根向量,k 为增益向量,不作图。
更详细的命令说明,可键入“help rlocus”在线帮助查阅。
例如:系统开环传递函数为)3)(1()(++=s s s k s G g方法一:根轨迹作图程序为k=1; %零极点模型的增益值z=[]; %零点p=[0,-1,-3]; %极点sys=zpk(z,p,k); %零点/极点/增益模型rlocus(sys)作出的根轨迹图如图1-2所示。
方法二:s=tf('s'); G1=1/(s*(s+1)*(s+3));rlocus(G1); 图1-2 gridK1=12;figure;step(feedback(G1*K1,1)) % 绘制K1=12的闭环单位反馈阶跃响应曲线闭合时域仿真simulink 模型:三、实验内容给定如下各系统的开环传递函数,作出它们的根轨迹图,并完成给定要求。
北京理工大学 控制理论基础实验报告
控制理论基础实验1.控制系统的模型建立2.控制系统的暂态特性分析3.根轨迹分析4.系统的频率特性分析一、实验目的实验一1.掌握利用MATLAB建立控制系统模型的方法。
2.掌握系统的各种模型表述及相互之间的转换关系。
3.学习和掌握系统模型连接的等效变化。
实验二1.学习和掌握利用MATLAB进行系统时域响应求解和仿真的方法。
2.考察二阶系统的时间响应,研究二阶系统参数对系统暂态特性的影响。
实验三1.学习和掌握利用MATLAB绘制根轨迹图的方法2.学习和掌握利用系统根轨迹图分析系统的性能。
实验四1.学习和掌握利用MATLAB绘制系统Nyquist图和Bode图的方法。
2.学习和掌握利用系统的频率特性分析系统的性能。
二、实验原理1)传递函数模型(TF)gtf=tf(num,den)2)零极点增益模型(ZPK)Gzpk=zpk(z,p,k)3)状态空间模型(SS)Gss=ss(a,b,c,d)4)三种模型之间的转换TF→ZPK:z pk(sys)TF→SS:ss(sys)ZPK→TF:t f(sys)ZPK→SS:s s(sys)SS→TF:tf(sys)SS→ZPK:z pk(sys)5)绘制系统零极点图Pzmap(gzpk);Grid on;6)系统模型的串联G(s)=G1(s)*G2(s)7)系统模型的并联G(s)=G1(s)+G2(s)8)系统模型的反馈连接T=feedback(G,H)T=feedback(G,H,sign)9)绘制阶跃响应step(sys)step(sys,T)10)线性时不变系统仿真工具ltiview11)绘制系统根轨迹图rlocus(sys)rlocus(sys,k)[r,k]=rlocus(sys)12)计算鼠标选择点处根轨迹增益值和闭环极点值[k,poles]=rlocfind(sys)13)在连续系统根轨迹或零极点图上绘制出栅格线sgrid(‘new’)sgrid(z,Wn)14)绘制系统的Nyquist图nyquist(SYS)nyquist(sys,w)15)绘制系统的Bode图bode(sys)bode(sys,w)16)从频率响应数据中计算幅度裕度,相位裕度及对应角频率margin(sys)[mag,phase]=bode(sys,w)三、实验结果实验一1)零极点图2)零极点图3)总串联函数Transfer function:10 s^6 + 170 s^5 + 1065 s^4 + 3150 s^3 + 4580 s^2 + 2980 s + 525---------------------------------------------------------------------------------------------------------------------- s^9 + 24 s^8 + 226 s^7 + 1084 s^6 + 2905 s^5 + 4516 s^4 + 4044 s^3 + 1936 s^2 + 384 s 4)闭环传递函数Transfer function:2.25 s^2 + 7.5 s + 6-------------------------------------------------------0.25 s^4 + 1.25 s^3 + 2 s^2 + 5.5 s + 65)闭环传递函数Transfer function:20 s^3 + 160 s^2 + 400 s + 320-------------------------------------------------------------------------s^6 + 10 s^5 + 35 s^4 + 44 s^3 + 82 s^2 + 116 s - 48%1num=[2 18 40]; den=[1 5 8 6]; gtf=tf(num,den) gzpk=zpk(gtf) gss=ss(gtf) pzmap(gzpk);grid on%2a=[0 1 0 00 0 1 00 0 0 1-1 -2 -3 -4];b=[0 0 0 0]’;c=[10 2 0 0];d=0;gss=ss(a,b,c,d); gtf=tf(gss); gzpk=zpk(gss); pzmap(gzpk)grid on%3g1a=[2 6 5]; g1b=[1 4 5 2];g2a=[1 4 1];g2b=[1 9 8 0];g3z=[-3 -7];g3p=[-1 -4 -6];g3k=5;g1tf=tf(g1a,g1b);g2tf=tf(g2a,g2b);g3zpk=zpk(g3z,g3p,g3k);g3tf=tf(g3zpk);g=g1tf*g2tf*g3tf%4g1=tf([1],[1 1]);g2=tf(1,[0.5 1]);g3=g2;g4=tf(3,[1 0]);g=feedback((g1+g2)*g4,g3)%5g1=tf(10,[1 1]);g2=tf(2,[1 1 0]);g3=tf([1 3],[1 2]);g4=tf([5 0],[1 6 8]);g=feedback(g1*(feedback(g2, g3,1)),g4)实验二12(1)t d=0.272 t r=0.371 t p=0.787 t s=1.19ϭ=9%(2)(3)(4)ξ变大,延迟时间,上升时间,峰值时间,调整时间均越来越长,超调量开始时减小,然后保持不变。
MATLAB 语言及其用 实验(答案)
《MATLAB 语言及其用》实验指导书目录实验一Matlab 使用方法和程序设计........................实验二控制系统的模型及其转换.............................实验三控制系统的时域、频域和根轨迹分析...........实验四动态仿真集成环境-Simulink.........................实验一Matlab使用方法和程序设计一、实验目的1、掌握Matlab软件使用的基本方法;2、熟悉Matlab的数据表示、基本运算和程序控制语句3、熟悉Matlab绘图命令及基本绘图控制4、熟悉Matlab程序设计的基本方法二、实验内容:1、帮助命令使用help命令,查找 sqrt(开方)函数的使用方法;在 CommandWindowL里输入help,接在在search里输入sqr即可。
sqrtSquare rootSyntaxB = sqrt(X)DescriptionB = sqrt(X) returns the square root of each element of the array X. For the elements of X that are negative or complex, sqrt(X) produces complex results.TipsSee sqrtm for the matrix square root. Examplessqrt((-2:2)')ans =0 + 1.4142i0 + 1.0000i1.00001.4142See Alsonthroot | realsqrt | sqrtm2、矩阵运算(1)矩阵的乘法已知A=[1 2;3 4]; B=[5 5;7 8];求A^2*BA=[1 2;3 4];B=[5 5;7 8];C=A^2*B>> format compactC =105 115229 251(2)矩阵除法已知 A=[1 2 3;4 5 6;7 8 9];B=[1 0 0;0 2 0;0 0 3];A\B,A/BA=[1 2 3;4 5 6;7 8 9];B=[1 0 0;0 2 0;0 0 3];C=A\B,D=A/BC =1.0e+016 *0.3152 -1.2609 0.9457-0.6304 2.5218 -1.89130.3152 -1.2609 0.9457D =1.0000 1.0000 1.00004.0000 2.5000 2.00007.0000 4.0000 3.0000(3)矩阵的转置及共轭转置已知A=[5+i,2-i,1;6*i,4,9-i];求A.', A'A=[5+i,2-i,1;6*i,4,9-i];B=A.', C=A'B =5.0000 + 1.0000i 0 +6.0000i2.0000 - 1.0000i 4.00001.0000 9.0000 - 1.0000iC =5.0000 - 1.0000i 0 -6.0000i2.0000 + 1.0000i 4.00001.0000 9.0000 + 1.0000i(4)使用冒号表达式选出指定元素已知: A=[1 2 3;4 5 6;7 8 9];求A中第3列前2个元素;A中所有列第2,3行的元素;A=[1 2 3;4 5 6;7 8 9];B1=A([1,2],[3])B2=A([2,3],:)B1 =36B2 =4 5 67 8 9方括号[]用magic函数生成一个4阶魔术矩阵,删除该矩阵的第四列A=magic(4)B=A(:,[1,2,3])或A=magic(4)A(:,4)=[]A =16 2 3 135 11 10 89 7 6 124 14 15 1B =16 2 35 11 109 7 64 14 153、多项式(1)求多项式4xxp的根=x)2(3--Y=[1 0 -2 -4];S=roots(Y)S =2.0000-1.0000 + 1.0000i-1.0000 - 1.0000i(2)已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ,求矩阵A的特征多项式;把矩阵A作为未知数代入到多项式中;A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4]P=poly(A)polyval(P,A)A =1.2000 3.0000 5.0000 0.90005.0000 1.7000 5.00006.00003.0000 9.0000 0 1.00001.00002.00003.00004.0000P =1.0000 -6.9000 -77.2600 -86.1300 604.5500ans =1.0e+003 *0.3801 -0.4545 -1.9951 0.4601-1.9951 0.2093 -1.9951 -2.8880-0.4545 -4.8978 0.6046 0.43530.4353 0.0840 -0.4545 -1.16174、基本绘图命令(1)绘制余弦曲线 y=cos(t),t∈[0,2π](2)在同一坐标系中绘制余弦曲线y=cos(t-0.25)和正弦曲线y=sin(t-0.5),t∈[0,2π](1)t=[0:0.05:2*pi];y=cos(t);plot(t,y)-1-0.8-0.6-0.4-0.20.20.40.60.81(2)t=[0:0.05:2*pi];y1=cos(t-0.25);y2=sin(t-0.5);plot(t,y1)hold onplot(t,y2)5、基本绘图控制绘制[0,4π]区间上的x1=10sint曲线,并要求:(1)线形为点划线、颜色为红色、数据点标记为加号;(2)坐标轴控制:显示范围、刻度线、比例、网络线(3)标注控制:坐标轴名称、标题、相应文本;t=[0:0.1:4*pi];x1=10*sin(t);plot(t,x1,'r-.+'); %画图,显示红色、点划线、标记加号;axis([0,15,-10,10]); %定义显示范围,横轴为[0,15],纵轴为[-10,10];title('曲线x1=10sint'); %显示标题;xlabel('T轴');ylabel('X1轴'); %显示坐标轴名称;set(gca,'xminortick','on');set(gca,'yminortick','on'); %显示刻度线;grid on %显示网络线T 轴X 1轴6、基本程序设计(1)编写命令文件:计算1+2+…+n<2000 时的最大n 值; (2)编写函数文件:分别用for 和while 循环结构编写程序,求2的0到n 次幂的和。
实验三 线性系统的根轨迹分析
实验三 线性系统的根轨迹分析09电信 任旭乐 20095042046一、 实验目的1.熟悉Matlab 的基本操作;2.掌握利用Matlab 函数实现系统根轨迹的绘制及设计的方法。
3.能够根据所得结果对系统进行性能分析。
二、 实验内容1、已知单位负反馈系统的开环传递函数为: (1)试画出K=0 →∞时的闭环系统根轨迹; (2)求出临界时的K 值及闭环极点; (3)求出使系统稳定的K 值的区间; (4)利用Matlab 函数将剩余的根求出。
程序: a=[1 0]; b=[0.05 1]; c=[0.05 0.2 1]; d=conv(a,b); e=conv(c,d); G=tf([1],e); figure(1); rlocus(G);[k,pole]=rlocfind(G);解:(1)根轨迹如图所示。
(2)临界时k=4.62;闭环极点p=0.336+4.34j (3)由图可知:0<k<4.62时系统稳定。
Root LocusReal AxisI m a g i n a r y A x i s-80-60-40-200204060-60-40-20204060System: G Gain: 4.62P ole: 0.336 + 4.34i Damping: -0.0772Overshoot (%): 128Frequency (rad/sec): 4.35System: G Gain: 0P ole: 0Damping: -1Overshoot (%): 0Frequency (rad/sec): 0System: GGain: 8.5P ole: -19.5Damping: 1Overshoot (%): 0Frequency (rad/sec): 19.5根轨迹2()(0.051)(0.050.21)KG s s s s s =+++2、已知单位负反馈系统的开环传递函数为:(1)试画出K=0 →∞时的闭环系统根轨迹;(2)找出ζ=0.707附近的点,绘制出其相应的单位阶跃响应曲线。
实验三 MATLAB用于根轨迹分析
实验三 MATLAB用于根轨迹分析
一、实验目的
通过使用MATLAB完成根轨迹绘制、部分分式展开以及根轨迹分析等工作。
二、实验原理
绘制根轨迹可用函数rlocus(num,den)或rlocus(num,den,k)。
其中num,den分别对应系统开环传递函数的分子系数和分母系数构成的数组。
如果参数k是指定的,将按照给定的参数绘制根轨迹图,否则k是自动确定的,k的变化范围为0到∞。
三、实验内容
用MATLAB绘制系统的根轨迹图。
四、实验代码
1、
num=[1];
den=[1 3 2 0];
rlocus(num,den)
2、
Gc=tf(1,[1 5]);
Go=tf([1 1],[1 8 0]);
H=tf(1,[1 2]);
rlocus(Gc*Go*H);
v=[-10 10 -10 13];
axis(v);
grid on
五、实验结果
1、
2、
六、实验总结
本次实验通过MATLAB实现了由系统结构图绘制根轨迹图。
七、实验心得
本次实验相对于前两次实验来说比较简单、较为容易实现,但是需要结合其它相关的函数比如说Gain、Pole、Damping等函数来加以理解。
自动控制原理实验报告_2
实验一matlab基本指令练习例1:num=[1,5];>> den=[1,2,3,4,5];>> G=tf(num,den)Transfer function:s + 5-----------------------------s^4 + 2 s^3 + 3 s^2 + 4 s + 5例2:num=6*[1,5];den=conv(conv([,3,1],[1,3,1]),[1,6]);>> tf(num,den)Transfer function:6 s + 30----------------------------------3 s^4 + 28 s^3 + 66 s^2 + 37 s + 6例3:Z=[-1.9294;-0.0353+0.9287j;-0.0353-0.9287j];P=[-0.9567+1.2272j;-0.9567-1.2272j;0.0433+0.6412j;0.0433-0.6412j]; G=zpk(Z,P,KGain)Zero/pole/gain:6 (s+1.929) (s^2 + 0.0706s + 0.8637)----------------------------------------------(s^2 - 0.0866s + 0.413) (s^2 + 1.913s + 2.421)例4:G1=tf(1,[1,2,1]);>> G2=tf(1,[1,1]);>> G=feedback(G1,G2)Transfer function:s + 1---------------------s^3 + 3 s^2 + 3 s + 2G1=tf(1,[1,2,1]);G2=tf(1,[1,1]);G=feedback(G1,G2,1)Transfer function:s + 1-----------------s^3 + 3 s^2 + 3 s例5G1=tf([1,7,24,24],[1,10,35,50,24]);G2=tf([10,5],[1,0]);H=tf([1],[0.01,1]);G_a=feedback(G1*G2,H)Transfer function:0.1 s^5 + 10.75 s^4 + 77.75 s^3 + 278.6 s^2 + 361.2 s + 120-------------------------------------------------------------------- 0.01 s^6 + 1.1 s^5 + 20.35 s^4 + 110.5 s^3 + 325.2 s^2 + 384 s + 120 例6:num=[6.8,61.2,95.2];>> den=[1,7.5,22,19.5,0];>> G=tf(num,den);>> G1=zpk(G)Zero/pole/gain:6.8 (s+7) (s+2)-------------------------s (s+1.5) (s^2 + 6s + 13)例7:Z=[-2,-7];>> P=[0,-3-2j,-3+2j,-1.5];>> K=6.8;>> G=zpk(Z,P,K);>> G1=tf(G)Transfer function:6.8 s^2 + 61.2 s + 95.2-------------------------------s^4 + 7.5 s^3 + 22 s^2 + 19.5 s例8:实验二应用MATLAB进行控制系统的根轨迹分析1.To get started, select MATLAB Help or Demos from the Help menu.>> b=[1 1];>> a1=[1 0];>> a2=[1 -1];>> a3=[1 4 16];>> a=conv(a1,a2);>> a=conv(a,a3);>> rlocus(b,a)>> p=1.5i;>> [k,poles]=rlocfind(b,a,p) k =22.5031poles =-1.5229 + 2.7454i-1.5229 - 2.7454i0.0229 + 1.5108i0.0229 - 1.5108i>> p=1.5108i;>> [k,poles]=rlocfind(b,a,p) k =22.6464poles =-1.5189 + 2.7382i -1.5189 - 2.7382i0.0189 + 1.5197i0.0189 - 1.5197i>> p=1.5197i;>> [k,poles]=rlocfind(b,a,p) k =22.7642poles =-1.5156 + 2.7323i-1.5156 - 2.7323i0.0156 + 1.5269i0.0156 - 1.5269i>> p=1.5269i;>> [k,poles]=rlocfind(b,a,p)k =22.8593poles =-1.5129 + 2.7275i-1.5129 - 2.7275i0.0129 + 1.5329i0.0129 - 1.5329i>> p=1.5329i;>> [k,poles]=rlocfind(b,a,p) k = 22.9385poles =-1.5107 + 2.7235i-1.5107 - 2.7235i0.0107 + 1.5378i0.0107 - 1.5378i>> p=1.5378i;>> [k,poles]=rlocfind(b,a,p) k =23.0031poles =-1.5088 + 2.7202i-1.5088 - 2.7202i0.0088 + 1.5418i0.0088 - 1.5418i>> p=1.5418i;>> [k,poles]=rlocfind(b,a,p) k =23.0558poles =-1.5073 + 2.7175i-1.5073 - 2.7175i0.0073 + 1.5451i0.0073 - 1.5451i>> p=1.5451i;>> [k,poles]=rlocfind(b,a,p) k =23.0992poles =-1.5061 + 2.7152i-1.5061 - 2.7152i0.0061 + 1.5479i0.0061 - 1.5479i>> p=1.5479i;>> [k,poles]=rlocfind(b,a,p) k =23.1361poles =-1.5051 + 2.7133i-1.5051 - 2.7133i0.0051 + 1.5502i0.0051 - 1.5502i >> p=1.5502i;>> [k,poles]=rlocfind(b,a,p) k = 23.1663poles =-1.5042 + 2.7118i-1.5042 - 2.7118i0.0042 + 1.5521i0.0042 - 1.5521i>> p=1.5521i;>> [k,poles]=rlocfind(b,a,p) k =23.1913poles =-1.5035 + 2.7105i-1.5035 - 2.7105i0.0035 + 1.5537i0.0035 - 1.5537i>> p=1.5537i;>> [k,poles]=rlocfind(b,a,p) k =23.2123poles =-1.5029 + 2.7094i-1.5029 - 2.7094i0.0029 + 1.5550i0.0029 - 1.5550i>> p=1.5550i;>> [k,poles]=rlocfind(b,a,p) k =23.2293poles =-1.5024 + 2.7085i-1.5024 - 2.7085i0.0024 + 1.5561i0.0024 - 1.5561i>> p=1.5561i;>> [k,poles]=rlocfind(b,a,p) k =23.2438poles =-1.5020 + 2.7077i-1.5020 - 2.7077i0.0020 + 1.5570i0.0020 - 1.5570i>> p=1.5570i;>> [k,poles]=rlocfind(b,a,p) k =23.2556poles =-1.5017 + 2.7071i-1.5017 - 2.7071i0.0017 + 1.5578i0.0017 - 1.5578i>> p=1.5578i;>> [k,poles]=rlocfind(b,a,p) k =23.2661poles =-1.5014 + 2.7066i-1.5014 - 2.7066i0.0014 + 1.5584i0.0014 - 1.5584i>> p=1.5584i;>> [k,poles]=rlocfind(b,a,p) k =23.2740 poles =-1.5012 + 2.7062i-1.5012 - 2.7062i0.0012 + 1.5589i0.0012 - 1.5589i>> p=1.5589i;>> [k,poles]=rlocfind(b,a,p) k =23.2805poles =-1.5010 + 2.7058i-1.5010 - 2.7058i0.0010 + 1.5593i0.0010 - 1.5593i>>实验三根轨迹(课本)例4-16:num=1;den=conv([1,0],conv([1,4],[1,4,20]));rlocus(num,den)例4-17:num=1;den=conv([1,1],[1,2,9]);sys=tf(num,den);rlocus(sys)[r,k]=rlocus(sys)r =1.0e+002 *Columns 1 through 6-0.0100 + 0.0283i -0.0095 + 0.0283i -0.0090 + 0.0283i -0.0082 +0.0285i -0.0067 + 0.0289i -0.0043 + 0.0300i-0.0100 - 0.0283i -0.0095 - 0.0283i -0.0090 - 0.0283i -0.0082 -0.0285i -0.0067 - 0.0289i -0.0043 - 0.0300i-0.0100 -0.0110 -0.0119 -0.0136 -0.0166 -0.0215Columns 7 through 12-0.0024 + 0.0312i -0.0009 + 0.0323i 0.0031 + 0.0363i 0.0079 +0.0420i 0.0136 + 0.0497i 0.0204 + 0.0598i-0.0024 - 0.0312i -0.0009 - 0.0323i 0.0031 - 0.0363i 0.0079 -0.0420i 0.0136 - 0.0497i 0.0204 - 0.0598i-0.0251 -0.0281 -0.0363 -0.0459 -0.0572 -0.0708Columns 13 through 150.0231 + 0.0640i 0.7526 + 1.3212i Inf0.0231 - 0.0640i 0.7526 - 1.3212i Inf-0.0763 -1.5352 Infk =1.0e+006 *Columns 1 through 110 0.0000 0.0000 0.0000 0.0000 0.00000.0000 0.0000 0.0000 0.0001 0.0001Columns 12 through 150.0003 0.0003 3.5494 Inf>> kg=spine(real(r(2,(9:12))),k(9:12),0)kg =24.0000例4-18:num=[1];>> den=conv([1,0],conv([1,1],conv([1,3.5],[1,6,13]))); >> rlocus(num,den);>> axis equal>> [kg,p]=rlocfind(num,den)Select a point in the graphics windowselected_point =-11.5789 +10.9846ikg =6.4164e+005p =-16.6028-6.5912 +13.7908i-6.5912 -13.7908i9.6426 + 8.5111i9.6426 - 8.5111i例4-19:n=[1,4];d=[1,1,0];>> sys=tf(n,d);>> rlocus(sys)>> axis equal>> [r,k]=rlocus(sys);>> ri=r(2,10:18)ri =Columns 1 through 6-0.9648 + 1.6697i -1.1358 + 1.9484i -1.5056 + 2.4038i -1.8755 + 2.7361i -2.2756 + 3.0044i -2.6756 + 3.2009iColumns 7 through 9-3.0757 + 3.3385i -3.4758 + 3.4242i -3.9086 + 3.4629i>> t=10:18;>> ma=min(angle(ri));>> ti=spline(angle(ri),t,mati =10>> hold on>> plot([0,2*real(r(2,10))],[0,2*imag(r(2,10))]);>> [wn,z]=damp(r(2,10))wn =1.9284z =0.5003>> mpmax=exp(-z*pi/sqrt(1-z*z))mpmax =0.1628实验四典型环节及阶跃响应测试一.比例环节二.积分环节三. 微分环节四. 惯性环节五.震荡环节。
自动控制理论实验报告
自动控制理论实验报告实验二控制系统的时域分析一、实验目的学习利用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.利用计算机完成控制系统的根轨迹作图2.了解控制系统根轨迹图的一般规律3.利用根轨迹图进行系统分析二实验要点1. 预习什么是系统根轨迹?2. 闭环系统根轨迹绘制规则。
三实验方法(一)方法:当系统中的开环增益k 从0到变化时,闭环特征方程的根在复平面上的一组曲线为根轨迹。
实验3 根轨迹分析
-20
-15
-10 Real Axis (seconds )
-1
-5
0
5
>> p = [-1 -3 -10];
Root Locus 30
20
Imaginary Axis (seconds-1)
10
0
-10
-20
-30 -12
-10
-8
-6
-4
-1
-2
0
2
Real Axis (seconds )
>> p = [-1 -3 -5];
10
0
-10
-20
-30 -14
-12
-10
-8
-6
-4
-2
0
2
Real Axis (seconds -1)
(2)
可知此时增益 K 为 3.38。 [a,pole]=rlocfind(sys) Select a point in the graphics window selected_point = -2.1213 + 0.0000i a = 3.3812 pole = -11.7575 -2.1213 + 0.0000i -2.1213 - 0.0000i (3) >> z = [-5]; >> p = [-1 -3 -40]; >> k = [1]; >> sys1 = zpk(z, p, k) sys1 = (s+5) -----------------(s+1) (s+3) (s+40) Continuous-time zero/pole/gain model. >>rlocus(sys1)
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
实验三 控制系统的根轨迹分析
一、实验目的
1.利用MATLAB 完成控制系统的根轨迹作图;
2.了解控制系统根轨迹图的一般规律;
3.利用根轨迹进行系统分析。
二、实验原理
与根轨迹相关的MATLAB 函数:
1.绘制根轨迹的函数为rlocus ,常用格式为:
rlocus(sys) sys 为系统开环传递函数名称;
rlocus(num,den,k) num,den 为开环传递函数分子分母多项式,k 为根轨迹增益。
k 的范围可以指定,若k 未给出,则默认k 从0→∞,绘制完整的根轨迹;
r= rlocus(num,den) 返回变量格式,不作图,计算所得的闭环根r ;
[r,k]= rlocus(num,den) 返回变量格式,不作图,计算所得的闭环根r 和开环增益k 。
2.利用函数rlocfind( )可以显示根轨迹上任意一点的相关数值,以此判断对应根轨迹增益下闭环系统的稳定性。
[k,r]=rlocfind(num,den) 运行后会有一个十字光标提示用户,在根轨迹上选择点,用鼠标单击选择后,在命令窗口就会显示此点的根轨迹增益及此时的所有闭环极点值。
例1 )
4)(1()(++=s s s k s G r k 在命令窗口输入:
k=1;
z=[];
p=[0,-1,-4];
[num,den]=zp2tf(z,p,k);
rlocus(num,den);
title(’G k 根轨迹’)
[k,r]=rlocfind(num,den)
3.当开环传递函数不是标准形式,无法直接求出零极点,可用pzmap( )绘制系统的零极点图。
pzmap(num,den) 在s 平面上作零极点图;
pzmap(num,den) 返回变量格式,不作图,计算零极点。
三、实验内容
给定如下各系统的开环传递函数,作出它们的根轨迹图,并完成给定要求。
1. )
2)(1()(1++=s s s k s G r k 要求:
(1) 准确记录根轨迹的起点、终点与根轨迹条数;
答:起点为0,-1,-2;终点为无穷处;共三条根轨迹。
(2) 确定根轨迹的分离点与相应的根轨迹增益;
答:分离点为(-0.463,0),k=0.382
(3) 确定临界稳定时的根轨迹增益k 1。
答:k1=6
z=[];
p=[0,-1,-2];
[num,den]=zp2tf(z,p,k);
rlocus(num,den);
title('gk1根轨迹');
[k,r]=rlocfind(num,den)
2. )
164)(1()1()(22++-+=s s s s s k s G r k 要求:
确定根轨迹与虚轴交点并确定系统稳定的根轨迹增益k r 范围。
答:与虚轴交点为(0,+-1.63)和(0,+-2.56);系统稳定的根轨迹增益24.1<k<35.9 num=[1,1];den=conv([1,0],conv([1,-1],[1,4,16]));
rlocus(num,den);
title('gk1根轨迹');
[k,r]=rlocfind(num,den)
3. )
2()3()(3++=s s s k s G r k 要求:
(1)确定系统具有最大超调量时的根轨迹增益;
答:k=0.456
(2)确定系统阶跃响应无超调时的根轨迹增益取值范围。
答:分离点为:(-1.27,0),无超调根轨迹范围:0~0.18
z=[-3];
p=[0,-2];
[num,den]=zp2tf(z,p,k);
rlocus(num,den);
title('gk3根轨迹');
[k,r]=rlocfind(num,den)
4.绘出课本中各题作业的根轨迹。
4.2 (1)
(2)
(3)
(4)
四、实验结论
根据绘制的各根轨迹,分析闭环极点在s平面上的位置与系统性能的关系。