说明1 自动控制原理的电子版内容以胡寿松教授主编的第五版“自动控制原理”为基础,以PowerPoint 2000和MATLAB6.5为工具,以帮助教师更好地讲好自控、帮助学生更好地学好自控为目的而制作的。 本课件大部分内容都是以点击鼠标的方式分步出现的,点击鼠标右键选择“定位”,然后再点击“幻灯片漫游”,可进入各章节学习。使用者在使用前应先看看各章说明,即可理解其含意。
说明2 课件3 ~6为第一章的内容。制作目的是节省画图时间,便于教师讲解。 课件6要强调串联并联反馈的特征,在此之前要交待相邻综合点与相邻引出点的等效变换。 课件7中的省略号部分是反过来说,如‘合并的综合点可以分开’等。最后一条特别要讲清楚,这是最容易出错的地方! 课件10先要讲清H1和H3的双重作用,再讲分解就很自然了。 课件11 、12 、13是直接在结构图上应用梅逊公式,制作者认为没必要将结构图变为信号流图后再用梅逊公式求传递函数。
说明3 课件17~30为第三章的内容。 课件17~19中的误差带均取为稳态值的5%,有超调的阶跃响应曲线的上升时间为第一次到达稳态值的时间。 课件20要讲清T的求法,T与性能指标的关系。 课件21要说明这是无零点的二阶系统。 课件22要交待Φ(s)的分母s2项的系数,且分子分母常数项相等。 课件28小结中的3个问题答案:1、系统稳定且 ; 2、非单位反馈输出端定义的误差 可通过等效变换后使用;3 、系统稳定。
说明4 课件32~42为第四章的内容。 课件32中的‘注意’应在观看‘rltool’后讲解。若不演示‘rltool’也可以。 课件33结论1和2与书中的相同,结论3分为n>m,n=m,n<m这3种情况介绍,其中n为开环极点数,m为开环零点数。 课件34根轨迹出现后,先介绍图上方的C(s)=6实际是K*=6,图中的3个小方块为K*=6所对应的3个闭环极点,然后验证模值条件和相角条件。 课件35要强调是1+,不能是1-,分子分母中的因子s的系数为1,不能为-1,K*不能为负。 课件41先回顾180o根轨迹的模值方程和相角方程,然后再介绍零度根轨迹的模值方程和相角方程。
说明5 课件44~63为第五章内容 课件44要说明几个问题:1.给一个稳定的系统输入一个正弦,其稳态输出才是正弦,幅值改变相角改变;2.不稳定的系统输出震荡发散,该振荡频率与输入正弦的频率有无关系?3.不稳定的系统输入改为阶跃时,其输出曲线类似,此时用运动模态来解释。 课件45中的省略号内容为:输入初始角不为零时如何处理,输入为余弦时没必要改为正弦。 课件57种的几点说明内容为:1. 增加k 值曲线上下平移,2. 取不同的值时,修正值不同,详细情况参考课件57。
第一章 自动控制的一般概念 1-1 自动控制的基本原理与方式 1-2 自动控制系统示例 1-3 自动控制系统的分类 1-4 对自动控制系统的基本要求
飞机示意图 反馈电位器 给定电位器
飞机方块图 扰动 给定装置 放大器 θ0 θc 舵机 飞机 反馈电位器 垂直陀螺仪 俯仰角控制系统方块图
液位控制系统 控制器 Q1 浮子 电位器 减速器 c 电动机 用水开关 SM Q2 if
第二章 控制系统的数学模型 2-1 时域数学模型 2-2 复域数学模型 2-3 结构图与信号流图
结构图三种基本形式 串 联 并 联 反 馈 G1 G2 G2 G1 G1 G2 G1 G2 1+ G1 G2 G1 G2
结构图等效变换方法 1 三种典型结构可直接用公式 2 相邻综合点可互换位置、可合并… 3 相邻引出点可互换位置、可合并… 注意事项: 1 三种典型结构可直接用公式 2 相邻综合点可互换位置、可合并… 3 相邻引出点可互换位置、可合并… 注意事项: 1 不是典型结构不可直接用公式 2 引出点综合点相邻,不可互换位置
H2 引出点移动 G1 G2 G3 G4 H3 H1 请你写出结果,行吗? G4 1 G1 G2 G3 G4 H3 H2 H1 a b
综合点移动 向同类移动 无用功 G1 G2 G3 H1 G2 错! G2 H1 G1 G3 G1
作用分解 G1 G4 H3 G2 G3 H1 G1 G4 G2 G3 H3 H1 H3 H1
Pk—从R(s)到C(s)的第k条前向通路传递函数 = ∑Pk△k △ : 梅逊公式介绍 R-C △称为系统特征式 △= - ∑La + ∑LbLc -∑LdLeLf+… 1 其中: — 所有单独回路增益之和 ∑La ∑LbLc—所有两两互不接触回路增益乘积之和 ∑LdLeLf—所有三个互不接触回路增益乘积之和 Pk—从R(s)到C(s)的第k条前向通路传递函数 △k称为第k条前向通路的余子式 去掉第k条前向通路后所求的△ △k求法: △k=1-∑LA+ ∑LBLC- ∑LDLELF+…
梅逊公式例R-C G4(s) 请你写出答案,行吗? G3(s) H3(s) H3(s) H3(s) H3(s) H3(s) H3(s) R(s) C(s) △1=1 △2=1+G1H1 G4(s) G3(s) C(s) R(s) =? 请你写出答案,行吗? G1(s) G2(s) G3(s) P1=G1G2G3 P2= G4G3 L1= –G1 H1 L2= – G3 H3 L3= – G1G2G3H3H1 L4= – G4G3 L5 = – G1G2G3 L1L2= (–G1H1) (–G3H3) = G1G3H1H3 L1L4=(–G1H1)(–G4G3)=G1G3G4H1
E(s)= 1 梅逊公式求E(s) P1=1 N(s) P2= - G3G2H3 △2= 1 P2△2=? P1= –G2H3 △1= 1 G3(s) G2(s) H3(s) E(S) R(s) N(s) G1(s) H1(s) H2(s) C(s) G3(s) G2(s) H3(s) R(s) E(S) G1(s) G3(s) H1(s) G2(s) H3(s) H2(s) R(s) C(s) N(s) E(S) G1(s) G3(s) H1(s) G2(s) H3(s) H2(s) C(s) N(s) R(s) E(S) G1(s) G3(s) H1(s) G2(s) H3(s) H2(s) R(s) C(s) N(s) E(S) R(s) E(S) P2= - G3G2H3 △2= 1 P2△2=? G1(s) H1(s) H2(s) C(s) P1= –G2H3 △1= 1 P1=1 △1=1+G2H2 P1△1= ? R(s)[ ] (1+G2H2) + (- G3G2H3) + (–G2H3) N(s) E(s)= 1 + G2H2 + G1G2H3 -G1H1G2 H2 - G1H1
信号流图 前向通路两条 四个单独回路,两个回路互不接触 (1 g) – b C(s) d a b c + e d = R(s) – a f h 前向通路两条 四个单独回路,两个回路互不接触 (1 g) – b C(s) d a b c + e d = R(s) – a f b g c h – 1 – – e f h g + a h f c
第三章 线性系统的时域分析法 3-1 时域性能指标 3-2 一阶系统时域分析 3-3 二阶系统时域分析 3-4 稳定性分析 3-1 时域性能指标 3-2 一阶系统时域分析 3-3 二阶系统时域分析 3-4 稳定性分析 3-6 稳态误差计算
动态性能指标定义1 h(t) t h(t) t h(t) t A B 超调量σ% = 100% 超调量σ% = A B 100% A B 时间tr 上 升 峰值时间tp A B 超调量σ% = 100% 调节时间ts h(t) t 超调量σ% = A B 100% A B 峰值时间tp 时间tr 上 升 调节时间ts
动态性能指标定义2 h(t) t 调节时间 ts 上升时间tr
动态性能指标定义3 h(t) t σ%= B A 100% A B tp tr ts
一阶系统时域分析 问 时间常数 无零点的一阶系统 Φ(s)= e- 单位斜坡响应 单位阶跃响应 ? k(t)= h(t)=1-e-t/T Ts+1 k , T 时间常数 (画图时取k=1,T=0.5) k(t)= T 1 e- t h(t)=1-e-t/T c(t)=t-T+Te-t/T r(t)= δ(t) r(t)= 1(t) r(t)= t 单 位 脉 冲 响 应 k(0)= T 1 单位斜坡响应 h’(0)=1/T 单位阶跃响应 K’(0)= T 1 2 h(T)=0.632h(∞) h(2T)=0.865h(∞) T h(3T)=0.95h(∞) h(4T)=0.982h(∞) 1 、3个图各如何求T? 2 、调节时间ts=? 问 ? 3 、r(t)=vt时,ess=? 4、求导关系
- ωn -cosωnt - ωn e- ω t 二阶系统单位 阶跃响应定性分析 + h(t)= ωn Φ(s)= s2+2 ωns+ωn2 j T1 1 T2 >1 =1 0< <1 =0 j ± √ 2 - 1 S1,2= - ωn ωn >1 j h(t)= 1 T2 t T1 1 e + S1,2= - ωn -ωn = =1 h(t)= 1 -(1+ωnt) e-ω t n 过阻尼 临界阻尼 j - ±j √1- 2 ωn S1,2= 0< <1 j sin(ωdt+β) e- ω t h(t)= √1- 2 1 n S1,2 = ±j ωn =0 h(t)= 1 -cosωnt 欠阻尼 零阻尼
e e sin( ) e-π sin( ) 欠阻尼二阶系统动态性能分析与计算 h(t)= 1- β - ωn h(t)= 1- t + β √1- 2 1 - ωnt sin( t+ ωd β ) Φ(s)= s2+2 ωns+ωn2 ωn 2 j ωd = ωn√1- 2 ωn π- β ωd 得 tr= 令h(t)=1取其解中的最小值, β 0 < <1时: - ωn 令h(t)一阶导数=0,取其解中的最小值, 得 tp= π ωd S1,2= - ωn ±j √1- 2 ωn 由σ%= h(∞) h(tp) -h(∞) 100% 得 σ% = e-π 100% h(t)= 1- √1- 2 1 e - ωnt sin( ωd t + β ) (0 ﹤ ≤ 0.8) 由包络线求调节时间
劳思表介绍 ε ε -8 ε ε ε s6+2s5+3s4+4s3+5s2+6s+7=0 (10-6)/2=2 劳 思 表 劳斯表特点 ε 设系统特征方程为: s6+2s5+3s4+4s3+5s2+6s+7=0 1 2 4 6 3 5 7 s6 s5 s0 s1 s2 s3 s4 1 2 4 6 3 5 7 (10-6)/2=2 (6-4)/2=1 4 1 2 (6-14)/1= -8 劳 思 表 劳斯表特点 1 2 7 -8 ε 1 1 2 7 ε -8 1 右移一位降两阶 2 +8 ε 2 劳思行列第一列不动 7 ε 3 次对角线减主对角线 -8(2 +8) - ε 7 2 4 每两行个数相等 7 ε 5 分母总是上一行第一个元素 6 一行可同乘以或同除以某正数 7 第一列出现零元素时, 用正无穷小量ε代替。
劳思判据 系统稳定的必要条件: 特征方程各项系数 均大于零! 系统稳定的充分条件: 劳思表第一列元素不变号! -s2-5s-6=0稳定吗? 有正有负一定不稳定! 缺项一定不稳定! -s2-5s-6=0稳定吗? 系统稳定的充分条件: 劳思表第一列元素不变号! 若变号系统不稳定! 变号的次数为特征根在s右半平面的个数!
劳思表出现零行 错啦!!! s4+5s3+7s2+5s+6=0 劳 思 表 劳斯表出现零行系统一定不稳定 s1,2=±j 设系统特征方程为: ① 有大小相等符号相反的 特征根时会出现零行 s4+5s3+7s2+5s+6=0 ② 由零行的上一行构成 辅助方程: 1 7 6 劳 思 表 1 1 5 5 s2+1=0 1 6 1 6 对其求导得零行系数: 2s1 继续计算劳斯表 2 第一列全大于零,所以系统稳定 1 错啦!!! 解辅助方程得对称根: s1,2=±j 劳斯表出现零行系统一定不稳定 1 劳斯表何时会出现零行? 由综合除法可得另两个根为s3,4= -2,-3 2 出现零行怎么办? 3 如何求对称的根?
误差定义 G(s) H(s) G(s) G(s) H(s) 总误差怎么求? ˊ ˊ 1 R(s) E(s) C(s) B(s) R(s) E(s)=R(s)-C(s) 输入端定义: E(s)=R(s)-B(s)=R(s)-C(s)H(s) G1(s) H(s) R(s) C(s) G2(s) N(s) 输出端定义: E(s)=C希-C实= -C(s) R(s) H(s) ˊ G(s) H(s) R(s) E(s) C(s) 1 ˊ En(s)=C希-C实= –Cn(s) 总误差怎么求?
kp kv ka 典型输入下的稳态误差与静态误差系数 ess= ess= ess= r(t)=R·1(t) R(s)=R/s G(s) H(s) R(s) E(s) C(s) ess= 1+ k sν R lim →0 s kp E(s)=R(s) 1+G(s)H(s) 1 r(t)=V·t R(s)=V/s2 ess= s· V lim →0 s k sν kv 若系统稳定, 则可用终值定理求ess r(t)=At2/2 R(s)=A/s3 ess= lim s 1+ k sν G0H0 R(s) →0 s ess= s2· A lim →0 s k sν ka
取不同的ν ∞ ∞ ∞ ∞ ∞ Kp=? 小结: Kv=? Ka=? k k k 1 ess= 2 ess= ess= 3 R·1(t) 稳态误差 静态误差系数 R·1(t) V·t At2/2 R·1(t) V·t At2/2 R 1+ k ∞ A k ∞ k k 0型 V k ∞ Ⅰ型 k ∞ ∞ Ⅱ型 1 ess= s· V lim →0 s k sν r(t)=R·1(t) Kp=? Kv=? Ka=? 啥时能用表格? r(t)=At2/2 小结: ess= 1+ k sν R lim →0 s 2 ess= s2· A lim →0 s k sν 非单位反馈怎么办? 3 表中误差为无穷时系统还稳定吗? r(t)=V·t
全 减小和消除误差的方法(1,2) essn= -limsC(s) =-lim Gn(s) k1 k2 s(T2s+1) 1 按扰动的全补偿 N(s) R(s) Gn(s) T1s+1 k1 s(T2s+1) k2 C(s) E(s) s (T1s+1)(T2s+1) + k1k2 (T1s+1)+ k1Gn(s) N(s) 令R(s)=0,En(s) = -C(s) = 令分子=0,得Gn(s) = - (T1s+1)/k1 这就是按扰动的全补偿 2 按扰动的稳态补偿 设系统稳定,N(s)=1/s ,则 t从0→∞全过程 全 k1k2 1+ k1Gn(s) essn= -limsC(s) =-lim s→0 ∴Gn(s)= -1/k1 各种干扰信号
减小和消除误差的方法(3,4) Gr(s) k1 s(T2s+1) k2 s (T1s+1)(T2s+1) + k1k2 N(s) R(s) Gr(s) T1s+1 k1 s(T2s+1) k2 C(s) E(s) 3 按输入的全补偿 s (T1s+1)(T2s+1) + k1k2 - k2 (T1s+1)Gr(s) R(s) 令N(s)=0, Er(s)= 令分子=0,得Gr(s)= s (T2s+1)/ k2 4 按输入的稳态补偿 设系统稳定,R(s)= 1/s2 则 1- k2 s Gr(s) k1k2 k2 s ∴Gr(s)= essr= limsEr(s)= lim s→0
第四章 线性系统的根轨迹法 4-1 根轨迹概念 4-2 绘制根轨迹的基本法则 4-3 广义根轨迹
注意: 根轨迹概念 k s(0.5s+1) 特征方程: 特征根:s1,2= -1±√1-2k K:0 ~ ∞ 0.5<k<∞时,s1,2=-1±j√2k-1 -2 -1 j 注意: 一组根对应同一个K; K一变,一组根变; K一停,一组根停; 演示rltool
G H 闭环零极点与开环零极点的关系 G ∏(s-pi ) ; ∏(s-zi H ∏(s-pj ) ∏(s-zj ∏(s-zi ) (s)= KG * ∏(s-pi q i=1 ) ; ∏(s-zi f H (s)= KH * ∏(s-pj h j=1 ) ∏(s-zj l ∏(s-zi f i=1 ) ∏(s-pj h j=1 * KG Φ(s )= ∏(s-pi q i=1 ) h j=1 ∏(s-pj ∏(s-zi f i=1 ) + KG * KH ∏(s-zj l j=1 结论:1 零点、 2 极点、3 根轨迹增益
模值条件与相 角条件的应用 -0.825 =0.466 ω n=2.34 -1.09+j2.07 78.8o 92.49o 127.53o 求模求角例题 模值条件与相 角条件的应用 -0.825 =0.466 ω n=2.34 -1.09+j2.07 78.8o 92.49o 2.61 2.26 127.53o 66.27o 2.072 2.11 s1=-0.825 s2,3= -1.09±j2.07 -1.5 -1 -2 0.5 2.26×2.11×2.61 2.072 K*= = 6.0068 92.49o- 66.27o- 78.8o- 127.53o= –180o
1 + K* = ∏ s zj ( - ) ∏ s pi ( - ) pi 根轨迹方程 特征方程 1+GH = 0 Zj 开环零点“○”,是常数! Zj j=1 m ∏ s zj ( - ) 1 + K* = i=1 n ∏ s pi ( - ) pi 开环极点“×”, 也是常数! 根轨迹增益K* ,不是定数,从0 ~ ∞变化 这种形式的特征方程就是根轨迹方程
K* = 1 1 + K* = K* = -1 ∏ ︱ s - zj pi ∏ ( s - zj pi ) ∏ ︱ s - zj pi 相角条件: 根轨迹的模值条件与相角条件 ∑∠(s-zj) -∑∠(s-pj) = (2k+1) π k=0, ±1, ±2, … j=1 i=1 m n 绘制根轨迹的充要条件 j=1 m n K* = 1 ∏ ︱ s - zj pi i=1 j=1 m n 1 + K* = ∏ ( s - zj pi ) i=1 K* = m n j=1 ∏ ︱ s - zj pi i=1 -1 模值条件: 确定根轨迹上某点对应的K*值
= 1 K* K* = 1 ∏ ︱ s - zj pi ∏ ︱ s - zj pi ∑ 绘制根轨迹的基本法则 n n , 1 根轨迹的条数 就是特征根的个数 7 与虚轴的交点 可由劳思表求出 或 令s=jω解出 2 根轨迹对称于 轴 实 8 起始角与终止角 3 根轨迹起始于 开环极点 ,终止于 开环零点 ( ) ∞ ( ) ∞ 4 ∣n-m∣条渐近线对称于实轴,均起于σa 点,方 j=1 m n = ∏ ︱ s - zj pi i=1 1 K* j=1 m n K* = 1 ∏ ︱ s - zj pi i=1 ∑pi-∑zj ∣n-m∣ i=1 j=1 n m σa = (n≠m?) 举例 φa= (2k+1)π n-m 向由φa确定: k= 0,1,2, … 5 实轴上某段右侧零、极点个数之和为奇数,则该段是根轨迹 实轴上的根轨迹 6 根轨迹的会合与分离 1 说明什么 2 d的推导 3 分离角定义 j=1 m ∑ i=1 n d-pi 1 d-zj = λL= (2k+1)π L , k= 0,1,2, … L为来会合的根轨迹条数 无零点时右边为零
根轨迹示例1 j j j j j j 同学们,头昏了吧? j j j j j
根轨迹示例2 n=[1 2];d=conv([1 2 5],[[1 6 10]);rlocus(n,d) j j j j j j n=[1 2];d=conv([1 2 5],[[1 6 10]);rlocus(n,d) n=1;d=conv([1 2 0],[1 2 2]);rlocus(n,d) j j j j j j
零度根轨迹 1– 1+ 请注意:G(s)H(s)的分子分母均首一 K*:0 ~ + K*:0 ~ – 特征方程为以下形式时,绘制零度根轨迹 1. K*:0 ~ + 1– 2. K*:0 ~ – 1+ 请注意:G(s)H(s)的分子分母均首一
K* = ∏ ︱ s - zj pi 2 kπ 零度 ∑∠(s-zj) -∑∠(s-pj) = (2k+1) π 模值条件: 相角条件: n 零度根轨迹的模值条件与相角条件 K* = m n j=1 ∏ ︱ s - zj pi i=1 模值条件: ∑∠(s-zj) -∑∠(s-pj) = (2k+1) π k=0, ±1, ±2, … j=1 i=1 m n 相角条件: 2 kπ
1 = K* ∏ ︱ s - zj pi ∑ 绘制零度根轨迹的基本法则 不变! 不变! 不变! 不变! 2kπ 偶 n 不变! , 不变! 根轨迹的条数 就是特征根的个数 不变! 不变! 2 根轨迹对称于 轴 实 不变! 3 根轨迹起始于 ,终止于 开环极点 开环零点 ( ) ∞ 4 ∣n-m∣条渐近线对称于实轴,起点 ∑pi-∑zj ∣n-m∣ i=1 j=1 n m σa = 不变! j=1 m n = ∏ ︱ s - zj pi i=1 1 K* 2kπ 渐近线方向: φa= (2k+1)π n-m k= 0,1,2, … 5 实轴上某段右侧零、极点个数之和为 奇 数,则该段是根轨迹 偶 6 根轨迹的分离点 j=1 m ∑ i=1 n d-pi 1 d-zj = k= 0,1,2, … λL= (2k+1)π L , 不变! 7 与虚轴的交点 不变! 8 起始角与终止角 变了
第五章 线性系统的频域分析法 5-1 频率判据 5-2 典型环节与开环频率特性 5-3 频域稳定判据 5-4 稳定裕度 5-5 闭环频域性能指标
频率特性的概念 不 40 设系统结构如图, 由劳思判据知系统稳定。 结论 给系统输入一个幅值不变频率不断增大的正弦, 曲线如下: 给稳定的系统输入一个正弦,其稳态输出是与输入 结论 同频率的正弦,幅值随ω而变,相角也是ω的函数。 Ar=1 ω=0.5 ω=1 ω=2 ω=2.5 ω=4
相角问题 , , … ① 稳态输出迟后于输入的角度为: A A B A 360o φ= ②该角度与ω有 关系 ∴为φ(ω) B B 关系 ∴为φ(ω) , B B ③该角度与初始 角度无关 ∴ , …
= = + + = + 设系统稳定,则正弦输入时输出为: ∏(s-si) ∏(s-zj) kΦ * s2+ω2 Arω Arω 频率特性 设系统稳定,则正弦输入时输出为: 频率特性 ∏(s-si) ∏(s-zj) kΦ * 1 n m s2+ω2 Arω = Φ(s) (s+jω)(s-jω) Arω C(s)=Φ(s)R(s)= ct(t)=∑ aies t i s-si ai ∑ 1 n = + s+jω B1 s-jω B2 ∵系统稳定,∴ ct(∞)=0 s+jω B1 + s-jω B2 Φ(jω) Ar 2j (s-jω) + = Φ(-jω) -2j(s+jω) Cs(s)= Φ(jω)ejωt Φ(-jω) e-jωt Ar 2j cs(t) Ar Φ(jω) ej∠Φ(jω) ejωt e-j∠Φ(jω) e-jωt 2j Φ(jω) = a(ω)+ j b(ω) c(ω)+ j d(ω) Φ(-jω) = c(ω)- j d(ω) a(ω)- j b(ω) ∠Φ(-jω) ∠Φ(jω) Φ(-jω) Φ(jω) Ar Φ(jω) sin(ωt+∠ Φ(jω))
对数坐标系
倒置的坐标系
积分环节L(ω) ① G(s)= ② G(s)= ③ G(s)= 1 s 10 s 1 5s 10 0.2 2 1 0.1 L(ω)dB ω 20 40 -40 -20 100 [-20] [-20] [-20]
微分环节L(ω) ② G(s)= ③ G(s)= ① G(s)= s 2s 0.1s 10 0.2 2 1 0.1 L(ω)dB ω 0dB 20 40 -40 -20 100 [+20] [+20] [+20]
惯性环节G(jω) G(s) = 0.5s+1 1 0.25 ω2+1 A(ω)= 1 φ(ω) = -tg-10.5 ω ω 0.5 1 0.5 1 2 4 5 8 20 φo(ω) A(ω) -14.5 -26.6 -45 -63.4 -68.2 -76 -84 1 0.97 0.89 0.71 0.45 0.37 0.24 0.05 j 1 Im[G(jω)] Re[G(jω)]
惯性环节L(ω) ① G(s)= ② G(s)= 1 100 0.5s+1 s+5 10 0.2 2 1 0.1 L(ω)dB ω 0dB 20 40 -40 -20 100 26dB [-20] 0o - 30o - 45o - 60o - 90o [-20]
一阶微分L(ω) (0.25s+0.1) ① G(s)= 0.5s+1 0.3 ② G(s)= [+20] [+20] L(ω)dB 10 40 -40 -20 100 0o +30o + 45o + 60o + 90o [+20] [+20]
振荡环节G(jω) (0< <1) (0< <0.707)
振荡环节G(jω)曲线 (Nyquist曲线) j 1
振荡环节L(ω) 10 0.2 2 1 0.1 L(ω)dB ω 0dB 20 40 -40 -20 100 [-40]
? 振荡环节再分析 w k G (s) = + w + w 2 ωr ωn ω = 友情提醒:φ (ωn)= - 90o S S r 2 n (0< <0.707) 0dB L(ω)dB ω ? 0< <0.5 = 0.5 0.5< <1 20lgk ωr ωn [-40] ω = r 友情提醒:φ (ωn)= - 90o
二阶微分 对数幅频渐近曲线 幅相曲线 ωn j L(ω)dB [+40] 1 ω 0dB 几点说明… 0< <0.707时有峰值: j 1 0dB L(ω)dB ω [+40] ωn 几点说明…
绘制L(ω)例题 绘制 的L(ω)曲线 10 0.2 2 1 0.1 L(ω)dB ω 0dB 20 40 -40 -20 100 [-20] [-40] [-20] [-40] 低频段: 时为38db 时为52db 转折频率:0.5 2 30 斜率: -20 +20 -20
解: 开环幅相曲线的绘制 求交点: 例题1:绘制 的幅相曲线。 令 . 6 4 , 5 )] j ( G Re[ = + w - 例题1:绘制 的幅相曲线。 解: 求交点: 令 . 6 4 , 5 )] j ( G Re[ 2 = + w - -25 Im[G(jω)] Re[G(jω)] 无实数解,与虚轴无交点 曲线如图所示: 开环幅相曲线的绘制
同时成立! 稳定裕度的定义 特点: G(jω)曲线过(-1, j0)点时, G(jω) =1 ∠ G(jω) = -180o 若z=p-2N中p=0,则G(jω)过(-1,j0)点时, 系统临界稳定,见下图: j 1 -1 G(jω) 特点: G(jω)曲线过(-1, j0)点时, G(jω) =1 同时成立! ∠ G(jω) = -180o
γ h= h 相角裕度 =1 ∠G(jωc) – γ = –180o 1 =180o +∠G(jωc) j G(jωx) G(jωx) -1 γ ωc 幅值裕度 h= G(jωx) 1 ∠G(jωc) G(jω) 相角裕度 =180o +∠G(jωc) γ 稳定裕度的定义续1
= –γ γ=180+ ∠ G(jωc) hdB=-20lg 幅值裕度: c ωx ωc –180o x 相角裕度: 20lg 稳定裕度的定义续2
第六章 线性系统的校正方法 6-1 系统的设计与校正 6-2 串联超前校正 6-3 串联滞后-超前校正
超前校正网络 关键思路:让 ω m= Lc(ωm)=10lga Gc(s)= 1+aTs 1+Ts a﹥1 低频段:1 (0dB) aT 0o 1/aT 1/T 低频段:1 (0dB) 1 aT 1 T 转折频率: 20lga 斜 率: 10lga [+20] [-20] ω=0 ω=∞ 0o +90o 0o -90o φm 0o dφc(ω) dω = 0 令 ωm ωm= 1 aT T = 1 T 关键思路:让 ω m= 得 φm=arcsin a-1 a+1 Lc(ωm)=10lga
例6-3 系统如图,试设计超前校正网络, 使r(t)=t 时
ω=10 Lc(ω)= 20lgb 迟后校正网络 Gc(s)= 1 时 bT c(ω) ≈ -5o~ -9o j 1+bTs 1+Ts ω=0 ω=∞ 低频段: 1 (0dB) 0o +90o 转折频率: 1 bT T 0o -90o 斜 率: [+20] [-20] 0o 1 bT ω=10 1 bT 时 c(ω) ≈ -5o~ -9o j Lc(ω)= 20lgb
例6-4 设计校正网络使图示系统 ω= 2.7时φo(2.7)= –133o OK
滞后-超前校正网络 [-20] [+20] -10lgα -20lgα φm
例6-5 设未校正系统开环传递函数如下,试设计校正网络使: 1)在最大指令速度为180/s时, 位置滞后误差不超过1o; 2) 相角裕度为 45o±3o; 3) 幅值裕度不低于10dB; 4)动态过程调节时间ts不超过3秒。
w w a a=50 w =2 3.5 = w ¢ 例6-5图1 b b ω 采用滞后超前校正 26.8 =100, 取 0dB 20 40 60 80 -20 -40 -60 -80 0.1 1 10 100 ω [-20] 采用滞后超前校正 [-60] 26.8 3.5 b w a =100, a=50 0.5s+1 0.01s+1 =58.25o, 3.5 取 =45o,ts=2.7s, 3.5 c = w ¢ ∴可取 a w =1 由(6-8) ~(6-10)求得 b w 取 =2 降阶 j0(3.5) = -180o L0(3.5)=26.8dB
√ h 例6-5图2 =27.7dB ts=1.65s 3.29 = w ¢ 嘿嘿ok! 42.8o g 180(s+1) G(s) = c = w ¢ g 42.8o h =27.7dB √ ts=1.65s G(s) = 180(s+1) s(s/6+1)(50s+1)(0.01s+1)
第七章 线性离散系统分析 7-1 信号的采样与保持 7-2 z变换 7-3 脉冲传递函数 7-4 离散系统性能
零阶保持器 T=0.2 T=0.4 T=0.8 T=3
[1(t)+t]*=[1(t)]*+[t]* Z域等效变换 [1(t)+t]*=[1(t)]*+[t]* R(s) E(s) E*(s) R(s) E*(s) R(s) B(s) E*(s) B(s) B(s)
采样信号的频谱 ωs ≥ 2ωh T≤π/ωh ωs = 2ωh 或: 滤波器的宽度满足什么 条件时能从 得到 ??! δT(t) = 连续信号的频谱为 ωs=2π/T为采样角频率, ??! 采样信号的频谱为 Cn是傅氏系数,其值为: ωh ωs ≥ 2ωh -ωh 或: T≤π/ωh δT(t) = -3ωs -2ωs -ωs -ωh ωh ωs 2ωs 3ωs ωh -ωh -ωs ωs ωh -ωh ωs 2ωs 3ωs -3ωs -2ωs -ωs ωs = 2ωh
脉冲响应
脉冲响应 2 r(t) 1 2 K(t) 0.03
脉冲响应
G(s) G(z) z Σ Σ 脉冲传递函数的意义 r*(t)=δ(t),c(t)=K(t) r*(t)=Σ r(nT) δ(t-nT) o r*(t)=δ(t-T),c(t)=K(t-T) R(z)= n=0 o Σ r(nT) z -n G(z)是加权脉冲序列的z变换 线性定常离散系统的位移不变性 r*(t)=r(nT)δ(t-nT),c(t)= r(nT)K(t-nT) e*(t)= Σ e(kT)δ(t-kT) k=0 o c*(t)=Σ c(nT) δ(t-nT) n=0 o 根据离散卷积定义得知, 下式右边的Z变换为R(z)K(z) o K[(k-n)T] δ(t-kT) r(nT) k=0 Σ n=0 = δ(t-nT) c*(t) o k=0 c(kT) δ(t-kT) c*(t)=Σ r(nT) K(kT-nT) δ(t-kT) o k=0 ∑ = o K[(k-n)T] δ(t-kT) r(nT) k=0 Σ = c*(t) C(z)=R(z)K(z)
采样拉氏变换的两个重要性质 G*(s)=G*(s+jnωs) [E*(s)G1(s) G2(s)]*=E*(s)[G1(s) G2(s)]* 1)采样函数的拉氏变换具有周期性 G*(s)=G*(s+jnωs) 2)离散信号可从离散符号中提出来 [E*(s)G1(s) G2(s)]*=E*(s)[G1(s) G2(s)]* 设G1(s)G2(s)=G (s), 则有: [E*(s)G(s)]*= ∵E*(s)与∑无关, 所以有: =E*(s)[G(s)]*
闭环实极点分布与相应的动态响应形式 Z平面 Im 1 Re
闭环复极点分布与相应的动态响应形式 Im –1 1 Re
第八章 非线性系统分析 8-1 非线性系统特点 8-2 非线性系统分析的描述函数法
根与相轨迹 λ1 j λ2 j λ2 λ1 节点 不稳定节点 j j 稳定焦点 不稳定节点 j λ1 λ2 j 鞍点 中心
非线性环节的正弦响应 y(t) ωt y(t) ωt ωt y(t) y(t) ωt
X(t) = Asin ωt y(t) ≈ Y1sin(ωt+φ1) 描述函数的定义 ∞ X(t) = Asin ωt y(t)= A0+∑(Ancosnωt+Bnsin nωt) n=1 ∞ y(t) ≈ Y1sin(ωt+φ1) = A0+∑Yn(sin nωt+φn) n=1 非线性环节可近似认为具有和线性环节相类似的频率响应形式 为此,定义正弦信号作用下,非线性环节的稳态输出中一次谐波 分量和输入信号的复数比为非线性环节的描述函数,用N(A)表示: 1 1 1 1 N(A) = N(A) ej∠N(A) = 若A0=0,且当n>1时,yn均很小,则可近似认为非线性环节的 正弦响应仅有一次谐波分量! y(t) A1cos t+B1sin t Y1sin(ωt+φ1) φ1= arctgA1/B1 ω ≈
{ = N(A)= = 死区特性的描述函数 x>△ x <△ x<-△ A> △ X(t)= Asinωt k –△ x(t) y(t) ωt ψ y(t) π - { y(t)= k(x- △) k(x+△) x>△ x <△ x<-△ ψ π - Ψ ω π 2 ψ=△/A, ψ=√1-(△/A)2 cos sin ωt π 2π A x(t) ψ A> △ π-ψ X(t)= Asinωt y(t) ≈ B1sinωt B1+jA1 B1 A = = x(t)=Asinωt N(A)= A 死区特性的描述函数