8.4 离散控制系统的数学模型_自动控制原理（下）-QQ阅读女生现言网

书名：自动控制原理（下）
作者名：宋永端主编
本章字数：4986字
更新时间：2022-01-14 17:21:30

8.4 离散控制系统的数学模型

为了研究离散系统的性能，需要建立离散系统的数学模型。与连续系统的数学模型类似，在经典控制理论中，线性离散系统的数学模型主要采用差分方程和脉冲传递函数。本节主要介绍线性离散系统的差分方程及其解法，脉冲传递函数的基本概念，以及开环脉冲传递函数、闭环脉冲传递函数的求法。

8.4.1 差分方程

对于连续控制系统，输入信号r（t）和输出信号c（t）之间的关系是用描述系统运动的微分方程来描述的，微分方程则是由输入信号r（t）和输出信号c（t）及其各阶导数构成的。在离散控制系统中，输入信号和输出信号都是离散信号，因此，只能用输入脉冲序列r（kT）和输出脉冲序列c（kT）及其各阶差分所组成的差分方程来描述系统输入信号及输出信号之间的动态特性。为简便起见，通常都省略掉采样周期T，将r（kT）和c（kT）直接写成r（k）和c（k）。

1.差分的定义

所谓差分，是指采样信号在两个相邻采样时刻的采样值之差。取差分的方式有前向差分和后向差分两种。如果当前时刻k的各阶差分全部依赖于当前时刻k和未来时刻k+1、k+2、…的采样值，则称为前向差分；如果当前时刻k的各阶差分全部依赖于当前时刻k和历史时刻k-1、k-2、…的采样值，则称为后向差分。

设连续信号x（t）经采样后的脉冲序列为x（kT），简记为x（k）。一阶前向差分定义为

Δx（k）=x（k+1）-x（k）

二阶前向差分定义为

Δ²x（k）=Δx（k+1）-Δx（k）=x（k+2）-2x（k+1）+x（k）

n阶前向差分定义为

Δⁿx（k）=Δⁿ^-1x（k+1）-Δⁿ^-1x（k）

同理，一阶后向差分定义为

Δx（k）=x（k）-x（k-1）

二阶后向差分定义为

Δ²x（k）=Δx（k）-Δx（k-1）=x（k）-2x（k-1）+x（k-2）

n阶后向差分定义为

Δⁿx（k）=Δⁿ^-1x（k）-Δⁿ^-1x（k-1）

2.差分方程

线性时不变连续系统的数学模型可用下列微分方程表示

式中，r（t）、c（t）分别表示系统的输入信号和输出信号。对于式（8-34）进行离散化，即可将其化为离散系统的差分方程。

设系统的采样周期为T，当T足够小时，函数r（t）在t=kT处的一阶导数近似为

同理，二阶导数近似为

用同样的方法，可以近似表示在t=kT处r（t）的其他各阶导数以及c（t）的各阶导数。因此，可得到用后向差分方程表示线性时不变离散系统的数学模型的一般表达式为

用前向差分方程表示线性时不变离散系统的数学模型的一般表达式为

式（8-35）和式（8-36）中，a_i（i=0，1，2，…，n）和b_j（j=0，1，2，…，m）均为常系数。式（8-35）和式（8-36）所表示的差分方程称为n阶线性常系数差分方程。

前向差分方程和后向差分方程并无本质区别，前向差分方程多用于描述非零初始条件下的离散系统，后向差分方程多用于描述零初始条件下的离散系统。若不考虑初始条件，就离散系统输入变量和输出变量之间的动态关系而言，两种差分方程形式完全等价。

3.差分方程的求解

线性常系数差分方程的求解通常采用迭代法和z变换法。

（1）迭代法

已知离散系统的差分方程，并且给定输出序列的初值，则可以递推计算出输出序列。

例8-10 已知下列二阶差分方程

c（k）-5c（k-1）+6c（k-2）=r（k）

并且给定输入序列r（k）=1，初始条件为c（0）=0，c（1）=1。试用迭代法求输出序列c（k）（k=0，1，2，…，10）。

解：由给定的差分方程可得递推关系

c（k）=r（k）+5c（k-1）-6c（k-2）

则根据初始条件及递推关系，求得

（2）z变换法

用z变换法求解差分方程，完全类似于用拉普拉斯变换法求解微分方程的方法。如果已知线性时不变差分方程，则具体方法就是对差分方程两端取z变换，并利用z变换的实数位移定理，将差分方程变成以z为变量的代数方程，再根据初始条件和给定输入信号z变换表达式，求得输出信号的z变换表达式，然后取z反变换，即可求得输出序列c（k）。

例8-11 已知下列二阶差分方程

c^∗（t+2T）+3c^∗（t+T）+2c^∗（t）=0

或 c（k+2）+3c（k+1）+2c（k）=0

初始条件为c（0）=0，c（1）=1。试用z变换法求解。

解：根据实数位移定理，对差分方程的每一项进行z变换，得

Z[c（k+2）]=z²C（z）-z²c（0）-zc（1）

