根据上课 ppt 做的笔记,不定期更新
自动控制指在没有人直接参与的情况下,利用外加的设备或装置(称控制装置或控制器),使机器、设备或生产过程(统称被控对象(plant))的某个工作状态或参数(即被控量)自动地按照预定的规律运行
以反馈控制原理为基础,以传递函数为工具,主要研究单输入、单输出反馈控制系统的理论问题,又常称为反馈控制理论
线性系统满足叠加原理和齐次性(均匀性)
线性系统的数学描述:
线性函数,如线性微分方程、线性差分方程、线性代数方程
在理想情况下,自控系统的被控量和给定值,在任何时候相等,没有误差,也不受干扰的影响,即 c(t)≡r(t)
一般来说整个调节过程分为两个阶段:
系统的稳定性表现为时域响应的收敛性
稳 指系统重新恢复平衡工作状态的能力(稳定性) 和动态过程的振荡倾向(平稳性)
快 指动态过程进行的时间长短(快速性)
准 系统过渡到新的平衡状态后或系统受扰重新恢复平衡之后(过渡过程结束),系统所保持的精度(准确性)
数学模型的类型:动态模型、静态模型
控制系统属于动态系统
时 域:微分方程(差分方程)
复数域:传递函数(脉冲传递函数)
频 域:频率特性
列写元件微分方程的步骤可归纳如下:
建立控制系统的微分方程时
一般先由系统原理线路图画出系统方块图,并分别列写组成系统各元件的微分方程;然后,消去中间变量便得到描述系统输出量与输入量之间关系的微分方程
列写系统各元件的微分方程时
一是应注意信号传送的单向性,即前一个元件的输出是后一个元件的输入,一级一级地单向传送
二是应注意前后连接的两个元件中,后级对前级的负载效应
微分方程的解 = 齐次微分方程的通解 + 非齐次微分方程的任一特解
用拉氏变换法求解线性定常微分方程的过程可归结如下:
将非线性微分方程线性化——小偏差法
实际工作点在某一平衡点 (x0,y0) 附近
非线性 y = f(x) 连续可导
x 在一个很小的范围内变化, Δ x \Delta{x} Δx 很小(小偏差)
可列写线性化增量方程
Δ y = K Δ x , K = d f ( x ) d x ∣ x = x 0 \Delta{y}=K\Delta{x},K=\frac{df(x)}{dx}|_{x=x_0} Δy=KΔx,K=dxdf(x)∣x=x0
线性定常系统传递函数的定义
零初始条件下,系统输出量的拉氏变换与输入量的拉氏变换之比
传递函数是复变量 s 的有理真分式函数,具有复变函数的所有性质;m≤n 且所有系数均为实数
传递函数是一种用系统参数表示输出量与输入量之间关系的表达式它只取决于系统或元件的结构和参数,而与输入量的形式无关,也不反映系统内部的任何信息
传递函数与微分方程具有相通性
传递函数是在零初始条件下定义的。控制系统的零初始条件有两方面的含义:
指输入量是在 t≥0 时才作用于系统,因此,输入量及其各阶导数在 t=0 时均为零;
指输入量加于系统之前,系统处于稳定的工作状态,即输出量及其各阶导数在t =0 时的值也为零,现实的工程控制系统多属此类情况。因此,传递函数可表征控制系统的动态性能,并用以求出在给定输入量时系统的零初始条件响应,即由拉氏变换的卷积定理,有
c ( t ) = L − 1 [ C ( s ) ] = L − 1 [ G ( s ) R ( s ) ] = ∫ 0 t r ( τ ) g ( t − τ ) d τ = ∫ 0 t r ( t − τ ) g ( τ ) d τ c(t)=L^{-1}[C(s)]=L^{-1}[G(s)R(s)]=\int_{0}^{t}r(\tau)g(t-\tau)d\tau =\int_{0}^{t}r(t-\tau)g(\tau)d\tau c(t)=L−1[C(s)]=L−1[G(s)R(s)]=∫0tr(τ)g(t−τ)dτ=∫0tr(t−τ)g(τ)dτ
式中, g ( t ) = L − 1 [ G ( s ) ] g(t)=L^{-1}[G(s)] g(t)=L−1[G(s)]是系统的脉冲响应
传递函数 G(s) 的拉氏反变换是脉冲响应 g(t)
脉冲响应(也称脉冲过渡函数)g(t)是系统在单位脉冲 δ ( t ) \delta(t) δ(t)输入时的输出响应
传递函数的表达形式
一般表达形式
G ( s ) = C ( s ) R ( s ) = b 0 s m + b 1 s m − 1 + . . . + b m − 1 s + b m a 0 s n + a 1 s n − 1 + . . . + a n − 1 s + a n = M ( s ) N ( s ) G(s)=\frac{C(s)}{R(s)}=\frac{b_0s^m+b_1s^{m-1}+...+b_{m-1}s+b_m}{a_0s^n+a_1s^{n-1}+...+a_{n-1}s+a_n}=\frac{M(s)}{N(s)} G(s)=R(s)C(s)=a0sn+a1sn−1+...+an−1s+anb0sm+b1sm−1+...+bm−1s+bm=N(s)M(s)
零-极点表达形式——最高阶次项系数为1
G ( s ) = b 0 ( s − z 1 ) ( s − z 2 ) . . . ( s − z m ) a 0 ( s − p 1 ) ( s − p 2 ) . . . ( s − p n ) = K ∗ ∏ i = 1 m ( s − z i ) ∏ j = 1 n ( s − p j ) G(s)=\frac{b_0(s-z_1)(s-z_2)...(s-z_m)}{a_0(s-p_1)(s-p_2)...(s-p_n)}=K^*\frac{\prod_{i=1}^{m}(s-z_i)}{\prod_{j=1}^{n}(s-p_j)} G(s)=a0(s−p1)(s−p2)...(s−pn)b0(s−z1)(s−z2)...(s−zm)=K∗∏j=1n(s−pj)∏i=1m(s−zi)
系数 K ∗ = b 0 / a 0 K^*=b_0/a_0 K∗=b0/a0称为传递系数或根轨迹增益。这种用零点和极点表示传递函数的方法在根轨迹法中使用较多
在复数平面上表示传递函数的零点和极点的图形,称为传递函数的零极点分布图。在图中一般用“o”表示零点,用“×”表示极点。传递函数零极点分布图可以更形象地反映系统的全面特性。
典型环节形式——常数项系数为1
G ( s ) = b m ( τ 1 s + 1 ) ( τ 2 s + 1 ) . . . ( τ m s + 1 ) a n ( T 1 s + 1 ) ( T 2 s + 1 ) . . . ( T n s + 1 ) = K ∏ i = 1 m ( τ i s + 1 ) ∏ j = 1 n ( T j s + 1 ) G(s)=\frac{b_m(\tau_1s+1)(\tau_2s+1)...(\tau_ms+1)}{a_n(T_1s+1)(T_2s+1)...(T_ns+1)}=K\frac{\prod_{i=1}^{m}(\tau_is+1)}{\prod_{j=1}^{n}(T_js+1)} G(s)=an(T1s+1)(T2s+1)...(Tns+1)bm(τ1s+1)(τ2s+1)...(τms+1)=K∏j=1n(Tjs+1)∏i=1m(τis+1)
K称传递系数或增益,传递函数这种表示形式在频率法中使用较多
传递函数的极点就是微分方程的特征根
极点/特征根决定了所描述系统自由运动(零输入响应)的模态
传递函数的零点并不形成自由运动的模态,但它们却影响各模态响应中所占的比重(系数),因而也影响响应曲线的形状
强迫运动(零初始条件响应)
自由运动(零输入响应),是系统“固有”的成分
典型环节
比例环节
积分环节
惯性环节(非周期环节)
振荡环节
G ( s ) = C ( s ) R ( s ) = w n 2 s 2 + 2 ζ w n s + w n 2 = 1 T 2 s 2 + 2 ζ T s + 1 G(s)=\frac{C(s)}{R(s)}=\frac{{w_n}^2}{s^2+2\zeta{w_n}s+{w_n}^2 }=\frac{1}{T^2s^2+2\zeta{T}s+1} G(s)=R(s)C(s)=s2+2ζwns+wn2wn2=T2s2+2ζTs+11
微分环节
延迟环节
无源网络:为了改善控制系统的性能,常在系统中引入无源网络作为校正元件。无源网络通常由电阻、电容和电感组成
控制系统的结构图是由许多对信号进行单向运算的方框和一些信号流向线组成,它包含四种基本单元:
系统结构图也是控制系统的一种数学模型
建议:第1个方程最好能包含输入,相加点,输出
信号流图与结构图
前向通路:信号从输入节点到输出节点传递时,每个节点只通过一次的通路,叫前向通路
前向通路总增益:前向通路上各支路增益之乘积,称前向通路总增益,一般用 p k p_k pk 表示
回路:起点和终点在同一节点,而且信号通过每一节点不多于一次的闭合通路称为单独回路,简称回路
回路增益:回路中所有支路增益之乘积叫回路增益,用 L a L_a La 表示
不接触回路:回路与回路之间没有公共节点时,这种回路叫做不接触回路。可以有两个或两个以上不接触的回路
不接触回路增益: L b L c L_bL_c LbLc 为所有互不接触的单独回路中,每次取其中2个回路的回路增益的乘积; L d L e L f L_dL_eL_f LdLeLf 为所有互不接触的单独回路中,每次取其中3个回路的回路增益的乘积
梅森增益公式
求取从任意源节点到任意阱节点之间传递函数的梅森增益公式记为:
P = 1 Δ ∑ k = 1 n p k Δ k P=\frac{1}{\Delta}\sum_{k=1}^{n}p_k\Delta_k P=Δ1k=1∑npkΔk
为了研究有用输入作用对系统输出C(s)的影响,需要求有用输入作用下的闭环传递函数C(s)/R(s)。同样,为了研究扰动作用N(s)对系统输出C(s)的影响,也需要求取扰动作用下的闭环传递函数C(s)/N(s)
此外,在控制系统的分析和设计中,还常用到在输入信号 R(s) 或扰动 N(s) 作用下,以误差信号 E(s) 作为输出量的闭环误差传递函数 E(s)/R(s) 或E(s)/N(s)
在一定条件下,系统的输出只取决于反馈通路传递函数 H(s) 及输入信号R(s),既与前向通路传递函数无关,也不受扰动作用的影响。特别是当H(s)=1,即单位反馈时,C(s) ≈ R(s),从而近似实现了对输入信号的完全复现,且对扰动具有较强的抑制能力。
绝不允许将各种闭环传递函数进行叠加后求其输出响应
线性系统的稳定性仅取决于系统自身的固有特性,而与外界条件无关。因此,设线性系统在初始条件为零时,作用一个理想单位脉冲,这时系统的输出增量为脉冲响应c(t)。这相当于系统在扰动信号作用下,输出信号偏离原平衡工作点的问题。若 t → ∞ t\to\infty t→∞ 时,脉冲响应:
lim t → ∞ c ( t ) = 0 \lim_{t\to\infty}c(t)=0 t→∞limc(t)=0
即输出增量收敛于原平衡工作点,则线性系统是稳定的
根据李雅普诺夫稳定性理论,线性控制系统的稳定性可叙述如下:
若线性控制系统在初始扰动的影响下,其动态过程随时间的推移逐渐衰减并趋于零(原平衡工作点),则称系统渐近稳定,简称稳定;反之,若在初始扰动影响下,系统的动态过程随时间的推移而发散,则称系统不稳定
线性系统稳定的充分必要条件:闭环系统特征方程的所有根均具有负实部;或者说,闭环传递函数的极点均严格位于左半s平面
若特征根中具有一个或一个以上零实部根,而其余的特征根均具有负实部,则脉冲响应c(t)趋于常数,或趋于等幅正弦振荡,按照稳定性定义,此时系统不是渐近稳定的,处于稳定和不稳定的临界状态,常称为临界稳定情况
在经典控制理论中,只有渐近稳定的系统才称为稳定系统;否则,称为不稳定系统
线性系统稳定的充分且必要条件是:劳思表中第一列各值为正
如果劳思表第一列中出现小于零的数值,系统不稳定,且第一列各系数符号的改变次数,代表特征方程的正实部根的数目
只有当系统稳定时,研究稳态误差才有意义
把在阶跃函数作用下没有原理性稳态误差的系统,称为无差系统;而把具有原理性稳态误差的系统,称为有差系统
非线性因素所引起的系统稳态误差,则称为附加稳态误差,或结构性稳态误差
稳态——时间t趋于无穷大时系统的输出状态
误差——系统实际输出与期望输出之差
稳态误差——对于稳定的系统,暂态响应随时间的增长而衰减,当衰减到可以忽略的程度时,则希望的稳态响应与实际的稳态响应之差称为稳态误差。稳态误差是系统达到稳态时,系统精度的度量
稳态误差定义
e s s ( ∞ ) = lim t → ∞ ε ( t ) = lim t → ∞ ( c 0 ( t ) − c ( t ) ) e_{ss}(\infty)=\lim_{t\to\infty}\varepsilon(t)=\lim_{t\to\infty}(c_0(t)-c(t)) ess(∞)=t→∞limε(t)=t→∞lim(c0(t)−c(t))
因系统的输入可分为给定和扰动输入,所以稳态误差也可分为给定稳态误差和扰动稳态误差
误差:是从系统输出端来定义,它定义为系统输出量的希望值与实际值之差:
ε ( s ) = C 0 ( s ) − C ( s ) \varepsilon(s)=C_0(s)-C(s) ε(s)=C0(s)−C(s)
偏差:是在系统输入端来定义,定义为系统给定量与反馈量之差:
E ( s ) = R ( s ) − B ( s ) = H ( s ) ε ( s ) E(s)=R(s)-B(s)=H(s)\varepsilon(s) E(s)=R(s)−B(s)=H(s)ε(s)
偏差值与误差值的差别为H(s)
对于单位反馈系统,偏差与误差的值是一样的
一般地用偏差代替误差
ε s s ( ∞ ) = e s s ( ∞ ) = lim t → ∞ e ( t ) = lim t → ∞ ( r ( t ) − b ( t ) ) \varepsilon_{ss}(\infty)=e_{ss}(\infty)=\lim_{t\to\infty}e(t)=\lim_{t\to\infty}(r(t)-b(t)) εss(∞)=ess(∞)=t→∞lime(t)=t→∞lim(r(t)−b(t))
如果有理函数sE(s)除在原点处有惟一的极点外,在s右半平面及虚轴上处处解析(即不存在极点),即sE(s)的极点均位于s左半平面(包括坐标原点),则可根据拉氏变换的终值定理,方便地求出系统的稳态误差:
e s s ( ∞ ) = lim t → ∞ e ( t ) = lim s → 0 s E ( s ) e_{ss}(\infty)=\lim_{t\to\infty}e(t)=\lim_{s\to0}sE(s) ess(∞)=t→∞lime(t)=s→0limsE(s)
上式算出的稳态误差是误差信号e(t)的稳态分量 e s s ( t ) e_{ss}(t) ess(t)在t趋于无穷时的数值,故有时称为终值误差或稳态误差终值,它不能反映e(t)随时间 t 的完整变化规律,具有一定的局限性
当输入信号形式一定时,系统是否存在稳态误差就取决于开环传递函数描述的系统结构
在一般情况下,分子阶次为m,分母阶次为n的开环传递函数 G k G_k Gk可表示为:
G ( s ) = G ( s ) H ( s ) = k ∏ i = 1 m ( τ i s + 1 ) s v ∏ j = 1 n ( T j s + 1 ) ( n ≥ m ) G(s)=G(s)H(s)=\frac{k\prod_{i=1}^{m}(\tau_is+1)}{s^v\prod_{j=1}^{n}(T_js+1)}(n\ge{m}) G(s)=G(s)H(s)=sv∏j=1n(Tjs+1)k∏i=1m(τis+1)(n≥m)
式中,K为开环增益; τ i \tau_i τi 和 T j T_j Tj 为时间常数; v v v为开环系统在s平面坐标原点上的极点的重数。现在的分类方法是以 v v v的数值来划分的: v v v =0,称为0型系统; v v v =1,称为I型系统; v v v =2,称为II型系统。当 v v v >2时,除复合控制系统外,使系统稳定是相当困难的。因此除航天控制系统外,III型及III型以上的系统几乎不采用。
这种以开环系统在 s 平面坐标原点上的极点数来分类的方法,其优点在于:可以根据已知的输入信号形式迅速判断系统是否存在原理性稳态误差及稳态误差的大小。它与按系统的阶次进行分类的方法不同,阶次m与n的大小与系统的型别无关,且不影响稳态误差的数值
系统的静态位置误差系数 K p K_p Kp: k p = lim s → 0 G k ( s ) k_p=\lim_{s\to0}G_k(s) kp=lims→0Gk(s)
系统的静态速度误差系数 K v K_v Kv: k v = lim s → 0 s G k ( s ) k_v=\lim_{s\to0}sG_k(s) kv=lims→0sGk(s)
系统的静态加速度误差系数 K a K_a Ka: k a = lim s → 0 s 2 G k ( s ) k_a=\lim_{s\to0}s^2G_k(s) ka=lims→0s2Gk(s)
如果系统为非单位反馈系统,其 H ( s ) = K h H(s)=K_h H(s)=Kh为常数,那么系统输出量的希望值为 R ’ ( s ) = R ( s ) / K h R’(s)=R(s)/K_h R’(s)=R(s)/Kh,系统输出端的稳态位置误差为: e s s ′ = e s s / K h e'_{ss}=e_{ss}/K_h ess′=ess/Kh
给定稳态误差级数 (动态误差系数)
利用动态误差系数法,可以研究输入信号几乎为任意时间函数时的系统稳态误差变化,因此动态误差系数又称广义误差系数。为了求取动态误差系数,写出稳态误差信号的拉氏变换式: E ( s ) = Φ e ( s ) R ( s ) E(s)=\Phi_e(s)R(s) E(s)=Φe(s)R(s)
将误差传递函数在 Φ e ( s ) = 0 \Phi_e(s)=0 Φe(s)=0的邻域内展成泰勒级数,得:
Φ e ( s ) = 1 1 + G ( s ) H ( s ) = Φ e ( 0 ) + Φ ˙ e ( 0 ) s + 1 2 ! Φ ¨ e ( 0 ) s 2 + ⋯ \Phi_e(s)=\frac{1}{1+G(s)H(s)}=\Phi_e(0)+\dot\Phi_e(0)s+\frac{1}{2!}\ddot\Phi_e(0)s^2+\cdots Φe(s)=1+G(s)H(s)1=Φe(0)+Φ˙e(0)s+2!1Φ¨e(0)s2+⋯
于是
E ( s ) = Φ e ( 0 ) R ( s ) + Φ ˙ e ( 0 ) s R ( s ) + 1 2 ! Φ ¨ e ( 0 ) s 2 R ( s ) + ⋯ + 1 l ! Φ e ( l ) ( 0 ) s l R ( s ) + ⋯ E(s)=\Phi_e(0)R(s)+\dot\Phi_e(0)sR(s)+\frac{1}{2!}\ddot\Phi_e(0)s^2R(s)+\cdots+\frac{1}{l!}\Phi_e^{(l)}(0)s^lR(s)+\cdots E(s)=Φe(0)R(s)+Φ˙e(0)sR(s)+2!1Φ¨e(0)s2R(s)+⋯+l!1Φe(l)(0)slR(s)+⋯
上述无穷级数收敛于s=0的邻域,称为给定输入作用下的稳态误差级数,相当于在时间域内 t → ∞ t\to\infty t→∞时成立。因此,当所有初始条件均为零时,对上式进行拉氏反变换,就得到作为时间函数的稳态误差表达式:
e s s ( t ) = ∑ i = 0 ∞ C i r ( i ) ( t ) C i = 1 i ! Φ e ( i ) ( 0 ) i = 0 , 1 , 2 , ⋯ e_{ss}(t)=\sum_{i=0}^{\infty}C_ir^{(i)}(t)\\ C_i=\frac{1}{i!}\Phi_e^{(i)}(0)\\ i=0,1,2,\cdots ess(t)=i=0∑∞Cir(i)(t)Ci=i!1Φe(i)(0)i=0,1,2,⋯
为了 减小或消除系统在输入信号和扰动作用下的稳态误差,可以采取以下措施:
增大系统开环增益或扰动作用点之前系统的前向通道增益
增大系统开环增益K以后,对于0型系统,可以减小系统在阶跃输入时的位置误差;对于I型系统,可以减小系统在斜坡输入时的速度误差;对于Ⅱ型系统,可以减小系统在加速度输入时的加速度误差
增大扰动点之后系统的前向通道增益,不能改变系统对扰动的稳态误差数值
在系统的前向通道或主反馈通道设置串联积分环节
扰动作用点之前的前向通道积分环节数与主反馈通道积分环节数之和决定系统响应扰动作用的型别,该型别与扰动作用点之后前向通道的积分环节数无关
反馈控制系统中,设置串联积分环节或增大开环增益以消除或减小稳态误差的措施,必然导致降低系统的稳定性,甚至造成系统不稳定,从而恶化系统的动态性能
采用串级控制抑制内回路扰动
当控制系统中存在多个扰动信号,且控制精度要求较高时,宜采用串级控制方式,可以显著抑制内回路的扰动影响
采用复合控制方法
复合控制系统是在系统的反馈控制回路中加入前馈通路,组成一个前馈控制与反馈控制相结合的系统
线性定常系统在任意输入信号作用下的时间响应,均可分解为零输入响应与零状态响应之和
y ( t ) = y z i + y z s y(t)=y_{zi}+y_{zs} y(t)=yzi+yzs
y z i y_{zi} yzi:零输入响应
只与系统的极点类型和分布有关,与系统的零点无关,与系统的初始条件有关
反映系统的稳定性
y z s y_{zs} yzs:零状态响应
微分方程 T d c ( t ) d t + c ( t ) = r ( t ) T\frac{dc(t)}{dt}+c(t)=r(t) Tdtdc(t)+c(t)=r(t)
开环传函 G k ( s ) = 1 T s G_k(s)=\frac{1}{Ts} Gk(s)=Ts1
闭环传函 Φ ( s ) = C ( s ) R ( s ) = 1 T s + 1 \Phi(s)=\frac{C(s)}{R(s)}=\frac{1}{Ts+1} Φ(s)=R(s)C(s)=Ts+11
特征方程 T s + 1 = 0 Ts+1=0 Ts+1=0
单位阶跃响应
t = 3 T , c ( t ) = 0.95 t=3T,c(t)=0.95 t=3T,c(t)=0.95 —— 5%误差带
t = 4 T , c ( t ) = 0.98 t=4T,c(t)=0.98 t=4T,c(t)=0.98 —— 2%误差带
可以用时间常数T去度量系统的输出量的数值
初始斜率为1/T
无超调
稳态误差 e s s = 0 e_{ss}=0 ess=0
单位脉冲响应
单位斜坡响应
出现稳态误差( e s s = T e_{ss}=T ess=T)
一阶系统不能实现对加速度输入函数的跟踪
T 2 d 2 c ( t ) d t 2 + 2 ζ T d c ( t ) d t + c ( t ) = r ( t ) d 2 c ( t ) d t 2 + 2 ζ w n d c ( t ) d t + w n 2 c ( t ) = w n 2 r ( t ) T^2\frac{d^2c(t)}{dt^2}+2\zeta{T}\frac{dc(t)}{dt}+c(t)=r(t)\\ \frac{d^2c(t)}{dt^2}+2\zeta{w_n}\frac{dc(t)}{dt}+w^2_nc(t)=w^2_nr(t) T2dt2d2c(t)+2ζTdtdc(t)+c(t)=r(t)dt2d2c(t)+2ζwndtdc(t)+wn2c(t)=wn2r(t)
G k ( s ) = 1 T 2 s 2 + 2 ζ T s = w n 2 s 2 + 2 ζ w n s G_k(s)=\frac{1}{T^2s^2+2\zeta{T}s}=\frac{w^2_n}{s^2+2\zeta{w_n}s} Gk(s)=T2s2+2ζTs1=s2+2ζwnswn2
Φ ( s ) = C ( s ) R ( s ) = w n 2 s 2 + 2 ζ w n s + w n 2 = 1 T 2 s 2 + 2 ζ T s + 1 \Phi(s)=\frac{C(s)}{R(s)}=\frac{{w_n}^2}{s^2+2\zeta{w_n}s+{w_n}^2 }=\frac{1}{T^2s^2+2\zeta{T}s+1} Φ(s)=R(s)C(s)=s2+2ζwns+wn2wn2=T2s2+2ζTs+11
s 1 , 2 = − ζ w n ∓ j w n 1 − ζ 2 s_{1,2}=-\zeta{w_n}\mp{jw_n\sqrt{1-\zeta^2}} s1,2=−ζwn∓jwn1−ζ2
ζ=0,无(零)阻尼
0<ζ<1,欠阻尼
ζ=1,临界阻尼(重极点)
ζ>1,过阻尼
由于阻尼比ξ为负,指数因子具有正幂指数,因此系统的动态过程为发散正弦振荡或单调发散的形式,从而表明 ζ < 0 \zeta<0 ζ<0的二阶系统是不稳定的
如果ξ=0,则特征方程有一对纯虚根, s 1 , 2 = ± j w n s_{1,2}=\pm{jw_n} s1,2=±jwn,对应于 s 平面虚轴上一对共轭极点,可以算出系统的阶跃响应为等幅振荡,此时系统相当于无阻尼情况
如果 0 < ζ < 1 0<\zeta<1 0<ζ<1,则特征方程有一对具有负实部的共轭复根, s 1 , 2 = − ζ w n ∓ j w n 1 − ζ 2 s_{1,2}=-\zeta{w_n}\mp{jw_n\sqrt{1-\zeta^2}} s1,2=−ζwn∓jwn1−ζ2 ,对应于s平面左半部的共轭复数极点,相应的阶跃响应为衰减振荡过程,此时系统处于欠阻尼情况
如果ξ=1,则特征方程具有两个相等的负实根, s 1 , 2 = − w n s_{1,2}=-w_n s1,2=−wn,对应于s平面负实轴上的两个相等实极点,相应的阶跃响应非周期地趋于稳态输出,此时系统处于临界阻尼情况。
如果ξ>1,则特征方程有两个不相等的负实根 s 1 , 2 = − ζ w n ± w n ζ 2 − 1 s_{1,2}=-\zeta{w_n}\pm{w_n\sqrt{\zeta^2-1}} s1,2=−ζwn±wnζ2−1 ,对应于s平面负实轴上的两个不等实极点,相应的单位阶跃响应也是非周期地趋于稳态输出,但响应速度比 1 临界阻尼情况缓慢,因此称为过阻尼情况
二阶控制系统的设计,一般取ξ=0.4~0.8
二阶系统的单位阶跃响应
根的实部反应在响应的指数部分(模态形状:是否衰减)
根的虚部反应在响应的正弦部分(模态形状:是否振荡)
ζ = 0 \zeta=0 ζ=0 临界稳定状态
ξ≈0.707 时调节时间最短(称为最佳阻尼比)
s 1 , 2 = − ζ w n ± j w n 1 − ζ 2 = − σ ± j w d s_{1,2}=-\zeta{w_n}\pm{jw_n\sqrt{1-\zeta^2}}=-\sigma\pm{jw_d} s1,2=−ζwn±jwn1−ζ2=−σ±jwd
衰减系数 σ \sigma σ(实部)是闭环极点到虚轴之间的距离
阻尼振荡频率 w d w_d wd (虚部)是闭环极点到实轴之间的距离
w d = w n ( 1 − ζ 2 ) w_d=w_n(\sqrt{1-\zeta^2}) wd=wn(1−ζ2)
自然频率 w n w_n wn是闭环极点到坐标原点之间的距离
w n w_n wn与负实轴夹角的余弦正好是阻尼比ξ,即 cos φ = ζ \cos\varphi=\zeta cosφ=ζ
结论:调节时间与闭环极点的实部数值成反比。闭环极点距虚轴的距离越远,系统的调节时间越短。由于阻尼比值主要根据对系统超调量的要求来确定,所以调节时间主要由自然频率决定。若能保持阻尼比值不变而加大自然频率值,则可在不改变超调量的情况下缩短调节时间
欠阻尼二阶系统阶跃响应
峰值时间 t p = π w d = π w n 1 − ζ 2 t_p=\frac{\pi}{w_d}=\frac{\pi}{w_n\sqrt{1-\zeta^2}} tp=wdπ=wn1−ζ2π 与根(极点)的虚部成反比
超调量 σ p = e − π ζ 1 − ζ 2 × 100 % \sigma_p=e^{-\frac{\pi\zeta}{\sqrt{1-\zeta^2}}}\times100\% σp=e−1−ζ2πζ×100%
上升时间 t r = π − φ w n 1 − ζ 2 t_r=\frac{\pi-\varphi}{w_n\sqrt{1-\zeta^2}} tr=wn1−ζ2π−φ 第一次达到稳态值时的时间
调节时间 t s = 3.5 ζ w n ( Δ = 5 % ) t_s=\frac{3.5}{\zeta{w_n}}(\Delta=5\%) ts=ζwn3.5(Δ=5%) t s = 4.4 ζ w n ( Δ = 2 % ) t_s=\frac{4.4}{\zeta{w_n}}(\Delta=2\%) ts=ζwn4.4(Δ=2%)
调节时间与闭环极点的实部数值成反比
延迟时间 t d ≈ 1 + 0.7 ζ w n t_d\approx\frac{1+0.7\zeta}{w_n} td≈wn1+0.7ζ 或 t d ≈ 1 + 0.6 ζ + 0.2 ζ 2 w n t_d\approx\frac{1+0.6\zeta+0.2\zeta^2}{w_n} td≈wn1+0.6ζ+0.2ζ2 第一次达到稳态值一半的时间
延迟时间与震荡频率成反比,与ξ成正比
由于阻尼比值主要根据对系统超调量的要求来确定,所以调节时间主要由自然频率决定
对于欠阻尼二阶系统,极点的阻尼角(阻尼比)决定响应的平稳性;阻尼比(阻尼角)一定时,极点与虚轴的距离( ζ w n \zeta{w_n} ζwn)决定响应的快速性
欠阻尼二阶系统单位阶跃响应稳态误差为零
添加零点对典型二阶系统暂态响应特性影响
添加闭环零点
添加开环零点增大阻尼比,不影响系统的稳定性,系统的动态性能可以得到改善,不改变系统的稳态精度
输出量的导数同样可以用来改善系统的性能。通过将输出的速度信号反馈到系统输入端,并与误差信号比较,其效果与比例一微分控制相似,可以增大系统阻尼,改善系统动态性能
比例-微分控制与测速反馈控制的比较
根轨迹简称根迹,它是开环系统某一参数从0变到 ∞ \infty ∞时,闭环系统特征方程式的根在 s 平面上变化的轨迹。
根轨迹的分类:
180度根轨迹 、参数根轨迹、零度根轨迹、根轨迹簇
当闭环系统没有零点与极点相消时,闭环特征方程式的根就是闭环传递函数的极点,我们常简称为闭环极点
从已知的开环零、极点位置及某一变化的参数来求取闭环极点的分布,实际上就是解决闭环特征方程式的求根问题
因为系统的稳定性由系统闭环极点唯一确定,而系统的稳态性能和动态性能又与闭环零、极点在 s 平面上的位置密切相关,所以根轨迹图不仅可以直接给出闭环系统时间响应的全部信息,而且可以指明开环零、极点应该怎样变化才能满足给定的闭环系统的性能指标要求
根轨迹法应用前提条件是受控系统的传递函数没有零极点相消
根轨迹是系统所有闭环极点的集合
根轨迹方程 K ∗ ∏ j = 1 m ( s − z j ) ∏ i = 1 n ( s − p i ) = − 1 \text{根轨迹方程 }K^*\frac{\prod_{j=1}^{m}(s-z_j)}{\prod_{i=1}^{n}(s-p_i)}=-1 根轨迹方程 K∗∏i=1n(s−pi)∏j=1m(s−zj)=−1
模值条件:
K ∗ = ∏ i = 1 n ∣ s − p i ∣ ∏ j = 1 m ∣ s − z j ∣ K^*=\frac{\prod_{i=1}^{n}|s-p_i|}{\prod_{j=1}^{m}|s-z_j|} K∗=∏j=1m∣s−zj∣∏i=1n∣s−pi∣
相角条件:180度相角差
由于闭环特征方程中的某些系数是根轨迹增益 K ∗ K^* K∗ 的函数,所以当 K ∗ K^* K∗ 从零到无穷大连续变化时,特征方程的某些系数也随之而连续变化,因而特征方程式根的变化也必然是连续的,故根轨迹具有连续性
根轨迹必对称于实轴的原因是显然的,因为闭环特征方程式的根只有实根和复根两种,实根位于实轴上,复根必共轭,而根轨迹是根的集合,因此根轨迹对称于实轴
根据对称性,只需做出上半 s 平面的根轨迹部分,然后利对称关系就可以画出下半 s 平面的根轨迹部分
一般情况下,常见的根轨迹分离点是位于实轴上的两条根轨迹分支的分离点
如果根轨迹位于实轴上两个相邻的开环极点之间,其中一个可以是无限极点,则在这两个极点之间至少存在一个分离点
同样,如果根轨迹位于实轴上两个相邻的开环零点之间,其中一个可以是无限零点,则在这两个零点之间也至少有一个分离点
对于特定 K ∗ K^* K∗值下的闭环极点,可用模值条件确定
在控制系统中,除根轨迹增益 K ∗ K^* K∗ 以外,其他情形下的根轨迹统称为广义根轨迹。包括:
通常,将负反馈系统中 K ∗ K^* K∗由零变到无穷大时的根轨迹叫做常规根轨迹(或180度根轨迹)
以非开环增益/根轨迹增益为可变参数绘制的根轨迹称为参数根轨迹,以区别于以开环增益K/根轨迹增益 K ∗ K^* K∗为可变参数的常规根轨迹
一般说来,零度根轨迹的来源有两个方面:
如果系统的所有开环极点和零点都位于s平面的左半平面,则称为最小相位系统。
若系统的开环传递函数在右半S平面有零点或极点,则该系统称为非最小相位系统
闭环零点对系统动态性能的影响,相当于减小闭环系统的阻尼,从而使系统的过渡过程有出现超调的趋势,并且这种作用将随闭环零点接近坐标原点的程度而加强
在开环系统中增加极点,会使根轨迹向右移动,从而降低系统的相对稳定性,增加系统响应的调节时间
如果闭环零、极点相距很近,那么这样的闭环零、极点常称为偶极子。偶极子有实数偶极子和复数偶极子之分,而复数偶极子必共轭出现
不难看出,只要偶极子不十分接近坐标原点,它们对系统动态性能的影响就甚微,从而可以忽略它们的存在
闭环零点的存在,将使系统的峰值时间提前,这相当于减小闭环系统的阻尼,从而使超调量加大,当闭环零点接近坐标原点时,这种作用尤甚
一般说来,闭环零点对调节时间的影响是不定的
闭环实数主导极点对系统性能的影响是:闭环实数主导极点的作用,相当于增大系统的阻尼,使峰值时间迟后,超调量下降
如果实数极点比共轭复数极点更接近坐标原点,甚至可以使振荡过程变为非振荡过程
控制系统中的信号可以表示为不同频率正弦信号的合成。控制系统频率特性反映正弦信号作用下系统响应的性能。应用频率特性研究线性系统的经典方法称为频域分析法。频域分析法具有以下特点:
在正弦信号作用下,输出稳态分量的幅值和相位的变化,称为幅值比 A ( w ) A(w) A(w) 和相位差 φ ( w ) \varphi(w) φ(w) , 且皆为输入正弦信号频率 w w w 的函数
{ A ( w ) = ∣ G ( j w ) ∣ φ ( w ) = ∠ G ( j w ) \left\{\begin{matrix} A(w)=\left | {G(jw)} \right | \\ \varphi(w)=\angle{G(jw)} \end{matrix}\right. {A(w)=∣G(jw)∣φ(w)=∠G(jw)
表明,对于稳定的线性定常系统,由谐波输入产生的输出稳态分量仍然是与输入同频率的谐波函数,而幅值和相位的变化是频率 w w w 的函数,且与系统数学模型相关。
为此,定义谐波输入下,输出响应中与输入同频率的谐波分量与谐波输入的幅值之比 A ( w ) A(w) A(w) 为幅频特性,相位之差 φ ( w ) \varphi(w) φ(w) 为相频特性,并称其指数表达形式:
G ( j w ) = A ( w ) e j φ ( w ) G(jw)=A(w)e^{j\varphi(w)} G(jw)=A(w)ejφ(w)
为系统的频率特性
频率特性的定义既可以适用于稳定系统,也可适用于不稳定系统。稳定系统的频率特性可以用实验方法确定,即在系统的输入端施加不同频率的正弦信号,然后测量系统输出的稳态响应,再根据幅值比和相位差作出系统的频率特性曲线。频率特性也是系统数学模型的一种表达形式
对于不稳定系统,输出响应稳态分量中含有由系统传递函数的不稳定极点产生的呈发散或振荡的分量,所以不稳定系统的频率特性不能通过实验方法确定
频率特性定义(物理意义)
稳定系统的频率特性为零初始条件下,输出和输入的傅氏变换之比
G ( j w ) = C ( j w ) R ( j w ) = G ( s ) ∣ s = j w G(jw)=\frac{C(jw)}{R(jw)}=G(s)|_{s=jw} G(jw)=R(jw)C(jw)=G(s)∣s=jw
在工程分析和设计中,通常把线性系统的频率特性画成曲线,再运用图解法进行研究。常用的频率特性曲线有以下三种:
幅相频率特性曲线
又称为幅相曲线或极坐标图或奈奎斯特曲线(Nyquist Plot)
以横轴为实轴、纵轴为虚轴,构成复数平面。对于任一给定的频率 w w w , 频率特性值为复数。若将频率特性表示为实数和虚数的形式,则实部为实轴坐标值,虚部为虚轴坐标值
若将频率特性表示为复指数形式,则为复平面上的向量,而向量的长度为频率特性的幅值,向量与实轴正方向的夹角等于频率特性的相位
由于幅频特性为 w w w 的偶数,相频特性为 w w w 的奇函数,则 w w w 从零变化至 + ∞ +\infty +∞ 和 w w w 从零变化至 − ∞ -\infty −∞ 的幅相曲线关于实轴对称,因此一般只绘制 w w w 从零变化至 + ∞ +\infty +∞ 的幅相曲线
在系统幅相曲线中,频率 w w w 为参变量,一般用小箭头表示 w w w 增大时幅相曲线的变化方向
对数频率特性曲线
又称为伯德曲线或伯德图(Bode Diagram)
对数频率特性曲线由对数幅频曲线和对数相频曲线组成
对数频率特性曲线的横坐标按 lg w \lg{w} lgw 分度,单位为弧度/秒(rad/s)
对数幅频曲线的纵坐标按:
对数幅频特性: L ( w ) = 20 lg ∣ G ( j w ) ∣ = 20 lg A ( w ) \text{对数幅频特性:}L(w)=20\lg{|G(jw)|}=20\lg{A(w)} 对数幅频特性:L(w)=20lg∣G(jw)∣=20lgA(w)
线性分度,单位是分贝(dB)
对数相频曲线的纵坐标按 φ ( w ) \varphi(w) φ(w) 线性分度,单位为度 ( ∘ ^\circ ∘)。由此构成的坐标系称为半对数坐标系
对数幅相曲线
又称尼科尔斯曲线或尼科尔斯图(Nichols Plot)
其特点是纵坐标为 L ( w ) L(w) L(w) ,单位为分贝 (dB);横坐标为 φ ( w ) \varphi(w) φ(w) ,单位为度 ( ∘ ^\circ ∘) ,均为线性分度,频率 w w w 为参变量
在尼科尔斯曲线对应的坐标系中,可以根据系统开环和闭环的关系,绘制关于闭环幅频特性的等 M 簇线和闭环相频特性的等 α \alpha α 簇线,因而根据频域指标要求确定校正网络,简化系统的设计过程
由于开环传递函数的分子和分母多项式的系数皆为实数,因此系统开环零极点或为实数或为共轭复数
根据开环零极点可将分子和分母多项式分解成因式,再将因式分类,即得典型环节
典型环节可分为两类:一类为最小相位环节;另一类为非最小相位环节
最小相位环节有下列七种:
比例环节 K (K >0)
对数幅频特性曲线
平行于频率轴的一条直线,值为 20 lg k 20\lg{k} 20lgk
对数相频特性曲线
为零度直线
惯性环节 1/(Ts+1) (T >0)
一阶微分环节 Ts+1 (T >0)
振荡环节 1 / ( s 2 / w n 2 + 2 ζ s / w n + 1 ) 1/(s^2/{w_n}^2+2\zeta{s}/w_n+1) 1/(s2/wn2+2ζs/wn+1) ( w n > 0 , 0 ≤ ζ < 1 ) (w_n>0,0\le\zeta<1) (wn>0,0≤ζ<1)
谐振频率: w r = w n 1 − 2 ζ 2 , 0 < ζ ≤ 2 2 w_r=w_n\sqrt{1-2{\zeta}^2},0<\zeta\le\frac{\sqrt{2}}{2} wr=wn1−2ζ2,0<ζ≤22
谐振峰值: M r = 1 2 ζ 1 − ζ 2 , 0 < ζ ≤ 2 2 M_r=\frac{1}{2\zeta\sqrt{1-\zeta^2}},0<\zeta\le\frac{\sqrt{2}}{2} Mr=2ζ1−ζ21,0<ζ≤22
由于 L a ( w ) L_a(w) La(w)与阻尼比 ζ \zeta ζ 无关,用渐近线近似表示对数幅频曲线存在误差,误差的大小不仅和 w w w 有关,而且也和 ζ \zeta ζ 有关,误差曲线为一曲线簇,根据误差曲线可以修正渐近特性曲线而获得准确曲线。
根据对数幅频特性定义还可知:
(1)非最小相位振荡环节与振荡环节的对数幅频渐近特性曲线相同
(2)二阶微分环节和非最小相位二阶微分环节与振荡环节的对数幅频渐近特性曲线关于0dB线对称
Nyquist 曲线
二阶微分环节 s 2 / w n 2 + 2 ζ s / w n + 1 s^2/{w_n}^2+2\zeta{s}/w_n+1 s2/wn2+2ζs/wn+1 ( w n > 0 , 0 ≤ ζ < 1 ) (w_n>0,0\le\zeta<1) (wn>0,0≤ζ<1)
积分环节 1/s
微分环节 s
延迟环节与延迟系统
非最小相位环节共有五种:
除了比例环节外,非最小相位环节和与之相对应的最小相位环节的区别在于开环零极点的位置
已知系统开环传递函数 G k ( s ) G_k(s) Gk(s),绘制开环辐相曲线
F ( s ) = ( s − z 1 ) ( s − z 2 ) … ( s − z m ) ( s − p 1 ) ( s − p 2 ) … ( s − p n ) = ∑ j = 1 m ( s − z j ) ∑ i = 1 n ( s − p i ) F(s)=\frac{(s-z_1)(s-z_2)\dots(s-z_m)}{(s-p_1)(s-p_2)\dots(s-p_n)}=\frac{\sum_{j=1}^{m}(s-z_j)}{\sum_{i=1}^{n}(s-p_i)} F(s)=(s−p1)(s−p2)…(s−pn)(s−z1)(s−z2)…(s−zm)=∑i=1n(s−pi)∑j=1m(s−zj)
如果 F(s) 是非奇异的(除有限个孤立奇点 p 1 p_1 p1 、 p 2 p_2 p2 … p n p_n pn 外处处解析),则有:
设s平面上不通过 F(s) 的任一零点和极点的封闭曲线 Γ \Gamma Γ 包围 F(s) 的P个极点和Z个零点,则 s 从封闭曲线 Γ \Gamma Γ上任一点 A 起,顺时针沿 Γ \Gamma Γ 移动一周,再回到 A 点时,
则相应地, 在F(s) 平面上,封闭曲线 Γ \Gamma Γ 的象或映射亦从点 F(A) 起,逆时针方向绕原点 R = P − Z R=P-Z R=P−Z 圈,回到 F(A) 点止,形成一条闭合曲线 Γ F \Gamma_F ΓF
封闭曲线 Γ \Gamma Γ 及其象或映射 Γ F \Gamma_F ΓF 的形状对于R值的确定没有影响,也即,R值只与封闭曲线 Γ \Gamma Γ 包围F(s)零点、极点的个数有关,而与封闭曲线 Γ \Gamma Γ 的形状无关
F(s)具有以下特点:
系统的闭环稳定性取决于系统闭环传递函数极点即F(s)的零点的位置
考虑选取 s 平面上可以围 住其整个右半平面的封闭曲线 Γ \Gamma Γ (称为N氏围线),若N氏围线 Γ \Gamma Γ 没有围住F(s)的零点,则闭环系统稳定;若围住了1个零点,系统就不稳定
设Z为F(s)零点在右半s平面中的个数,则Z =0 表示F(s)在右半s平面无零点,闭环特征方程F(s)=0的根全部位于s平面的左半平面,系统稳定;Z≠0,系统不稳定
设P为F(s)极点在右半s平面中的个数, 即为系统开环极点在右半s平面的个数;若已知 G k ( s ) G_k(s) Gk(s),则P已知
原始的Nyquist稳定判据步骤:
已知: G k ( s ) G_k(s) Gk(s)
则P已知(P为 G k ( s ) G_k(s) Gk(s) 在右半s平面上开环极点的个数)
在F(s) 平面上画出 F ( s ) = 1 + G k ( s ) F(s)=1+G_k(s) F(s)=1+Gk(s) 曲线
则可知R(当s在s 平面N氏围线 Γ \Gamma Γ 上顺时针方向移动一周时,F(s)平面上 F ( s ) = 1 + G k ( s ) F(s)=1+G_k(s) F(s)=1+Gk(s) 曲线逆时针绕原点的圈数)
在自变量s平面(N氏围线 Γ \Gamma Γ) →因变量F(s) 平面(象或映射 Γ F \Gamma_F ΓF )上直接应用幅角原理,由R=P-Z,可得:Z=P-R
若 Z =0 ,则系统稳定
若 Z≠0,则系统不稳定
F(s)平面上的原点,就是G(s)H(s)平面上的(-1,j0)点
封闭曲线 Γ F \Gamma_F ΓF 包围 F(s) 平面原点的圈数,就等于封闭曲线 Γ G ( s ) H ( s ) \Gamma_{G(s)H(s)} ΓG(s)H(s) 包围 F(s) 平面点 (-1, j0) 的圈数
0型系统:当自变量s从s平面的原点出发,顺时针沿N氏围线 Γ \Gamma Γ 移动一周时,G(s) H(s)平面上曲线 Γ G ( s ) H ( s ) \Gamma_{G(s)H(s)} ΓG(s)H(s) 映射到G(jω)H(jω)平面上的曲线 Γ G ( j w ) H ( j w ) \Gamma_{G(jw)H(jw)} ΓG(jw)H(jw) ,即为当
ω:-∞→+∞的Nyquist曲线(全闭合曲线)
I型或I型以上系统:当自变量s从s平面的原点出发,顺时针沿N氏围线 Γ \Gamma Γ 移动一周时,G(s)H(s) 平面上曲线 Γ G H \Gamma_{GH} ΓGH 映射到 G(jω)H(jω) 平面上的曲线 Γ G ( j w ) H ( j w ) \Gamma_{G(jw)H(jw)} ΓG(jw)H(jw) ,即为当 w: 0+ →+∞、 -∞→ 0-、0- → 0+的Nyquist曲线(全闭合曲线)
具有等幅振荡环节的系统:当自变量 s 从 s 平面的原点出发,顺时针沿N氏围线 Γ \Gamma Γ 移动一周时,G(s)H(s) 平面上曲线 Γ G ( s ) H ( s ) \Gamma_{G(s)H(s)} ΓG(s)H(s) 映射到G(jω)H(jω)平面上的曲线 Γ G ( j w ) H ( j w ) \Gamma_{G(jw)H(jw)} ΓG(jw)H(jw),即为当 w: 0 → w n − w_{n-} wn− 、 w n − w_{n-} wn−→ w n + w_{n+} wn+ 、 w n + w_{n+} wn+→+∞、 -∞→ w n + w_{n+} wn+、
− w n + -w_{n+} −wn+→ − w n − -w_{n-} −wn−、 − w n − -w_{n-} −wn− → 0 的 Nyquist 曲线(全闭合曲线)
设 s 平面上不通过 F(s) 的任一零点和极点的N氏围线 Γ \Gamma Γ (围住整个右半s平面)包围 F(s) 平面的 P 个极点和 Z 个零点,则自变量 s 从原点顺时针沿 Γ \Gamma Γ 移动一周时, 在 G(jω)H(jω) 平面上, N 氏围线 Γ \Gamma Γ 的象或映射 Γ G ( j w ) H ( j w ) \Gamma_{G(jw)H(jw)} ΓG(jw)H(jw) 曲线(全闭合 Nyquist 曲线w: 0+ →+∞、 -∞→ 0-、0- → 0+) ,将以逆时针方向绕 (-1,j0) 点 R=P-Z 圈
实用的Nyquist稳定判据(全闭合Nyquist曲线)步骤:
已知: G k ( s ) G_k(s) Gk(s)
则v、P已知(P为 G k ( s ) G_k(s) Gk(s) 在右半s平面上开环极点的个数)
画出w: 0 →+∞、 -∞→ 0 开环(全闭合)Nyquist曲线
则可根据画好的全闭合Nyquist曲线,从图上获知Nyquist曲线逆时针绕(-1,j0)点的圈数,即可得R(当s在s 平面N氏围线 Γ \Gamma Γ 上顺时针方向移动一周时)
在s平面上N氏围线 Γ \Gamma Γ →G(jω)H(jω)平面上的 Γ G ( j w ) H ( j w ) \Gamma_{G(jw)H(jw)} ΓG(jw)H(jw) 曲线(全闭合Nyquist曲线),应用幅角原理,由R=P-Z(+逆,-顺)可得:Z=P-R
若Z =0 ,则系统稳定
若Z≠0,则系统不稳定
半闭合Nyquist曲线
设s平面上不通过F(s) 的任一零点和极点的N氏围线 Γ \Gamma Γ (围住整个右半s平面)包围F(s)平面的P个极点和Z个零点,则自变量s 从原点顺时针沿 Γ \Gamma Γ 移动一周时, 在G(jω)H(jω)平面上,N氏围线 Γ \Gamma Γ 的象或映射 Γ G ( j w ) H ( j w ) \Gamma_{G(jw)H(jw)} ΓG(jw)H(jw) 曲线(半闭合Nyquist曲线 w: 0→ + ∞ +\infty +∞),将以逆时针方向绕(-1,j0) 点 R=P-Z 圈( R = 2 ( N + − N − ) = 2 N R=2(N^+-N^-)=2N R=2(N+−N−)=2N)
扩展的Nyquist稳定判据(半闭合Nyquist曲线)步骤:(常用)
已知: G k ( s ) G_k(s) Gk(s);
则v、P已知(P为 G k ( s ) G_k(s) Gk(s) 在右半s平面上开环极点的个数)
画出w: 0 → + ∞ +\infty +∞开环(半闭合)Nyquist曲线;
**0型系统:**画w: 0 → + ∞ +\infty +∞ Nyquist曲线;
I型及以上系统:
具有 s = ± j w n s=±jw_n s=±jwn系统:
则可根据画好的半闭合Nyquist曲线获知,位于(-1,j0)点左侧,从上向下穿越(II → III)次数和N+(逆),以及从下向上穿越(III → II)次数和N-(顺),即可得 R = 2 ( N + − N − ) = 2 N R=2(N^+-N^-)=2N R=2(N+−N−)=2N(当s在s 平面N氏围线 Γ \Gamma Γ上顺时针方向移动一周时)
在s平面上N氏围线 Γ \Gamma Γ →G(jω)H(jω)平面上的 Γ G ( j w ) H ( j w ) \Gamma_{G(jw)H(jw)} ΓG(jw)H(jw) 曲线(半闭合Nyquist曲线),应用幅角原理,由R=P-Z,可得: Z=P-R
若Z =0 ,则系统稳定
若Z≠0,则系统不稳定
奈氏判据基于复平面的半闭合曲线 Γ G H \Gamma_{GH} ΓGH 判定系统的闭环稳定性,由于半闭合曲线 Γ G H \Gamma_{GH} ΓGH 可以转换为半对数坐标下的曲线,因此可以推广运用奈氏判据,其关键问题是需要根据半对数坐标下的 Γ G H \Gamma_{GH} ΓGH 曲线确定穿越次数 N 或 N+和N-
复平面 Γ G H \Gamma_{GH} ΓGH 曲线一般由两部分组成:开环幅相曲线和开环系统存在积分环节和等幅振荡环节时所补作的半径为无穷大的虚圆弧。而N 的确定取决于 A(w)>1 时 Γ G H \Gamma_{GH} ΓGH 穿越负实轴的次数,因此应建立和明确以下对应关系:
穿越点确定
{ A ( w c ) = ∣ G ( j w c ) H ( j w c ) ∣ = 1 L ( w c ) = 20 lg A ( w C ) = 0 \left\{\begin{matrix} A(w_c)=|G(jw_c)H(jw_c)|=1 \\ L(w_c)=20\lg A(w_C)=0 \end{matrix}\right. {A(wc)=∣G(jwc)H(jwc)∣=1L(wc)=20lgA(wC)=0
w c w_c wc 为截止频率
对于复平面的负实轴和开环对数相频特性,当取频率为穿越频率 w x w_x wx 时:
φ ( w x ) = ( 2 k + 1 ) π , k = 0 , ± 1 , … \varphi(w_x)=(2k+1)\pi,k=0,\pm1,\dots φ(wx)=(2k+1)π,k=0,±1,…
设半对数坐标下 Γ G H \Gamma_{GH} ΓGH 的对数幅频曲线和对数相频曲线分别为 Γ L \Gamma_{L} ΓL 和 Γ φ \Gamma_{\varphi} Γφ 。由于 Γ L \Gamma_{L} ΓL 等于L(w)曲线,则 Γ G H \Gamma_{GH} ΓGH 在A(w)>1时,穿越负实轴的点等于 Γ G H \Gamma_{GH} ΓGH 在半对数坐标下,对数幅频特性L(w)>0时对数相频特性曲线 Γ φ \Gamma_{\varphi} Γφ 与 ( 2 k + 1 ) π ; k = 0 , ± 1 , … (2k+1)\pi; k=0, \pm1, \dots (2k+1)π;k=0,±1,…, 平行线的交点
Γ φ \Gamma_{\varphi} Γφ 确定
开环系统无虚轴上极点时, Γ φ \Gamma_{\varphi} Γφ 等于 φ ( w ) \varphi(w) φ(w) 曲线
开环系统存在积分环节 1 s v , ( v > 0 ) \frac{1}{s^v},(v>0) sv1,(v>0) 时,复数平面的 Γ G H \Gamma_{GH} ΓGH 曲线,需从w=0+的开环幅相曲线的对应点G( j0+) H( j0+) 起,逆时针补作 v × 9 0 ∘ v\times90^{\circ} v×90∘ 半径为无穷大的虚圆弧。对应地,需从对数相频特性曲线 w 较小且 L(w)>0 的点处向上补作 v × 9 0 ∘ v\times90^{\circ} v×90∘ 的虚直线, φ ( w ) \varphi(w) φ(w) 曲线和补作的虚直线构成 Γ φ \Gamma_{\varphi} Γφ
开环系统存在等幅振荡环节: 1 ( s 2 + w n 2 ) v 1 , v 1 > 0 \frac{1}{{(s^2+w_n^2)}^{v_1}},v_1>0 (s2+wn2)v11,v1>0
时,复数平面的 Γ G H \Gamma_{GH} ΓGH 曲线,需从 w = w n − w = w_{n-} w=wn−的开环幅相曲线的对应点 G ( j w n − ) H ( j w n − ) G( jw_{n-})H(jw_{n-}) G(jwn−)H(jwn−) 起,逆时针补作 v 1 × 18 0 ∘ v_1\times180^{\circ} v1×180∘ 半径为无穷大的虚圆弧至 w = w n + w = w_{n+} w=wn+的对应点 G ( j w n + ) H ( j w n + ) G( jw_{n+})H(jw_{n+}) G(jwn+)H(jwn+) 处。对应地,需从对数相频特性曲线 φ ( w n − ) \varphi(w_{n-}) φ(wn−)点起向上补作 v 1 × 18 0 ∘ v_1\times180^{\circ} v1×180∘ 的虚直线至 φ ( w n + ) \varphi(w_{n+}) φ(wn+) 处, φ ( w ) \varphi(w) φ(w)曲线和补作的虚直线构成 Γ φ \Gamma_{\varphi} Γφ
穿越次数计算
正穿越一次: Γ G H \Gamma_{GH} ΓGH 由上向下穿越 (-1, j0) 点左侧的负实轴一次,等价于:在L(w)>0时, Γ φ \Gamma_{\varphi} Γφ 由下向上穿越 ( 2 k + 1 ) π (2k+1)\pi (2k+1)π 线一次。
负穿越一次: Γ G H \Gamma_{GH} ΓGH 由下向上穿越 (-1, j0) 点左侧的负实轴一次,等价于:在L(w)>0时, Γ φ \Gamma_{\varphi} Γφ 由上向下穿越 ( 2 k + 1 ) π (2k+1)\pi (2k+1)π 线一次
正穿越半次: Γ G H \Gamma_{GH} ΓGH 由上向下止于或由上向下起于 (-1, j0) 点左侧的负实轴,等价于:在L(w)>0时, Γ φ \Gamma_{\varphi} Γφ 由下向上止于或由下向上起于 ( 2 k + 1 ) π (2k+1)\pi (2k+1)π 线一次
负穿越半次: Γ G H \Gamma_{GH} ΓGH 由下向上止于或由下向上起于(-1, j0)点左侧的负实轴,等价于:在 L(w)>0 时, Γ φ \Gamma_{\varphi} Γφ 由上向下止于或由上向下起于 ( 2 k + 1 ) π (2k+1)\pi (2k+1)π 线一次
应该指出的是,补作的虚直线所产生的穿越皆为负穿越
Bode图(半对数坐标 Γ G H \Gamma_{GH} ΓGH曲线)上的Nyquist稳定判据步骤:
已知: G k ( s ) G_k(s) Gk(s)
则v、P已知(P为 G k ( s ) G_k(s) Gk(s) 在右半s平面上开环极点的个数)
画出 w : 0 → + ∞ w: 0 →+\infty w:0→+∞开环半对数坐标 Γ G H \Gamma_{GH} ΓGH 曲线(Bode图)
从画好的对数幅频特性L(w)>0时,对数相频特性曲线 Γ φ \Gamma{\varphi} Γφ与 ( 2 k + 1 ) π ; k = 0 , ± 1 , … (2k+1)\pi; k=0, \pm1, \dots (2k+1)π;k=0,±1,…, 平行线的交点,可得R=2(N±N-)=2N
在s平面上N氏围线 Γ \Gamma Γ →G(jω)H(jω)平面上的 Γ G ( j w ) H ( j w ) \Gamma_{G(jw)H(jw)} ΓG(jw)H(jw) 曲线(半闭合Nyquist曲线),应用幅角原理,由R=P-Z,可得: Z=P-R
若Z =0 ,则系统稳定;若Z≠0,则系统不稳定
当开环传递函数的某些系数发生变化时, Γ G H \Gamma_{GH} ΓGH 包围 (-1, j0) 点的情况亦随之改变
在稳定性研究中,称 (-1, j0)点为临界点,而闭合曲线 Γ G H \Gamma_{GH} ΓGH 相对于临界点的位置即偏离临界点的程度,反映系统的相对稳定性
频域的相对稳定性即稳定裕度常用相角裕度 γ \gamma γ 和幅值裕度 h h h 来度量
设 w c w_c wc 为系统的截止频率(固定模为1,比较相角离 (-1, j0)点的远近):
A ( w c ) = ∣ G ( j w c ) H ( j w c ) ∣ = 1 A(w_c)=|G(jw_c)H(jw_c)|=1 A(wc)=∣G(jwc)H(jwc)∣=1
定义相角裕度为:
γ = 18 0 ∘ + ∠ G ( j w c ) H ( j w c ) \gamma=180^{\circ}+\angle{G(jw_c)H(jw_c)} γ=180∘+∠G(jwc)H(jwc)
相角裕度 γ \gamma γ 的含义是,对于闭环稳定系统,如果系统开环相频特性再滞后 γ \gamma γ 度,则系统将处于临界稳定状态
设 w x w_x wx 为系统的穿越频率(固定相角为-π,比较模离 (-1, j0)点的远近) :
φ ( w x ) = ∠ G ( j w x ) H ( j w x ) = ( 2 k + 1 ) π ; k = 0 , ± 1 , … \varphi(w_x)=\angle{G(jw_x)H(jw_x)}=(2k+1)\pi;k=0,\pm1,\dots φ(wx)=∠G(jwx)H(jwx)=(2k+1)π;k=0,±1,…
定义幅值裕度为:
h = 1 ∣ G ( j w x ) H ( j w x ) ∣ h=\frac{1}{|G(jw_x)H(jw_x)|} h=∣G(jwx)H(jwx)∣1
幅值裕度 h h h 的含义是,对于闭环稳定系统,如果系统开环幅频特性再增大 h 倍,则系统将处于临界稳定状态
对数坐标下,幅值裕度按下式定义:
h = − 20 lg ∣ G ( j w x ) H ( j w x ) ∣ ( d B ) h=-20\lg{|G(jw_x)H(jw_x)|}(dB) h=−20lg∣G(jwx)H(jwx)∣(dB)
仅用相角裕度或幅值裕度,都不足以反映系统的稳定程度。 对非最小相位系统并不都适用
开环频域性能指标
稳:相角裕度、幅值裕度
快:幅值穿越频率(截止频率、剪切频率)、相角穿越频率
准:稳态误差
在开环频率特性曲线上Nyquist曲线、Bode图
闭环频域性能指标
稳:谐振峰值
快:带宽频率、谐振频率
准:零频幅值
在闭环频率特性曲线上Nyquist曲线、Bode图、闭环幅频曲线、等M/N圆、尼科尔斯图
零频幅值 M 0 M_0 M0
w = 0 M ( 0 ) = M 0 = C ( 0 ) R ( 0 ) M ( 0 ) = { k 1 + k , v = 0 1 , v > 1 w=0\\ M(0)=M_0=\frac{C(0)}{R(0)}\\ M(0)=\begin{cases} \frac{k}{1+k},v=0 \\ 1,v>1 \end{cases} w=0M(0)=M0=R(0)C(0)M(0)={1+kk,v=01,v>1
M ( 0 ) = 1 M(0)=1 M(0)=1 表示系统阶跃响应的终值等于输入,静差为零
M ( 0 ) ≠ 1 M(0)\ne1 M(0)=1表示系统有静差
M 0 M_0 M0 与1相差的大小,反映了系统的稳态精度, M 0 M_0 M0 越接近1,系统精度越高
谐振频率 w r w_r wr
闭环幅频特性出现峰值时的频率称为谐振频率 w r w_r wr,它在一定程度上反映了系统的快速性
w r w_r wr 越大,系统瞬态响应越快
谐振峰值 M r M_r Mr
闭环幅频特性最大值 M r = M ( w r ) M_r=M(w_r) Mr=M(wr) 称为谐振峰值 M r M_r Mr
峰值大,表明系统对某个频率的正弦输入信号反映强烈,有共振的倾向,即意味着系统的平稳性较差,阶跃响应将有过大的超调量
一般要求: M r < 1.5 M 0 M_r<1.5M_0 Mr<1.5M0
带宽频率 w b w_b wb 与系统带宽 0 ~ w b w_b wb
w b w_b wb :幅频特性 M ( w ) M(w) M(w) 的数值从 M ( 0 ) M(0) M(0) 衰减到0.707 M ( 0 ) M(0) M(0) 时所对应的频率
w b w_b wb 高: M ( w ) M(w) M(w) 从 M ( 0 ) M(0) M(0) 到 0.707 M ( 0 ) M(0) M(0) 所占据的频率区间0~ w b w_b wb 较宽,表
明系统复现快速变化的信号能力强,失真小,也即快速性好,阶跃响应的调节时间短
w b w_b wb 低: M ( w ) M(w) M(w) 从 M ( 0 ) M(0) M(0) 到 0.707 M ( 0 ) M(0) M(0) 所占据的频率区间0~ w b w_b wb 较窄,带 宽窄,系统反映迟钝,失真大,快速性差(复现快速信号能力差)
系统带宽与信号频谱的关系
( 0 , w b ) (0,w_b) (0,wb) 称为系统的带宽
带宽定义表明,对于高于带宽频率的正弦输入信号,系统输出将呈现较大的衰减
当控制输入为矩形波时,系统的跟踪能力取决于带宽覆盖矩形波频谱的范围,带宽大则处于较高频率范围的谱线衰减小,故失真小
当控制输入信号频谱特性具有收敛形式时,可按跟踪要求选定常数 ε \varepsilon ε ,若 w > w b w>w_b w>wb 时, ∣ F ( j w ) ∣ < ε |F(jw)|<\varepsilon ∣F(jw)∣<ε,从而确定系统所需要的带宽 w b w_b wb
系统分析应区分输入信号的性质、位置,根据其频谱或谱密度、相应的传递函数及噪声频段选择带宽
将控制系统中的各个变量看成一些信号而这些信号又是由许多不同频率的正弦信号合成的,各个变量的运动就是系统对各个不同频率的信号的响应的综合
建模
建立系统(1阶、2阶、任意阶次)数学模型:微分方程/传递函数/频率特性
画系统频率特性曲线
系统开环频率特性曲线(Nyquist曲线/波德图)
系统闭环频率特性(Nyquist曲线/波德图/幅频曲线)
系统分析(在开/闭环频率特性曲线上)
稳定性分析
性能分析(求系统的性能指标)
动态响应性能分析(稳、快)
稳态响应性能分析
所谓校正,就是在系统中加入一些其参数可以根据需要而改变的机构或装置,使系统整个特性发生变化,从而满足给定的各项性能指标
在控制系统设计中,采用的设计方法一般依据性能指标的形式而定:
目前,工程技术界多习惯采用频率法,故通常通过近似公式进行两种指标的互换
时域性能指标 | 频域开环性能指标 | 频域闭环性能指标 | |
---|---|---|---|
稳 | σ % \sigma\% σ% | γ \gamma γ | M r M_r Mr |
快 | t s , t p , t r , t d t_s,t_p,t_r,t_d ts,tp,tr,td | w c w_c wc | w r , w b w_r,w_b wr,wb |
准 | K , v ( e s s ) K,v(e_{ss}) K,v(ess) | K , v K,v K,v | M ( 0 ) , v M(0),v M(0),v |
无论采用哪种校正方式,都要求校正后的系统既能以所需精度跟踪输入信号,又能抑制噪声扰动信号
在控制系统实际运行中,输入信号一般是低频信号,而噪声信号则一般是高频信号
为了使系统具有较高的稳定裕度,希望系统开环对数幅频特性在截止频率 w c w_c wc 处的斜率为 -20dB/dec,但从要求系统具有较强的从噪声中辩识信号的能力来考虑, 却又希望 w c w_c wc 处的斜率小于 40dB/dec
通常, 一个设计良好的实际运行系统, 其相角裕度具有45º左右的数值:
要实现 45º 左右的相角裕度要求,开环对数幅频特性在中频区的斜率应为-20dB/dec,同时要求中频区占据一定的频率范围,以保证在系统参数变化时,相角裕度变化不大
过此中频区后,要求系统幅频特性迅速衰减,以削弱噪声对系统的影响。这是选择系统带宽应该考虑的一个方面。另一方面,进入系统输入端的信号,既有输入信号r(t),又有噪声信号 n(t),如果输入信号的带宽为 0 ∼ w M 0\sim{w_M} 0∼wM,噪声信号集中起作用的频带为 w 1 ∼ w n w_1\sim{w_n} w1∼wn,则控制系统的带宽频率通常取为:
w b = ( 5 ∼ 10 ) w M w_b=(5\sim10)w_M wb=(5∼10)wM
且使 w 1 ∼ w n w_1\sim{w_n} w1∼wn 处于 0 ∼ w b 0\sim{w_b} 0∼wb 范围之外
根据稳、快、准的要求,选择时注意:
按照校正装置在系统中的连接方式,控制系统校正方式可分为四种:
串联校正装置一般接在系统误差测量点之后和放大器之前,串接于系统前向通道之中
反馈校正装置一般接在系统局部反馈通路之中
前馈校正又称顺馈校正,是在系统主反馈回路之外采用的校正方式
前馈校正装置接在系统给定值(或指令、参考输入信号)之后及主反馈作用点之前的前向通道上。这种校正方式的作用相当于对给定值信号进行整形或滤波后,再送入反馈系统
另一种前馈校正装置接在系统可测扰动作用点与误差测量点之间,对扰动信号进行直接或间接测量,并经变换后接入系统,形成一条附加的对扰动影响进行补偿的通道
前馈校正可以单独作用于开环控制系统,也可以作为反馈控制系统的附加校正而组成复合控制系统
复合校正方式是在反馈控制回路中,加入前馈校正通路,组成一个有机整体;
复合校正装置与前馈校正一致,有 2 种接入方式:
在控制系统设计中,常用的校正方式为串联校正和反馈校正两种
一般来说,串联校正设计比反馈校正设计简单,也比较容易对信号进行各种必要形式的变换
串联校正装置又分无源和有源两类
在实际控制系统中,还广泛采用反馈校正装置。一般来说,反馈校正所需元件数目比串联校正多。由于反馈信号通常由系统输出端或放大器输出级供给,信号是从高功率点传向低功率点,因此反馈校正一般无需附加放大器
此外,反馈校正尚可消除系统原有部分参数波动对系统性能的影响。在性能指标要求较高的控制系统设计中,常常兼用串联校正与反馈校正两种方式
确定校正装置的具体形式时,应先了解校正装置所需提供的控制规律,以便选择相应的元件
包含校正装置在内的控制器,常常采用比例P、微分D、积分I等基本控制规律,或者采用这些基本控制规律的某些组合,如比例微分PD、比例积分PI、比例-积分-微分PID等组合控制规律,以实现对被控对象的有效控制
串联滞后-超前网络校正装置的特性(作用)
z 变换法是离散系统理论的数学基础。利用 z 变换法研究离散系统,可以把连续系统中的许多概念和方法,推广应用于线性离散系统
线性离散系统分析方法的基本思路
建模
建立线性离散系统数学模型:(前向或后向)差分方程、闭环脉冲传递函数
求时域响应
已知条件:
系统数学模型:脉冲传递函数(或:差分方程→脉冲传递函数)
输入(典型输入):单位阶跃信号(单位脉冲,单位斜坡,单位加速度)
求: c ∗ ( t ) / c ( n T ) / c ( k ) c^*(t)/c(nT)/c(k) c∗(t)/c(nT)/c(k)
系统分析
稳定性分析
性能分析(求系统的性能指标)
通常,把系统中的离散信号是脉冲序列形式的离散系统,称为采样控制系统或脉冲控制系统
一般说来,采样系统是对来自传感器的连续信息在某些规定的时间瞬时上取值。如果在有规律的间隔上,系统取到了离散信息,则这种采样称为周期采样
反之,如果信息之间的间隔是时变的,或随机的,则称为非周期采样,或随机采样
为了使两种信号在系统中能相互传递,在连续信号和脉冲序列之间要用采样器,而在脉冲序列和连续信号之间用保持器,以实现两种信号的转换
采样器和保持器,是采样控制系统中的两个特殊环节
在采样控制系统中,把连续信号转变为脉冲序列的过程称为采样过程,简称采样。实现采样的装置称为采样器或称采样开关。
用 T 表示采样周期,单位为s; f s = 1 / T f_s=1/T fs=1/T 表示采样频率,单位为1/s; w s = 2 π f s = 2 π / T w_s=2\pi{f_s}=2\pi/T ws=2πfs=2π/T 表示采样角频率,单位为rad/s
在实际应用中,采样开关多为电子开关,闭合时间极短,采样持续时间 τ \tau τ 远小于采样周期 T, 也远小于系统连续部分的最大时间常数。为了简化系统的分析,可认为 τ \tau τ 趋于零,即把采样器的输出近似看成一串强度等于矩形脉冲面积的理想脉冲 e ∗ ( t ) e^*(t) e∗(t)
在采样控制系统中,把脉冲序列转变为连续信号的过程称为信号复现过程。实现复现过程的装置称为保持器。
把数字序列形式的离散系统,称为数字控制系统或计算机控制系统
数字控制系统,是一种以数字计算机为控制器去控制具有连续工作状态的被控对象的闭环控制系统。因此,数字控制系统包括工作于离散状态下的数字计算机和工作于连续状态下的被控对象两大部分
计算机作为系统的控制器,其输入和输出只能是二进制编码的数字信号,即在时间上和幅值上都离散的信号,而系统中被控对象和测量元件的输入和输出是连续信号,所以在计算机控制系统中,需要应用A/D (模/数) 和 D/A (数/模) 转换器,以实现两种信号的转换
数字计算机在对系统进行实时控制时,每隔 T 秒进行一次控制修正,T 为采样周期
在每个采样周期中,控制器要完成对于连续信号的采样编码 (即A/D 过程) 和按控制律进行的数码运算,然后将计算结果由输出寄存器经解码网络将数码转换成连续信号(即D/A过程)
因此,A/D 转换器和 D/A 转换器是计算机控制系统中的两个特殊环节
线性定理
Z [ e 1 ( t ) ± e 2 ( t ) ] = E 1 ( z ) ± E 2 ( z ) Z [ a e ( t ) ] = a E ( z ) \mathscr{Z}[e_1(t){\pm}e_2(t)]=E_1(z){\pm}E_2(z)\\ \mathscr{Z}[ae(t)]=aE(z) Z[e1(t)±e2(t)]=E1(z)±E2(z)Z[ae(t)]=aE(z)
实数位移定理
实数位移定理又称平移定理。实数位移的含意,是指整个采样序列在时间轴上左右平移若干采样周期,其中向左平移为超前,向右平移为滞后
Z [ e ( t − k T ) ] = z − k E ( z ) Z [ e ( t + k T ) ] = z k [ E ( z ) − ∑ n = 0 k − 1 e ( n T ) z − n ] \mathscr{Z}[e(t-kT)]=z^{-k}E(z)\\ \mathscr{Z}[e(t+kT)]=z^k[E(z)-\sum^{k-1}_{n=0}e(nT)z^{-n}] Z[e(t−kT)]=z−kE(z)Z[e(t+kT)]=zk[E(z)−n=0∑k−1e(nT)z−n]
复数位移定理
e ∓ a t e ( t ) = E ( z e ∓ a T ) \mathscr{e^{{\mp}at}e(t)}=E(ze^{{\mp}aT}) e∓ate(t)=E(ze∓aT)
终值定理
如果函数e(t)的 z 变换为 E(z),函数序列 e(nT) 为有限值 (n=0,1,2,…),且极限存在,则函数序列的终值:
lim n → ∞ e ( n T ) = lim z → 1 ( z − 1 ) E ( z ) 也可表示为 e ( ∞ ) = lim n → ∞ e ( n T ) = lim z → 1 ( 1 − z − 1 ) E ( z ) \lim_{n\to\infty}e(nT)=\lim_{z\to1}(z-1)E(z)\\ \text{也可表示为 }e(\infty)=\lim_{n\to\infty}e(nT)=\lim_{z\to1}(1-z^{-1})E(z) n→∞lime(nT)=z→1lim(z−1)E(z)也可表示为 e(∞)=n→∞lime(nT)=z→1lim(1−z−1)E(z)
卷积定理
g ( n T ) = x ( n T ) ∗ y ( n T ) G ( z ) = X ( z ) . Y ( z ) g(nT)=x(nT)*y(nT)\\ G(z)=X(z).Y(z) g(nT)=x(nT)∗y(nT)G(z)=X(z).Y(z)
K p = lim z → 1 [ 1 + G ( z ) ] K_p=\lim_{z\to1}[1+G(z)] Kp=z→1lim[1+G(z)]
称为静态位置误差系数
K v = lim z → 1 ( z − 1 ) G ( z ) K_v=\lim_{z\to1}(z-1)G(z) Kv=z→1lim(z−1)G(z)
称为静态速度误差系数。0 型系统的 K v = 0 K_v=0 Kv=0 ,I 型系统的 K v K_v Kv 为有限值,Ⅱ型和 II 型以上系统的 K v = ∞ K_v=\infty Kv=∞
K a = lim z → 1 ( z − 1 ) 2 G ( z ) K_a=\lim_{z\to1}(z-1)^2G(z) Ka=z→1lim(z−1)2G(z)
称为静态加速度误差系数。0 型及 I 型系统的 K a = 0 K_a=0 Ka=0 ,Ⅱ型系统的 K a K_a Ka 为常值,III 型及 III 型以上系统的 K a = ∞ K_a=\infty Ka=∞