该系列博客主要讲述Matlab软件在自动控制方面的应用,如无自动控制理论基础,请先学习自动控制系列博文,该系列博客不再详细讲解自动控制理论知识。
自动控制理论基础相关链接:https://blog.csdn.net/qq_39032096/category_10287468.html?spm=1001.2014.3001.5482
博客参考书籍:《MATLAB/Simulink与控制系统仿真》。
1.离散控制系统基本概念
- 如果控制系统中的所有信号都是时间变量的连续函数,即这些信号在全部时间上都是已知的,则这样的系统称为连续时间系统,亦称连续系统;
- 如果控制系统中有一处或几处信号是一串脉冲或数码,即这些信号仅定义在离散时间上,则这样的系统称为离散时间系统,亦称离散系统;
- 把系统中的离散信号是脉冲序列形式的离散系统,称为采样控制系统或脉冲控制系统;把数字序列形式的离散系统,称为数字控制系统或计算机控制系统;
1.1 采样控制系统和数字控制系统
1.1.1 采样控制系统
-
采样系统是对来自传感器的连续信息在某些规定的时间瞬时上取值;
-
如果在有规律的间隔上,系统取到离散信息,则这种采样称为周期采样;如果信息之间的间隔是时变的,或随机的,则称为非周期采样;
-
在采样系统中不仅有模拟部件,还有脉冲部件;通常,测量元件、执行元件和被控对象是模拟元件,其输入和输出是连续信号,即时间上和幅值上都连续的信号,称为模拟信号;控制器中的脉冲元件,其输入和输出为脉冲序列,即时间上离散而幅值上连续的信号,称为离散模拟信号;
-
为了使模拟信号和离散模拟信号在系统中能相互传递,在连续信号和脉冲序列间采用采样器,在脉冲序列和连续信号间用保持器,实现两种信号的转换;
-
信号采样和复现
- 在采样控制系统中,把连续信号转变为脉冲序列的过程称为采样过程,简称采样,实现采样的装置称为采样器,或称采样开关;用 T T T表示采样周期,单位为 s {\rm s} s; f s = 1 / T f_s=1/T fs=1/T表示采样频率,单位为 s − 1 {\rm s}^{-1} s−1; ω s = 2 π f s = 2 π / T \omega_s=2\pi{f_s}=2\pi/T ωs=2πfs=2π/T表示采样角频率,单位为 r a d / s {\rm rad/s} rad/s;
- 在采样控制系统中,把脉冲序列转变为连续信号的过程称为信号复现过程,实现复现过程的装置称为保持器;
-
采样系统的典型结构
- 如果采样器位于系统闭环回路之外,或系统本身不存在闭合回路,则称为开环采样系统;如果采样器位于系统闭合回路内,则称为闭环采样系统;
- 上图中, S S S为理想采样开关,其采样瞬时的脉冲幅值,等于相应采样瞬时误差信号 e ( t ) e(t) e(t)的幅值,且采样持续时间 τ \tau τ趋于零; G h ( s ) G_h(s) Gh(s)为保持器的传递函数; G 0 ( s ) G_0(s) G0(s)为被控对象的传递函数; H ( s ) H(s) H(s)为测量变送反馈元件的传递函数;
1.1.2 数字控制系统
-
数字控制系统是一种以数字计算机为控制器去控制具有连续工作状态的被控对象的闭环控制系统;
-
数字控制系统包括工作于离散状态下的数字计算机和工作于连续状态下的被控对象两部分;
-
计算机作为系统的控制器,其输入和输出只能是二进制编码的数字信号,即在时间上和幅值上都离散的信号,系统中被控对象和测量元件的输入和输出是连续信号,在计算机控制系统中,应用 A / D {\rm A/D} A/D(模/数)和 D / A {\rm D/A} D/A(数/模)转换器,以实现这种信号的转换;
扫描二维码关注公众号,回复: 14867143 查看本文章 -
计算机控制系统的典型原理图:
-
A/D转换器
- A / D {\rm A/D} A/D转换器是把连续的模拟信号转换为离散数字信号的装置;
- A / D {\rm A/D} A/D转换包括两个过程:一是采样过程,即每隔 T T T秒对连续信号 e ( t ) e(t) e(t)进行一次采样,得到采样后的离散信号 e ∗ ( t ) e^*(t) e∗(t);二是量化过程,在计算机中,任何数值的离散信号必须表示为最小二进制的整数倍,成为数字信号,才能运算;采样信号 e ∗ ( t ) e^*(t) e∗(t)经量化后变成数字信号 e ‾ ∗ ( t ) \overline{e}^*(t) e∗(t)的过程,称为编码过程;
-
D/A转换器
- D / A {\rm D/A} D/A转换器是把离散的数字信号转换为连续模拟信号的装置;
- D / A {\rm D/A} D/A转换包括两个过程:一是解码过程,把离散数字信号转换为离散的模拟信号;二是复现过程,离散的模拟信号无法直接控制连续的被控对象,需要经过保持器将离散模拟信号复现为连续的模拟信号;
-
数字控制系统的典型结构图
1.2 信号的采样和保持
1.2.1 采样过程
-
把连续信号变换为脉冲序列的装置称为采样器,亦称采样开关;
-
假设采样器每个 T T T秒闭合一次,闭合的持续时间为 τ \tau τ;采样器的输入 e ( t ) e(t) e(t)为连续信号,输出 e ∗ ( t ) e^*(t) e∗(t)为宽度等于 τ \tau τ的调幅脉冲序列,在采样瞬时 n T ( n = 0 , 1 , 2 , … , ∞ ) nT(n=0,1,2,\dots,\infty) nT(n=0,1,2,…,∞)时出现;
-
在 t = 0 t=0 t=0时,采样器闭合 τ \tau τ秒,此时 e ∗ ( t ) = e ( t ) e^*(t)=e(t) e∗(t)=e(t), t = τ t=\tau t=τ后,采样器打开,输出 e ∗ ( t ) = 0 e^*(t)=0 e∗(t)=0,以后每隔 T T T秒重复一次此过程,采样过程要丢失采样间隔之间的信息;
-
采样过程如下图所示:
-
一般在分析时,可以认为 τ = 0 \tau=0 τ=0,采样器可以用一个理想采样器代替,采样过程看成是一个幅值调制过程;理想采样器看成一个载波为 δ T ( t ) \delta_{T}(t) δT(t)的幅值调制器;
-
理想采样过程如下图所示:
用数学形式描述上图调制过程,有:
e ∗ ( t ) = e ( t ) δ T ( t ) e^*(t)=e(t)\delta_T(t) e∗(t)=e(t)δT(t)
理想单位脉冲序列 δ T ( t ) \delta_T(t) δT(t)可以表示为:
δ T ( t ) = ∑ n = 0 ∞ δ ( t − n T ) \delta_T(t)=\sum_{n=0}^{\infty}\delta(t-nT) δT(t)=n=0∑∞δ(t−nT)
其中: δ ( t − n T ) \delta(t-nT) δ(t−nT)是出现在时刻 t = n T t=nT t=nT、强度为 1 1 1的单位脉冲;因此,有:
e ∗ ( t ) = e ( t ) ∑ n = 0 ∞ δ ( t − n T ) e^*(t)=e(t)\sum_{n=0}^{\infty}\delta(t-nT) e∗(t)=e(t)n=0∑∞δ(t−nT)
因为 e ( t ) e(t) e(t)数值仅在采样瞬时有意义,因此有:
e ∗ ( t ) = ∑ n = 0 ∞ e ( n T ) δ ( t − n T ) e^*(t)=\sum_{n=0}^{\infty}e(nT)\delta(t-nT) e∗(t)=n=0∑∞e(nT)δ(t−nT)
1.2.2 采样过程的数学描述
-
采样信号的拉氏变换
对采样信号 e ∗ ( t ) e^*(t) e∗(t)进行拉氏变换,有:
E ∗ ( s ) = L [ e ∗ ( t ) ] = L [ ∑ n = 0 ∞ e ( n T ) δ ( t − n T ) ] E^*(s)=L\left[e^*(t)\right]=L\left[\sum_{n=0}^{\infty}e(nT)\delta(t-nT)\right] E∗(s)=L[e∗(t)]=L[n=0∑∞e(nT)δ(t−nT)]
采样信号的拉氏变换为:
E ∗ ( s ) = ∑ n = 0 ∞ e ( n T ) e − n T s E^*(s)=\sum_{n=0}^{\infty}e(nT){\rm e}^{-nTs} E∗(s)=n=0∑∞e(nT)e−nTs -
香农采样定理
在设计离散系统时,香农采样定理是必须严格遵守的一条准则,香农定理指明了从采样信号中不失真地复现原连续信号所必需的理论上的最小采样周期 T T T;
香农采样定理指出:如果采样器的输入信号 e ( t ) e(t) e(t)具有有限带宽,并且直到 ω h \omega_h ωh的频率分量,则使信号 e ( t ) e(t) e(t)圆满地从采样信号 e ∗ ( t ) e^*(t) e∗(t)中恢复过来的采样周期 T T T,满足下列条件:
T ≤ 2 π 2 ω h ,即 ω s ≥ 2 ω h T≤\frac{2\pi}{2\omega_h},即\omega_s≥2\omega_h T≤2ωh2π,即ωs≥2ωh
其中: ω h \omega_h ωh为连续频谱的最大角频率; ω s \omega_s ωs为采样角频率; -
采样周期的选取
控制过程 采样周期 T / s T/{\rm s} T/s 控制过程 采样周期 T / s T/{\rm s} T/s 流量 1 1 1 温度 20 20 20 压力 5 5 5 成分 20 20 20 液面 5 5 5 从频域性能指标看,控制系统的闭环频率响应通常具有低通滤波特性,当随动系统的输入信号的频率高于其闭环幅频特性的谐振频率 ω r \omega_r ωr时,信号通过系统将会很快衰减,因此可认为通过系统的控制信号的最高频率分量为 ω r \omega_r ωr;
随动系统的采样角频率可近似取为:
ω s = 10 ω c 或 T = 1 40 t s \omega_s=10\omega_c或T=\frac{1}{40}t_s ωs=10ωc或T=401ts
1.2.3 信号保持
-
保持器的数学描述
-
保持器是具有外推功能的元件,保持器的外推作用表现为现在时刻的输出信号取决于过去时刻离散信号的外推;
-
采用如下多项式外推公式描述保持器:
e ( n T + Δ t ) = a 0 + a 1 Δ t + a 2 ( Δ t ) 2 + ⋯ + a m ( Δ t ) m e(nT+\Delta{t})=a_0+a_1\Delta{t}+a_2(\Delta{t})^2+\dots+a_m(\Delta{t})^m e(nT+Δt)=a0+a1Δt+a2(Δt)2+⋯+am(Δt)m
其中: Δ t \Delta{t} Δt是以 n T nT nT时刻为原点的坐标; -
上式表示:现在时刻的输出 e ( n T + Δ t ) e(nT+\Delta{t}) e(nT+Δt)值,取决于 Δ t = 0 , − T , − 2 T , … , − m T \Delta{t}=0,-T,-2T,\dots,-mT Δt=0,−T,−2T,…,−mT各过去时刻的离散信号 e ∗ ( n T ) , e ∗ [ ( n − 1 ) T ] , e ∗ [ ( n − 2 ) T ] , … , e ∗ [ ( n − m ) T ] e^*(nT),e^*[(n-1)T],e^*[(n-2)T],\dots,e^*[(n-m)T] e∗(nT),e∗[(n−1)T],e∗[(n−2)T],…,e∗[(n−m)T]的 ( m + 1 ) (m+1) (m+1)个值;
-
外推公式中 ( m + 1 ) (m+1) (m+1)个待定系数 a i ( i = 0 , 1 , 2 , … , m ) a_i(i=0,1,2,\dots,m) ai(i=0,1,2,…,m),唯一由过去各采样时刻 ( m + 1 ) (m+1) (m+1)个离散信号值 e ∗ [ ( n − i ) T ] ( i = 0 , 1 , 2 , … , m ) e^*[(n-i)T](i=0,1,2,\dots,m) e∗[(n−i)T](i=0,1,2,…,m)来确定,因此系数 a i a_i ai有唯一解;这样保持器称为 m m m阶保持器;若 m = 0 m=0 m=0,称为零阶保持器; m = 1 m=1 m=1,称为一阶保持器;
-
-
零阶保持器
零阶保持器的数学表达式为:
e ( n T + Δ t ) = e ( n T ) , 0 ≤ Δ t < T e(nT+\Delta{t})=e(nT),0≤\Delta{t}<T e(nT+Δt)=e(nT),0≤Δt<T
零阶保持器是一种按常值外推的保持器,把前一采样时刻 n T nT nT的采样值 e ( n T ) e(nT) e(nT)一直保持到下一时刻 ( n + 1 ) T (n+1)T (n+1)T到来之前,从而使采样信号 e ∗ ( t ) e^*(t) e∗(t)变成阶梯信号 e h ( t ) e_h(t) eh(t);零阶保持器输出特性如下图所示:
把阶梯信号 e h ( t ) e_h(t) eh(t)的中点连接,可以得到与连续信号 e ( t ) e(t) e(t)形状一致但在时间上落后 T / 2 T/2 T/2的响应 e [ t − ( T / 2 ) ] e[t-(T/2)] e[t−(T/2)];
零阶保持器的传递函数:
G h ( s ) = 1 s − e − T s s = 1 − e − T s s G_h(s)=\frac{1}{s}-\frac{ {\rm e}^{-Ts}}{s}=\frac{1-{\rm e}^{-Ts}}{s} Gh(s)=s1−se−Ts=s1−e−Ts
零阶保持器特性:- 低通特性。由于幅频特性的幅值随频率值增大而迅速衰减,即零阶保持器基本上是一个低通滤波器;在 ω = ω s / 2 \omega=\omega_s/2 ω=ωs/2时,幅值只有初值的 63.7 % 63.7\% 63.7%,且截止频率不止一个,因此,零阶保持器除允许主要频谱分量通过外,还允许部分高频频谱分量通过,从而造成数字控制系统的输出中存在纹波;
- 相角滞后特性。零阶保持器要产生相角滞后,且随 ω \omega ω增大而加大,在 ω = ω s \omega=\omega_s ω=ωs处,相角滞后可达 − 180 ° -180° −180°,从而使闭环系统的稳定性变差;
- 时间滞后特性。零阶保持器的输出为阶梯信号 e h ( t ) e_h(t) eh(t),其平均响应为: e [ t − ( T / 2 ) ] e[t-(T/2)] e[t−(T/2)],表明其输出比输入在时间上滞后 T / 2 T/2 T/2,相当于给系统增加了一个延迟时间为 T / 2 T/2 T/2的延迟环节,使系统总的相角滞后增大,对系统的稳定性不利;零阶保持器的阶梯输出同时增加了系统输出中的纹波;
- 一阶保持器复现原信号准确度高,但幅频特性普遍较大,允许通过的信号高频分量较多,更容易造成纹波;一阶保持器的相角滞后比零阶保持器大,在 ω = ω s \omega=\omega_s ω=ωs时,可达 − 280 ° -280° −280°,对系统稳定性更加不利;
1.3. z z z变换理论
1.3.1 z z z变换定义
设连续函数 e ( t ) e(t) e(t)是可拉氏变换的,则拉氏变换定义为:
E ( s ) = ∫ 0 ∞ e ( t ) e − s t d t E(s)=\int_0^{\infty}e(t){\rm e}^{-st}{\rm d}t E(s)=∫0∞e(t)e−stdt
采样信号表达式为:
e ∗ ( t ) = ∑ n = 0 ∞ e ( n T ) δ ( t − n T ) e^*(t)=\sum_{n=0}^{\infty}e(nT)\delta(t-nT) e∗(t)=n=0∑∞e(nT)δ(t−nT)
采样信号 e ∗ ( t ) e^*(t) e∗(t)的拉氏变换为:
E ∗ ( s ) = ∑ n = 0 ∞ e ( n T ) e − n s T E^*(s)=\sum_{n=0}^{\infty}e(nT){\rm e}^{-nsT} E∗(s)=n=0∑∞e(nT)e−nsT
令 z = e s T z={\rm e}^{sT} z=esT, T T T为采样周期, z z z是在复数平面上定义的一个复变量,称为 z z z变换算子;
采样信号 e ∗ ( t ) e^*(t) e∗(t)的 z z z变换定义为:
E ( z ) = E ∗ ( s ) ∣ s = 1 T ln z = ∑ n = 0 ∞ e ( n T ) z − n E(z)=\left.E^*(s)\right|_{s=\frac{1}{T}\ln{z}}=\sum_{n=0}^{\infty}e(nT)z^{-n} E(z)=E∗(s)∣s=T1lnz=n=0∑∞e(nT)z−n
1.3.2 z z z变换方法
-
级数求和法
级数求和法直接根据 z z z变换的定义,将 E ( z ) E(z) E(z)写成展开形式:
E ( z ) = e ( 0 ) + e ( T ) z − 1 + e ( 2 T ) z − 2 + ⋯ + e ( n T ) z − n + ⋯ + E(z)=e(0)+e(T)z^{-1}+e(2T)z^{-2}+\dots+e(nT)z^{-n}+\dots+ E(z)=e(0)+e(T)z−1+e(2T)z−2+⋯+e(nT)z−n+⋯+
Example1: 求单位阶跃函数 1 ( t ) 1(t) 1(t)的 z z z变换。解:
依题意可得, e ( n T ) = 1 ( n = 0 , 1 , 2 , … , ∞ ) e(nT)=1(n=0,1,2,\dots,\infty) e(nT)=1(n=0,1,2,…,∞),可得:
E ( z ) = 1 + z − 1 + z − 2 + ⋯ + z − n + ⋯ + E(z)=1+z^{-1}+z^{-2}+\dots+z^{-n}+\dots+ E(z)=1+z−1+z−2+⋯+z−n+⋯+
在上式中,若 ∣ z − 1 ∣ < 1 |z^{-1}|<1 ∣z−1∣<1,则无穷级数是收敛的,利用等比级数求和公式,可得 1 ( t ) 1(t) 1(t)的 z z z变换的闭合形式:
E ( z ) = 1 1 − z − 1 = z z − 1 E(z)=\frac{1}{1-z^{-1}}=\frac{z}{z-1} E(z)=1−z−11=z−1z
Example2: 设 e ( t ) = δ T ( t ) = ∑ n = 0 ∞ δ ( t − n T ) e(t)=\delta_T(t)=\displaystyle\sum_{n=0}^{\infty}\delta(t-nT) e(t)=δT(t)=n=0∑∞δ(t−nT),求理想脉冲序列 δ T ( t ) \delta_T(t) δT(t)的 z z z变换;解:
依题意可得:
e ∗ ( t ) = δ T ( t ) = ∑ n = 0 ∞ δ ( t − n T ) e^*(t)=\delta_T(t)=\sum_{n=0}^{\infty}\delta(t-nT) e∗(t)=δT(t)=n=0∑∞δ(t−nT)
由拉氏变换可知:
E ∗ ( s ) = ∑ n = 0 ∞ e − n s T E^*(s)=\sum_{n=0}^{\infty}{\rm e}^{-nsT} E∗(s)=n=0∑∞e−nsT
因此:
E ( z ) = ∑ n = 0 ∞ z − n = 1 + z − 1 + z − 2 + ⋯ + E(z)=\sum_{n=0}^{\infty}z^{-n}=1+z^{-1}+z^{-2}+\dots+ E(z)=n=0∑∞z−n=1+z−1+z−2+⋯+
将上式写出闭合形式,可得 δ T ( t ) \delta_T(t) δT(t)的 z z z变换为:
E ( z ) = 1 1 − z − 1 = z z − 1 , ∣ z − 1 ∣ < 1 E(z)=\frac{1}{1-z^{-1}}=\frac{z}{z-1},\left|z^{-1}\right|<1 E(z)=1−z−11=z−1z, z−1 <1
由实例可知,相同的 z z z变换 E ( z ) E(z) E(z)对应于相同的采样函数 e ∗ ( t ) e^*(t) e∗(t),但不一定对应于相同的连续函数 e ( t ) e(t) e(t); -
部分分式法
先求出已知连续时间函数 e ( t ) e(t) e(t)的拉氏变换 E ( s ) E(s) E(s),再将有理分式函数 E ( s ) E(s) E(s)展成部分分式之和的形式,使每一部分分式对应简单的时间函数,然后可得各部分对应的 z z z变换,即可求出 E ( s ) E(s) E(s)对应的 z z z变换 E ( z ) E(z) E(z);
Example3: 设 e ( t ) = sin ω t e(t)=\sin\omega{t} e(t)=sinωt,求 E ( z ) E(z) E(z);
解:
对 e ( t ) = sin ω t e(t)=\sin\omega{t} e(t)=sinωt取拉氏变换,可得:
E ( s ) = ω s 2 + ω 2 E(s)=\frac{\omega}{s^2+\omega^2} E(s)=s2+ω2ω
将上式展开为部分分式:
E ( s ) = 1 2 j ( 1 s − j ω − 1 s + j ω ) E(s)=\frac{1}{2{\rm j}}\left(\frac{1}{s-{\rm j}\omega}-\frac{1}{s+{\rm j}\omega}\right) E(s)=2j1(s−jω1−s+jω1)
可得:
E ( z ) = 1 2 j ( z z − e j ω T − z z − e − j ω T ) = 1 2 j [ z ( e j ω T − e − j ω T ) z 2 − z ( e j ω T + e − j ω T ) + 1 ] E(z)=\frac{1}{2{\rm j}}\left(\frac{z}{z-{\rm e}^{ {\rm j}\omega{T}}}-\frac{z}{z-{\rm e}^{-{\rm j}\omega{T}}}\right)=\frac{1}{2{\rm j}}\left[\frac{z({\rm e}^{ {\rm j}\omega{T}}-{\rm e}^{-{\rm j}\omega{T}})}{z^2-z({\rm e}^{ {\rm j}\omega{T}}+{\rm e}^{-{\rm j}\omega{T}})+1}\right] E(z)=2j1(z−ejωTz−z−e−jωTz)=2j1[z2−z(ejωT+e−jωT)+1z(ejωT−e−jωT)]
化简:
E ( z ) = z sin ω T z 2 − 2 z cos ω T + 1 E(z)=\frac{z\sin\omega{T}}{z^2-2z\cos\omega{T}+1} E(z)=z2−2zcosωT+1zsinωT
1.3.3 z z z变换表
序号 | 时间函数 e ( t ) e(t) e(t) | 拉普拉斯变换 E ( s ) E(s) E(s) | z z z变换 E ( z ) E(z) E(z) |
---|---|---|---|
1 1 1 | δ ( t − n T ) \delta(t-nT) δ(t−nT) | e − s n T {\rm e}^{-snT} e−snT | z − n z^{-n} z−n |
2 2 2 | δ ( t ) \delta(t) δ(t) | 1 1 1 | 1 1 1 |
3 3 3 | 1 ( t ) 1(t) 1(t) | 1 s \displaystyle\frac{1}{s} s1 | z z − 1 \displaystyle\frac{z}{z-1} z−1z |
4 4 4 | t t t | 1 s 2 \displaystyle\frac{1}{s^2} s21 | T z ( z − 1 ) 2 \displaystyle\frac{Tz}{(z-1)^2} (z−1)2Tz |
5 5 5 | t 2 2 ! \displaystyle\frac{t^2}{2!} 2!t2 | 1 s 3 \displaystyle\frac{1}{s^3} s31 | T 2 z ( z + 1 ) 2 ( z − 1 ) 3 \displaystyle\frac{T^2z(z+1)}{2(z-1)^3} 2(z−1)3T2z(z+1) |
6 6 6 | a t / T a^{t/T} at/T | 1 s − ( 1 / T ) ln a \displaystyle\frac{1}{s-(1/T)\ln{a}} s−(1/T)lna1 | z z − a \displaystyle\frac{z}{z-a} z−az |
7 7 7 | e − a t {\rm e}^{-at} e−at | 1 s + a \displaystyle\frac{1}{s+a} s+a1 | z z − e − a T \displaystyle\frac{z}{z-{\rm e}^{-aT}} z−e−aTz |
8 8 8 | t e − a t t{\rm e}^{-at} te−at | 1 ( s + a ) 2 \displaystyle\frac{1}{(s+a)^2} (s+a)21 | T z e − a T ( z − e − a T ) 2 \displaystyle\frac{Tz{\rm e}^{-aT}}{(z-{\rm e}^{-aT})^2} (z−e−aT)2Tze−aT |
9 9 9 | 1 2 t 2 e − a t \displaystyle\frac{1}{2}t^2{\rm e}^{-at} 21t2e−at | 1 ( s + a ) 3 \displaystyle\frac{1}{(s+a)^3} (s+a)31 | T 2 z e − a T 2 ( z − e − a T ) 2 + T 2 z e − 2 a T ( z − e − a T ) 3 \displaystyle\frac{T^2z{\rm e}^{-aT}}{2(z-{\rm e}^{-aT})^2}+\frac{T^2z{\rm e}^{-2aT}}{(z-{\rm e}^{-aT})^3} 2(z−e−aT)2T2ze−aT+(z−e−aT)3T2ze−2aT |
10 10 10 | 1 − e − a t 1-{\rm e}^{-at} 1−e−at | a s ( s + a ) \displaystyle\frac{a}{s(s+a)} s(s+a)a | ( 1 − e − a T ) z ( z − 1 ) ( z − e − a T ) \displaystyle\frac{(1-{\rm e}^{-aT})z}{(z-1)(z-{\rm e}^{-aT})} (z−1)(z−e−aT)(1−e−aT)z |
11 11 11 | sin ω t \sin\omega{t} sinωt | ω s 2 + ω 2 \displaystyle\frac{\omega}{s^2+\omega^2} s2+ω2ω | z sin ω T z 2 − 2 z cos ω T + 1 \displaystyle\frac{z\sin\omega{T}}{z^2-2z\cos\omega{T}+1} z2−2zcosωT+1zsinωT |
12 12 12 | cos ω t \cos\omega{t} cosωt | s s 2 + ω 2 \displaystyle\frac{s}{s^2+\omega^2} s2+ω2s | z ( z − cos ω T ) z 2 − 2 z cos ω T + 1 \displaystyle\frac{z(z-\cos\omega{T})}{z^2-2z\cos\omega{T}+1} z2−2zcosωT+1z(z−cosωT) |
13 13 13 | 1 − cos ω t 1-\cos\omega{t} 1−cosωt | ω 2 s ( s 2 + ω 2 ) \displaystyle\frac{\omega^2}{s(s^2+\omega^2)} s(s2+ω2)ω2 | z z − 1 − z ( z − cos ω T ) z 2 − 2 z cos ω T + 1 \displaystyle\frac{z}{z-1}-\frac{z(z-\cos\omega{T})}{z^2-2z\cos\omega{T}+1} z−1z−z2−2zcosωT+1z(z−cosωT) |
14 14 14 | e − a t sin ω t {\rm e}^{-at}\sin\omega{t} e−atsinωt | ω ( s + a ) 2 + ω 2 \displaystyle\frac{\omega}{(s+a)^2+\omega^2} (s+a)2+ω2ω | z e − a t sin ω T z 2 − 2 z e − a t cos ω T + e − 2 a T \displaystyle\frac{z{\rm e}^{-at}\sin\omega{T}}{z^2-2z{\rm e}^{-at}\cos\omega{T}+{\rm e}^{-2aT}} z2−2ze−atcosωT+e−2aTze−atsinωT |
15 15 15 | e − a t cos ω t {\rm e}^{-at}\cos\omega{t} e−atcosωt | s + a ( s + a ) 2 + ω 2 \displaystyle\frac{s+a}{(s+a)^2+\omega^2} (s+a)2+ω2s+a | z 2 − z e − a T cos ω T z 2 − 2 z e − a T cos ω T + e − 2 a T \displaystyle\frac{z^2-z{\rm e}^{-aT}\cos\omega{T}}{z^2-2z{\rm e}^{-aT}\cos\omega{T}+{\rm e}^{-2aT}} z2−2ze−aTcosωT+e−2aTz2−ze−aTcosωT |
1.3.4 z z z变换性质
-
线性定理
若 E 1 ( z ) = Z [ e 1 ( t ) ] , E 2 ( z ) = Z [ e 2 ( t ) ] , a E_1(z)=Z[e_1(t)],E_2(z)=Z[e_2(t)],a E1(z)=Z[e1(t)],E2(z)=Z[e2(t)],a为常数,则:
Z [ e 1 ( t ) ± e 2 ( t ) ] = E 1 ( z ) ± E 2 ( z ) ; Z [ a e ( t ) ] = a E ( z ) ,式中 : E ( z ) = Z [ e ( t ) ] Z[e_1(t)±e_2(t)]=E_1(z)±E_2(z);Z[ae(t)]=aE(z),式中:E(z)=Z[e(t)] Z[e1(t)±e2(t)]=E1(z)±E2(z);Z[ae(t)]=aE(z),式中:E(z)=Z[e(t)] -
实数位移定理(平移定理)
实数位移含义:整个采样序列在时间轴上左右平移若干个采样周期,其中向左平移为超前,向右平移为滞后;如果函数 e ( t ) e(t) e(t)是可拉普拉斯变换的,其 z z z变换为 E ( z ) E(z) E(z),则有:
Z [ e ( t − k T ) ] = z − k E ( z ) ; Z [ e ( t + k T ) ] = z k [ E ( z ) − ∑ n = 1 k − 1 e ( n T ) z − n ] Z[e(t-kT)]=z^{-k}E(z);Z[e(t+kT)]=z^k[E(z)-\sum_{n=1}^{k-1}e(nT)z^{-n}] Z[e(t−kT)]=z−kE(z);Z[e(t+kT)]=zk[E(z)−n=1∑k−1e(nT)z−n] -
复数位移定理
如果函数 e ( t ) e(t) e(t)是可拉普拉斯变换的,其 z z z变换为 E ( z ) E(z) E(z),则有:
Z [ e ∓ a t e ( t ) ] = E ( z e ± a T ) Z[{\rm e}^{\mp{at}}e(t)]=E(z{\rm e}^{\pm{aT}}) Z[e∓ate(t)]=E(ze±aT) -
终值定理
如果函数 e ( t ) e(t) e(t)的 z z z变换为 E ( z ) E(z) E(z),函数序列 e ( n T ) e(nT) e(nT)为有限值 ( n = 0 , 1 , 2 , … ) (n=0,1,2,\dots) (n=0,1,2,…),且极限 lim n → ∞ e ( n T ) \displaystyle\lim_{n\rightarrow\infty}e(nT) n→∞lime(nT)存在,则函数序列的终值:
lim n → ∞ e ( n T ) = lim z → 1 ( z − 1 ) E ( z ) \lim_{n\rightarrow\infty}e(nT)=\lim_{z\rightarrow1}(z-1)E(z) n→∞lime(nT)=z→1lim(z−1)E(z) -
卷积定理
设 x ( n T ) x(nT) x(nT)和 y ( n T ) y(nT) y(nT)为两个采样函数,其离散卷积积分定义为:
x ( n T ) ∗ y ( n T ) = ∑ k = 0 ∞ x ( k T ) y [ ( n − k ) T ] x(nT)*y(nT)=\sum_{k=0}^{\infty}x(kT)y[(n-k)T] x(nT)∗y(nT)=k=0∑∞x(kT)y[(n−k)T]
则卷积定理如下:若
g ( n T ) = x ( n T ) ∗ y ( n T ) g(nT)=x(nT)*y(nT) g(nT)=x(nT)∗y(nT)
则有:
G ( z ) = X ( z ) ⋅ Y ( z ) G(z)=X(z)·Y(z) G(z)=X(z)⋅Y(z)
其中:
X ( z ) = ∑ k = 0 ∞ x ( k T ) z − k , Y ( z ) = ∑ n = 0 ∞ y ( n T ) z − n , G ( z ) = Z [ g ( n T ) ] = Z [ x ( n T ) ∗ y ( n T ) ] X(z)=\sum_{k=0}^{\infty}x(kT)z^{-k},Y(z)=\sum_{n=0}^{\infty}y(nT)z^{-n},G(z)=Z[g(nT)]=Z[x(nT)*y(nT)] X(z)=k=0∑∞x(kT)z−k,Y(z)=n=0∑∞y(nT)z−n,G(z)=Z[g(nT)]=Z[x(nT)∗y(nT)]
1.3.5 z z z反变换
z z z反变换:已知 z z z变换表达式 E ( z ) E(z) E(z)求相应离散序列 e ( n T ) e(nT) e(nT)的过程,记为:
e ( n T ) = Z − 1 [ E ( z ) ] e(nT)=Z^{-1}[E(z)] e(nT)=Z−1[E(z)]
-
部分分式法(查表法)
设已知的 z z z变换函数 E ( z ) E(z) E(z)无重极点,求出 E ( z ) E(z) E(z)的极点为 z 1 , z 2 , … , z n z_1,z_2,\dots,z_n z1,z2,…,zn,再将 E ( z ) / z E(z)/z E(z)/z展成:
E ( z ) z = ∑ i = 1 n A i z − z i \frac{E(z)}{z}=\sum_{i=1}^n\frac{A_i}{z-z_i} zE(z)=i=1∑nz−ziAi
其中: A i A_i Ai为 E ( z ) / z E(z)/z E(z)/z在极点 z i z_i zi处的留数;由上式写出 E ( z ) E(z) E(z)的部分分式展开式:
E ( z ) = A i z z − z i E(z)=\frac{A_iz}{z-z_i} E(z)=z−ziAiz
然后逐项查表,得到:
e i ( n T ) = Z − 1 [ A i z z − z i ] , i = 1 , 2 , 3 , … , n e_i(nT)=Z^{-1}\left[\frac{A_iz}{z-z_i}\right],i=1,2,3,\dots,n ei(nT)=Z−1[z−ziAiz],i=1,2,3,…,n
最后写出 E ( z ) E(z) E(z)对应的采样函数:
e ∗ ( t ) = ∑ n = 0 ∞ ∑ i = 1 n e i ( n T ) δ ( t − n T ) e^*(t)=\sum_{n=0}^{\infty}\sum_{i=1}^ne_i(nT)\delta(t-nT) e∗(t)=n=0∑∞i=1∑nei(nT)δ(t−nT) -
幂级数法(综合除法)
z z z变换函数 E ( z ) E(z) E(z)可以表示为:
E ( z ) = b 0 + b 1 z − 1 + b 2 z − 2 + ⋯ + b m z − m 1 + a 1 z − 1 + a 2 z − 2 + ⋯ + a n z − n , m ≤ n E(z)=\frac{b_0+b_1z^{-1}+b_2z^{-2}+\dots+b_mz^{-m}}{1+a_1z^{-1}+a_2z^{-2}+\dots+a_nz^{-n}},m≤n E(z)=1+a1z−1+a2z−2+⋯+anz−nb0+b1z−1+b2z−2+⋯+bmz−m,m≤n
其中: a i ( i = 1 , 2 , … , n ) a_i(i=1,2,\dots,n) ai(i=1,2,…,n)和 b j ( j = 0 , 1 , 2 , … , m ) b_j(j=0,1,2,\dots,m) bj(j=0,1,2,…,m)均为常系数;对上式做综合除法,得到 z − 1 z^{-1} z−1升幂排列的幂级数展开式:
E ( z ) = c 0 + c 1 z − 1 + c 2 z − 2 + ⋯ + c n z − n + ⋯ = ∑ n = 0 ∞ c n z − n E(z)=c_0+c_1z^{-1}+c_2z^{-2}+\dots+c_nz^{-n}+\dots=\sum_{n=0}^{\infty}c_nz^{-n} E(z)=c0+c1z−1+c2z−2+⋯+cnz−n+⋯=n=0∑∞cnz−n
如果得到的无穷幂级数是收敛的,则由 z z z变化定义可知,幂级数展开式中的系数 c n ( n = 0 , 1 , 2 , … , n ) c_n(n=0,1,2,\dots,n) cn(n=0,1,2,…,n)就是采样脉冲序列 e ∗ ( t ) e^*(t) e∗(t)的脉冲强度 e ( n T ) e(nT) e(nT),则 E ( z ) E(z) E(z)对应的采样函数为:
e ∗ ( t ) = ∑ n = 0 ∞ c n δ ( t − n T ) e^*(t)=\sum_{n=0}^{\infty}c_n\delta(t-nT) e∗(t)=n=0∑∞cnδ(t−nT) -
反演积分法(留数法)
E ( z ) E(z) E(z)的幂级数展开式为:
E ( z ) = ∑ n = 0 ∞ e ( n T ) z − n = e ( 0 ) + e ( T ) z − 1 + e ( 2 T ) z − 2 + ⋯ + E(z)=\sum_{n=0}^{\infty}e(nT)z^{-n}=e(0)+e(T)z^{-1}+e(2T)z^{-2}+\dots+ E(z)=n=0∑∞e(nT)z−n=e(0)+e(T)z−1+e(2T)z−2+⋯+
用 z n − 1 z^{n-1} zn−1乘以幂级数展开式两端:
E ( z ) z n − 1 = e ( 0 ) z n − 1 + e ( T ) z n − 2 + ⋯ + e ( n T ) z − 1 + … E(z)z^{n-1}=e(0)z^{n-1}+e(T)z^{n-2}+\dots+e(nT)z^{-1}+\dots E(z)zn−1=e(0)zn−1+e(T)zn−2+⋯+e(nT)z−1+…
根据柯西留数定理,设函数 E ( z ) z n − 1 E(z)z^{n-1} E(z)zn−1除有限极点 z 1 , z 2 , … , z k z_1,z_2,\dots,z_k z1,z2,…,zk外,在域 G G G上是解析的;如果有闭合路径 Γ \Gamma Γ包含了这些极点,则有:
e ( n T ) = 1 2 π j ∮ Γ E ( z ) z n − 1 d z = ∑ i = 1 k R e s [ E ( z ) z n − 1 ] z → z i e(nT)=\frac{1}{2\pi{ {\rm j}}}\oint_{\Gamma}E(z)z^{n-1}{\rm d}z=\sum_{i=1}^k{\rm Res}[E(z)z^{n-1}]_{z\rightarrow{z_i}} e(nT)=2πj1∮ΓE(z)zn−1dz=i=1∑kRes[E(z)zn−1]z→zi
其中: R e s [ E ( z ) z n − 1 ] z → z i {\rm Res}[E(z)z^{n-1}]_{z\rightarrow{z_i}} Res[E(z)zn−1]z→zi表示函数 E ( z ) z n − 1 E(z)z^{n-1} E(z)zn−1在极点 z i z_i zi处的留数;因此, E ( z ) E(z) E(z)对应的采样函数为:
e ∗ ( t ) = ∑ n = 0 ∞ e ( n T ) δ ( t − n T ) e^*(t)=\sum_{n=0}^{\infty}e(nT)\delta(t-nT) e∗(t)=n=0∑∞e(nT)δ(t−nT)
关于函数 E ( z ) z n − 1 E(z)z^{n-1} E(z)zn−1在极点处留数计算方法:若 z i ( i = 1 , 2 , … , k ) z_i(i=1,2,\dots,k) zi(i=1,2,…,k)为单极点,则:
R e s [ E ( z ) z n − 1 ] z → z i = lim z → z i [ ( z − z i ) E ( z ) z n − 1 ] {\rm Res}\left[E(z)z^{n-1}\right]_{z\rightarrow{z_i}}=\lim_{z\rightarrow{z_i}}\left[(z-z_i)E(z)z^{n-1}\right] Res[E(z)zn−1]z→zi=z→zilim[(z−zi)E(z)zn−1]
若 E ( z ) z n − 1 E(z)z^{n-1} E(z)zn−1有 n n n阶重极点 z i z_i zi,则:
R e s [ E ( z ) z n − 1 ] z → z i = 1 ( n − 1 ) ! lim z → z i d n − 1 [ ( z − z i ) n E ( z ) z n − 1 ] d z n − 1 {\rm Res}\left[E(z)z^{n-1}\right]_{z\rightarrow{z_i}}=\frac{1}{(n-1)!}\lim_{z\rightarrow{z_i}}\frac{ {\rm d}^{n-1}\left[(z-z_i)^nE(z)z^{n-1}\right]}{ {\rm d}z^{n-1}} Res[E(z)zn−1]z→zi=(n−1)!1z→zilimdzn−1dn−1[(z−zi)nE(z)zn−1]
1.4 MATLAB中离散控制系统相关的函数
【描述模型特性相关的函数】
函数功能 | 离散控制系统中的函数 | 连续控制系统中的函数 |
---|---|---|
阻尼系数和固有频率 | d d a m p {\rm ddamp} ddamp | d a m p {\rm damp} damp |
稳态增益 | d d c g a i n {\rm ddcgain} ddcgain | d c g a i n {\rm dcgain} dcgain |
相对于白噪声的协方差响应 | d c o v a r {\rm dcovar} dcovar | c o v a r {\rm covar} covar |
可控性和可观性 | d g r a m {\rm dgram} dgram | g r a m {\rm gram} gram |
【时域响应相关的函数】
函数功能 | 离散控制系统中的函数 | 连续控制系统中的函数 |
---|---|---|
脉冲响应 | d i m p u l s e {\rm dimpulse} dimpulse | i m p u l s e {\rm impulse} impulse |
零输入响应 | d i n i t i a l {\rm dinitial} dinitial | i n i t i a l {\rm initial} initial |
任意输入下的响应 | d l s i m {\rm dlsim} dlsim | l s i m {\rm lsim} lsim |
阶跃响应 | d s t e p {\rm dstep} dstep | s t e p {\rm step} step |
在网格上绘制根轨迹 | z g r i d {\rm zgrid} zgrid | s g r i d {\rm sgrid} sgrid |
【频域响应相关的函数】
函数功能 | 离散控制系统中的函数 | 连续控制系统中的函数 |
---|---|---|
绘制 B o d e {\rm Bode} Bode图 | d b o d e {\rm dbode} dbode | b o d e {\rm bode} bode |
绘制 N i c h o l s {\rm Nichols} Nichols图 | d n i c h o l s {\rm dnichols} dnichols | n i c h o l s {\rm nichols} nichols |
绘制 N y q u i s t {\rm Nyquist} Nyquist图 | d n y q u i s t {\rm dnyquist} dnyquist | n y q u i s t {\rm nyquist} nyquist |