Discrete-time systems(离散时间系统)

一、分类

---

1）连续时间系统与离散时间系统
连续时间系统的数学模型用微分方程描述
离散时间系统的数学模型用差分方程描述
2）即时系统与动态系统
即时系统：当前输出只有当前输入有关，与过去的输入无关。
动态系统：当前输出不仅与当前输入有关，还与过去的输入有关。
3）线性系统与非线性系统
线性系统：具有叠加性与均匀性（齐次性）的系统
非线性系统：不满足叠加性或均匀性的系统。
4）时变系统与时不变系统
时变系统：系统参数随时间变化的系统。
时不变系统：系统参数不随时间变化的系统。
5）因果系统与非因果系统
因果系统是指某时刻的输出只取决于此时刻和此时刻之前输入的系统。如果系统现在的输出还取决与未来的输入，则不符合因果关系，是非因果系统，是不实际的系统。
因果系统的充要条件是：单位冲激响应h(n) = 0, n<0 。 即系统的输出与n<0时刻无关。

二、离散时间系统

---

2.1、离散时间线性系统

1、可加性
若 $y_1(n)=T[x_1(n)],y_2(n)=T[x_2(n)]$
则 $y_1(n)+y_2(n)=T[x_1(n)]+T[x_2(n)]=T[x_1(n)+x_2(n)]$
2、比例性
$$a_1{y}_1(n)=a_{1}T[x_1(n)]=T[a_1{x}_1(n)]$$
$$a_2{y}_2(n)=a_{2}T[x_2(n)]=T[a_2{x}_2(n)]$$
$$a_1{y}_1(n)+a_2{y}_2(n)=a_{1}T[x_1(n)]+a_{2}T[x_2(n)]=T[a_1{x}_1(n)]+T[a_2{x}_2(n)]$$
注：
1、线性方程不一定是线性系统，如y(n)=4x(n)+6,它是增量线性系统。
2、系统是线性，那么在全部时间上，零输入等于零输出。

2.2、离散时间移不变系统

1、定义
系统参数是不随时间而变化的，即系统响应与激励加于系统的时刻无关。
若 T[x(n)]=y(n)
则 $T[x(n-n_0)]=y(n-n_0)$
2、若系统有一个移变的增益，如y(n)=nx(n)，则系统一定是移变系统。
3、若系统在时间轴上有压缩和扩展，如y(n)=x(Mn)，则系统一定是移变系统。

三、离散时间线性移不变系统（LSI系统）表示法

---

3.1 时域描述（时域分析）

3.1.1 单位抽样响应表示法

1、单位抽样响应表示法：当系统输入为单位抽样序列$\delta (n)$时，LSI系统输出序列为h(n)，即：
$$h(n)=T[\delta (n)]$$
2、利用卷积和关系描述
LSI系统输入为x(n),输出为y(n),得：
$$y(n)=x(n)\ast h(n)=\sum _{m=-\infty }^{\infty }x(m)h(n-m)=\sum _{m=-\infty }^{\infty }h(m)x(n-m)$$
理解：当输入序列x(n)与h(n)的卷积得到输出响应y(n)，所以h(n)描述了系统应有的特性，即系统得到正确的表示（描述）。
3、递推法求解

3.1.2 常系数线性差分方程

模型如下：
$$\sum _{k=0}^N{a}_{k}y(n-k)=\sum _{m=0}^M{b}_{m}x(n-m)$$
或
$$y(n)=\sum _{m=0}^M{b}_{m}x(n-m)-\sum _{k=1}^N{a}_{k}y(n-k)$$
所谓线性是指各输出y(n-k)项及各输入x(n-m)项都只有一次幂且不存在它们的相乘项，否则就是非线性的。
求解：
方法一：时域经典解法(求其次解与特解)
方法二：迭代法(递推法)
方法三：变换域方法(z变换)

3.2 变换域描述（频域分析）

3.2.1 系统函数H(z)

