MATLAB标准系统根轨迹和频域分析实验

1、实验三MATLAB系统根轨迹和频域分析实验一、实验目的1. 学习使用MATLAB求特征多项式的根，分析系统稳定性;2. 学习使用MATLAB由传递函数求零点和极点；3. 学习使用MATLAB绘制根轨迹；4. 掌握由根轨迹分析系统性能的方法；5. 学习使用MATLAB绘制Bode图和Nyquist图；6. 掌握使用Bode图和Nyquist图分析系统性能的方法。二、实验仪器计算机三、实验内容3.1特征多项式求解3.1.1直接求特征多项式的根设P为特征多项式的系数矢量，用MATLAB函数roots()可直接 求出方程P=0在复数范围内的解，该函数的调用格式为：v=roots(p)例二十三己知系统的

2、特征多项式为：，+ 3/+2，+卄1特征方程的解可由下面的MATLAB命令得出：p=lA3,2,l,lv=roots(p)结果显示：0.3202+1.7042i0.3202-1.7042i-0.72090.0402+0.6780i0.0402-0.6780i利用多项式求根函数roots(),可方便的求出系统的零点和极点， 然后根据零极点分析系统稳定性和其他性能。3.1.2由根创建多项式如果己知多项式的因式分解式或特征根，可由MATLAB函数 poly()直接得岀特征多项式系数矢量，其调用格式为：p=poly(v)o如上题中：v =0.3202+1.7042i； 0.3202-1.7042i；

3、-0.7209； 0.0402+0.6780i； 0.0402-0.6780i；p=poly(v)结果显示：p= 1.0000 -0.0000 3.0000 2.0000 1.0000 1.0000由此可见，函数roots()与函数poly()互为逆运算。3.1.3多项式求值在MALAB中通过函数polyval()可求得多项式在给定点的值，该 函数的调用格式为：polyval(p,v )o对于上题中的p值，求取多项式在x点的值，可输入如下命令：p=l,0,3,2,l,l；x=lpolyval(p,x)结果显示:ans=83.1.5由传递函数求零点和极点在MATLAB控制系统工具箱中，给出了由传

4、递函数对象G求系 统零点和极点的函数，其调用格式分别为：Z=tzero(G)P=pole(G)注意：上式中要求的G必须是零极点模型对象。_6.2? +61.2s+ 95.2例二十四 己知传递函数为=+22/+19.孔输入如下命令：num=6.&61.2,95.2；den=l,7.5,22,19.5,0；G=tf(num,den)；Gl=zpk(G)；Z=tzero(G)；P= pole(G)；结果为：Z= -7-2P=01/1-3.0000 - 2.0000i -1.50003.1.6零极点分布图在MATLAB中，可利用pzmap()函数绘制连续系统的零、极点图,从而分析系统的稳定性，该函数调

5、用格式为：pzmap(num,den)o妇+ 2+W+4s + 6例二十五 给定传递函数：+ 2利用下列命令可自动打开一个图形窗口，显示该系统的零、极点分布 图。用鼠标点击图中零、极点可自动显示其坐标值。num=3,2,5,4,6dcn=l,3,4,2,7,2pzmap(num,den)title(4Pole-Zero Map5) % 图形标题ILLX.o0X.rrr-1-0.5Real Axis0R)le-Zero Map5 o.o 5osxv 一 BeuJ-20.53.2根轨迹法控制系统的稳定性，由其闭环极点唯一确定，而系统过渡过程的基 本特性，则与闭环零极点在s平面的位置有关。根轨迹法就

6、是在已知 控制系统开环传递函数零极点分部的基础上，研究某些参数变化时控 制系统闭环传递函数零极点分布影响的一种图解方法。利用根轨迹法, 能够分析系统的瞬态响应特性以及参数变化对瞬态响应特性的影响。 也可以根据对瞬态响应的要求去确定可变参数或调整零极点的位置和 个数。因此，根轨迹法可以用于解决线性系统的分析和综合问题。3.2.1求系统根轨迹rlocus命令可求得系统的根轨迹格式：r, k = rlocus (num, den)r, k = rlocus (num, den, k)不带输出变量时则绘出系统的根轨迹图，带输出变量时给出一组r, k的对应数据。若给定了 k的取值范围，则该命令将按要求绘

