用MATLAB 分析闭环系统的频率特性 1、等M 圆图与等N 圆图原理 1.1设有单位系统如图1示。其闭环频率特性G B (j )与开环频率特性G K (j )的关系为 )(j G 1)(j G )(j X )(j X )(j G K K i 0B ωωωωω+== (1) 图 1 可将其开环频率特性G K (j )写成 G K (j )=U ()+jV() (2) 则闭环频率特性为 )(j B )e M(jV U 1jV U )G (j 1)G (j )(j G ωαωωωω=+++=+= (3) 式中 M()——闭环的幅频特性 ()——闭环的相频特性 闭环的幅频特性为 2 12222V )U (1V U |jV U 1||jV U |M ??????++++++= (4) 所以 222 22 V U)(1V U M +++= (5) 则有 2 22 2222 1)-(M M V )1-M M (U =++ (6) 显然,式(6)是一个元的方程,他表明了开环的实频U 、虚频V 和闭环的幅频M 之间 G K (j ) X i (j ) X 0(j )
的的关系,该圆方程的圆心坐标为(1M M 22--,j0),半径为|1-M M |2。当M 取不同的值时,便可以得到一簇圆,如图1,该图称为等M 圆图(邮称为等幅值轨迹图)。 有闭环的相频特性为 )V U U V (tg )U 1V (tg )U V (tg )jV U 1jV U (221-1-1-++=+=+++∠=-α (7) 令22V U U V tg N ++==α,上式可改为 22224N 1N )2N 1(V )21(U +=+++ (8) 可见式(8)也是一个圆方程,他表明了U 、V 与N 之间的关系。该圆方程的圆心坐标为 |。-,半径为|-1N )2N 1j ,21(2当N 取不同的值时,可画出一簇圆,如图2所示。该 方法复杂,也不准确,我们用一个具体的力来说明一下用MATLAB 解决这类问题的方
实验三 典型环节频率特性分析 一.实验目的 1. 学习频率特性分析仪的使用; 2. 掌握频率特性测试方法; 3. 掌握由对象频率特性求传递函数的方法。 二.实验设备及简介 1. 实验设备 TD4011A 频率特性分析仪,微计算机,打印机。 2. TD4011A 频率特性分析仪简介 数字键区 信号发生器输出 图2 TD4011A 频率特性分析仪面板图 TD4011A 分析仪如图1所示,由信号发生器和分析器组成。其面板图如图2所示。主要按键功能: ⑴.上档键 — DELAY — 延迟时间。分0.1s 、1s 、10s 三档。每按一次,循环改变一次。 CYCLE — 积分周数。分 1、10、100、1000三档。每按一次,循环改变一次。积分周数大精度高。 AMPL — 信号发生器输出电压值。 FREQ — 信号发生器输出频率值。 F MAX — 扫频(即频率按顺序变化)频率上限。 F MIN — 扫频频率下限。 D LOG — 对数扫频增量(每倍频程扫频步数) D LIN — 线性扫频增量(单位:Hz ) PROGRAM — 前后面板输入选择。0为前面板输入,1为后面板输入。用数字键区 ※ 以上功能设定,均由 图1 TD4011A 频率特性分析仪
⑵.下档键— 下档功能中AUTO、30mV、300mV、3V、30V、300V为输入量程选择; ; ※下档键功能均为灯亮有效。 ⑶.中档键— RECYCLE —发生器输出连续扫频信号; SINGLE —发生器输出单步扫频信号; STOP —测量停止。只有此键灯亮时才能对面板状态进行设定; HOLD —将发生器信号保持在扫频范围内的某一频率上; LOG↑—对数上扫(即发生器信号频率按对数规律由F MIN至F MAX变化); LIN↑—线性上扫;LOG↓—对数下扫;LIN↓—线性下扫; OFF —关断扫频; LOCAL —与计算机进行通讯; PRINT —打印,实验中此功能不用;PROGRAM —信号源停止时的相位设置。 ※按键位于中档键标识处,特别提醒上档和下档功能设定时,均要按在中档位置。 ⑷.数字键区 Select —选择键; RESET —复位键,强迫进入初始状态; CLEAR —发生器显示清零,不清内存; ENTER —确认键。 三.实验内容及步骤 1. 已知环节1,测试其频率特性。 ⑴实验前准备内容 ①环节1网络如图3所示,K R200 1 =,K R2 2 =,F Cμ1 =。求 出环节传递函数及频率特性。 ②绘制环节的对数幅频渐近线,并进行修正;绘制其对数相频特性; 由对数频率特性确定实验频率f的最佳范围(注意:π ω2 = f仪器允许频率范围0.001 ~1000Hz)。 ⑵实验步骤 ①打开计算机,进入WIN98操作系统,将软件狗插入USB接口。
实验四 线性系统的频率特性 一、实验目的: 1. 测量线性系统的幅频特性 2. 复习巩固周期信号的频谱测量 二、实验原理: 我们讨论的确定性输入信号作用下的集总参数线性非时变系统,又简称线性系统。线性系统的基本特性是齐次性与叠加性、时不变性、微分性以及因果性。对线性系统的分析,系统的数学模型的求解,可分为时间域方法和变换域方法。这里主要讨论以频率特性为主要研究对象,通过傅里叶变换以频率为独立变量。 设输入信号)(t v in ,其频谱)(ωj V in ;系统的单位冲激响应)(t h ,系统的频率特性 )(ωj H ;输出信号)(t v out ,其频谱)(ωj V out ,则 时间域中输入与输出的关系 )()()(t h t v t v in out *= 频率域中输入与输出的关系 )()()(ωωωj H j V j V in out ?= 时间域方法和变换域方法并没有本质区别,两种方法都是将输入信号分解为某种基本单元,在这些基本单元的作用下求得系统的响应,然后再叠加。变换域方法可以将时域分析中的微分、积分运算转化为代数运算,将卷积积分变换为乘法;在信号处理时,将输入时间信号用一组变换系数(谱线)来表示,根据信号占有的频带与系统通带间的关系来分析信号传输,判别信号中带有特征性的分量,比时域法简便和直观。 三、实验方法: 1. 输入信号的选取 这里输入信号选取周期矩形信号,并且要求 τ T 不为整数。这是因为周期矩形信号具有丰富的谐波分量,通过观察系统的输入、输出波形的谐波的变化,分析系统滤波特性。周期矩形信号可以分解为直流分量和许多谐波分量;由于测量频率点的数目有限,因此需要排除谐波幅度为零的频率点,周期矩形信号谐波幅度为零的频率点是 Ω KT ,其中1=K 、2、3、… 。 图11.1 输入的周期矩形信号时域波形 t
第三章 系统频率特性 系统的时域分析是分析系统的直接方法,比较直观,但离开计算机仿真,分析高阶系统是困难的。系统频域分析是工程广为应用的系统分析和综合的间接方法。频率分析不仅可以了解系统频率特性,如截止频率、谐振频率等,而且可以间接了解系统时域特性,如快速性,稳定性等,为分析和设计系统提供更简便更可靠的方法。 本章首先阐明频率响应的特点,给出计算频率响应的方法,接着介绍Nyquist 图和Bode 图的绘制方法、系统的稳定裕度及系统时域性能指标计算。 3.1 频率响应和频率特性 3.1.1 一般概念 频率响应是指系统对正弦输入的稳态响应。考虑传递函数为G(s)的线性系统,若输入正弦信号 t X t x i i ωsin )(= (3.1-1) 根据微分方程解的理论,系统的稳态输出仍然为与输入信号同频率的正弦信号,只是其幅值和相位发生了变化。输出幅值正比于输入的幅值i X ,而且是输入正弦频率ω的函数。输出的相位与i X 无关,只与输入信号产生一个相位差?,且也是输入信号频率ω的函数。即线性系统的稳态输出为 )](sin[)()(00ω?ωω+=t X t x (3.1-2)
由此可知,输出信号与输入信号的幅值比是ω的函数,称为系统的幅频特性,记为)(ωA 。输出信号与输入信号相位差也是ω的函数,称为系统的相频特性,记为)(ω?。 幅频特性: )()()(0ωωωi X X A = (3.1-3) 相频特性: )()()(0ω?ω?ω?i -= (3.1-4) 频率特性是指系统在正弦信号作用下,稳态输出与输入之比对频率的关系特性,可表示为: )()()(0ωωωj X j X j G i = (3.1-5) 频率特性)(ωj G 是传递函数)(s G 的一种特殊形式。任何线性连续时间系统的频率特性都可由系统传递函数中的s 以ωj 代替而求得。 )(ωj G 有三种表示方法: )()()(ω?ωωj e A j G = (3.1-6) )()()(ωωωjV U j G += (3.1-7) )(sin )()cos()()(ω?ωωωωjA A j G += (3.1-8) 式中,实频特性: )(cos )()(ω?ωωA U = 虚频特性:
开环系统频率特性曲线的绘制方法 (一) 已知系统开环传递函数G k (s ),绘制Nyquist 曲线(开环幅相曲线) 一、ω:0+→+∞ 1、由已知的G k (s )求()()k k s j G j G s ωω==,A (ω),φ(ω) ,P (ω),Q (ω); 11211222 1 1 2 2 1 2 1 1 2 2 1 2 1121 12221 1221 2 1 1 2 2 1 2 22222 2 2 2(1)[(1)2](1)[(1)2]()()(1)[(1)2](1)[(1)2] m m m m j k j k k k j k j k k k k v n n n n i l i l l l i l i l l l j T j j T j k G j j j T j j T j ωωωωωξωξωωωωωωωωωωωξωξωωω ω+-+---= +-+---∏∏∏∏∏∏∏∏ (1) 式中:分子多项式中最小相位环节的阶次和为111212m m m =+, 分子多项式中非最小相位环节的阶次和为212222m m m =+, 分母多项式中最小相位环节的阶次和为111212n n n v =++, 分母多项式中非最小相位环节的阶次和为212222n n n =+, 分子多项式阶次之和为12m m m =+,分母多项式阶次之和为12n n n =+。 注:式中仅包含教材p192所列5种非最小相位环节,不包含形如1Ts -、 11Ts -、2 2 121 n n s s ξωω+-、22 21n n s s ξωω+-等非最小相位环节。 2、求N 氏曲线的起点 当ω→0+时,(1)式可近似为: 0lim ()()k v k G j j ωωω+ →→ (2) 于是,N 氏曲线的起点取决于开环放大系数k 和系统的型v 。 ① 当0v =时,N 氏曲线起始于实轴上的一点(k ,0)或(-k ,0); ② 当0v >时,N 氏曲线起始于无穷远点: 0k >时,沿着角度()2 v π?ω=-?起始于无穷远点; 0k <时,沿着角度()2 v π?ωπ=--?起始于无穷远点。 ③ 当0v <时,N 氏曲线起始于原点: 0k >时,沿着角度()2 v π?ω=?起始于原点; 0k <时,沿着角度()2 v π?ωπ=-+?起始于原点。 3、求N 氏曲线的终点 当ω→+∞时,(1)式中各环节的相角分别为:
课程名称:_________控制理论(甲)实验_______指导老师:_____ ____成绩:__________________ 实验名称:___典型环节和系统频率特性的测量___实验类型:________________同组学生:__________ 一、实验目的 二、实验原理 三、实验接线图 四、实验设备 五、实验步骤 六、实验数据记录 七、实验数据分析 八、实验结果或结论 一、实验目的 1.了解典型环节和系统的频率特性曲线的测试方法; 2.根据实验求得的频率特性曲线求取传递函数。 二、实验原理 1.系统(环节)的频率特性 设G(S)为一最小相位系统(环节)的传递函数。如在它的输入端施加一幅值为X m 、频率为ω的正弦信号,则系统的稳态输出为 )sin()()sin(?ωω?ω+=+=t j G Xm t Y y m 由式①得出系统输出,输入信号的幅值比相位差 )() (ωωj G Xm j G Xm Xm Ym == (幅频特性) )()(ωωφj G ∠= (相频特性) 式中)(ωj G 和)(ωφ都是输入信号ω的函数。 2.频率特性的测试方法 2.1 沙育图形法测试 2.1.1幅频特性的测试 由于 m m m m X Y X Y j G 22)(= = ω 改变输入信号的频率,即可测出相应的幅值比,并计算 m m X Y A L 22log 20)(log 20)(==ωω (d B ) 其测试框图如下所示:
图5-1 幅频特性的测试图(沙育图形法) 注:示波器同一时刻只输入一个通道,即系统(环节)的输入或输出。 2.1.2相频特性的测试 图5-2 相频特性的测试图(沙育图形法) 令系统(环节)的输入信号为:t X t X m ωsin )(= (5-1) 则其输出为 )sin()(φω+=t Y t Y m (5-2) 对应的沙育图形如图5-2所示。若以t 为参变量,则)(t X 与)(t Y 所确定点的轨迹将在示波器的屏幕上形成一条封闭的曲线(通常为椭圆),当t=0时,0)0(=X 由式(5-2)得 )sin()0(φm Y Y = 于是有 m m Y Y Y Y 2) 0(2sin )0(sin )(1 1--==ωφ (5-3) 同理可得 m X X 2) 0(2sin )(1 -=ωφ (5-4) 其中: )0(2Y 为椭圆与Y 轴相交点间的长度; )0(2X 为椭圆与X 轴相交点间的长度。 式(5-3)、(5-4)适用于椭圆的长轴在一、三象限;当椭圆的长轴在二、四时相位φ的计算公式变为 m Y Y 2) 0(2sin 180)(1 0--=ωφ 或 m X X 2)0(2sin 180)(10--=ωφ
5.3 系统开环频率特性的绘制 对自动控制系统进行频域分析时,通常是根据开环系统的频率特性来判断闭环系统的稳定性和估算闭环系统时域响应的各项性能指标,或者根据开环系统的频率特性绘制闭环系统的频率特性,然后再分析及估算时域性能指标。因此,掌握开环系统的频率特性曲线的绘制和特点是十分重要的。 5.3.1 开环幅相曲线的绘制 开环系统的幅相频率特性曲线简称为开环幅相曲线。准确的开环幅相曲线可以根据系统的开环幅频特性和相频特性的表达式,用解析计算法绘制。显然,这种方法比较麻烦。在一般情况下,只需要绘制概略开环幅相曲线,概略开环幅相曲线的绘制方法比较简单,但是概略曲线应保持准确曲线的重要特征,并且在要研究的点附近有足够的准确性。 下面首先介绍幅相频率特性曲线的一般规律与特点,然后举例说明概略绘制开环幅相曲线的方法。 设系统开环传递函数的一般形式为 ) 1()1()()(11 ++= ∏∏-==s T s s K s H s G j v n j v m i i τ )(m n ≥ (5-49) 式中,K 为开环增益;v 为系统中积分环节的个数。 则系统的开环频率特性为 ) 1() ()1()()(1 1∏∏-==++= v n j j v m i i T j j j K j H j G ωωωτωω (5-50) 1.开环幅相曲线的起点 在低频段当0→ω时,由式(5-50)可得 )90(0 lim ) (lim )()(lim ??-→→→==v j v v e K j 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ω时,曲线渐近于与虚轴的平行的直线,其横坐标
典型环节(或系统)的频率特性测量 一·实验目的 1.学习和掌握测量典型环节(或系统)频率特性曲线的方法和技能。 2.学习根据实验所得频率特性曲线求取传递函数的方法。 二·实验要求 1.用实验方法完成一阶惯性环节的频率特性曲线测试。 2.用实验方法完成典型二阶系统开环频率特性曲线的测试。 3.根据测得的频率特性曲线求取各自的传递函数。 4.用软件仿真方法求取一阶惯性环节频率特性和典型二阶系统开环频率特性,并与实验所得结果比较。 三·实验原理 掌握改变正弦波信号幅值和频率的方法。利用实验箱上的模拟电路单元,参考本实验附录设计并连接“一阶惯性环节”模拟电路(如用U9+U8连成)或“两个一阶惯性环节串联”的模拟电路(如用U9+U11连成)。 四·实验所用仪器 PC微机(含实验系统上位机软件)、ACT-I实验箱、USB2.0通讯线 五·实验步骤和方法 1.用实验方法完成一阶惯性环节的频率特性曲线测试。 2.用实验方法完成典型二阶系统开环频率特性曲线的测试。 3.根据测得的频率特性曲线求取各自的传递函数。 4.用软件仿真方法求取一阶惯性环节频率特性和典型二阶系统开环频率特性,并与实验所得结果比较。 具体步骤: 1.熟悉实验箱上的信号源,掌握改变正弦波信号幅值和频率的方法。利用实验箱上的模拟电路单元,参考本实验附录设计并连接“一阶惯性环节”模拟电路(如用U9+U8连成)或“两个一阶惯性环节串联”的模拟电路(如用U9+U11连成)。 2.利用实验设备完成一阶惯性环节的频率特性曲线测试。 无上位机时,利用实验箱上的信号源单元U2所输出的正弦波信号作为环节输入,即连接箱上U2的“正弦波”与环节的输入端(例如对一阶惯性环节即图1.5.2的Ui)。然后用示波器观测该环节的输入与输出(例如对一阶惯性环节即测试图1.5.2的Ui和Uo)。注意调节U2的正弦波信号的“频率”电位器RP5与“幅值”电位器RP6,测取不同频率时环节输出的增益和相移(测相移可用“李沙育”图形),从而画出环节的频率特性。 有上位机时,必须在熟悉上位机界面操作的基础上,充分利用上位机提供的虚拟示波器与信号发生器功能。为了利用上位机提供的虚拟示波器与信号发生器功能,接线方式将不同于上述无上位机情况。仍以一阶惯性环节为例,此时将Ui连到实验箱 U3单元的O1或O2(D/A通道的输出端,这个是通过上位机选择其中的一路输出),将Uo连到实验箱 U3单元的I1(A/D通道的输入端),然后再将你选择的D/A输出通道测试信号O1(如果选择的是O1)连接到这组A/D输入的另一采集输入端I2,然后连接设备与上位机的USB通信线。接线完成,
实验六 控制系统的频率特性分析 1.已知系统传递函数为:1 2.01)(+=s s G ,要求: (1) 使用simulink 进行仿真,改变正弦输入信号的频率,用示波器观察输 出信号,记录不同频率下输出信号与输入信号的幅值比和相位差,即 可得到系统的幅相频率特性。 F=10时 输入: 输出:
F=50时 输入:输出: (2)使用Matlab函数bode()绘制系统的对数频率特性曲线(即bode图)。 提示:a)函数bode()用来绘制系统的bode图,调用格式为: bode(sys) 其中sys为系统开环传递函数模型。 参考程序: s=tf(‘s’); %用符号表示法表示s G=1/(0.2*s+1); %定义系统开环传递函数 bode(G) %绘制系统开环对数频率特性曲线(bode图)
实验七连续系统串联校正 一.实验目的 1.加深理解串联校正装置对系统动态性能的校正作用。 2. 对给定系统进行串联校正设计,并通过matlab实验检验设计的正确性。二.实验内容 1.串联超前校正 系统设计要求见课本例题6-3,要求设计合理的超前校正环节,并完成以下内容用matlab画出系统校正前后的阶跃相应,并记录系统校正前后的超调量及调节时间 num=10; 1)figure(1) 2)hold on
3)figure(1) 4)den1=[1 1 0]; 5)Gs1=tf(num,den1); 6)G1=feedback(Gs1,1,-1); 7)Step(G1) 8) 9)k=10; 10)figure(2) 11)GO=tf([10],[1,1,0]); 12)Gc=tf([0.456,1],[1,00114]); 13)G=series(G0,Gc); 14)G1=feedback(G,1); 15)step(G1);grid
【实验目的】 1. 掌握测量典型一阶系统和二阶系统的频率特性曲线的方法; 2. 掌握软件仿真求取一、二阶系统的开环频率特性的方法; 3. 学会用Nyquist 判据判定系统的稳定性。 【实验设备与软件】 1. labACT 实验台与虚拟示波器 2. MATLAB 软件 【实验原理】 1.系统的频率特性测试方法 对于现行定常系统,当输入端加入一个正弦信号)sin()(t X t X m ωω=时,其稳态输出是一个与输入信号频率相同,但幅值和相位都不同的正弦信号 )si n ()()si n ()(ψωωψω+=+=t j G X t Y s Y m m 。 幅频特性:m m X Y j G /)(=ω,即输入与输出信号的幅度比值,通常转换成 )(lg 20ωj G 形式。 相频特性:)(arg )(ωω?j G =,可以直接基于虚拟示波器读取,也可以用“李沙育图行”法得到。 可以将用Bode 图或Nyquist 图表示幅频特性和相频特。 在labACT 试验台采用的测试结构图如下: 被测定稳 定系统对于实验就是有源放大电路模拟的一、二阶稳定系统。 2.系统的频率测试硬件原理 1)正弦信号源的产生方法 频率特性测试时,一系列不同频率输入正弦信号可以通过下图示的原理产生。按
照某种频率不断变化的数字信号输入到DAC0832,转换成模拟信号,经一级运放将其转换为模拟电压信号,再经过一个运放就可以实现双极性电压输出。 根据数模转换原理,知 R V N V 8012 - = (1) 再根据反相加法器运算方法,得 R R R V N V N V R R V R R V 1281282282201210--=??? ??+-?-=???? ??+-= (2) 由表达式可以看出输出时双极性的:当N 大于128时,输出为正;反之则为负;当输入为128时,输出为0. 在labACT 实验箱上使用的参考电压时5V 的,内部程序可以产生频率范围是对一阶系统是0.5 H Z ~64H Z 、对二阶系统是0.5 H Z ~16 H Z 的信号,并由B2单元的OUT2输出。 2)被测对象输出信号的采样方法 对被测对象的输出信号夏阳,首先将其通过LM324与基准电压进行比较嵌位,再通过CD14538进行脉冲整形,一保证有足够的IRQ 采样时间,最后将信号送到处理器的IRQ6脚,向处理器申请中断,在中断中对模拟量V y 进行采样并模数转换,进而进行处理与计算幅值与相位。途中采用ADC089采集模拟量,以单极性方式使用,所以在出现振荡的情况下需要加入一个二极管,将V y 出现负值时将其直接拉倒0。
第三章 系统频率特性 系统的时域分析是分析系统的直接方法,比较直观,但离开计算机仿真,分析高阶系统是困难的。系统频域分析是工程广为应用的系统分析和综合的间接方法。频率分析不仅可以了解系统频率特性,如截止频率、谐振频率等,而且可以间接了解系统时域特性,如快速性,稳定性等,为分析和设计系统提供更简便更可靠的方法。 本章首先阐明频率响应的特点,给出计算频率响应的方法,接着介绍Nyquist 图和Bode 图的绘制方法、系统的稳定裕度及系统时域性能指标计算。 3.1 频率响应和频率特性 3.1.1 一般概念 频率响应是指系统对正弦输入的稳态响应。考虑传递函数为G(s)的线性系统,若输入正弦信号 t X t x i i ωsin )(= (3.1-1) 根据微分方程解的理论,系统的稳态输出仍然为与输入信号同频率的正弦信号,只是其幅值和相位发生了变化。输出幅值正比于输入的幅值i X ,而且是输入正弦频率ω的函数。输出的相位与i X 无关,只与输入信号产生一个相位差?,且也是输入信号频率ω的函数。即线性系统的稳态输出为 )](sin[)()(00ω?ωω+=t X t x (3.1-2)
由此可知,输出信号与输入信号的幅值比是ω的函数,称为系统的幅频特性,记为)(ωA 。输出信号与输入信号相位差也是ω的函数,称为系统的相频特性,记为)(ω?。 幅频特性: )()()(0ωωωi X X A = (3.1-3) 相频特性: )()()(0ω?ω?ω?i -= (3.1-4) 频率特性是指系统在正弦信号作用下,稳态输出与输入之比对频率的关系特性,可表示为: )()()(0ωωωj X j X j G i = (3.1-5) 频率特性)(ωj G 是传递函数)(s G 的一种特殊形式。任何线性连续时间系统的频率特性都可由系统传递函数中的s 以ωj 代替而求得。 )(ωj G 有三种表示方法: )()()(ω?ωωj e A j G = (3.1-6) )()()(ωωωjV U j G += (3.1-7) )(sin )()cos()()(ω?ωωωωjA A j G += (3.1-8) 式中,实频特性: )(cos )()(ω?ωωA U = 虚频特性:
实验 4 系统的频率特性分析 一、实验目的 (1)为学习和掌握利用MATLAB 绘制系统Nyquist 图和Bode 图的方法。 (2)为学习和掌握利用系统的频率特性分析系统的性能。 二、实验原理 系统的频率特性是一种图解方法,运用系统的开环频率特性曲线,分析闭环系统的性 能,如系统的稳态性能、暂态性能。常用的频率特性曲线有Nyquist 图和Bode 图。在MATLAB 中,提供了绘制Nyquist 图和Bode 图的专门函数。 1. Nyquist 图 nyquist 函数可以用于计算或绘制连续时间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 图数据的幅度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。
5.6 利用开环频率特性分析系统的性能 在频域中对系统进行分析、设计时,通常是以频域指标作为依据的,但是不如时域指标来得直接、准确。因此,须进一步探讨频域指标与时域指标之间的关系。考虑到对数频率特性在控制工程中应用的广泛性,本节将以Bode 图为基点,首先讨论开环对数幅频特性)(ωL 的形状与性能指标的关系,然后根据频域指标与时域指标的关系估算出系统的时域响应性能。 实际系统的开环对数幅频特性)(ωL 一般都符合如图5-49所示的特征:左端(频率较低的部分)高;右端(频率较高的部分)低。将)(ωL 人为地分为三个频段:低频段、中频段和高频段。低频段主要指第一个转折点以前的频段;中频段是指截止频率c ω附近的频段;高频段指频率远大于c ω的频段。这三个频段包含了闭环系统性能不同方面的信息,需要分别进行讨论。 需要指出,开环对数频率特性三频段的划分是相对的,各频段之间没有严格的界限。一般控制系统的频段范围在Hz 100~01.0之间。这里所述的“高频段”与无线电学科里的“超高频”、“甚高频”不是一个概念。 5.6.1 )(ωL 低频渐近线与系统稳态误差的关系 系统开环传递函数中含积分环节的数目(系统型别)确定了开环对数幅频特性低频渐近线的斜率,而低频渐近线的高度则取决于开环增益的大小。因此, )(ωL 低频段渐近线集中反映了系统跟踪控制信号的稳态精度信息。根据)(ωL 低 频段可以确定系统型别υ和开环增益K ,利用第3章中介绍的静态误差系数法可以确定系统在给定输入下的稳态误差。 图5-49 对数频率特性三频段的划分
5.6.2 )(ωL 中频段特性与系统动态性能的关系 开环对数幅频特性的中频段是指截止频率c ω附近的频段。设开环部分纯粹由积分环节构成,图5-50(a )所示的对数幅频特性对应一个积分环节,斜率为 dec dB /20-,相角 90)(-=ω?,因而相角裕度 90=γ;图5-50(b )的对数幅频特性对应两个积分环节,斜率为dec dB /40-,相角 180)(-=ω?,因而相角裕度 0=γ。 图5-50 )(ωL 中频段对稳定性的影响 一般情况下,系统开环对数幅频特性的斜率在整个频率范围内是变化的,故截止频率c ω处的相角裕度γ应由整个对数幅频特性中各段的斜率所共同确定。在 c ω处,)(ωL 曲线的斜率对相角裕度γ的影响最大,远离c ω的对数幅频特性,其斜率对γ的影响就很小。为了保证系统有满意的动态性能,希望)(ωL 曲线以dec dB /20-的斜率穿过dB 0线, 并保持较宽的频段。截止频率c ω和相角裕度γ是系统开环频域指标,主要由中频段决定,它与系统动态性能指标之间存在着密切关系,因而频域指标是表征系统动态性能的间接指标。 1 二阶系统 典型二阶系统的结构图可用图5-51表示。其中开环传递函数为 2 ()(01)(2) n n G s s s ωξξω=<<+ 相应的闭环传递函数为 2 22 2)(n n n s s s ωξωω++=Φ (1)γ和%σ的关系: 系统开环频率特性为 图5-51 典型二阶系统结构图
实验五 典型环节和系统频率特性的测量 一、实验目的 1. 了解典型环节和系统的频率特性曲线的测试方法; 2. 根据实验求得的频率特性曲线求取传递函数。 二、实验设备 同实验一。 三、实验内容 1. 惯性环节的频率特性测试; 2. 二阶系统频率特性测试; 3. 无源滞后—超前校正网络的频率特性测试; 4. 由实验测得的频率特性曲线,求取相应的传递函数; 5. 用软件仿真的方法,求取惯性环节和二阶系统的频率特性。 四、实验原理 1. 系统(环节)的频率特性 设G(S)为一最小相位系统(环节)的传递函数。如在它的输入端施加一幅值为Xm 、频率为ω的正弦信号,则系统的稳态输出为 )sin()()sin(?ωω?ω+=+=t j G Xm t Y y m ① 由式①得出系统输出,输入信号的幅值比和相位差 )()(ωωj G Xm j G Xm Xm Ym == (幅频特性) )()(ωφω=∠j G (相频特性) 式中)(ωj G 和)(ωφ都是输入信号ω的函数。 2. 频率特性的测试方法 2.1 李沙育图形法测试 2.1.1幅频特性的测试 由于 m m m m X Y X Y j G 22)(== ω 改变输入信号的频率,即可测出相应的幅值比,并计算 m m X Y A L 22log 20)(log 20)(==ωω (dB ) 其测试框图如下所示:
图5-1 幅频特性的测试图(李沙育图形法) 注:示波器同一时刻只输入一个通道,即系统(环节)的输入或输出。 2.1.2相频特性的测试 图5-2 相频特性的测试图(李沙育图形法) 令系统(环节)的输入信号为:t X t X m ωsin )(= (5.1) 则其输出为 )sin()(φω+=t Y t Y m (5.2) 对应的李沙育图形如图5-2所示。若以t 为参变量,则)(t X 与)(t Y 所确定点的轨迹将在示波器的屏幕上形成一条封闭的曲线(通常为椭圆),当t=0时,0)0(=X 由式(5.2)得 )sin()0(φm Y Y = 于是有 m m Y Y Y Y 2) 0(2sin )0(sin )(1 1--==ωφ (5.3) 同理可得 m X X 2) 0(2sin )(1 -=ωφ (5.4) 其中 )0(2Y 为椭圆与Y 轴相交点间的长度; )0(2X 为椭圆与X 轴相交点间的长度。 式(5.3)、(5.4)适用于椭圆的长轴在一、三象限;当椭圆的长轴在二、四时相位φ的计算公式变为 m Y Y 2) 0(2sin 180)(1 0--=ωφ 或 m X X 2)0(2sin 180)(10--=ωφ
实验四二阶开环系统的频率特性曲线
实验报告 课程名称控制工程基础 实验项目实验四二阶开环系统的频率特性曲线 专业电子科学与技术班级一 姓名学号 指导教师实验成绩 2014年5月29日
实验四 二阶开环系统的频率特性曲线 一、实验目的 1.研究表征系统稳定程度的相位裕度γ和幅值穿越频率c ω对系统的影响。 2. 了解和掌握二阶开环系统中对数幅频特性L(w )和相频特性)(ω?,实频特性Re (w )和虚频特性Im (w )的计算。 3.了解和掌握欠阻尼二阶开环系统中的相位裕度γ和幅值穿越频率c ω的计算。 4.观察和分析欠阻尼二阶开环系统波德图中的相位裕度γ和幅值穿越频率ωc ,与计算值作比对。 二、实验仪器 PC 机一台,实验箱 三、实验内容及操作步骤 本实验用于观察和分析二阶开环系统的频率特性曲线。 由于Ⅰ型系统含有一个积分环节,它在开环时响应曲线是发散的,因此欲获得其开环频率特性时,还是需构建成闭环系统,测试其闭环频率特性,然后通过公式换算,获得其开环频率特性。 自然频率:TiT K =n ω 阻尼比:KT Ti 2 1=ξ (3-2-1) 谐振频率:221ξωω-=n r 谐振峰值:2121 lg 20)(ξξω-=r L (3- 2-2) 计算欠阻尼二阶闭环系统中的幅值穿越频率ωc 、相位裕度γ: 幅值穿越频率: 24241ξξωω-+?=n c (3-2-3)
相位裕度: 4 24 1 2 2 arctan ) ( 180 ξ ξ ξ ω ? γ + + - = + = c (3-2-4)γ值越小,Mp%越大,振荡越厉害;γ值越大,Mp%小,调节时间ts越长,因此为使二阶闭环系统不致于振荡太厉害及调节时间太长,一般希望: 30°≤γ≤70°(3-2-5)本实验所构成的二阶系统符合式(3-2-5)要求。 被测系统模拟电路图的构成如图1所示。 图1 实验电路 本实验将数/模转换器(B2)单元作为信号发生器,自动产生的超低频正弦信号的频率从低到高变化(0.5Hz~16Hz),OUT2输出施加于被测系统的输入端r(t),然后分别测量被测系统的输出信号的开环对数幅值和相位,数据经相关运算后在虚拟示波器中显示。 实验步骤: (1)将数/模转换器(B2)输出OUT2作为被测系统的输入。 (2)构造模拟电路:安置短路套及测孔联线表同笫3.2.2 节《二阶闭环系统的频率特性曲线测试》。 (3)运行、观察、记录: ①将数/模转换器(B2)输出OUT2作为被测系统的输入,运行LABACT 程序,在界面的自动控制菜单下的线性控制系统的频率响应分析-实验项目,选择二阶系统,就会弹出虚拟示波器的界面,点击开始,实验开始后,实验机将自动产生0.5Hz~16H等多种频率信号,等待将近十分钟,测试结束后,观察闭环对数幅频、相频曲线和幅相曲线。
姓名:陈,H 学号:XXXXXXXX 班级:电气 实验五 典型环节和系统频率特性的测量 一、实验目的 1.了解典型环节和系统的频率特性曲线的测试方法; 2.根据实验求得的频率特性曲线求取传递函数。 二、实验设备 1.THBDC-1型 控制理论·计算机控制技术实验平台; 2.PC 机一台(含“THBDC-1”软件)、USB 数据采集卡、37针通信线1根、16芯数据排线、USB 接口线。 三、实验内容 1.惯性环节的频率特性测试; 2.二阶系统频率特性测试; 3.无源滞后—超前校正网络的频率特性测试; 4.由实验测得的频率特性曲线,求取相应的传递函数; 5.用软件仿真的方法,求取惯性环节和二阶系统的频率特性。 四、实验原理 1.系统(环节)的频率特性 设G(S)为一最小相位系统(环节)的传递函数。如在它的输入端施加一幅值为X m 、频率为ω的正弦信号,则系统的稳态输出为 )sin()()sin(?ωω?ω+=+=t j G Xm t Y y m 由式①得出系统输出,输入信号的幅值比相位差 ) () (ωωj G Xm j G Xm Xm Ym == (幅频特性) )()(ωωφj G ∠= (相频特性) 式中)(ωj G 和)(ωφ都是输入信号ω的函数。 2.频率特性的测试方法 2.1 李沙育图形法测试 2.1.1幅频特性的测试 由于 m m m m X Y X Y j G 22)(== ω 改变输入信号的频率,即可测出相应的幅值比,并计算 m m X Y A L 22l o g 20)(log 20)(==ωω (dB ) 其测试框图如下所示:
图5-1 幅频特性的测试图(李沙育图形法) 注:示波器同一时刻只输入一个通道,即系统(环节)的输入或输出。 2.1.2相频特性的测试 图5-2 相频特性的测试图(李沙育图形法) 令系统(环节)的输入信号为:t X t X m ωsin )(= (5-1) 则其输出为 )s i n ()(φω+=t Y t Y m (5-2) 对应的李沙育图形如图5-2所示。若以t 为参变量,则)(t X 与)(t Y 所确定点的轨迹将在示波器的屏幕上形成一条封闭的曲线(通常为椭圆),当t=0时,0)0(=X 由式(5-2)得 )s i n ()0(φm Y Y = 于是有 m m Y Y Y Y 2)0(2sin )0(sin )(1 1 --==ωφ (5-3) 同理可得 m X X 2)0(2s i n )(1 -=ωφ (5-4) 其中: )0(2Y 为椭圆与Y 轴相交点间的长度; )0(2X 为椭圆与X 轴相交点间的长度。 式(5-3)、(5-4)适用于椭圆的长轴在一、三象限;当椭圆的长轴在二、四时相位φ的计算公式变为 m Y Y 2)0(2sin 180 )(1 --=ωφ 或 m X X 2)0(2s i n 180 )(1 --=ωφ 下表列出了超前与滞后时相位的计算公式和光点的转向。
广西大学实验报告纸 姓名: 指导老师:胡老师 成绩: 学院:电气工程学院 专业:自动化 班级:121 实验内容:零、极点对限性控制系统的影响 2014年 11月 16 日 【实验时间】2014年11月14日 【实验地点】宿舍 【实验目的】 1. 掌握测量典型一阶系统和二阶系统的频率特性曲线的方法; 2. 掌握软件仿真求取一、二阶系统的开环频率特性的方法; 3. 学会用Nyquist 判据判定系统的稳定性。 【实验设备与软件】 1. labACT 实验台与虚拟示波器 2. MATLAB 软件 【实验原理】 1.系统的频率特性测试方法 对于现行定常系统,当输入端加入一个正弦信号)sin()(t X t X m ωω=时,其稳态输出是一个与输入信号频率相同,但幅值和相位都不同的正弦信号 )sin()()sin()(ψωωψω+=+=t j G X t Y s Y m m 。 幅频特性:m m X Y j G /)(=ω,即输入与输出信号的幅度比值,通常转换成)(lg 20ωj G 形式。 相频特性:)(arg )(ωω?j G =,可以直接基于虚拟示波器读取,也可以用“李沙育图行”法得到。 可以将用Bode 图或Nyquist 图表示幅频特性和相频特。 在labACT 试验台采用的测试结构图如下:
被测定稳定系统对于实验就是有源放大电路模拟的一、二阶稳定系统。 2.系统的频率测试硬件原理 1)正弦信号源的产生方法 频率特性测试时,一系列不同频率输入正弦信号可以通过下图示的原理产生。按照某种频率不断变化的数字信号输入到DAC0832,转换成模拟信号,经一级运放将其转换为模拟电压信号,再经过一个运放就可以实现双极性电压输出。 根据数模转换原理,知 R V N V 8 012- = (1) 再根据反相加法器运算方法,得 R R R V N V N V R R V R R V 1281282282201210--=??? ??+-?-=??? ? ??+-= (2) 由表达式可以看出输出时双极性的:当N 大于128时,输出为正;反之则为负;当输入为128时,输出为0. 在labACT 实验箱上使用的参考电压时5V 的,内部程序可以产生频率范围是对一阶系统是0.5 H Z ~64H Z 、对二阶系统是0.5 H Z ~16 H Z 的信号,并由B2单元的OUT2输出。
第五章线性系统的频域分析法 一、频率特性四、稳定裕度 二、开环系统的典型环节分解 五、闭环系统的频域性能指标 和开环频率特性曲线的绘制 三、频率域稳定判据 本章主要内容: 1 控制系统的频带宽度 2 系统带宽的选择 3 确定闭环频率特性的图解方法 4 闭环系统频域指标和时域指标的转换 五、闭环系统的频域性能指标
1 控制系统的频带宽度 1 频带宽度 当闭环幅频特性下降到频率为零时的分贝值以下3分贝时,对应的频率称为带宽频率,记为ωb。即当ω>ωb 而频率范围(0,ωb)称为系统带宽。 根据带宽定义,对高于带宽频率的正弦输入信号,系统输出将呈现较大的衰减,因此选取适当的带宽,可以抑制高频噪声的影响。但带宽过窄又会影响系统正弦输入信号的能力,降低瞬态响应的速度。因此在设计系统时,对于频率宽度的确定必须兼顾到系统的响应速度和抗高频干扰的要求。 2、I型和II型系统的带宽 2、系统带宽的选择 由于系统会受多种非线性因素的影响,系统的输入和输出端不可避免的存在确定性扰动和随机噪声,因此控制系统的带宽的选择需综合考虑各种输入信号的频率范围及其对系统性能的影响,即应使系统对输入信号具有良好的跟踪能力和对扰动信号具有较强的抑制能力。 总而言之,系统的分析应区分输入信号的性质、位置,根据其频谱或谱密度以及相应的传递函数选择合适带宽,而系统设计主要是围绕带宽来进行的。 3、确定闭环频率特性的图解方法
1、尼科尔斯图线 设开环和闭环频率特性为 4、闭环系统频域指标和时域指标的转换 工程中常用根据相角裕度γ和截止频率ω估算时域指标的两种方法。 相角裕度γ表明系统的稳定程度,而系统的稳定程度直接影响时域指标σ%、ts。 1、系统闭环和开环频域指标的关系 系统开环指标截止频率ωc与闭环带宽ωb有着密切的关系。对于两个稳定程度相仿的系统,ωc大的系统,ωb也大;ωc小的系统,ωb也小。 因此ωc和系统响应速度存在正比关系,ωc可用来衡量系统的响应速度。又由于闭环振荡性指标谐振Mr和开环指标相角裕度γ都能表征系统的稳定程度。 系统开环相频特性可表示为