Z[3c（k+1）]=3zC（z）-3zc（0）

Z[2c（k）]=2C（z）

将以上各式及初始条件代入给定的差分方程中，得到如下代数方程

（z²+3z+2）C（z）=z

解出

查z变换表，得C（z）的z反变换为

或写成 c（k）=（-1）^k-（-2）^k （k=0，1，2，…）

8.4.2 脉冲传递函数

在离散控制系统中，采用z变换，可以求解线性常系数差分方程，来研究离散控制系统的动态性能。但是，如果把z变换的作用仅仅理解为求解线性常系数差分方程，显然是不够的。z变换更为重要的意义在于导出线性离散系统的脉冲传递函数，这给线性离散系统的分析和综合带来极大的方便。

1.脉冲传递函数的定义

设线性时不变离散系统的差分方程的一般表达式为

a₀c（k）+a₁c（k-1）+a₂c（k-2）+…+a_n_-1c（k-n+1）+a_nc（k-n）=b₀r（k）+b₁r（k-1）+…+b_m_-1r（k-m+1）+b_mr（k-m）

如果当t<0时，输入脉冲序列各采样值r（-T）、r（-2T）、…以及输出脉冲序列各采样值c（-T）、c（-2T）、…均为零，即零初始条件，则在零初始条件下对上式两边取z变换，经整理后可以得到

称G（z）为线性时不变离散系统的脉冲传递函数，或叫作z传递函数。

基于以上的讨论，线性时不变离散系统的脉冲传递函数的定义为：在零初始条件下，系统输出脉冲序列c（k）的z变换C（z）与输入脉冲序列r（k）的z变换R（z）之比，即

脉冲传递函数在离散系统的结构图上的表示如图8-10所示，图中，T为采样周期。然而，对大多数实际系统来说，其输出往往是连续信号c（t），而不是采样信号c^∗（t），如图8-11所示。此时，可以在系统输出端虚设一个理想采样开关，如图中虚线所示，它与输入信号的采样开关同步工作，并具有相同的采样周期。如果系统的实际输出c（t）比较平滑，且采样频率较高，则可用c^∗（t）近似描述c（t）。必须指出，虚设的采样开关是不存在的，它表明了脉冲传递函数只能描述输出连续函数c（t）在采样时刻上的离散值c^∗（t）。

图8-10 离散系统的结构图

图8-11 实际离散系统的结构图

2.脉冲传递函数的物理意义

线性时不变离散系统如图8-11所示，当输入信号为单位脉冲函数δ（t）时，其输出即为系统的单位脉冲响应g（t），或称冲击响应，又称脉冲过渡函数；如果输入信号为δ（t-a），则系统的输出为g（t-a）。现假设输入信号为r（t），经采样后为一脉冲序列，即

式中，r（0），r（T），r（2T），… 对应各采样时刻t=0，T，2T，… 的脉冲强度。根据叠加原理，输出量c（t）为一系列脉冲响应之和，即

由于当t<0时，g（t）=0，所以当n>k时，g（kT-nT）=0。于是，当t=kT时，输出脉冲序列为

根据上式及z变换的定义，可得输出量c（t）的z变换C（z）为

在上式的最后一个等式右端的两项中，第一项对应输入信号r^∗（t）的z变换R（z），第二项对应单位脉冲响应函数g^∗（t）的z变换G（z）。因此，上式可以写成

C（z）=R（z）G（z）

脉冲传递函数为

由式（8-39）可知，脉冲传递函数的物理意义为：脉冲传递函数G（z）是系统脉冲过渡函数g（t）经采样后g^∗（t）的z变换。

3.脉冲传递函数的求法

根据脉冲传递函数的定义或脉冲传递函数的物理意义，可得到求取脉冲传递函数的两种方法：①由差分方程求脉冲传递函数；②由传递函数G（s）求脉冲传递函数G（z）。

例8-12 已知离散系统的差分方程为

c（k+2）-2c（k+1）+c（k）=Tr（k+1）

试求脉冲传递函数G（z）。

解：令c（1）=c（0）=0，r（0）=0，利用实数位移定理，对差分方程两端取z变换，得

（z²-2z+1）C（z）=TzR（z）

则有

例8-13 已知开环离散系统连续部分的传递函数为

试求对应的脉冲传递函数G（z）。

解：将G（s）展开为部分分式

查z变换表，得G（s）的z变换为

8.4.3 离散控制系统的动态结构图

离散系统的结构图与连续系统的绘制方法基本相同，其差别仅在于某些位置增加了采样开关。由于脉冲传递函数的定义和传递函数的定义在形式上完全相同，因此在进行结构图的简化变换时，所遵循的等效原则是一致的，即变换前后信号要完全等效。但由于系统中连续信号和离散信号并存，简化法则不再与连续系统相一致。由于采样开关的数目和位置不同，化简后求出的脉冲传递函数也会截然不同。

1.开环系统的脉冲传递函数