$$H(z)=ϝ[h(n)]=\sum _{n=-\infty }^{\infty }h(n)z^{-n}$$
此时，在z域输入输出关系为
$$Y(z)=X(z)H(z)$$
当系统的起始状态为零时，将3.1.2节常系数差分方程等式两端取z变换可得：
$$H(z)=\frac{Y(z)}{X(z)}=\frac{{\sum }_{m=0}^M{b}_m{z}^{-m}}{1+{\sum }_{k=1}^N{a}_k{z}^{-k}}=\frac{{\sum }_{m=0}^M{b}_m{z}^{-m}}{{\sum }_{k=0}^N{a}_k{z}^{-k}}$$
除了由各个$a_k$、$b_m$决定特性外，还必须给定收敛域范围，才能唯一的确定一个LSI系统。
通过分子分母多项式因式分解可以改写为：
$$H(z)=G\frac{{\prod }_{r=1}^M(1-z_r{z}^{-1})}{{\prod }_{k=1}^N(1-p_k{z}^{-1})}$$
其中$z_r$是H(z)的零点，$p_k$是H(z)的极点。

3.2.2 频域响应表征

若系统函数在z平面单位圆内收敛，则当$z=e^{jw}$时$H(z){|}_{z=e^{jw}}=H{e}^{jw}$存在，称之为系统的频率响应。
用差分方程表征：
$$H(e^{jw})=\frac{Y(e^{jw})}{X(e^{jw})}=\sum _{n=-\infty }^{\infty }h(n)e^{ejwn}$$
用差分方程各系数$a_k$、$b_m$表征：
$$H(e^{jw})=\frac{Y(e^{jw})}{X(e^{jw})}=\frac{{\sum }_{m=0}^M{b}_m{e}^{-jwm}}{1+{\sum }_{k=1}^N{a}_k{e}^{-jwk}}$$
当起始状态为零时，LSI系统的频率响应时由系统本身的冲激响应h(n)或各系数$a_k$、$b_m$决定的，与输入输出无关。
其中$e^{jwn}$称为LSI系统的特征函数，而把$H(e^{jw})$称为特征值。
$$H(e^{jw})=|H(e^{jw})|e^{jarg[H(e^{jw})]}=Re[H(e^{jw})]+jIm[H(e^{jw})]$$
1、$|H(e^{jw})|$表示幅度特性（幅频特性），是系统对输入序列$e^{jwn}$的增益。
2、$arg[H(e^{jw})]$表示相位特性（相频特性），是系统对输入序列$e^{jwn}$的相角改变量。
3、$Re[H(e^{jw})]$和$Im[H(e^{jw})]$用来表示频率响应$H(e^{jw})$的实部与虚部。可以用来观察零极点的分布。

3.3 IIR系统与FIR系统

1、从系统的单位抽样响应h(n)来看，由于设计的方法不同，故划分为两种系统：
(1) 无限长单位抽样响应（IIR）系统：冲激响应无限长。
(2) 有限长单位抽样响应（FIR）系统：冲激响应有限长。
2、从系统函数H(z)的零点、极点来看。
(1) 全零点（MA）系统（滑动平均系统），在z平面中没有极点（全部极点在z=0处），只有零点。
(2) 全极点（AR）系统（自回归系统），在z平面中没有零点，只有极点。
(3) 零极点系统（ARMA）系统（自回归滑动平均），在z平面中既有零点，也有极点。
3、从结构类型来看，可以由差分方程看出，
(1) 递归型结构，IIR系统
(2) 非递归型结构，FIR系统输出只和输入有关，只有零点，没有极点。
FIR系统用零点极点抵消的办法，则也可以采用含有递归结构的系统。

四、离散时间系统的状态空间

---

4.1 状态定义

状态 系统在任何时间$t\ge t_0$的行为
状态变量 是确定状态系统的最小一组变量，用来完全描述动态系统的行为。即一旦给出$t\ge t_0$的输出量和$t=t_0$的初始状态，根据这些状态变量，就可以完全确定系统的未来状态。
状态向量 如果完全描述一个给定系统的行为需要n个状态变量，n个状态变量都可以看做状态向量的n个分量。
状态空间 设$x_1,x_2,...,x_n$为状态变量，那么由$x_1$轴,$x_2$轴,…,$x_n$轴所组成的n为空间称为状态空间。任何状态都可以状态空间中的一点来表示。