7、出图形或数组或者输 出指定增益k所对应的r值。每条根轨迹都以不同的颜色区别。 例二十八某系统开环传递函数为：小、2kGG)= ? +要绘制系统的根轨迹，则输入：n = 2d=13 2 0rlocus(n ,d)执行后得到下面图形。Ftoot Locus-4-6-5-4-321Real Axis若要得到指定增益k值对应的r值则输入：n = 2d=13 2 0r, k = rlocus (n, d, 5)结果如下：r= -3.30890.1545 + 1.7316i0.1545 - 1.7316ik= 53.2.2求根轨迹增益rlocfind命令可求得给定根的根轨迹增益。格式：k, poles =

8、 rlocfind (n, d)k, poles = rlocfind (num, den, p)当代有输出变量时，可得到所有极点的座标数据和增益值。不带输 出时只得到所选点的座标和增益值。注意：在执行这条命令前最好先执行一次根轨迹的绘图命令，这样 就可直接在根轨迹图上选取我们感兴趣的点。其中的p是系统的根，由此可得到对应的增益值。3.2.3 绘制和wn格sgrid命令是在图形中绘制出阻尼系数和自然频率栅格，其阻尼系数从01,步长为0.1。 命令格式： sgrid sgrid (z, wn)例二十七绘制系统带栅格的根轨迹图5 + 1+ 2$ + 3则执行:n=11%加入栅格d = 12rloc

9、us(n,d) sgrid当该命令带有指定的z (g), w ()时，则将按指定的参数绘制 有关图形。例二十八 在上题中绘制g=0.8,曾=2的根轨迹图执行：n = 1, 1 ；d = 1, 2, 3；rlocus (n, d)z = 0.8；w = 2；%加入指定的栅格sgrid (z, w)得到所需图形。Ftoot Locus1.5sxv Ajeu 一 Beta-1.51cccc1i1435325-2-1.51050Real Axis3.3 频域法频域分析法是利用频率特性研究控制系统的一种方法。频率特性 是指系统或环节在正弦信号作用下，稳态输出与输入之比对于频率的 关系。在控制系统的频域分

10、析法中常用到的坐标系统是极坐标系和对数 坐标系。在分析方法中常用的有三种：Bode图、Nyquist曲线和Nichols 图。3.3.1波特图法bode命令可获得连续系统的波特图或有关数据组。命令格式：mag, phase, w = bode (num, den)mag, phase, w = bode (num, den, w)当不带输出变量时则直接绘出图形。而带有输出变量时则得到一组 相关数据。其中的W是频率的取值范围，若缺省则该项由函数自动确绘图时的横坐标是以对数分度的。为了指定频率的范围，可采用 以下命令格式：Iogspace(dl,d2)或 ogspace(dl,d2,n)式是在指定

11、频率范围内按对数距离分成50等份的，即在两个十 进制数wl=10dl和w2=10d2之间产生一个由50个点组成的分量，矢量 中的点数50是一个默认值。_例如要在wl =0.1rad/s与w2= 100rad/s之间的频区画伯德图，则 输入命令时，dl=loglO(wl),d2=loglO(w2),在此频区按对数距离等分成 50个频率点，返回到工作空间中，即：w=logspace(-l,2)要对计算点数进行人工设定，则采用公式。例如要在wl = l与w2=1000之间产生100个对数等分点，可输入 以下命令：w=logspace(0,3,100)利用波特图我们可以分析系统的幅、相裕度、带宽、稳定

12、性、扰动 抑制能力等问题。利用波特图分析系统稳定性的方法如下： a：相角裕量0,幅值裕量k0,系统是稳定的。b： r = 0； k = 0,系统为临界稳定。c： r0； k0, k0,因此闭环后系统是稳定的。3.3.2求增益和相位裕度从前面的例题中可以看出要求系统的相位和增益裕度的准确值直 接调用bode命令是不太容易的。而使用margin命令则可以较容易的得 到所需值。margin可求出开环系统的幅值裕度和相角裕度，其格式为：margin (num, den)gm, pm, wcg, wcp margin (mag, phase, w)margin (num, den)可计算系统的相角裕度和