当开环离散系统由几个环节串联组成时，其脉冲传递函数的求法与连续系统情况不完全相同。即使两个开环离散系统的组成环节完全相同，但由于采样开关的数目和位置不同，求出的开环脉冲传递函数也不相同。

（1）串联环节的脉冲传递函数

1）环节间有采样开关隔开的情况。设开环离散系统如图8-12a所示，在两个串联连续环节G₁（s）和G₂（s）之间有采样开关隔开。根据脉冲传递函数的定义，考虑到离散信号d^∗（t）的存在，由图8-12a可得

式（8-40）表明，有采样开关隔开的两个线性连续环节串联时的脉冲传递函数，等于这两个环节各自的脉冲传递函数之积。这一结论，可以推广到类似的n个环节相串联时的情况。

2）环节间无采样开关隔开的情况。设开环离散系统如图8-12b所示，在两个串联连续环节G₁（s）和G₂（s）之间没有采样开关隔开。两个串联连续环节G₁（s）和G₂（s）可以简化为一个连续环节G₁（s）G₂（s），于是开环系统的脉冲传递函数为

式（8-41）表明，没有采样开关隔开的两个线性连续环节串联时的脉冲传递函数，等于这两个环节传递函数相乘后的相应z变换。这一结论也可以推广到类似的n个环节相串联时的情况。

图8-12 环节串联时的开环离散系统

例8-14 设开环离散系统如图8-12a、b所示，其中G₁（s）=1/（s+a），G₂（s）=1/（s+b），试求系统的开环脉冲传递函数G（z）。

解：如图8-12a所示，环节间有采样开关隔开时，

如图8-12b所示，环节间没有采样开关隔开时，

显然，在串联环节之间有无采样开关隔离时，其总的脉冲传递函数是不相同的。但是，不同之处仅表现在其零点不同，极点仍然一样。这也是离散系统特有的现象。

（2）有零阶保持器时的开环脉冲传递函数

设有零阶保持器的开环离散系统如图8-13所示。图中，G_h（s）为零阶保持器的传递函数，G_p（s）为连续部分传递函数，两个串联环节之间无采样开关隔离。由于G_h（s）不是s的有理分式函数，因此不便于直接用求串联环节的脉冲传递函数的方法求开环系统脉冲传递函数。

图8-13 带零阶保持器的开环离散系统

由图8-13，将零阶保持器与连续部分相串联的传递函数写成如下形式

式中，W（s）=G_p（s）/s，并注意到e^-^Ts为延迟一个采样周期的延迟环节，则根据实数位移定理可得

由式（8-42）可得到求取有零阶保持器时开环脉冲传递函数的一个很有用的结论：若W（s）所对应的z变换为W（z），则（1-e^-^Ts）W（s）所对应的z变换为（1-z^-1）W（z）。

（3）连续信号进入连续环节时开环离散系统的输出表达式

设开环离散系统如图8-14所示。当开环离散系统的输入端无采样开关时，连续的输入信号r（t）就直接进入连续环节G₁（s），将求不出开环脉冲传递函数，而只能求得系统的输出表达式C（z）。

由图8-14可得

D（z）=Z[D（s）]=Z[G₁（s）R（s）]=G₁R（z）

则有

图8-14 连续信号进入连续环节时的开环离散系统

2.闭环系统的脉冲传递函数

闭环系统的脉冲传递函数定义为：闭环离散控制系统输出信号的z变换C（z）与输入信号的z变换R（z）之比，即

应当注意，当连续的输入信号直接进入连续环节时，将求不出闭环脉冲传递函数，只能求得系统的输出表达式C（z）。

在离散系统中，由于采样开关在系统中设置的不同，结构形式就不一样，因此，闭环离散系统没有唯一的典型结构图形式，系统的闭环脉冲传递函数就没有一般的计算公式，只能根据系统的实际结构具体地求取。

为了便于求出闭环脉冲传递函数，需要了解采样函数拉普拉斯变换的基本关系式。

假设X（s）和Y（s）表示连续信号x（t）和y（t）的拉普拉斯变换，X^∗（s）和Y^∗（s）表示采样信号x^∗（t）和y^∗（t）的拉普拉斯变换，则

由式（8-44）可知，若采样函数的拉普拉斯变换Y^∗（s）与连续函数的拉普拉斯变换X（s）相乘后再采样，则Y^∗（s）可以从采样符号中提出来。

证明：根据式（8-4），有

而

于是

若令m=k+n，则由上式可得

式（8-46）表明，Y^∗（s）是以采样角频率ω_s为周期的周期函数。将式（8-46）代入式（8-45），可得

求闭环脉冲传递函数的具体方法：选择系统输入变量和输出变量，并取采样开关输入端的变量为中间变量。用s域象函数列写方程组，然后对方程组中的各变量进行采样后取z变换，消去中间变量，得到闭环脉冲传递函数或输出表达式。注意：在列写中间变量的s域象函数方程时，避免出现输出连续函数的象函数C（s），以免无法得到闭环脉冲传递函数或输出表达式。

例8-15 设闭环离散系统如图8-15所示，试求系统的闭环脉冲传递函数。