4.2 状态方程的建立

1、在状态空间分析中使用三种类型的变量，它们是输入变量、状态变量和输出变量。
2、对于给定的系统，其状态空间的表达式不上唯一的，但是其状态变量的数量是相同的。
3、动态系统必定包含记忆元件，它在$t\ge t_1$时能够记忆输入量的值。积分器可以作为记忆装置，积分器的输出量可以作为系统的状态变量，能够完全描述系统的动态特性。积分器的数目和状态变量的数目相同。
在系统中包含n个积分器，积分器的n个输出量为状态变量$x_1(t),x_t(t),...,x_n(t)$。
状态变量：$\mathbf{x}(t)=\left[\begin{array}{c}{x}_1(t)\\ x_2(t)\\ .\\ .\\ .\\ x_n(t)\end{array}\right]$ 输入变量：$\mathbf{u}(t)=\left[\begin{array}{c}{u}_1(t)\\ u_2(t)\\ .\\ .\\ .\\ u_n(t)\end{array}\right]$ 输出变量：$\mathbf{y}(t)=\left[\begin{array}{c}{y}_1(t)\\ y_2(t)\\ .\\ .\\ .\\ y_n(t)\end{array}\right]$
得到状态方程：
$\dot{\mathbf{x}}(t)=\mathbf{f}(\mathbf{x},\mathbf{u},\mathbf{t})=\left[\begin{array}{c}{f}_1(x_1,x_2,...,x_n;u_1,u_2,...,u_r;t)\\ f_2(x_1,x_2,...,x_n;u_1,u_2,...,u_r;t)\\ .\\ .\\ .\\ f_n(x_1,x_2,...,x_n;u_1,u_2,...,u_r;t)\end{array}\right]$
得到输出方程：
$\mathbf{y}(t)=\mathbf{g}(\mathbf{x},\mathbf{u},\mathbf{t})=\left[\begin{array}{c}{g}_1(x_1,x_2,...,x_n;u_1,u_2,...,u_r;t)\\ g_2(x_1,x_2,...,x_n;u_1,u_2,...,u_r;t)\\ .\\ .\\ .\\ g_n(x_1,x_2,...,x_n;u_1,u_2,...,u_r;t)\end{array}\right]$
将其线性化：
状态方程：$\dot{\mathbf{x}}(t)=\mathbf{A}(t)\mathbf{x}(t)+\mathbf{B}(t)\mathbf{u}(t)$
输出方程：$\mathbf{y}(t)=\mathbf{C}(t)\mathbf{x}(t)+\mathbf{D}(t)\mathbf{u}(t)$
$\mathbf{A}(t)$为状态矩阵，$\mathbf{B}(t)$为输入矩阵，$\mathbf{C}(t)$为输出矩阵，$\mathbf{D}(t)$为直接输出矩阵。
如果系统为时不变系统，那么$\mathbf{f}和\mathbf{g}$不含时间t，则称系统为定常系统，那么$\mathbf{A},\mathbf{B},\mathbf{C},\mathbf{D}$为常数，系统简化为：
状态方程：$\dot{\mathbf{x}}(t)=\mathbf{A}\mathbf{x}(t)+\mathbf{B}\mathbf{u}(t)$
输出方程：$\mathbf{y}(t)=\mathbf{C}\mathbf{x}(t)+\mathbf{D}\mathbf{u}(t)$

4.3 状态方程的求解

