自动控制原理实验六 线性系统的频域分析
长安大学：自动控制原理第五章 线性系统的频域分析
A () 1 0 T
() 0
() 90
V() A() sin ()
长安大学信息工程学院
自动控制理论
第五章
二、研究频率特性的意义 1、频率特性是控制系统在频域中的一种数学模型，是研究自 动控制系统的另一种工程方法。 2、根据系统的频率性能间接地揭示系统的动态特性和稳态特 性，可以简单迅速地判断某些环节或参数对系统性能的影响， 指出系统改进的方向。 3、频率特性可以由实验确定，这对于难以建立动态模型的系 统来说，很有用处。 三、频率特性的求取方法 1、已知系统的系统方程，输入正弦函数求其稳态解，取输 出稳态分量和输入正弦的复数比； 2、根椐传递函数来求取； 3、通过实验测得。
设
x c (t) ae jt ae jt b1es1t b2es2t ... b1esn t
A AG( j) ( s j ) | s j s 2 2 2j
( t 0)
对于稳定的系统, -s1,s2,…,sn 其有负实部
x c (t) ae jt ae jt
a G(s)
a G (s)
CHANG’AN UNIVERSITY
A AG( j) ( s j ) | s j s 2 2 2j
长安大学信息工程学院
自动控制理论
第五章
a
AG( j) 2j
AG( j) a 2j
G( j) | G( j) | e jG( j) | G( j) | e jG( j)
幅频特性 相频特性 实频特性 虚频特性
CHANG’AN UNIVERSITY
A() | G ( j) | U 2 () V 2 () 1 V() () G( j) tg U () U() A() cos()
《自动控制原理》实验3.线性系统的频域分析
《自动控制原理》实验3.线性系统的频域分析实验三线性系统的频域分析一、实验目的1．掌握用MATLAB语句绘制各种频域曲线。
2．掌握控制系统的频域分析方法。
二、基础知识及MATLAB函数频域分析法是应用频域特性研究控制系统的一种经典方法。
它是通过研究系统对正弦信号下的稳态和动态响应特性来分析系统的。
采用这种方法可直观的表达出系统的频率特性，分析方法比较简单，物理概念明确。
1．频率曲线主要包括三种：Nyquist图、Bode图和Nichols图。
1）Nyquist图的绘制与分析MATLAB中绘制系统Nyquist图的函数调用格式为：nyquist(num,den) 频率响应w的范围由软件自动设定 nyquist(num,den,w) 频率响应w的范围由人工设定[Re,Im]= nyquist(num,den) 返回奈氏曲线的实部和虚部向量，不作图2s?6例4-1：已知系统的开环传递函数为G(s)?3，试绘制Nyquists?2s2?5s?2图，并判断系统的稳定性。
num=[2 6]; den=[1 2 5 2]; nyquist(num,den)极点的显示结果及绘制的Nyquist图如图4-1所示。
由于系统的开环右根数P=0，系统的Nyquist曲线没有逆时针包围（-1，j0）点，所以闭环系统稳定。
p =-0.7666 + 1.9227i -0.7666 - 1.9227i -0.4668图4-1 开环极点的显示结果及Nyquist图若上例要求绘制??(10?2,103)间的Nyquist图，则对应的MATLAB语句为：num=[2 6]; den=[1 2 5 2];w=logspace(-1,1,100); 即在10-1和101之间，产生100个等距离的点nyquist(num,den,w)2）Bode图的绘制与分析系统的Bode图又称为系统频率特性的对数坐标图。
Bode图有两张图，分别绘制开环频率特性的幅值和相位与角频率?的关系曲线，称为对数幅频特性曲线和对数相频特性曲线。
自动控制原理实验报告四
自动控制原理实验报告实验时间：201Ｘ年Ｘ 月Ｘ 日 地点：ＸＸＸＸ 实验报告人（签名）：倪马 同组实验人（签名）：1 实验名称：线性系统的频域分析2 实验目的：(1)掌握二阶开环系统的对数频率特性、幅相频率特性、实频特性和虚频特性；(2)掌握欠阻尼二阶闭环系统中的自然频率、阻尼比对开环参数幅值穿越频率、相位裕度的影响，以及幅值穿越频率和相位裕度的计算；(3)掌握二阶开环系统对数频率特性曲线和幅相频率特性曲线的特点及绘制方法。
3 实验内容：（1）根据实验讲义上模拟电路图和接线要求，在LabACTn 自控／计控原理实验机的对应接口上连接好线路；（2）根据实验讲义的具体要求进行“运行”操作，并观察实验曲线，根据曲线计算对应参数——一阶惯性环节的转折频率、二阶闭环系统的谐振频率r ω&谐振峰值)(r L ω，改变被测系统的各项电路参数，画出其系统模拟电路图。
4 实验步骤 4.1 实验操作4.1.1 一阶惯性环节的频率特性曲线惯性环节的频率特性测试电路见图3-2-1，改变被测系统的各项电路参数，画出其系统模拟电路图，及频率特性曲线，並计算和测量其转折频率，填入实验报告。
一阶惯性环节的转折频率：T /1=ω图3-2-1 惯性环节的频率特性测试电路图3-2-1电路的增益K=1，惯性时间常数 T=0.1，转折频率：s /1rad .0/1==T ω 实验内容及步骤（1）构造模拟电路：按图3-2-1安置短路套及插孔连线。
（2）运行、观察、记录：① 选择系统的频域分析／一阶惯性环节频率特性曲线，将弹出频率特性扫描点设置表，用户可在…频率特性扫描点‟设置表中根据需要填入各个扫描点角频率，设置完后，点击《确认》后，将弹出…频率特性曲线‟实验界面，点击《开始》，即可按表中规定的角频率值，按序自动产生多种频率信号，画出频率特性曲线。
② 测试结束后（约五分钟），将显示被测系统的对数幅频、相频特性曲线（伯德图）和幅相曲线（奈奎斯特图），界面“显示选择”选择了“伯德图”。
线性系统的频域分析实验心得
线性系统的频域分析实验心得
1·熟练掌握用 MATLA语句绘制频域曲线。
2·掌握控制系统频域范围内的分析校正方法。
3掌握用频率特性法进行串联校正设计的思路和步骤
某单位负反馈控制系统的开环传递函数4为，试设计一超前校正装置,G(s)1、' s(s 1)K. 20s 150使校正后系统的静态速度误差系数，相位裕量，增益裕量20lgK10dB
绘制伯德图程序，以及计算穿越频率，相位裕量ans =相位 Inf 9.0406频率Inf 3.1425>e=5; r=50; rO=9; >>[gm1,pm1,wcg1,wcp1]=marg in(num 0,de nO);phic=(r-rO+e)*pi/180;
[gm1,pm1,wcg1,wcp1]=margi n(num 0,de nO);>>alpha=(1+s in (phic))/(1-si n(phic))[gm1,pm1,wcg1,wcp1]=margin(num 0,de n0); alpha =6.1261 [gm1,pm1,wcg1,wcp1]=marg in(num 0,de n0);lgm1,pm1,wcg1,wcp1]
通过MATLAB寸系统进行校正，可以清晰明了的显示矫正过程，以及矫正结果，方便快捷。
这种基于MATLAB的方法对于系统的设计非常实用。
值得以后再学习过程中认真领悟学习!! ! ! !。
自动控制原理课件：线性系统的频域分析
包围坐标原点 − 周。
m
F (s)
K1 ( s z j )
j 1
n
i 1
( s pi )
24
• 02
基本概念
m
1 G ( s) H ( s) F ( s)
K1 ( s z j )
j 1
在 平面上的映射曲线 F 1 G ( j ) H ( j )将按逆时针方向
围绕坐标原点旋转 = − 周。
如果在s平面上，s沿着奈奎斯特回线顺时针方向移动一周时，
在 平面上的映射曲线围绕坐标原点按逆时针方向旋转 =
周，则系统为稳定的。
26
根据
( 1, j 0)
L( ) 20 lg K 20 lg 1 12 2 20 lg 1 22 2
( ) arctg 1 arctg 2
τ2
20dB / dec 1
2
L3 ( )
L2 ( )
40dB / dec
( )
0
L( )
90
A( ) 1, ( )
L ( ) 20 lg A( ) 0
L( )
jQ( )
L( ) 0
0
( )
1
0
1
P( )
1
0
30
60
16
5.3
系统开环频率特性图
设开环系统由n个典型环节串联组成
G(s ) G 1(s )G 2(s ) G n(s )
这意味着 的映射曲线 F 围绕原点运动的情况，相当于
自动控制原理实验报告
3.1 线性系统的时域分析3.1.1 典型环节的模拟研究一. 实验目的1． 了解和掌握各典型环节模拟电路的构成方法、传递函数表达式及输出时域函数表达式2． 观察和分析各典型环节的阶跃响应曲线，了解各项电路参数对典型环节动态特性的影响二．实验内容及步骤观察和分析各典型环节的阶跃响应曲线，了解各项电路参数对典型环节动态特性的影响.。
改变被测环节的各项电路参数，画出模拟电路图，阶跃响应曲线，观测结果，填入实验报告。
1)．观察比例环节的阶跃响应曲线典型比例环节模拟电路如图3-1-1所示。
图3-1-1 典型比例环节模拟电路传递函数：01(S)(S)(S)R R K KU U G i O ===；单位阶跃响应：K )t (U = 实验步骤：注：‘S ST ’用短路套短接！ （1）将函数发生器（B5）所产生的周期性矩形波信号（OUT ），作为系统的信号输入（Ui ）；该信号为零输出时，将自动对模拟电路锁零。
（2）构造模拟电路：按图3-1-1安置短路套及测孔联线，表如下：（a ）安置短路套（b ）测孔联线（3打开虚拟示波器的界面，点击开始，按下信号发生器（B1）阶跃信号按钮（0→+4V 阶跃），观测A5B 输出端（Uo ）的实际响应曲线。
示波器的截图详见虚拟示波器的使用。
实验数据纪录：2)．观察惯性环节的阶跃响应曲线典型惯性环节模拟电路如图3-1-2所示。
图3-1-2 典型惯性环节模拟电路传递函数：C R T R R K TSKU U G i O 1011(S)(S)(S)==+==单位阶跃响应：)1()(0Tte K t U --=实验步骤：注：‘S ST ’用短路套短接！ （1）将函数发生器（B5）所产生的周期性矩形波信号（OUT ），作为系统的信号输入（Ui ）； （1）该信号为零输出时，将自动对模拟电路锁零。
（2）构造模拟电路：按图3-1-4安置短路套及测孔联线，表如下： （a（3）运行、观察、记录：打开虚拟示波器的界面，点击开始，按下信号发生器（B1）阶跃信号按钮时（0→+4V 阶跃），观测A5B输出端（Uo）响应曲线，等待完整波形出来后，移动虚拟示波器横游标到输出稳态值×0.632处，得到与输出曲线的交点，再移动虚拟示波器两根纵游标，从阶跃开始到输出曲线的交点，量得惯性环节模拟电路时间常数T。
《自动控制原理》Matlab求解控制系统频域分析实验
《自动控制原理》Matlab求解控制系统频域分析实验
一、实验目的
1、加深了解系统频率特性的概念。
2、学习使用Matlab软件绘制Nyquist图、Bode图的基本方法。
3、掌握典型环节的频率特性。
二、实验仪器
Matlab2014b版
三、实验原理
1.奈奎斯特图（幅相频率特性图）
MATLAB为用户提供了专门用于绘制奈奎斯特图的函数nyquist
axis([-2,0.4,-1.5,1.5]);
k=500;
num=[1,10];
den=conv([1,0],conv([1,1],conv([1,20],[h,50])))；
w=logspace(-1,3,200)
bode(k*num,den,w);
grid;
五、实验原始数据记录与数据处理
六、实验结果与分析讨论
范围是自动确定的。当需要指定幅值范围和相角范围时，则需用下面的功能指令：
[mag,phase,w]=bode(num,den,w)
四、实验内容及步骤
z=[]:
p=[0,-1,-2]:
k=5;
g=zpk(z,p,k):
nyquist(g);
w=0.5：0.1:10:
figure(2):
nyquist(g:w);
自动控制原理--第5章 频域分析法
L() 20lg | G( j) | 20lg 2T 2 1
arctanT
当=0时，L()=0dB， =0， 曲线起始于坐标原点；当=1/T时， L()=-3dB， =-45；
自动控制原理
30
5-4 频域稳定性判据
一、映射定理
闭环特征函数 F(s)=1+G(s)H(s)
T
如果τ>T，则∠G(j)>0°，极坐标曲线在第Ⅰ象限变化；如果τ<T， 则∠G(j)<0°，极坐标曲线在第Ⅳ象限变化，如图所示。
自动控制原理
16
5.3.2 对数坐标图
通过半对数坐标分别表示幅频特性和相频特性的图形， 称为对数坐称图或波德(Bode)图。
1.对数坐标 对数频率特性曲线由对数幅频特性和相频特性两部分
系统的传递函数为 C(s) G(s)
R(s)
假定输入信号r(t)为
r(t) Asint
R(s) L[ Asint] A
A
s 2 2 (s j)(s j)
自动控制原理
7
G(s)
K (s z1 )(s z2 )(s zm ) (s s1 )(s s2 )(s sn )
nm
2j
AG( j) sin(t )
B sin(t )
G( j ) G( j ) e jG( j) G( j) e j
即
G( j) G(s) s j
这里的结论同RC网络讨论的结果是一致的。
自动控制原理
10
5.3 频率特性的图示方法
频率特性的图示方法主要有三种，即极坐标图、对数坐 标图和对数幅相图，现分述如下。
所以K=10。因此，所求开环传递函数
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
G=
0.375
-----------------
s^3 + 3 s^2 + 2 s
Continuous-time transfer
（ 1 ） 低 频 段 增 益 -20db/dec ， 高 频 段 -60db/dec,低频段渐进相位角为-90°，高频 段为-270°，增益裕度 Gm=1.5000,相位裕度 Pm= 11.4304°
二. 实验原理及内容
1、频率特性函数 G( j ) 。
频率特性函数为：
G(
jw)
b0 ( j ) m b1 ( j ) m1 bm1 ( j ) bm a0 ( j ) n a1 ( j ) n an1 ( j ) an
由下面的 MATLAB 语句可直接求出 G(jw)。
[mag,phase,w]=bode(num,den)
[mag,phase,w]=bode(num,den,w)
由于伯德图是半对数坐标图且幅频图和相频图要同时在一个绘图窗口中绘制，因此，要用到半
对数坐标绘图函数和子图命令。
（1） 对数坐标绘图函数 利用工作空间中的向量 x，y 绘图，要调用 plot 函数，若要绘制对数或半对数坐标图， 只需要用相应函数名取代 plot 即可，其余参数应用与 plot 完全一致。 （2） 子图命令
i=sqrt(-1) % 求取-1 的平方根
GW=polyval(num，i*w)./polyval(den，i*w)
2、用 MATLAB 作奈魁斯特图。
控制系统工具箱中提供了一个 MATLAB 函数 nyquist( )，该函数可以用来直接求解
Nyquist 阵列或绘制奈氏图。当命令中不包含左端返回变量时，nyquist（）函数仅在屏幕上
s(s 1)(s 2)
计算系统的增益裕度，相位裕度，在 Bode 图上标注低频段斜率，高频段斜率及低频段、高
频段的渐近相位角。(2) 如果希望增益裕度为 16dB，求出响应的 k 值，并验证
伯德图
结果分析
（2）计算-90°-tan-1(w)-tan-1(w/2)=-180°， 计算得到 wc=1.414，代入式子得到 K=0.375 验证：num=[0.375],den=[1 3 2 0],
此外，控制系统工具箱中提供了 margin()函数来求取给定线性系统幅值裕量和相位裕 量，该函数可以由下面格式来调用：
margin（num , den）;给定开环系统的数学模型，作 Bode 图，并在图上方标注幅值裕度 Gm 和对应频率ωg，相位裕度 Pm 和对应的频率ωc。
[Gm, Pm, Wcg, Wcp]=margin(G); 如果已知系统的频率响应数据，我们还可以由下面的格式调用此函数。
Gm=1 Pm =9.5374e-06 Wcg =4.4721 Wcp =4.4721
Pm 系统响应在 K=12 时产
生等幅振 荡，此时临
界稳定
K1 2
K=5
K1 2
K=15
5 ---------------------0.05 s^3 + 0.6 s^2 + s
Gm=2.40 Pm = 19.9079 Wcg =4.4721 Wcp =2.7992
function. Gm =16.0000 Pm = 74.3477 Wcg = 1.4142 Wcp =0.1836
6、系统对数频率稳定性分析
【自我实践 6-2】系统开环传递函数 G(s)
k
，试分析系统的稳定性。
s(0.5s 1)(0.1s 1)
K
响应曲线
运行结果
分析
12
Gm
K=12
---------------------0.05 s^3 + 0.6 s^2 + s
实验六 线性系统的频域分析
一. 实验目的
（1）熟练掌握使用 MATLAB 命令绘制控制系统 Nyquist 图的方法； （2）能够分析控制系统 Nyquist 图的基本规律； （3）加深理解控制系统乃奎斯特稳定性判据的实际应用； （4）学会利用奈氏图设计控制系统； （5）熟练掌握运用 MATLAB 命令绘制控制系统伯德图的方法； （6）了解系统伯德图的一般规律及其频域指标的获取方法； （7）熟练掌握运用伯德图分析控制系统稳定性的方法； （8）设计超前校正环节并绘制 Bode 图； （9）设计滞后校正环节并绘制 Bode 图。
【自我实践 6-3】某单位负反馈系统的开环传递函数 G(s)
31.6
，求(1)绘制
产生奈氏图，命令调用格式为：
nyquist(num,den) ; 作 Nyquist 图， nyquist(num,den,w); 作开环系统的奈氏曲线， 3、奈奎斯特稳定性判据（又称奈氏判据） 反馈控制系统稳定的充分必要条件是当ω从-∞变到∞时，开环系统的奈氏曲线不穿过 点（-1，j0）且逆时针包围临界点（-1，j0）点的圈数 R 等于开环传递函数的正实部极点数。 4、用 MATLAB 作伯德图 控制系统工具箱里提供的 bode()函数可以直接求取、绘制给定线性系统的伯德图。 命令的调用格式为：
15 ---------------------0.05 s^3 + 0.6 s^2 + s
Gm=0.800 Pm = -4.6920 Wcg = 4.4721 Wcp = 4.9916
Gm Pm 系统响应在 K 12 时产 生减幅振 荡，此时系 统稳定
Gm Pm 系统响应在 K 12 时产 生发散振 荡，此时系 统不稳定
[Gm, Pm, Wcg, Wcp]=margin(mag, phase, w); 其中(mag, phase, w)分别为频率响应的幅值、相位与频率向量。
如果系统的相角裕量 γ>45o,我们一般称该系统有较好的相角裕量。
【自我实践 6-1】某单位负反馈系统的开环传递函数 G(s)
k
，求(1) 当 k=4 时，
MATLAB 允许将一个图形窗口分成多个子窗口，分别显示多个图形，这就要用到 subplot()函数，其调用格式为：subplot(m，n，k) 5、 用 MATLEB 求取稳定裕量
同前面介绍的求时域响应性能指标类似，由 MATLAB 里 bode()函数绘制的伯德图也可 以采用游动鼠标法求取系统的幅值裕量和相位裕量。