13、幅值裕度，并绘制出 Bode 图。margin (mag, phase, w)可以由幅值裕度和相角裕度绘制岀Bode图, 其中，mag phase和w是由bode得到的幅值裕度、相角裕度和频率。当带输出变量引用函数时，仅计算幅值裕度、相角裕度及幅值穿 越频率wcg和相角穿越频率wcp,不绘制Bode图。例三十：求例二十九中系统开环传递函数的相对稳定裕度。执行：n = 11 ；d = 4 3 2 0；margin （n, d）执行后得到相应图形和有关数据k= 15.6db r = 46.46040200-204022510 2Bode DigramGm= 15.6 dB (at 1.41 rad

14、.sec), Rn = 46.4 deg (at 0.615 rad/sec)10 110Frequency (rad/sec)101sp) apnpuees3.3.3奈奎斯特法奈奎斯特法是利极坐标图对系统进行分析的一种方法。频率特性的 极坐标图是当w由零变化到无穷大时，表示在极坐标上的幅与相角的 关系图，采用极坐标图，可以在一张图纸上描绘出整个频域的频率响 应。奈奎斯特稳定数据如下：（1）开环系统稳定时，如果曲线不包围（-1,川）点，则闭环系统 是稳定的，否则为不稳定的。（2）开环系统不稳定时，如果曲线反时针方向环绕（-1,丿0）点的 次数N等于右半平而内的极点数p,那么闭环系统是稳定的，否

15、则是 不稳定的。nyquist命令可以求得连续系统的奈奎斯特曲线。命令格式：re, im, w = nyquist (num, den)re, im, w = nyquist (num, den, w)当带有输出变量时，可得到相应的一组数据，不带输出变量时，贝9 绘出奈奎斯特曲线。也可用指定频率向量w指定所要绘制的曲线范围。 例三十一：系统开环传递函数为G(s)H($)=s + 2s +s* + 05绘出系统的奈奎斯特图并判断系统闭环后是否稳定。输入：n = 1 ；d = 1210.5 ；nyquist (n, d)执行后可得到所需图形。-31.5lL-厂-1 一一一 、III-1V111rr

16、11 1rrr00.5Real Axis1Nyquist Diagram22sixv从图形上可以看出由于曲线不包围(-1,)0 )点，因此闭环后系 统是稳定的。在某些场合我们需要在奈奎斯特曲线上加上单位圆帮助我们了解 相位，幅值裕量的粗值。下面的方法可以使我们获得带单位圆的奈奎 斯特图，同时还可获得不同增益下的奈奎斯特图。例三十二：系统开环传递函数为（5+ 2）（?+25+ 5）绘制k= 10, 26, 50时的带单位圆的奈奎斯特图，并估算系统的增僧 裕量。单位圆的绘制是通过绘沪的实部与虚部的轨迹而获得。输入：n = 10od = conv(l 2,1 2 5)。 w = 0:0.01:10o

17、e = exp(j *w)。r = real(e)oi = imag(e)oa,b = nyquist(n,d,w)o 果nl = 26odl = dod2 = do%取k= 10时的值%确定频率范I韦1%给岀指数函数ejw%求指数函数的实部，结果不显示%求函数cjw的虚部，结果不显示 %求指定频率范围内的奈氏值，不显示结%取 k = 26%保留原分母矢量al,bl = nyquist(nl,dl,w)o n2 = 50o%求k = 0.5时的奈氏值，结果不显示% k = 50%分母保留a2,b2 = nyquist(n2,d2,w)oplot(rj,a9b,a 1 ,b 1 ,a2,b2),

18、grid %绘出：r, i； a, b； bi； a?, b?； 的对应图形并加上栅格。执行以上程序后可在上图上得k = 10, 26, 50并加有单位圆的奈奎斯特图。一般来说由于此列的关系显示的图形不是一个正规的圆。从图形上我们可以看出开环增益K对闭环系统稳定性的影响。当K值变 化时，幅频特性成比例变化，而相频特性不受影响。因此取K = 26时,曲线恰好通过（-1,）0）点，这是临界稳定状态；当Kv26时，幅相曲 线G（将从（-1,）0 ）点的右方穿过负实轴，不再包围（-1,丿0 ）点, 这时闭环系统是稳定的；而K26时，开环频率特性5问随着。从0变 化到+ 00时，顺时针方向围绕（-1打0）点一圈，即N = 一1,可求得闭 环系统在右半s平而的极点数为：Z = P-2N =0-2x（-1） = 2 ,所以