状态方程的求解包括时域和变换域两种方法，对于连续系统则通过拉普拉斯变换，对于离散系统则通过z变换。
对于连续系统，则通过拉普拉斯变换，得到：
状态方程：$s\mathbf{X}(s)-\mathbf{x}(0)=\mathbf{A}\mathbf{X}(s)+\mathbf{B}\mathbf{U}(s)$
输出方程：$\mathbf{Y}(s)=\mathbf{C}\mathbf{X}(s)+\mathbf{D}\mathbf{U}(s)$
对于离散系统，则通过z变换，得到：
状态方程：$z\mathbf{X}(z)-\mathbf{x}(0)=\mathbf{A}\mathbf{X}(z)+\mathbf{B}\mathbf{U}(z)$
输出方程：$\mathbf{Y}(z)=\mathbf{C}\mathbf{X}(z)+\mathbf{D}\mathbf{U}(z)$
推导过程省略，得系统的转移函数（系统函数）：
$$H(z)=\mathbf{C}(z\mathbf{I}-\mathbf{A}{)}^{-1}\mathbf{B}+\mathbf{D}$$

4.4 状态矢量的线性变换

对同一系统而言，一各种形式表示的状态矢量之间存在着线性变换关系。
$\mathbf{\hat{A}} = \mathbf{PAP^{-1}} $
$\hat{\mathbf{B}}=\mathbf{P}\mathbf{B}$
$\hat{\mathbf{C}}=\mathbf{C}{\mathbf{P}}^{-1}$
$\hat{\mathbf{D}}=\mathbf{D}$
$\mathbf{P}$ 为新状态方程的变换矩阵，具体变换过程请查询相关书籍。

五、系统分析

---

5.1 稳定性

5.1.1 稳定条件

任何系统只要满足有界输入产生有界输出（BIBO）条件，就是稳定系统，即
$$|x(n)|\le M<\infty$$
$$|y(n)|\le P<\infty$$
LSI系统稳定的充要条件是其单位抽样相应h(n)决对可和，即
$$\sum _{n=-\infty }^{\infty }h(n)=P<\infty$$

5.1.2 根轨迹

其他篇章介绍

5.1.3 奈奎斯特稳定性判据

其他篇章介绍

5.1 可控制性

控制作用是否必然可使系统在有限时间内，从起始状态指引到所要求状态，这就是可控性问题。
所谓可控性被定义为：当系统用状态方程描述时，给定系统的任意初始状态，可以找到容许的输入量(即控制矢量)，在有限时间之内把系统的所有状态引向状态空间的原点（即零状态），如果可以做到这一点，则称系统实施完全可控制的。如果只有对部分状态变量可以做到这一点，则称系统是不完全可控制的。
可控判别矩阵（可控阵）：
$$\mathbf{M}=[\mathbf{B},\mathbf{A}\mathbf{B},{\mathbf{A}}^2\mathbf{B},...,{\mathbf{A}}^{k-1}\mathbf{B}]$$
其中$\mathbf{A}$为状态矩阵，$\mathbf{B}$为输入矩阵，如果矩阵$\mathbf{M}$满秩（即矩阵的每一行都非零），则系统完全可控。

5.2 可观测性

是否可以做到，通过观测有限时间内的输出量，而识别出系统的起始状态。因为一旦能够根据输出量识别出起始状态，则任一状态也就唯一确定。
系统的可观性定义为：如果系统用状态方程来描述时，在给定控制后，能在有限时间间隔内（$0<t<t_1$）根据系统输出唯一地确定系统的所有起始状态，则称系统完全可观；若只能确定部分起始状态，则称系统不完全可观。
可观性判别矩阵（可观阵）：
$$\mathbf{N}=\left[\begin{array}{c}\mathbf{C}\\ \mathbf{C}\mathbf{A}\\ .\\ .\\ .\\ \mathbf{C}{\mathbf{A}}^{\mathbf{k}-1}\end{array}\right]$$
其中$\mathbf{A}$为状态矩阵，$\mathbf{C}$为输出矩阵。如果矩阵$\mathbf{N}$满秩（即矩阵的每一行都非零），则系统完全可观。

六、参考文献

---

1、《经典版 数字信号处理教程 第四版》，程佩青编著，清华大学出版社。
2、《信号与系统 第三版》上下册，郑君里、应启珩、杨为理编著，高等教育出版社。