该系列博客主要讲述Matlab软件在自动控制方面的应用,如无自动控制理论基础,请先学习自动控制系列博文,该系列博客不再详细讲解自动控制理论知识。
自动控制理论基础相关链接:https://blog.csdn.net/qq_39032096/category_10287468.html?spm=1001.2014.3001.5482
博客参考书籍:《MATLAB/Simulink与控制系统仿真》。
2.离散系统数学模型与稳定性分析
2.1 离散系统数学模型
2.1.1 线性常系数差分方程
对于一般的线性定常离散系统, k k k时刻的输出 c ( k ) c(k) c(k)不但与 k k k时刻的输入 r ( k ) r(k) r(k)有关,且与 k k k时刻以前的输入 r ( k − 1 ) , r ( k − 2 ) , … , r(k-1),r(k-2),\dots, r(k−1),r(k−2),…,有关,同时还与 k k k以前的输出 c ( k − 1 ) , c ( k − 2 ) , … , c(k-1),c(k-2),\dots, c(k−1),c(k−2),…,有关;这种关系可用 n n n阶后向差分方程描述:
c ( k ) + a 1 c ( k − 1 ) + a 2 c ( k − 2 ) + ⋯ + a n − 1 c ( k − n + 1 ) + a n c ( k − n ) = b 0 r ( k ) + b 1 r ( k − 1 ) + ⋯ + b m − 1 r ( k − m + 1 ) + b m r ( k − m ) \begin{aligned} &c(k)+a_1c(k-1)+a_2c(k-2)+\dots+a_{n-1}c(k-n+1)+a_nc(k-n)\\\\ =&b_0r(k)+b_1r(k-1)+\dots+b_{m-1}r(k-m+1)+b_mr(k-m) \end{aligned} =c(k)+a1c(k−1)+a2c(k−2)+⋯+an−1c(k−n+1)+anc(k−n)b0r(k)+b1r(k−1)+⋯+bm−1r(k−m+1)+bmr(k−m)
亦可表示为:
c ( k ) = − ∑ i = 1 n a i c ( k − i ) + ∑ j = 0 m b j r ( k − j ) c(k)=-\sum_{i=1}^na_ic(k-i)+\sum_{j=0}^mb_jr(k-j) c(k)=−i=1∑naic(k−i)+j=0∑mbjr(k−j)
其中: a i ( i = 1 , 2 , … , n ) 、 b j ( j = 0 , 1 , 2 , … , m ) a_i(i=1,2,\dots,n)、b_j(j=0,1,2,\dots,m) ai(i=1,2,…,n)、bj(j=0,1,2,…,m)为常系数, m ≤ n m≤n m≤n;
上式称为 n n n阶线性常系数差分方程;
线性定常离散系统亦可采用如下 n n n阶前向差分方程描述:
c ( k + n ) + a 1 c ( k + n − 1 ) + a 2 c ( k + n − 2 ) + ⋯ + a n − 1 c ( k + 1 ) + a n c ( k ) = b 0 r ( k + m ) + b 1 r ( k + m − 1 ) + b 2 r ( k + m − 2 ) + ⋯ + b m − 1 r ( k + 1 ) + b m r ( k ) \begin{aligned} &c(k+n)+a_1c(k+n-1)+a_2c(k+n-2)+\dots+a_{n-1}c(k+1)+a_nc(k)\\\\ =&b_0r(k+m)+b_1r(k+m-1)+b_2r(k+m-2)+\dots+b_{m-1}r(k+1)+b_mr(k) \end{aligned} =c(k+n)+a1c(k+n−1)+a2c(k+n−2)+⋯+an−1c(k+1)+anc(k)b0r(k+m)+b1r(k+m−1)+b2r(k+m−2)+⋯+bm−1r(k+1)+bmr(k)
亦可表示为:
c ( k + n ) = − ∑ i = 1 n a i c ( k + n − i ) + ∑ j = 0 m b j r ( k + m − j ) c(k+n)=-\sum_{i=1}^na_ic(k+n-i)+\sum_{j=0}^mb_jr(k+m-j) c(k+n)=−i=1∑naic(k+n−i)+j=0∑mbjr(k+m−j)
2.1.2 脉冲传递函数
-
脉冲传递函数定义
设开环离散系统如下图所示:
如果系统初始条件为零,输入信号为 r ( t ) r(t) r(t),采样后 r ∗ ( t ) r^*(t) r∗(t)的 z z z变换函数为 R ( z ) R(z) R(z),系统连续部分的输出为 c ( t ) c(t) c(t),采样后 c ∗ ( t ) c^*(t) c∗(t)的 z z z变换函数为 C ( z ) C(z) C(z),则线性定常离散系统的脉冲传递函数定义为:系统输出采样信号的 z z z变换与输入采样信号的 z z z变换之比,记为:
G ( z ) = C ( z ) R ( z ) = ∑ n = 0 ∞ c ( n T ) z − n ∑ n = 0 ∞ r ( n T ) z − n G(z)=\frac{C(z)}{R(z)}=\frac{\displaystyle\sum_{n=0}^{\infty}c(nT)z^{-n}}{\displaystyle\sum_{n=0}^{\infty}r(nT)z^{-n}} G(z)=R(z)C(z)=n=0∑∞r(nT)z−nn=0∑∞c(nT)z−n如果已知 R ( z ) 、 G ( z ) R(z)、G(z) R(z)、G(z),则在零初始条件下,线性定常离散系统的输出采样信号为:
c ∗ ( t ) = Z − 1 [ C ( z ) ] = Z − 1 [ G ( z ) R ( z ) ] c^*(t)=Z^{-1}[C(z)]=Z^{-1}[G(z)R(z)] c∗(t)=Z−1[C(z)]=Z−1[G(z)R(z)]
-
脉冲传递函数的求法
Example1: 设开环系统如下图所示,求相应的脉冲传递函数 G ( z ) G(z) G(z);
解:
将 G ( s ) G(s) G(s)展成部分分式:
G ( s ) = 1 s − 1 s + a G(s)=\frac{1}{s}-\frac{1}{s+a} G(s)=s1−s+a1
可得:
G ( z ) = z z − 1 − z z − e − a T = z ( 1 − e − a T ) ( z − 1 ) ( z − e − a T ) G(z)=\frac{z}{z-1}-\frac{z}{z-{\rm e}^{-aT}}=\frac{z(1-{\rm e}^{-aT})}{(z-1)(z-{\rm e}^{-aT})} G(z)=z−1z−z−e−aTz=(z−1)(z−e−aT)z(1−e−aT)
2.1.3 开环系统脉冲传递函数
-
采样拉氏变换两个重要性质
-
采样函数的拉氏变换具有周期性,即:
G ∗ ( s ) = G ∗ ( s + j k ω s ) G^*(s)=G^*(s+{\rm j}k\omega_s) G∗(s)=G∗(s+jkωs)
其中: ω s \omega_s ωs为采样角频率; -
若采样函数的拉氏变换 E ∗ ( s ) E^*(s) E∗(s)与连续函数的拉氏变换 G ( s ) G(s) G(s)相乘后再离散化,则 E ∗ ( s ) E^*(s) E∗(s)可以从离散符号中提出来,即:
[ G ( s ) E ∗ ( s ) ] ∗ = G ∗ ( s ) E ∗ ( s ) [G(s)E^*(s)]^*=G^*(s)E^*(s) [G(s)E∗(s)]∗=G∗(s)E∗(s)
-
-
有串联环节的开环系统脉冲传递函数
-
串联环节间有采样开关
设开环离散系统如下图所示,在两个串联连接环节 G 1 ( s ) 、 G 2 ( s ) G_1(s)、G_2(s) G1(s)、G2(s)间,有理想采样开关隔开;
串联环节间有采样开关开环系统脉冲传递函数为:
G ( z ) = C ( z ) R ( z ) = G 1 ( z ) G 2 ( z ) G(z)=\frac{C(z)}{R(z)}=G_1(z)G_2(z) G(z)=R(z)C(z)=G1(z)G2(z)
有理想采样开关隔开的两个线性连续环节串联时的脉冲传递函数,等于这两个环节各自的脉冲传递函数之积,可推广到 n n n个环节的串联; -
串联环节间无采样开关
设开环离散系统如下图所示,在两个串联连接环节 G 1 ( s ) 、 G 2 ( s ) G_1(s)、G_2(s) G1(s)、G2(s)间,无理想采样开关隔开;
串联环节间无采样开关开环系统脉冲传递函数为:
G ( z ) = C ( z ) R ( z ) = G 1 G 2 ( z ) G(z)=\frac{C(z)}{R(z)}=G_1G_2(z) G(z)=R(z)C(z)=G1G2(z)
没有理想采样开关隔开的两个线性连续环节串联时的脉冲传递函数,等于这两个环节传递函数乘积后的相应 z z z变换,可推广到 n n n个环节的串联;Example2: 设开环离散系统如下图所示,输入信号 r ( t ) = 1 ( t ) r(t)=1(t) r(t)=1(t),求系统 ( a ) {\rm (a)} (a)和 ( b ) {\rm (b)} (b)的脉冲传递函数 G ( z ) G(z) G(z)和输出的 z z z变换 C ( z ) C(z) C(z);
解:
输入 r ( t ) = 1 ( t ) r(t)=1(t) r(t)=1(t)的 z z z变换为:
R ( z ) = z z − 1 R(z)=\frac{z}{z-1} R(z)=z−1z
对于系统 ( a ) {\rm (a)} (a):
G 1 ( z ) = Z [ 1 s ] = z z − 1 , G 2 ( z ) = Z [ a s + a ] = a z z − e − a T G_1(z)=Z\left[\frac{1}{s}\right]=\frac{z}{z-1},G_2(z)=Z\left[\frac{a}{s+a}\right]=\frac{az}{z-{\rm e}^{-aT}} G1(z)=Z[s1]=z−1z,G2(z)=Z[s+aa]=z−e−aTaz
有:
G ( z ) = G 1 ( z ) G 2 ( z ) = a z 2 ( z − 1 ) ( z − e − a T ) G(z)=G_1(z)G_2(z)=\frac{az^2}{(z-1)(z-{\rm e}^{-aT})} G(z)=G1(z)G2(z)=(z−1)(z−e−aT)az2C ( z ) = G ( z ) R ( z ) = a z 3 ( z − 1 ) 2 ( z − e − a T ) C(z)=G(z)R(z)=\frac{az^3}{(z-1)^2(z-{\rm e}^{-aT})} C(z)=G(z)R(z)=(z−1)2(z−e−aT)az3
对于系统 ( b ) {\rm (b)} (b):
G 1 ( s ) G 2 ( s ) = a s ( s + a ) G_1(s)G_2(s)=\frac{a}{s(s+a)} G1(s)G2(s)=s(s+a)a
有:
G ( z ) = G 1 G 2 ( z ) = Z [ a s ( s + a ) ] = z ( 1 − e − a T ) ( z − 1 ) ( z − e − a T ) G(z)=G_1G_2(z)=Z\left[\frac{a}{s(s+a)}\right]=\frac{z(1-{\rm e}^{-aT})}{(z-1)(z-{\rm e}^{-aT})} G(z)=G1G2(z)=Z[s(s+a)a]=(z−1)(z−e−aT)z(1−e−aT)C ( z ) = G ( z ) R ( z ) = z 2 ( 1 − e − a T ) ( z − 1 ) ( z − e − a T ) C(z)=G(z)R(z)=\frac{z^2(1-{\rm e}^{-aT})}{(z-1)(z-{\rm e}^{-aT})} C(z)=G(z)R(z)=(z−1)(z−e−aT)z2(1−e−aT)
-
有零阶保持器时的开环系统脉冲传递函数
设有零阶保持器的开环离散系统如下图所示:
图中: G h ( s ) G_h(s) Gh(s)为零阶保持器传递函数, G 0 ( s ) G_0(s) G0(s)为连续部分传递函数,两个串联环节间无同步采样开关隔离;
有零阶保持器时,开环系统脉冲传递函数:
G ( z ) = C ( z ) R ( z ) = ( 1 − z − 1 ) Z [ G 0 ( s ) s ] G(z)=\frac{C(z)}{R(z)}=(1-z^{-1})Z\left[\frac{G_0(s)}{s}\right] G(z)=R(z)C(z)=(1−z−1)Z[sG0(s)]
-
2.1.4 闭环系统脉冲传递函数
一种常见的误差采样闭环离散系统结构图分析:
由上图可知,连续输出信号和误差信号拉氏变换为:
C ( s ) = G ( s ) E ∗ ( s ) , E ( s ) = R ( s ) − H ( s ) C ( s ) C(s)=G(s)E^*(s),E(s)=R(s)-H(s)C(s) C(s)=G(s)E∗(s),E(s)=R(s)−H(s)C(s)
因此有:
E ( s ) = R ( s ) − H ( s ) G ( s ) E ∗ ( s ) E(s)=R(s)-H(s)G(s)E^*(s) E(s)=R(s)−H(s)G(s)E∗(s)
误差采样信号 e ∗ ( t ) e^*(t) e∗(t)的拉氏变换为:
E ∗ ( s ) = R ∗ ( s ) − H G ∗ ( s ) E ∗ ( s ) E^*(s)=R^*(s)-HG^*(s)E^*(s) E∗(s)=R∗(s)−HG∗(s)E∗(s)
整理可得:
E ∗ ( s ) = R ∗ ( s ) 1 + H G ∗ ( s ) E^*(s)=\frac{R^*(s)}{1+HG^*(s)} E∗(s)=1+HG∗(s)R∗(s)
由于:
C ∗ ( s ) = [ G ∗ ( s ) E ∗ ( s ) ] ∗ = G ∗ ( s ) E ∗ ( s ) = G ∗ ( s ) 1 + H G ∗ ( s ) R ∗ ( s ) C^*(s)=\left[G^*(s)E^*(s)\right]^*=G^*(s)E^*(s)=\frac{G^*(s)}{1+HG^*(s)}R^*(s) C∗(s)=[G∗(s)E∗(s)]∗=G∗(s)E∗(s)=1+HG∗(s)G∗(s)R∗(s)
因此有:
E ( z ) = 1 1 + H G ( z ) R ( z ) , C ( z ) = G ( z ) 1 + H G ( z ) R ( z ) E(z)=\frac{1}{1+HG(z)}R(z),C(z)=\frac{G(z)}{1+HG(z)}R(z) E(z)=1+HG(z)1R(z),C(z)=1+HG(z)G(z)R(z)
定义闭环离散系统对于输入量的误差脉冲传递函数:
Φ e ( z ) = E ( z ) R ( z ) = 1 1 + H G ( z ) \Phi_e(z)=\frac{E(z)}{R(z)}=\frac{1}{1+HG(z)} Φe(z)=R(z)E(z)=1+HG(z)1
定义闭环离散系统对于输入量的脉冲传递函数:
Φ ( z ) = C ( z ) R ( z ) = G ( z ) 1 + H G ( z ) \Phi(z)=\frac{C(z)}{R(z)}=\frac{G(z)}{1+HG(z)} Φ(z)=R(z)C(z)=1+HG(z)G(z)
闭环离散系统特征方程:
D ( z ) = 1 + H G ( z ) = 0 D(z)=1+HG(z)=0 D(z)=1+HG(z)=0
注:只要误差信号 e ( t ) e(t) e(t)处没有采样开关,输入采样信号 r ∗ ( t ) r^*(t) r∗(t)便不存在,此时不能求出闭环离散系统对于输入量的脉冲传递函数,只能求出输出采样信号的 z z z变换函数 C ( z ) C(z) C(z);
2.1.5 典型闭环离散系统及输出 z z z变换函数
-
前四种典型离散系统如下图:
上图 C ( z ) C(z) C(z)求解依次为:
C ( z ) = G ( z ) R ( z ) 1 + G H ( z ) C ( z ) = R G 1 ( z ) G 2 ( z ) 1 + G 2 H G 1 ( z ) C ( z ) = G ( z ) R ( z ) 1 + G ( z ) H ( z ) C ( z ) = G 1 ( z ) G 2 ( z ) R ( z ) 1 + G 1 ( z ) G 2 H ( z ) \begin{aligned} &C(z)=\frac{G(z)R(z)}{1+GH(z)}\\\\ &C(z)=\frac{RG_1(z)G_2(z)}{1+G_2HG_1(z)}\\\\ &C(z)=\frac{G(z)R(z)}{1+G(z)H(z)}\\\\ &C(z)=\frac{G_1(z)G_2(z)R(z)}{1+G_1(z)G_2H(z)} \end{aligned} C(z)=1+GH(z)G(z)R(z)C(z)=1+G2HG1(z)RG1(z)G2(z)C(z)=1+G(z)H(z)G(z)R(z)C(z)=1+G1(z)G2H(z)G1(z)G2(z)R(z) -
后四种典型离散系统如下图:
上图 C ( z ) C(z) C(z)求解依次为:
C ( z ) = R G 1 ( z ) G 2 ( z ) G 3 ( z ) 1 + G 2 ( z ) G 1 G 3 H ( z ) C ( z ) = R G ( z ) 1 + H G ( z ) C ( z ) = G 1 ( z ) G 2 ( z ) R ( z ) 1 + G 1 ( z ) G 2 ( z ) H ( z ) C ( z ) = R ( z ) G ( z ) 1 + G ( z ) H ( z ) \begin{aligned} &C(z)=\frac{RG_1(z)G_2(z)G_3(z)}{1+G_2(z)G_1G_3H(z)}\\\\ &C(z)=\frac{RG(z)}{1+HG(z)}\\\\ &C(z)=\frac{G_1(z)G_2(z)R(z)}{1+G_1(z)G_2(z)H(z)}\\\\ &C(z)=\frac{R(z)G(z)}{1+G(z)H(z)} \end{aligned} C(z)=1+G2(z)G1G3H(z)RG1(z)G2(z)G3(z)C(z)=1+HG(z)RG(z)C(z)=1+G1(z)G2(z)H(z)G1(z)G2(z)R(z)C(z)=1+G(z)H(z)R(z)G(z)
2.2 离散系统稳定性分析
2.2.1 s s s域到 z z z域的映射
在 z z z变换定义中, z = e s T z={\rm e}^{sT} z=esT给出了 s s s域到 z z z域的关系, s s s域中的任意点可表示为 s = σ + j ω s=\sigma+{\rm j}\omega s=σ+jω,映射到 z z z域为:
z = e ( σ + j ω ) T = e σ T e j ω T z={\rm e}^{(\sigma+{\rm j}\omega)T}={\rm e}^{\sigma{T}}{\rm e}^{
{\rm j}\omega{T}} z=e(σ+jω)T=eσTejωT
s s s域到 z z z域的基本映射关系为:
∣ z ∣ = e σ T , ∠ z = ω T |z|={\rm e}^{\sigma{T}},\angle{z}=\omega{T} ∣z∣=eσT,∠z=ωT
令 σ = 0 \sigma=0 σ=0,相当于取 s s s平面的虚轴,当 ω → − ∞ \omega\rightarrow-\infty ω→−∞变到 ∞ \infty ∞时,映射到 z z z平面的轨迹是以原点为圆心的单位圆;
-
等 σ \sigma σ线映射
s s s平面上的等 σ \sigma σ垂线,映射到 z z z平面上的轨迹,是以原定为圆心,以 ∣ z ∣ = e σ T |z|={\rm e}^{\sigma{T}} ∣z∣=eσT为半径的圆,其中: T T T为采样周期;
s s s平面上的虚轴映射为 z z z平面上的单位圆,因此,左半 s s s平面上的等 σ \sigma σ线映射为 z z z平面上的同心圆,在单位圆内;右半 s s s平面上的等 σ \sigma σ线映射为 z z z平面上的同心圆,在单位圆外;
-
等 ω \omega ω线映射
在特定采样周期 T T T情况下, s s s平面上的等 ω \omega ω水平线映射到 z z z平面上的轨迹,是一簇从原点出发的射线,其相角 ∠ z = ω T \angle{z}=\omega{T} ∠z=ωT从正实轴计量,如下图所示:
s s s平面上 ω = ω s / 2 \omega=\omega_s/2 ω=ωs/2水平线,在 z z z平面上正好映射为负实轴;
-
等 ζ \zeta ζ线映射
s s s平面上的等 ζ \zeta ζ线用下式描述:
s = − ω tan β + j ω s=-\omega\tan\beta+{\rm j}\omega s=−ωtanβ+jω
其中: β \beta β为 ζ \zeta ζ线与虚轴之间的夹角;
z = e s T = e ( ω tan β + j ω ) T z={\rm e}^{sT}={\rm e}^{(\omega\tan\beta+{\rm j}\omega)T} z=esT=e(ωtanβ+jω)T
等 ζ \zeta ζ线从 s s s域到 z z z域的映射关系式为:
∣ z ∣ = e − ( 2 π / ω s ) ω tan β , ∠ z = 2 π ω ω s |z|={\rm e}^{-(2\pi/\omega_s)\omega\tan\beta},\angle{z}=\frac{2\pi\omega}{\omega_s} ∣z∣=e−(2π/ωs)ωtanβ,∠z=ωs2πω
除 β = 0 ° \beta=0° β=0°和 β = 90 ° \beta=90° β=90°外,当 β \beta β为常数时,左半 s s s平面上的等 ζ \zeta ζ线,映射为 z z z平面上单位圆内一簇收敛的对数螺旋线,其起点为 z z z平面上正实轴的 1 1 1处,终点为 z z z平面的原点;设 s s s平面上的主要带如下图所示,通过 z = e s T z={\rm e}^{sT} z=esT变换,映射为 z z z平面上的单位圆及单位圆内的负实轴,如下图所示:
s s s平面上所有的次要带,在 z z z平面上映射为相同的单位圆及单位圆内的负实轴;
2.2.2 离散系统稳定的充分必要条件
在线性定常连续系统中,系统稳定的充分必要条件是:系统齐次微分方程的解是收敛的,或系统特征方程式的根均具有负实部,或系统传递函数的极点均位于左半 s s s平面;
-
时域中离散系统稳定的充分必要条件
设线性定常差分方程为:
c ( k ) = − ∑ i = 1 n a i c ( k − i ) + ∑ j = 0 m b j r ( k − j ) c(k)=-\sum_{i=1}^na_ic(k-i)+\sum_{j=0}^mb_jr(k-j) c(k)=−i=1∑naic(k−i)+j=0∑mbjr(k−j)
其齐次差分方程为:
c ( k ) + ∑ i = 1 n a i c ( k − i ) = 0 c(k)+\sum_{i=1}^na_ic(k-i)=0 c(k)+i=1∑naic(k−i)=0
差分方程的特征方程为:
α n + a 1 α n − 1 + a 2 α n − 2 + ⋯ + a n = 0 \alpha^{n}+a_1\alpha^{n-1}+a_2\alpha^{n-2}+\dots+a_n=0 αn+a1αn−1+a2αn−2+⋯+an=0
当特征方程的根 ∣ α i ∣ < 1 ( i = 1 , 2 , … , n ) |\alpha_i|<1(i=1,2,\dots,n) ∣αi∣<1(i=1,2,…,n)时,必有: lim k → ∞ c ( k ) = 0 \displaystyle\lim_{k\rightarrow\infty}c(k)=0 k→∞limc(k)=0,因此,系统稳定的充分必要条件是:当且仅当差分方程所有特征根的模 ∣ α i ∣ < 1 ( i = 1 , 2 , … , n ) |\alpha_i|<1(i=1,2,\dots,n) ∣αi∣<1(i=1,2,…,n),相应的线性定常离散系统是稳定的; -
z z z域中离散系统稳定的充分必要条件
设离散系统特征方程为:
D ( z ) = 1 + G H ( z ) = 0 D(z)=1+GH(z)=0 D(z)=1+GH(z)=0
z z z域中,线性定常离散系统稳定的充分必要条件:当且仅当离散系统特征方程的全部特征根均分布在 z z z平面上的单位圆内,或所有特征根的模均小于1,即 ∣ z i ∣ < 1 ( i = 1 , 2 , … , n ) |z_i|<1(i=1,2,\dots,n) ∣zi∣<1(i=1,2,…,n),相应的线性定常离散系统是稳定的;
2.2.3 离散系统的稳定性判据
-
w w w变换与劳斯稳定判据
令
z = w + 1 w − 1 ⇒ w = z + 1 z − 1 z=\frac{w+1}{w-1}\Rightarrow{w=\frac{z+1}{z-1}} z=w−1w+1⇒w=z−1z+1
复变量 z z z与 w w w互为线性变换, w w w变换称为双线性变换;z z z平面与 w w w平面的对应关系如下图所示:
w w w平面的虚轴对应 z z z平面上的单位,圆周左半 w w w平面对应于 z z z平面上单位圆内的区域,右半 w w w平面对应 z z z平面上单位圆外的区域;
离散系统稳定的充分必要条件:由特征方程 1 + G H ( z ) = 0 1+GH(z)=0 1+GH(z)=0的所有根位于 z z z平面上的单位圆内转换为特征方程 1 + G H ( w ) = 0 1+GH(w)=0 1+GH(w)=0的所有根位于左半 w w w平面;
-
朱利稳定判据
朱利判据是根据离散系统的闭环特征方程 D ( z ) = 0 D(z)=0 D(z)=0的系数,判别其根是否位于 z z z平面上的单位圆内,从而判断该离散系统是否稳定。
设离散系统 n n n阶闭环特征方程为:
D ( z ) = a 0 + a 1 z + a 2 z 2 + ⋯ + a n z n = 0 , a n > 0 D(z)=a_0+a_1z+a_2z^2+\dots+a_nz^n=0,a_n>0 D(z)=a0+a1z+a2z2+⋯+anzn=0,an>0
按照下述方法构造 ( 2 n − 3 ) (2n-3) (2n−3)行、 ( n + 1 ) (n+1) (n+1)列朱利阵列。行数 z 0 z^0 z0 z 1 z^1 z1 z 2 z^2 z2 z 3 z^3 z3 … \dots … z n − k z^{n-k} zn−k … \dots … z n − 1 z^{n-1} zn−1 z n z^n zn 1 1 1 a 0 a_0 a0 a 1 a_1 a1 a 2 a_2 a2 a 3 a_3 a3 … \dots … a n − k a_{n-k} an−k … \dots … a n − 1 a_{n-1} an−1 a n a_n an 2 2 2 a n a_n an a n − 1 a_{n-1} an−1 a n − 2 a_{n-2} an−2 a n − 3 a_{n-3} an−3 … \dots … a k a_k ak … \dots … a 1 a_1 a1 a 0 a_0 a0 3 3 3 b 0 b_0 b0 b 1 b_1 b1 b 2 b_2 b2 b 3 b_3 b3 … \dots … b n − k b_{n-k} bn−k … \dots … b n − 1 b_{n-1} bn−1 4 4 4 b n − 1 b_{n-1} bn−1 b n − 2 b_{n-2} bn−2 b n − 3 b_{n-3} bn−3 b n − 4 b_{n-4} bn−4 … \dots … b k − 1 b_{k-1} bk−1 … \dots … b 0 b_0 b0 5 5 5 c 0 c_0 c0 c 1 c_1 c1 c 2 c_2 c2 c 3 c_3 c3 … \dots … c n − 2 c_{n-2} cn−2 6 6 6 c n − 2 c_{n-2} cn−2 c n − 3 c_{n-3} cn−3 c n − 4 c_{n-4} cn−4 c n − 5 c_{n-5} cn−5 … \dots … c 0 c_0 c0 ⋮ \vdots ⋮ ⋮ \vdots ⋮ ⋮ \vdots ⋮ ⋮ \vdots ⋮ ⋮ \vdots ⋮ 2 n − 5 2n-5 2n−5 p 0 p_0 p0 p 1 p_1 p1 p 2 p_2 p2 p 3 p_3 p3 2 n − 4 2n-4 2n−4 p 3 p_3 p3 p 2 p_2 p2 p 1 p_1 p1 p 0 p_0 p0 2 n − 3 2n-3 2n−3 q 0 q_0 q0 q 1 q_1 q1 q 2 q_2 q2 其中:
b k = ∣ a 0 a n − k a n a k ∣ , k = 0 , 1 , 2 , … , n − 1 c k = ∣ b 0 b n − k − 1 b n − 1 b k ∣ , k = 0 , 1 , 2 , … , n − 2 d k = ∣ c 0 c n − k − 2 c n − 2 c k ∣ , k = 0 , 1 , 2 , … , n − 3 ⋮ q 0 = ∣ p 0 p 3 p 3 p 0 ∣ , q 1 = ∣ p 0 p 2 p 3 p 1 ∣ , q 2 = ∣ p 0 p 1 p 3 p 2 ∣ , \begin{aligned} &b_k=\begin{vmatrix} a_0 & a_{n-k}\\ a_n & a_k \end{vmatrix},&k=0,1,2,\dots,n-1\\\\ &c_k= \begin{vmatrix} b_0 & b_{n-k-1}\\ b_{n-1} & b_k \end{vmatrix},&k=0,1,2,\dots,n-2\\\\ &d_k= \begin{vmatrix} c_0 & c_{n-k-2}\\ c_{n-2} & c_k \end{vmatrix},&k=0,1,2,\dots,n-3\\\\ &\space\vdots\\\\ &q_0= \begin{vmatrix} p_0 & p_{3}\\ p_3 & p_0 \end{vmatrix}, q_1= \begin{vmatrix} p_0 & p_{2}\\ p_3 & p_1 \end{vmatrix}, q_2= \begin{vmatrix} p_0 & p_{1}\\ p_3 & p_2 \end{vmatrix}, \end{aligned} bk= a0anan−kak ,ck= b0bn−1bn−k−1bk ,dk= c0cn−2cn−k−2ck , ⋮q0= p0p3p3p0 ,q1= p0p3p2p1 ,q2= p0p3p1p2 ,k=0,1,2,…,n−1k=0,1,2,…,n−2k=0,1,2,…,n−3
朱利稳定判据:特征方程 D ( z ) = 0 D(z)=0 D(z)=0的根,全部位于 z z z平面上单位圆内的充分必要条件是:
D ( 1 ) > 0 , D ( − 1 ) { > 0 ,当 n 为偶数时 < 0 ,当 n 为奇数时 D(1)>0,D(-1)\begin{cases}>0,当n为偶数时\\<0,当n为奇数时\end{cases} D(1)>0,D(−1){ >0,当n为偶数时<0,当n为奇数时
及以下 n − 1 n-1 n−1个约束条件成立:
∣ a 0 ∣ < a n , ∣ b 0 ∣ > ∣ b n − 1 ∣ , ∣ c 0 ∣ > ∣ c n − 2 ∣ , ∣ d 0 ∣ > ∣ d n − 3 ∣ , … , ∣ q 0 ∣ > ∣ q 2 ∣ |a_0|<a_n,|b_0|>|b_{n-1}|,|c_0|>|c_{n-2}|,|d_0|>|d_{n-3}|,\dots,|q_0|>|q_2| ∣a0∣<an,∣b0∣>∣bn−1∣,∣c0∣>∣cn−2∣,∣d0∣>∣dn−3∣,…,∣q0∣>∣q2∣
2.2.4 离散系统的稳态误差
设单位反馈误差采样系统如下图所示:
其中: G ( s ) G(s) G(s)为连续部分的传递函数, e ( t ) e(t) e(t)为系统连续误差信号, e ∗ ( t ) e^*(t) e∗(t)为系统采样误差信号,其 z z z变换为:
E ( z ) = R ( z ) − C ( z ) = [ 1 − Φ ( z ) ] R ( z ) = Φ e ( z ) R ( z ) E(z)=R(z)-C(z)=[1-\Phi(z)]R(z)=\Phi_e(z)R(z) E(z)=R(z)−C(z)=[1−Φ(z)]R(z)=Φe(z)R(z)
其中:
Φ e ( z ) = E ( z ) R ( z ) = 1 1 + G ( z ) \Phi_e(z)=\frac{E(z)}{R(z)}=\frac{1}{1+G(z)} Φe(z)=R(z)E(z)=1+G(z)1
如果 Φ e ( z ) \Phi_e(z) Φe(z)的极点全部位于 z z z平面上的单位圆内,即若离散系统是稳定的,可用 z z z变换的终值定理求出采样瞬时的稳态误差:
e s s ( ∞ ) = lim t → ∞ e ∗ ( t ) = lim z → 1 ( 1 − z − 1 ) E ( z ) = lim z → 1 ( z − 1 ) R ( z ) z [ 1 + G ( z ) ] e_{ss}(\infty)=\lim_{t\rightarrow\infty}e^*(t)=\lim_{z\rightarrow1}(1-z^{-1})E(z)=\lim_{z\rightarrow1}\frac{(z-1)R(z)}{z[1+G(z)]} ess(∞)=t→∞lime∗(t)=z→1lim(1−z−1)E(z)=z→1limz[1+G(z)](z−1)R(z)
2.2.5 离散系统的型别与静态误差系数
在离散系统中,把开环脉冲传递函数 G ( z ) G(z) G(z)具有 z = 1 z=1 z=1的极点数 ν \nu ν作为划分离散系统型别的标准,把 G ( z ) G(z) G(z)中 ν = 0 , 1 , 2 , … \nu=0,1,2,\dots ν=0,1,2,…的系统,称为 0 0 0型、Ⅰ型和Ⅱ型离散系统等;
-
单位阶跃输入时的稳态误差
当系统输入为单位阶跃函数 r ( t ) = 1 ( t ) r(t)=1(t) r(t)=1(t)时, z z z变换函数为:
R ( z ) = z z − 1 R(z)=\frac{z}{z-1} R(z)=z−1z
稳态误差为:
e s s ( ∞ ) = lim z → 1 1 1 + G ( z ) = 1 lim z → 1 [ 1 + G ( z ) ] = 1 K p e_{ss}(\infty)=\lim_{z\rightarrow1}\frac{1}{1+G(z)}=\frac{1}{\displaystyle\lim_{z\rightarrow1}[1+G(z)]}=\frac{1}{K_p} ess(∞)=z→1lim1+G(z)1=z→1lim[1+G(z)]1=Kp1
其中静态位置误差系数为:
K p = lim z → 1 [ 1 + G ( z ) ] K_p=\lim_{z\rightarrow1}[1+G(z)] Kp=z→1lim[1+G(z)]
在单位阶跃函数作用下, 0 0 0型离散系统在采样瞬时存在位置误差,Ⅰ型及Ⅰ型以上的离散系统,在采样瞬时没有位置误差; -
单位斜坡输入时的稳态误差
当系统输入为单位斜坡函数 r ( t ) = t r(t)=t r(t)=t时, z z z变换函数为:
R ( z ) = T z ( z − 1 ) 2 R(z)=\frac{Tz}{(z-1)^2} R(z)=(z−1)2Tz
稳态误差为:
e s s ( ∞ ) = lim z → 1 T ( z − 1 ) [ 1 + G ( z ) ] = T lim z → 1 ( z − 1 ) G ( z ) = T K v e_{ss}(\infty)=\lim_{z\rightarrow1}\frac{T}{(z-1)[1+G(z)]}=\frac{T}{\displaystyle\lim_{z\rightarrow1}(z-1)G(z)}=\frac{T}{K_v} ess(∞)=z→1lim(z−1)[1+G(z)]T=z→1lim(z−1)G(z)T=KvT
其中静态速度误差系数为:
K v = lim z → 1 ( z − 1 ) G ( z ) K_v=\lim_{z\rightarrow1}(z-1)G(z) Kv=z→1lim(z−1)G(z)
0 0 0型离散系统不能承受单位斜坡函数作用,Ⅰ型离散系统在单位斜坡函数作用下存在速度误差,Ⅱ型及Ⅱ型以上离散系统在单位斜坡函数作用下不存在稳态误差; -
单位加速度输入时的稳态误差
当系统输入为单位加速度函数 r ( t ) = t 2 / 2 r(t)=t^2/2 r(t)=t2/2时, z z z变换函数为:
R ( z ) = T 2 z ( z + 1 ) 2 ( z − 1 ) 3 R(z)=\frac{T^2z(z+1)}{2(z-1)^3} R(z)=2(z−1)3T2z(z+1)
稳态误差为:
e s s ( ∞ ) = lim z → 1 T 2 ( z + 1 ) 2 ( z − 1 ) 2 [ 1 + G ( z ) ] = T 2 lim z → 1 ( z − 1 ) 2 G ( z ) = T 2 K a e_{ss}(\infty)=\lim_{z\rightarrow1}\frac{T^2(z+1)}{2(z-1)^2[1+G(z)]}=\frac{T^2}{\displaystyle\lim_{z\rightarrow1}(z-1)^2G(z)}=\frac{T^2}{K_a} ess(∞)=z→1lim2(z−1)2[1+G(z)]T2(z+1)=z→1lim(z−1)2G(z)T2=KaT2
其中静态加速度误差系数为:
K a = lim z → 1 ( z − 1 ) 2 G ( z ) K_a=\lim_{z\rightarrow1}(z-1)^2G(z) Ka=z→1lim(z−1)2G(z)
0 0 0型和Ⅰ型离散系统不能承受单位加速度函数作用,Ⅱ型离散系统在单位加速度函数作用下存在加速度误差,Ⅲ型及Ⅲ型以上的离散系统在单位加速度函数作用下,不存在采样瞬时的稳态误差;单位反馈离散系统的稳态误差小结:
系统类型 位置误差 r ( t ) = 1 ( t ) r(t)=1(t) r(t)=1(t) 速度误差 r ( t ) = t r(t)=t r(t)=t 加速度误差 r ( t ) = 1 2 t 2 r(t)=\frac{1}{2}t^2 r(t)=21t2 0 0 0型 1 K p \displaystyle\frac{1}{K_p} Kp1 ∞ \infty ∞ ∞ \infty ∞ Ⅰ型 0 0 0 T K v \displaystyle\frac{T}{K_v} KvT ∞ \infty ∞ Ⅱ型 0 0 0 0 0 0 T 2 K a \displaystyle\frac{T^2}{K_a} KaT2 Ⅲ型 0 0 0 0 0 0 0 0 0
2.3 离散系统的数字校正
2.3.1 数字控制器的脉冲传递函数
设离散系统如下图所示:
其中: D ( z ) D(z) D(z)为数字控制器的脉冲传递函数, G ( z ) G(z) G(z)为保持器与被控对象的传递函数, H ( s ) H(s) H(s)为反馈测量装置的传递函数;
设 H ( s ) = 1 H(s)=1 H(s)=1, G ( s ) G(s) G(s)的 z z z变换为 G ( z ) G(z) G(z),可得系统的闭环脉冲传递函数:
Φ ( z ) = D ( z ) G ( z ) 1 + D ( z ) G ( z ) = C ( z ) R ( z ) \Phi(z)=\frac{D(z)G(z)}{1+D(z)G(z)}=\frac{C(z)}{R(z)} Φ(z)=1+D(z)G(z)D(z)G(z)=R(z)C(z)
误差脉冲传递函数:
Φ e ( z ) = 1 1 + D ( z ) G ( z ) = E ( z ) R ( z ) \Phi_e(z)=\frac{1}{1+D(z)G(z)}=\frac{E(z)}{R(z)} Φe(z)=1+D(z)G(z)1=R(z)E(z)
可得数字控制器的脉冲传递函数为:
D ( z ) = Φ ( z ) G ( z ) [ 1 − Φ ( z ) ] = 1 − Φ e ( z ) G ( z ) Φ e ( z ) , Φ e ( z ) = 1 − Φ ( z ) D(z)=\frac{\Phi(z)}{G(z)[1-\Phi(z)]}=\frac{1-\Phi_e(z)}{G(z)\Phi_e(z)},\Phi_e(z)=1-\Phi(z) D(z)=G(z)[1−Φ(z)]Φ(z)=G(z)Φe(z)1−Φe(z),Φe(z)=1−Φ(z)
2.3.2 最少拍系统设计
最少拍系统:指在典型输入作用下,能以有限拍结束响应过程,且在采样时刻上无稳态误差的离散系统;
常见的典型输入:单位阶跃函数、单位速度函数、单位加速度函数;
Z [ 1 ( t ) ] = z z − 1 = 1 1 − z − 1 Z [ t ] = T z ( z − 1 ) 2 = T z − 1 ( 1 − z − 1 ) 2 Z [ 1 2 t 2 ] = T 2 z ( z + 1 ) 2 ( z − 1 ) 3 = 1 2 T 2 z − 1 ( 1 + z − 1 ) ( 1 − z − 1 ) 3 \begin{aligned} &Z\left[1(t)\right]=\frac{z}{z-1}=\frac{1}{1-z^{-1}}\\\\ &Z\left[t\right]=\frac{Tz}{(z-1)^2}=\frac{Tz^{-1}}{(1-z^{-1})^2}\\\\ &Z\left[\frac{1}{2}t^2\right]=\frac{T^2z(z+1)}{2(z-1)^3}=\frac{\displaystyle\frac{1}{2}T^2z^{-1}(1+z^{-1})}{(1-z^{-1})^3} \end{aligned} Z[1(t)]=z−1z=1−z−11Z[t]=(z−1)2Tz=(1−z−1)2Tz−1Z[21t2]=2(z−1)3T2z(z+1)=(1−z−1)321T2z−1(1+z−1)
典型输入可表示为一般形式:
R ( z ) = A ( z ) ( 1 − z − 1 ) m R(z)=\frac{A(z)}{(1-z^{-1})^m} R(z)=(1−z−1)mA(z)
其中: A ( z ) A(z) A(z)是不含 ( 1 − z − 1 ) (1-z^{-1}) (1−z−1)因子的 z − 1 z^{-1} z−1多项式;
最少拍系统的设计原则:若系统广义被控对象 G ( z ) G(z) G(z)无延迟且在 z z z平面单位圆上及单位圆外无零极点,要求选择闭环脉冲传递函数 Φ ( z ) \Phi(z) Φ(z),使系统在典型输入作用下,经最少采样周期后能使输出序列在各采样时刻的稳态误差为零,达到完全跟踪的目的,从而确定所需要的数字控制器的脉冲传递函数 D ( z ) D(z) D(z);
误差信号 e ( t ) e(t) e(t)的 z z z变换为:
E ( z ) = Φ e ( z ) R ( z ) = Φ e ( z ) A ( z ) ( 1 − z − 1 ) m E(z)=\Phi_e(z)R(z)=\frac{\Phi_e(z)A(z)}{(1-z^{-1})^m} E(z)=Φe(z)R(z)=(1−z−1)mΦe(z)A(z)
离散系统的稳态误差为:
e s s ( ∞ ) = lim z → 1 ( 1 − z − 1 ) E ( z ) = lim z → 1 ( 1 − z − 1 ) A ( z ) ( 1 − z − 1 ) m Φ e ( z ) e_{ss}(\infty)=\lim_{z\rightarrow1}(1-z^{-1})E(z)=\lim_{z\rightarrow1}(1-z^{-1})\frac{A(z)}{(1-z^{-1})^m}\Phi_e(z) ess(∞)=z→1lim(1−z−1)E(z)=z→1lim(1−z−1)(1−z−1)mA(z)Φe(z)
使 e s s ( ∞ ) e_{ss}(\infty) ess(∞)为零的条件是: Φ e ( z ) \Phi_e(z) Φe(z)中包含 ( 1 − z − 1 ) m (1-z^{-1})^m (1−z−1)m的因子,即:
Φ e ( z ) = ( 1 − z − 1 ) m F ( z ) \Phi_e(z)=(1-z^{-1})^mF(z) Φe(z)=(1−z−1)mF(z)
其中: F ( z ) F(z) F(z)不含 ( 1 − z − 1 ) (1-z^{-1}) (1−z−1)因子的多项式,为了 D ( z ) D(z) D(z)简单,阶数最低,取 F ( z ) = 1 F(z)=1 F(z)=1。
-
单位阶跃输入
由 r ( t ) = 1 ( t ) r(t)=1(t) r(t)=1(t)时有 m = 1 , A ( z ) = 1 m=1,A(z)=1 m=1,A(z)=1可得:
Φ e ( z ) = 1 − z − 1 , Φ ( z ) = z − 1 \Phi_e(z)=1-z^{-1},\Phi(z)=z^{-1} Φe(z)=1−z−1,Φ(z)=z−1
数字控制器脉冲传递函数为:
D ( z ) = z − 1 ( 1 − z − 1 ) G ( z ) , E ( z ) = A ( z ) ( 1 − z − 1 ) m Φ e ( z ) = 1 D(z)=\frac{z^{-1}}{(1-z^{-1})G(z)},E(z)=\frac{A(z)}{(1-z^{-1})^m}\Phi_e(z)=1 D(z)=(1−z−1)G(z)z−1,E(z)=(1−z−1)mA(z)Φe(z)=1
e ( 0 ) = 1 , e ( T ) = e ( 2 T ) = ⋯ = 0 e(0)=1,e(T)=e(2T)=\dots=0 e(0)=1,e(T)=e(2T)=⋯=0,可见,最少拍系统经过一拍即可完全跟踪输入 r ( t ) = 1 ( t ) r(t)=1(t) r(t)=1(t),这样的离散系统称为一拍系统, t s = T t_s=T ts=T; -
单位斜坡输入
由 r ( t ) = t r(t)=t r(t)=t时,有 m = 2 , A ( z ) = T z − 1 m=2,A(z)=Tz^{-1} m=2,A(z)=Tz−1,可得:
Φ e ( z ) = ( 1 − z − 1 ) m F ( z ) = ( 1 − z − 1 ) 2 , Φ ( z ) = 1 − Φ e ( z ) = 2 z − 1 − z − 2 \Phi_e(z)=(1-z^{-1})^mF(z)=(1-z^{-1})^{2},\Phi(z)=1-\Phi_e(z)=2z^{-1}-z^{-2} Φe(z)=(1−z−1)mF(z)=(1−z−1)2,Φ(z)=1−Φe(z)=2z−1−z−2
数字控制器脉冲传递函数为:
D ( z ) = Φ ( z ) G ( z ) Φ e ( z ) = z − 1 ( 2 − z − 1 ) ( 1 − z − 1 ) 2 G ( z ) , E ( z ) = A ( z ) ( 1 − z − 1 ) m Φ e ( z ) = T z − 1 D(z)=\frac{\Phi(z)}{G(z)\Phi_e(z)}=\frac{z^{-1}(2-z^{-1})}{(1-z^{-1})^2G(z)},E(z)=\frac{A(z)}{(1-z^{-1})^m}\Phi_e(z)=Tz^{-1} D(z)=G(z)Φe(z)Φ(z)=(1−z−1)2G(z)z−1(2−z−1),E(z)=(1−z−1)mA(z)Φe(z)=Tz−1
e ( 0 ) = 0 , e ( T ) = T , e ( 2 T ) = e ( 3 T ) = ⋯ = 0 e(0)=0,e(T)=T,e(2T)=e(3T)=\dots=0 e(0)=0,e(T)=T,e(2T)=e(3T)=⋯=0,最少拍系统经过二拍即可完全跟踪输入 r ( t ) = t r(t)=t r(t)=t,这样的离散系统称为二拍系统,其调节时间 t s = 2 T t_s=2T ts=2T; -
单位加速度输入
由于 r ( t ) = t 2 / 2 r(t)=t^2/2 r(t)=t2/2时,有 m = 3 , A ( z ) = 1 2 T 2 z − 1 ( 1 + z − 1 ) m=3,A(z)=\displaystyle\frac{1}{2}T^2z^{-1}(1+z^{-1}) m=3,A(z)=21T2z−1(1+z−1),可得:
Φ e ( z ) = ( 1 − z − 1 ) 3 , Φ ( z ) = 3 z − 1 − 3 z − 2 + z − 3 \Phi_e(z)=(1-z^{-1})^3,\Phi(z)=3z^{-1}-3z^{-2}+z^{-3} Φe(z)=(1−z−1)3,Φ(z)=3z−1−3z−2+z−3
数字控制脉冲传递函数为:
D ( z ) = z − 1 ( 3 − 3 z − 1 + z − 2 ) ( 1 − z − 1 ) 3 G ( z ) , E ( z ) = A ( z ) = 1 2 T 2 z − 1 + 1 2 T 2 z − 2 D(z)=\frac{z^{-1}(3-3z^{-1}+z^{-2})}{(1-z^{-1})^3G(z)},E(z)=A(z)=\frac{1}{2}T^2z^{-1}+\frac{1}{2}T^2z^{-2} D(z)=(1−z−1)3G(z)z−1(3−3z−1+z−2),E(z)=A(z)=21T2z−1+21T2z−2
输出脉冲序列为:
C ( z ) = Φ ( z ) R ( z ) = 3 2 T 2 z − 2 + 9 2 T 2 z − 3 + ⋯ + n 2 2 T 2 z − n + ⋯ + C(z)=\Phi(z)R(z)=\frac{3}{2}T^2z^{-2}+\frac{9}{2}T^2z^{-3}+\dots+\frac{n^2}{2}T^2z^{-n}+\dots+ C(z)=Φ(z)R(z)=23T2z−2+29T2z−3+⋯+2n2T2z−n+⋯+
有:
e ( 0 ) = 0 , e ( T ) = 1 2 T 2 , e ( 2 T ) = 1 2 T 2 , e ( 3 T ) = e ( 4 T ) = ⋯ = 0 c ( 0 ) = c ( T ) = 0 , c ( 2 T ) = 1.5 T 2 , c ( 3 T ) = 4.5 T 2 , ⋯ \begin{aligned} &e(0)=0,e(T)=\frac{1}{2}T^2,e(2T)=\frac{1}{2}T^2,e(3T)=e(4T)=\dots=0\\\\ &c(0)=c(T)=0,c(2T)=1.5T^2,c(3T)=4.5T^2,\cdots \end{aligned} e(0)=0,e(T)=21T2,e(2T)=21T2,e(3T)=e(4T)=⋯=0c(0)=c(T)=0,c(2T)=1.5T2,c(3T)=4.5T2,⋯
最少拍系统经过三拍即可完全跟踪输入 r ( t ) = t 2 / 2 r(t)=t^2/2 r(t)=t2/2,这样的离散系统称为三拍系统,调节时间为 t s = 3 T t_s=3T ts=3T; -
最少拍系统的设计结果
小结:
- 从快速性而言,按单位斜坡输入设计的最少拍系统,在各种典型输入作用下,其动态过程均为二拍;
- 从准确性而言,系统对单位阶跃输入和单位斜坡输入,在采样时刻均无稳态误差,对单位加速度输入,采样时刻上的稳态误差为常量 T T T;
- 从动态性能而言,系统对单位斜坡输入下的响应性能较好,因为系统本身就是针对此设计的,但系统对单位阶跃输入响应性能较差,有 100 % 100\% 100%的超调量,因此按某种典型输入设计的最少拍系统,适应性较差;
- 从平稳性而言,在各种典型输入作用下系统进入稳态后,在非采样时刻一般均存在纹波,从而增加系统的机械磨损;
-
实例分析
Example8: 设单位反馈线性定常离散系统的连续部分和零阶保持器的传递函数分别为:
G 0 ( s ) = 10 s ( s + 1 ) , G h ( s ) = 1 − e − s T s G_0(s)=\frac{10}{s(s+1)},G_h(s)=\frac{1-{\rm e}^{-sT}}{s} G0(s)=s(s+1)10,Gh(s)=s1−e−sT
其中采样周期为 T = 1 s T=1s T=1s;要求系统在单位斜坡输入时实现最少拍控制,求数字控制器脉冲传递函数 D ( z ) D(z) D(z)。解:
系统开环传递函数:
G ( s ) = G 0 ( s ) G h ( s ) = 10 ( 1 − e − s T ) s 2 ( s + 1 ) G(s)=G_0(s)G_h(s)=\frac{10(1-e^{-sT})}{s^2(s+1)} G(s)=G0(s)Gh(s)=s2(s+1)10(1−e−sT)
因为:
Z [ 1 s 2 ( s + 1 ) ] = T z ( z − 1 ) 2 − ( 1 − e − T ) z ( z − 1 ) ( z − e − T ) Z\left[\frac{1}{s^2(s+1)}\right]=\frac{Tz}{(z-1)^2}-\frac{(1-{\rm e}^{-T})z}{(z-1)(z-{\rm e}^{-T})} Z[s2(s+1)1]=(z−1)2Tz−(z−1)(z−e−T)(1−e−T)z
因此有:
G ( z ) = 10 ( 1 − z − 1 ) [ T z ( z − 1 ) 2 − ( 1 − e − T ) z ( z − 1 ) ( z − e − T ) ] = 3.68 z − 1 ( 1 + 0.717 z − 1 ) ( 1 − z − 1 ) ( 1 − 0.368 z − 1 ) G(z)=10(1-z^{-1})\left[\frac{Tz}{(z-1)^2}-\frac{(1-{\rm e}^{-T})z}{(z-1)(z-{\rm e}^{-T})}\right]=\frac{3.68z^{-1}(1+0.717z^{-1})}{(1-z^{-1})(1-0.368z^{-1})} G(z)=10(1−z−1)[(z−1)2Tz−(z−1)(z−e−T)(1−e−T)z]=(1−z−1)(1−0.368z−1)3.68z−1(1+0.717z−1)
其中, r ( t ) = t r(t)=t r(t)=t,闭环脉冲传递函数和误差脉冲传递函数:
Φ ( z ) = 2 z − 1 ( 1 − 0.5 z − 1 ) , Φ e ( z ) = ( 1 − z − 1 ) 2 \Phi(z)=2z^{-1}(1-0.5z^{-1}),\Phi_e(z)=(1-z^{-1})^2 Φ(z)=2z−1(1−0.5z−1),Φe(z)=(1−z−1)2
可得:
D ( z ) = Φ ( z ) G ( z ) Φ e ( z ) = 0.543 ( 1 − 0.368 z − 1 ) ( 1 − 0.5 z − 1 ) ( 1 − z − 1 ) ( 1 + 0.717 z − 1 ) D(z)=\frac{\Phi(z)}{G(z)\Phi_e(z)}=\frac{0.543(1-0.368z^{-1})(1-0.5z^{-1})}{(1-z^{-1})(1+0.717z^{-1})} D(z)=G(z)Φe(z)Φ(z)=(1−z−1)(1+0.717z−1)0.543(1−0.368z−1)(1−0.5z−1)
2.3.3 无纹波最少拍系统设计
无纹波最少拍系统的设计要求:在某一种典型输入作用下设计的系统,其输出响应经过尽可能少的采样周期后,不仅在采样时刻输出可以完全跟踪输入,且在非采样时刻不存在纹波;
-
无纹波最少拍系统的必要条件
为了在稳态过程中获得无纹波的平滑输出 c ∗ ( t ) c^*(t) c∗(t),被控对象 G 0 ( s ) G_0(s) G0(s)必须有能力给出与输入 r ( t ) r(t) r(t)相同的平滑输出 c ( t ) c(t) c(t);
若针对单位斜坡输入 r ( t ) = t r(t)=t r(t)=t设计最少拍系统,则 G 0 ( s ) G_0(s) G0(s)的稳态输出也必须是斜坡函数,因此, G 0 ( s ) G_0(s) G0(s)必须至少有一个积分环节,使被控对象在零阶保持器常值输出信号作用下,稳态输出为等速变化量;若针对单位加速度输入 r ( t ) = t 2 / 2 r(t)=t^2/2 r(t)=t2/2设计最少拍系统,则 G 0 ( s ) G_0(s) G0(s)至少应包括两个积分环节;
若输入信号为:
r ( t ) = R 0 + R 1 t + 1 2 t 2 + ⋯ + 1 ( q − 1 ) ! R q − 1 t q − 1 r(t)=R_0+R_1t+\frac{1}{2}t^2+\dots+\frac{1}{(q-1)!}R_{q-1}t^{q-1} r(t)=R0+R1t+21t2+⋯+(q−1)!1Rq−1tq−1
无纹波最少拍系统的必要条件:被控对象传递函数 G 0 ( s ) G_0(s) G0(s)中,至少包含 ( q − 1 ) (q-1) (q−1)个积分环节; -
无纹波最少拍系统的附加条件
D ( z ) = Φ ( z ) G ( z ) Φ e ( z ) ⇒ D ( z ) Φ e ( z ) = Φ ( z ) G ( z ) D(z)=\frac{\Phi(z)}{G(z)\Phi_e(z)}\Rightarrow{D(z)\Phi_e(z)=\frac{\Phi(z)}{G(z)}} D(z)=G(z)Φe(z)Φ(z)⇒D(z)Φe(z)=G(z)Φ(z)
设广义对象脉冲传递函数为:
G ( z ) = P ( z ) Q ( z ) G(z)=\frac{P(z)}{Q(z)} G(z)=Q(z)P(z)
其中: P ( z ) P(z) P(z)为 G ( z ) G(z) G(z)的零点多项式, Q ( z ) Q(z) Q(z)为 G ( z ) G(z) G(z)的极点多项式;则有:
D ( z ) Φ e ( z ) = Φ ( z ) Q ( z ) P ( z ) D(z)\Phi_e(z)=\frac{\Phi(z)Q(z)}{P(z)} D(z)Φe(z)=P(z)Φ(z)Q(z)
D ( z ) Φ e ( z ) D(z)\Phi_e(z) D(z)Φe(z)成为 z − 1 z^{-1} z−1有限多项式的条件是: Φ ( z ) \Phi(z) Φ(z)的零点应抵消 G ( z ) G(z) G(z)的全部零点,即有:
Φ ( z ) = P ( z ) M ( z ) \Phi(z)=P(z)M(z) Φ(z)=P(z)M(z)
其中: M ( z ) M(z) M(z)为待定 z − 1 z^{-1} z−1多项式;上式为无纹波最少拍系统的附加条件。 -
无纹波最少拍系统设计
Example9: 设单位反馈线性定常离散系统的连续部分和零阶保持器的传递函数分别为:
G 0 ( s ) = 10 s ( s + 1 ) , G h ( s ) = 1 − e − s T s G_0(s)=\frac{10}{s(s+1)},G_h(s)=\frac{1-{\rm e}^{-sT}}{s} G0(s)=s(s+1)10,Gh(s)=s1−e−sT
其中采样周期为 T = 1 s T=1s T=1s;要求系统在单位斜坡输入时实现无纹波最少拍控制,求数字控制器脉冲传递函数 D ( z ) D(z) D(z)。解:
系统开环传递函数:
G ( s ) = G 0 ( s ) G h ( s ) = 10 ( 1 − e − s T ) s 2 ( s + 1 ) G(s)=G_0(s)G_h(s)=\frac{10(1-{\rm e}^{-sT})}{s^2(s+1)} G(s)=G0(s)Gh(s)=s2(s+1)10(1−e−sT)
因为:
Z [ 1 s 2 ( s + 1 ) ] = T z ( z − 1 ) 2 − ( 1 − e − T ) z ( z − 1 ) ( z − e − T ) Z\left[\frac{1}{s^2(s+1)}\right]=\frac{Tz}{(z-1)^2}-\frac{(1-{\rm e}^{-T})z}{(z-1)(z-{\rm e}^{-T})} Z[s2(s+1)1]=(z−1)2Tz−(z−1)(z−e−T)(1−e−T)z
因此有:
G ( z ) = 10 ( 1 − z − 1 ) [ T z ( z − 1 ) 2 − ( 1 − e − T ) z ( z − 1 ) ( z − e − T ) ] = 3.68 z − 1 ( 1 + 0.717 z − 1 ) ( 1 − z − 1 ) ( 1 − 0.368 z − 1 ) G(z)=10(1-z^{-1})\left[\frac{Tz}{(z-1)^2}-\frac{(1-e^{-T})z}{(z-1)(z-e^{-T})}\right]=\frac{3.68z^{-1}(1+0.717z^{-1})}{(1-z^{-1})(1-0.368z^{-1})} G(z)=10(1−z−1)[(z−1)2Tz−(z−1)(z−e−T)(1−e−T)z]=(1−z−1)(1−0.368z−1)3.68z−1(1+0.717z−1)
可见, G ( z ) G(z) G(z)有一个零点 z = − 0.717 z=-0.717 z=−0.717,有一个延迟因子 z − 1 z^{-1} z−1,且在单位圆上有一个极点 z = 1 z=1 z=1;零点补偿,根据无纹波附加条件, G ( z ) G(z) G(z)中 z = − 0.717 z=-0.717 z=−0.717零点应被 Φ ( z ) \Phi(z) Φ(z)零点对消;
因此,令 M ( z ) = a + b z − 1 M(z)=a+bz^{-1} M(z)=a+bz−1,其中 a 、 b a、b a、b待定,选择:
Φ ( z ) = z − 1 ( 1 + 0.717 z − 1 ) ( a + b z − 1 ) \Phi(z)=z^{-1}(1+0.717z^{-1})(a+bz^{-1}) Φ(z)=z−1(1+0.717z−1)(a+bz−1)
由最少拍条件下,在单位斜坡输入下,
Φ e ( z ) = ( 1 − z − 1 ) 2 , Φ ( z ) = 2 z − 1 ( 1 − 0.5 z − 1 ) \Phi_e(z)=(1-z^{-1})^2,\Phi(z)=2z^{-1}(1-0.5z^{-1}) Φe(z)=(1−z−1)2,Φ(z)=2z−1(1−0.5z−1)
因无纹波时,要求 Φ ( z ) \Phi(z) Φ(z)比有纹波时增加一阶,选择:
Φ e ( z ) = ( 1 − z − 1 ) 2 ( 1 + c z − 1 ) \Phi_e(z)=(1-z^{-1})^2(1+cz^{-1}) Φe(z)=(1−z−1)2(1+cz−1)
其中: c c c待定;可得:
1 − Φ e ( z ) = ( 2 − c ) z − 1 + ( 2 c − 1 ) z − 2 − c z − 3 1-\Phi_e(z)=(2-c)z^{-1}+(2c-1)z^{-2}-cz^{-3} 1−Φe(z)=(2−c)z−1+(2c−1)z−2−cz−3Φ ( z ) = z − 1 ( 1 + 0.717 z − 1 ) ( a + b z − 1 ) = a z − 1 + ( b + 0.717 a ) z − 2 + 0.717 b z − 3 \Phi(z)=z^{-1}(1+0.717z^{-1})(a+bz^{-1})=az^{-1}+(b+0.717a)z^{-2}+0.717bz^{-3} Φ(z)=z−1(1+0.717z−1)(a+bz−1)=az−1+(b+0.717a)z−2+0.717bz−3
对应系数相等,可得:
a = 1.408 , b = − 0.826 , c = 0.592 a=1.408,b=-0.826,c=0.592 a=1.408,b=−0.826,c=0.592
因此:
Φ ( z ) = 1.408 z − 1 ( 1 + 0.717 z − 1 ) ( 1 − 0.587 z − 1 ) \Phi(z)=1.408z^{-1}(1+0.717z^{-1})(1-0.587z^{-1}) Φ(z)=1.408z−1(1+0.717z−1)(1−0.587z−1)Φ e ( z ) = ( 1 − z − 1 ) 2 ( 1 + 0.592 z − 1 ) \Phi_e(z)=(1-z^{-1})^2(1+0.592z^{-1}) Φe(z)=(1−z−1)2(1+0.592z−1)
可得:
D ( z ) = Φ ( z ) G ( z ) Φ e ( z ) = 0.383 ( 1 − 0.368 z − 1 ) ( 1 − 0.587 z − 1 ) ( 1 − z − 1 ) ( 1 + 0.592 z − 1 ) D(z)=\frac{\Phi(z)}{G(z)\Phi_e(z)}=\frac{0.383(1-0.368z^{-1})(1-0.587z^{-1})}{(1-z^{-1})(1+0.592z^{-1})} D(z)=G(z)Φe(z)Φ(z)=(1−z−1)(1+0.592z−1)0.383(1−0.368z−1)(1−0.587z−1)
2.3.4 PID数字控制器的实现
P I D {\rm PID} PID控制器的传递函数为:
D ( s ) = U ( s ) X ( s ) = K 1 + K 2 s + K 3 s D(s)=\frac{U(s)}{X(s)}=K_1+\frac{K_2}{s}+K_3s D(s)=X(s)U(s)=K1+sK2+K3s
将其中的微分项和积分项进行离散化处理,可以得到 P I D {\rm PID} PID控制器的数字实现;
u ( k T ) = d x d t ∣ t = k T = 1 T { x ( k T ) − x [ ( k − 1 ) T ] } u(kT)=\left.\frac{
{\rm d}x}{
{\rm d}t}\right|_{t=kT}=\frac{1}{T}\left\{x(kT)-x[(k-1)T]\right\} u(kT)=dtdx
t=kT=T1{
x(kT)−x[(k−1)T]}
即:
U ( z ) = ( 1 − z − 1 ) T X ( z ) = z − 1 T z X ( z ) U(z)=\frac{(1-z^{-1})}{T}X(z)=\frac{z-1}{Tz}X(z) U(z)=T(1−z−1)X(z)=Tzz−1X(z)
u ( k T ) = u [ ( k − 1 ) T ] + T x ( k T ) u(kT)=u[(k-1)T]+Tx(kT) u(kT)=u[(k−1)T]+Tx(kT)
即:
U ( z ) = z − 1 U ( z ) + T X ( z ) U(z)=z^{-1}U(z)+TX(z) U(z)=z−1U(z)+TX(z)
整理可得:
U ( z ) = T z z − 1 X ( z ) U(z)=\frac{Tz}{z-1}X(z) U(z)=z−1TzX(z)
P I D {\rm PID} PID控制器在 z z z域的传递函数为:
D ( z ) = U ( z ) X ( z ) = K 1 + K 2 T z z − 1 + K 3 z − 1 T z D(z)=\frac{U(z)}{X(z)}=K_1+K_2\frac{Tz}{z-1}+K_3\frac{z-1}{Tz} D(z)=X(z)U(z)=K1+K2z−1Tz+K3Tzz−1
记 x ( k T ) = x ( k ) x(kT)=x(k) x(kT)=x(k),可得 P I D {\rm PID} PID控制器的差分方程:
u ( k ) = K 1 x ( k ) + K 2 [ u ( k − 1 ) + T x ( k ) ] + K 3 T [ x ( k ) − x ( k − 1 ) ] = [ K 1 + K 2 T + K 3 T ] x ( k ) − K 3 T x ( k − 1 ) + K 2 u ( k − 1 ) \begin{aligned} u(k)&=K_1x(k)+K_2[u(k-1)+Tx(k)]+\frac{K_3}{T}[x(k)-x(k-1)]\\\\ &=\left[K_1+K_2T+\frac{K_3}{T}\right]x(k)-\frac{K_3}{T}x(k-1)+K_2u(k-1) \end{aligned} u(k)=K1x(k)+K2[u(k−1)+Tx(k)]+TK3[x(k)−x(k−1)]=[K1+K2T+TK3]x(k)−TK3x(k−1)+K2u(k−1)