系统开环频率特性的绘制
- 格式:doc
- 大小:504.50 KB
- 文档页数:9
《自动控制原理》实践报告实验三系统频率特性曲线的绘制及系统分析熟悉利用计算机绘制系统伯德图、乃奎斯特曲线的方法,并利用所绘制图形分析系统性能。
一、实验目的1.熟练掌握使用MATLAB软件绘制Bode图及Nyquist曲线的方法;2.进一步加深对Bode图及Nyquist曲线的了解;3.利用所绘制Bode图及Nyquist曲线分析系统性能。
二、主要实验设备及仪器实验设备:每人一台计算机奔腾系列以上计算机,配置硬盘≥2G,内存≥64M。
实验软件:WINDOWS操作系统(WINDOWS XP 或WINDOWS 2000),并安装MATLAB 语言编程环境。
三、实验内容已知系统开环传递函数分别为如下形式, (1))2)(5(50)(++=s s s G (2))15)(5(250)(++=s s s s G(3)210()(21)s G s s s s +=++ (4))12.0)(12(8)(++=s s s s G (5)23221()0.21s s G s s s s ++=+++ (6))]105.0)125.0)[(12()15.0(4)(2++++=s s s s s s G 1.绘制其Nyquist 曲线和Bode 图,记录或拷贝所绘制系统的各种图形; 1、 程序代码: num=[50];den=conv([1 5],[1 2]); bode(num,den)num=[50];den=conv([1 5],[1 2]); nyquist(num,den)-80-60-40-20020M a g n i t u d e (d B)10-210-110101102103-180-135-90-450P h a s e (d e g )Bode DiagramFrequency (rad/sec)-1012345-4-3-2-11234Nyquist DiagramReal AxisI m a g i n a r y A x i s2、 程序代码: num=[250];den=conv(conv([1 0],[1 5]),[1 15]); bode(num,den)num=[250];den=conv(conv([1 0],[1 5]),[1 15]);-150-100-5050M a g n i t u d e (d B )10-110101102103-270-225-180-135-90P h a s e (d e g )Bode DiagramFrequency (rad/sec)nyquist(num,den)3、 程序代码: num=[1 10];den=conv([1 0],[2 1 1]); bode(num,den)-150-100-50050100M a g n i t u d e (d B)10-210-110101102103-270-225-180-135-90P h a s e (d e g )Bode DiagramFrequency (rad/sec)-1-0.9-0.8-0.7-0.6-0.5-0.4-0.3-0.2-0.10-15-10-551015System: sys Real: -0.132Imag: -0.0124Frequency (rad/sec): -10.3Nyquist DiagramReal AxisI m a g i n a r y A x i snum=[1 10];den=conv([1 0],[2 1 1]); nyquist(num,den)-25-20-15-10-5-200-150-100-5050100150200Nyquist DiagramReal AxisI m a g i n a r y A x i s-100-5050100M a g n i t u d e (d B )10-210-110101102-270-225-180-135-90P h a s e (d e g )Bode DiagramFrequency (rad/sec)4、 程序代码: num=[8];den=conv(conv([1 0],[2 1]),[0.2 1]); bode(num,den)-18-16-14-12-10-8-6-4-20-250-200-150-100-50050100150200250Nyquist DiagramReal AxisI m a g i n a r y A x i snum=[8];den=conv(conv([1 0],[2 1]),[0.2 1]); nyquist(num,den)5、 程序代码: num=[1 2 1]; den=[1 0.2 1 1]; bode(num,den)num=[1 2 1];den=[1 0.2 1 1]; nyquist(num,den)-40-30-20-10010M a g n i t u d e (d B )10-210-110101102-360-270-180-90P h a s e (d e g )Bode DiagramFrequency (rad/sec)-2.5-2-1.5-1-0.500.51 1.5-3-2-1123Nyquist DiagramReal AxisI m a g i n a r y A x i s-100-5050100M a g n i t u d e (d B )10-210-110101102-270-225-180-135-90P h a s e (d e g )Bode DiagramFrequency (rad/sec)6、 num=[2 4];den=conv(conv([1 0],[2 1]),[0.015625 0.05 1]); bode(num,den)num=[2 4];den=conv(conv([1 0],[2 1]),[0.015625 0.05 1]); nyquist(num,den)2.利用所绘制出的Nyquist 曲线及Bode 图对系统的性能进行分析:(1)利用以上任意一种方法绘制的图形判断系统的稳定性; 由Nyquist 曲线判断系统的稳定性,Z=P-2N 。
绘制开环对数幅频特性曲线
在控制系统的分析与设计中,开环传递函数是一个重要的概念。
开环传递函数可以反映系统的频率响应特性,为系统的调节与控制提供重要的依据。
本文将以某系统的开环传递函数为例,介绍如何绘制开环对数幅频特性曲线。
开环传递函数是指系统的输出与输入之间的函数关系。
在频率域中,它可以表示为H(jω),其中ω 是角频率。
对数幅频特性曲线是将开环传递函数的幅值以对数形式绘制的曲线,表示系统在不同频率下的幅值响应。
绘制对数幅频特性曲线的步骤如下:
计算开环传递函数的幅值:H(jω) = |H(jω)|
设定横坐标的范围和刻度,并按照要求选择合适的频率单位。
将开环传递函数的幅值用对数形式表示,即y = log |H(jω)|。
绘制曲线,并观察系统在不同频率下的幅值响应。
绘制完成后,我们就可以通过观察对数幅频特性曲线来了解系统的频率响应特性。
例如,当频率单位为赫兹时,如果对数幅频特性曲线呈现出典型的高通特性,即随着频率的增大,幅值也随之增大,那么这意味着系统在高频时更加灵敏。
如果对数幅频特性曲线呈现出典
型的低通特性,即随着频率的增大,幅值减小,那么这意味着系统在低频时更加稳定。
绘制开环对数幅频特性曲线是了解系统频率响应特性的重要方法。
通过观察对数幅频特性曲线,我们可以了解系统在不同频率下的幅值响应,为系统的调节与控制提供重要的依据。
5.3 系统开环频率特性的绘制对自动控制系统进行频域分析时,通常是根据开环系统的频率特性来判断闭环系统的稳定性和估算闭环系统时域响应的各项性能指标,或者根据开环系统的频率特性绘制闭环系统的频率特性,然后再分析及估算时域性能指标。
因此,掌握开环系统的频率特性曲线的绘制和特点是十分重要的。
5.3.1 开环幅相曲线的绘制开环系统的幅相频率特性曲线简称为开环幅相曲线。
准确的开环幅相曲线可以根据系统的开环幅频特性和相频特性的表达式,用解析计算法绘制。
显然,这种方法比较麻烦。
在一般情况下,只需要绘制概略开环幅相曲线,概略开环幅相曲线的绘制方法比较简单,但是概略曲线应保持准确曲线的重要特征,并且在要研究的点附近有足够的准确性。
下面首先介绍幅相频率特性曲线的一般规律与特点,然后举例说明概略绘制开环幅相曲线的方法。
设系统开环传递函数的一般形式为)1()1()()(11++=∏∏-==s T s s K s H s G j vn j v mi i τ )(m n ≥ (5-49)式中,K 为开环增益;v 为系统中积分环节的个数。
则系统的开环频率特性为)1()()1()()(11∏∏-==++=v n j jvmi i Tj j j K j H j G ωωωτωω (5-50)1.开环幅相曲线的起点在低频段当0→ω时,由式(5-50)可得 )90(0lim)(lim)()(lim ︒⋅-→→→==v j vve Kj K j H j G ωωωωωωω (5-51)由式(5-51)可知,当0→ω时,开环幅相曲线的起点取决于开环传递函数中积分环节的个数v 和开环增益K ,参见图5-23(a )。
0型(v =0)系统,开环幅相曲线起始于实轴上的)0,(j K 点。
Ⅰ型(v =1)系统,开环幅相曲线起始于相角为︒-90的无穷远处。
当+→0ω时,曲线渐近于与虚轴的平行的直线,其横坐标[])()(Re lim 0ωωωj H j G V x +→= (5-52)Ⅱ型(v =2)系统,开环幅相曲线起始于相角为︒-180的无穷远处。
当+→0ω时,曲线渐近于负实轴。
2.开环幅相曲线的终点在高频段∞→ω时,由于系统一般有m n ≥,故开环幅相曲线总是以顺时针方向趋于∞=ω点。
由式(5-50)可得︒--∞→=90)(0)()(lim m n j ej H j G ωωω (5-53)即开环幅相曲线以︒⨯--90)(m n 方向终止于坐标原点,如图5-23(b )所示。
3. 开环幅相曲线与实轴的交点开环幅相曲线与实轴的交点频率x ω可由下式求出,即令式(5-50)的虚部为零 []0)()(Im =ωωj H j G (5-54) 将求出的交点频率x ω代入式(5-50)的实部,即[])()(Re x x j H j G ωω (5-55) 由式(5-55)可计算出开环幅相曲线与实轴的交点坐标值。
4. 开环幅相曲线的变化范围(象限、单调性)在式(5-50)中,如果),,2,1(0m i i Λ==τ,即不存在一阶微分环节时,则当ω从∞→0变化过程中,开环幅相曲线的相角将单调减小,曲线平滑地变化;若式(5-50)中有一阶微分环节,则视这些环节时间常数的数值大小不同,开环幅相曲线的相角可能不是以同一方向单调地变化,这时曲线上将会出现凹凸现象。
下面举例说明概略开环幅相曲线的绘制。
图5-23 不同类型系统的幅相频率特性例5.1已知某单位反馈系统,其开环传递函数为 )1)(1()(21++=s T s T s Ks G ,试绘制概略开环幅相曲线。
解 系统的开环频率特性 )1)(1()(21++=T j T j j Kj G ωωωω422212222122121)(1)1)(()(ωωωωT T T T T T Kj T T K +++--+-=①曲线的起点:该系统是Ⅰ型系统,由式(5-51)可知,系统的幅相曲线起始于相角为︒-90的无穷远处。
②曲线的终点:该系统的3=n 、0=m ,由式(5-53)可知,系统的幅相曲线以︒-=︒--27090)03(方向终止于坐标原点。
③曲线的变化范围:该系统不存在一阶微分环节,因此,系统幅相曲线的相角将由︒-90单调减小到︒-270,曲线平滑地变化。
④低频渐近线:低频渐近线坐标为 [])(Re lim 0ωωj G V x +→=4222122221210)(1)(lim ωωωT T T T T T K ++++-=+→)(21T T K +-=⑤曲线与实轴的交点:令[]0)(Im =ωj G ,可求出开环幅相曲线与实轴交点处的频率为211T T x =ω将频率x ω代入[])(Re ωj G ,可得出开环幅相曲线与实轴交点的坐标为[]2121)(Re T T T KT j G x +-=ω系统概略开环幅相曲线如图5-24所示。
例5.2 已知系统的开环传递函数为 )1()1()()(122++=s T s s T K s H s G ,试绘制概略开环幅相曲线。
解 系统的开环频率特性图5-24 例5.1开环幅相曲线)1()1()()(122+-+=T j T j K j H j G ωωωωω)arctan 180(arctan 21222121)(1)(ωωωωωT T j e T T K -︒-++=该系统是Ⅱ型(v =2)系统。
开环幅相曲线的起点:∞=+)0(A ,︒-=+180)0(ϕ 开环幅相曲线的终点:0)(=∞A ,︒-=∞180)(ϕ系统的开环幅相曲线的形状视时间常数1T 和2T 的数值大小不同而不同,下面讨论两种典型情况。
①12T T >:由于12T T >,因此∞<<+ω0时,12arctan arctan T T ωω>,开环幅相曲线位于第三象限,见图5-25(a )所示。
②12T T <:由于12T T <,因此∞<<+ω0时,12arctan arctan T T ωω<,开环幅相曲线位于第二象限,见图5-25(b )所示。
5.3.2 开环对数频率特性曲线的绘制开环系统的对数频率特性曲线简称为开环对数频率特性曲线。
根据典型环节的对数频率特性曲线,能够方便地绘制出开环对数频率特性曲线。
设系统的开环传递函数由n 个典型环节串联组成,即 ∏===ni in s G s G s G s G s H s G 121)()()()()()(Λ系统的开环频率特性为图5-25 例5.2的开环幅相曲线)()(11)()()()()(1ωϕωϕωωωωωj jn i in i ie A eA j G j H j G ni i =∑===∏∏==故系统的开环对数幅频特性和开环对数相频特性分别为∑===ni iA A L 1)(lg 20)(lg 20)(ωωω (5-56)∑==ni i 1)()(ωϕωϕ (5-57)由式(5-56)和(5-57)表明,若系统开环传递函数由n 个典型环节串联组成,其对数幅频特性曲线和对数相频特性曲线可由各典型环节的对数频率特性曲线叠加而得。
例 5.3 已知单位反馈系统的开环传递函数 )1()(+=Ts s Ks G ,试绘制系统的开环对数频率特性曲线。
解 系统的开环传递函数是由三个典型环节的组成:比例环节K 、积分环节s1和惯性环节11+Ts 。
分别做出各典型环节对数频率特性曲线,如图5-26所示。
图5-26中1L 、2L 、3L 分别为比例环节、积分环节和惯性环节的对数幅频特性曲线;1ϕ、2ϕ、3ϕ分别为比例环节、积分环节和惯性环节的对数相频特性曲线。
将各典型环节的对数幅频特性曲线叠加,即得系统开环对数幅频特性曲线,见图5-26中的)(ωL 。
在交接频率附近加以修正可得到精确曲线。
将各典型环节的对数相频特性曲线叠加,即得系统开环对数相频特性曲线,见图5-26中的)(ωϕ。
鉴于系统开环对数幅频近似特性曲线在控制系统的分析与设计中具有十分重要的作用,以下着重介绍开环对数幅频近似特性曲线的绘制方法。
图5-26 例5.3伯德图分析图5-26中系统的对数幅频近似特性曲线,可以看出开环对数幅频特性曲线有如下特点: (1) 低频段的斜率为v 20-dB/dec ,其中v 为开环系统中所包含的串联积分环节的个数。
(2) 在1=ω时,低频段或其延长线(当1<ω的频率范围内有交接频率时)的分贝值是K lg 20,如图5-27所示。
低频段或其延长线与零分贝线的交点频率为v K =0ω。
(3)在典型环节的交接频率处,对数幅频近似特性曲线的斜率要发生变化,变化的情况取决于典型环节的类型。
若遇到1)1()(±+=Ts s G 的环节,在交接频率处斜率改变20±dB/dec ;当遇到122)12()(±++=Ts s T s G ζ的环节时,在交接频率处斜率改变40± dB/dec 。
了解了以上特点,可以根据开环传递函数直接绘制近似对数幅频特性曲线。
绘制对数幅频特性曲线的步骤归纳如下:(1)将系统的开环传递函数化成典型环节串联组成的标准形式; (2)根据系统的开环增益K ,计算K lg 20的分贝值;(3)在1=ω处,标出K L lg 20)1(=点,过(K lg 20,1)点绘制斜率为v 20-dB/dec 的低频段;(4)根据交接频率绘制出相应的线段:(5)若有必要,可以利用误差修正曲线,对交接频率附近的曲线进行修正,则可以得到精确的特性曲线。
例5.4 已知某系统的开环传递函数为)252)(15()1(50)()(2++++=s s s s s s H s G ,试绘制其开环对数幅频特性曲线。
解 此系统是由一个比例环节、一个一阶微分环节、一个积分环节、一个惯性环节、一个振荡环节组成。
①先将开环传递函数为)()(s H s G 化成典型环节串联组成的标准形式,即图5-271=ω时,近似对数幅频特性曲线低频段或其延长线的纵坐标为K lg 20)152.025)(15()1(2)()(22+⨯+++=s s s s s s H s G②由2=K ,可得02.6lg 20=K dB 。
③在图中1=ω处标出02.6)1(=L dB ,过(6.02,1)点画一条斜率为20-dB/dec 的直线,它就是低频段的渐近线。
④在横坐标上标出各典型环节的交接频率,即惯性环节的交接频率2.01=ω,一阶微分环节的交接频率12=ω,振荡环节的交接频率53=ω。
在各交接频率处依次改变斜率,直接绘制开环对数幅频特性曲线的渐近线。
在2.0=ω处,曲线斜率由20-dB/dec 变为40-dB/dec ;在1=ω处,曲线斜率由40-dB/dec 变为20-dB/dec ;在5=ω处,曲线斜率由20-dB/dec 变为60-dB/dec 。