LTI（复杂/离散）系统的零状态响应求解实例分析【离散傅里叶 or Z变换】

对于求解这类离散时间LTI系统零状态响应的题目，主要步骤如下：

1. 建立频率响应函数：将差分方程变换到频域，求解系统的频率响应函数 $$H(e^{j\Omega })$$
2. 输入信号的频域表示：将输入信号变换到频域，得到其离散时间傅里叶变换（DTFT） $$X(e^{j\Omega })$$
3. 求解输出的频域表示：通过 $$Y_{ZS}(e^{j\Omega })=H(e^{j\Omega })\cdot X(e^{j\Omega })$$ 计算输出信号在频域中的表示
4. 部分分式分解：对 $$Y_{ZS}(e^{j\Omega })$$ 进行部分分式分解，以便后续逆变换。
5. 离散时间傅里叶逆变换（IDTFT）：对分解后的结果进行IDTFT，得到系统的零状态响应 $$y_{ZS}[k]$$

→ 题目

已知描述某离散时间LTI系统的差分方程为：
$$y[k]-0.75y[k-1]+0.125y[k-2]=4x[k]+3x[k-1]$$

输入 $$x[k]=(0.75{)}^{k}u[k]$$，求该系统的零状态响应 $$y_{ZS}[k]$$

解答步骤

1. 对差分方程两边进行DTFT：

   差分方程在频域表示为：
   $$(1-0.75e^{-j\Omega }+0.125e^{-j2\Omega })Y_{ZS}(e^{j\Omega })=(4+3e^{-j\Omega })X(e^{j\Omega })$$

2. 根据离散系统频率响应的定义：

   $$H(e^{j\Omega })=\frac{Y_{ZS}(e^{j\Omega })}{X(e^{j\Omega })}=\frac{4+3e^{-j\Omega }}{1-0.75e^{-j\Omega }+0.125e^{-j2\Omega }}$$

3. 求输入的频域表示：

   输入 $$x[k]=(0.75{)}^{k}u[k]$$ 的DTFT为：
   $$X(e^{j\Omega })=\sum _{k=0}^{\infty }(0.75{)}^k{e}^{-j\Omega k}=\frac{1}{1-0.75e^{-j\Omega }}$$

4. 求系统输出的频域表示：

   $$Y_{ZS}(e^{j\Omega })=H(e^{j\Omega })X(e^{j\Omega })$$

   代入之前求得的表达式：
   $$Y_{ZS}(e^{j\Omega })=\frac{4+3e^{-j\Omega }}{1-0.75e^{-j\Omega }+0.125e^{-j2\Omega }}\cdot \frac{1}{1-0.75e^{-j\Omega }}$$

5. 进行部分分式分解：

   $$\frac{4+3e^{-j\Omega }}{(1-0.75e^{-j\Omega }+0.125e^{-j2\Omega })(1-0.75e^{-j\Omega })}=\frac{8}{1-0.25e^{-j\Omega }}-\frac{40}{1-0.5e^{-j\Omega }}+\frac{36}{1-0.75e^{-j\Omega }}$$

6. 进行IDTFT（逆离散时间傅里叶变换）：

   $$y_{ZS}[k]=8(0.25{)}^{k}u[k]-40(0.5{)}^{k}u[k]+36(0.75{)}^{k}u[k]$$

因此，系统的零状态响应为：

$$y_{ZS}[k]=8(0.25{)}^{k}u[k]-40(0.5{)}^{k}u[k]+36(0.75{)}^{k}u[k]$$

使用Z变换来解决

这道题的求解过程包括以下几个步骤：

1. 将差分方程转换为Z变换域：
   将差分方程$$y[k]-0.75y[k-1]+0.125y[k-2]=4x[k]+3x[k-1]$$应用Z变换

2. 求出系统的传递函数：
   通过Z变换，求出系统的传递函数$$H(z)$$

3. 求出零状态响应：
   利用输入信号$$x[k]=(0.75{)}^{k}u[k]$$的Z变换，求解零状态响应$$Y_{zs}[k]$$

具体步骤如下：

1. 差分方程

已知差分方程为：
$$y[k]-0.75y[k-1]+0.125y[k-2]=4x[k]+3x[k-1]$$

输入信号：
$$x[k]=(0.75{)}^{k}u[k]$$

2. 对差分方程进行Z变换

应用Z变换，并记住初始条件为零（零状态响应），可以得到：
$$Y(z)-0.75z^{-1}Y(z)+0.125z^{-2}Y(z)=4X(z)+3z^{-1}X(z)$$

3. 代入输入信号的Z变换

输入信号$$x[k]=(0.75{)}^{k}u[k]$$的Z变换为：
$$X(z)=\frac{z}{z-0.75}$$

4. 求解系统传递函数

将X(z)代入差分方程的Z变换中：
$$Y(z)\left(1-0.75z^{-1}+0.125z^{-2}\right)=X(z)\left(4+3z^{-1}\right)$$

传递函数H(z)为：
$$H(z)=\frac{Y(z)}{X(z)}=\frac{4+3z^{-1}}{1-0.75z^{-1}+0.125z^{-2}}$$

将分母和分子同时乘以$$z^2$$以简化：
$$H(z)=\frac{4z+3}{z^2-0.75z+0.125}$$

5. 求零状态响应Y(z)

用$$X(z)=\frac{z}{z-0.75}$$代入：
$$Y_{zs}(z)=H(z)\cdot X(z)=\frac{4z+3}{z^2-0.75z+0.125}\cdot \frac{z}{z-0.75}$$

简化得到：
$$Y_{zs}(z)=\frac{4z^2+3z}{(z^2-0.75z+0.125)(z-0.75)}$$

6. 部分分式展开

将$$Y_{zs}(z)$$部分分式展开：
$$Y_{zs}(z)=\frac{A}{z-0.25}+\frac{B}{z-0.5}+\frac{C}{z-1}$$

求解A, B, C的值并进行逆Z变换。计算得到：
$$A=1,B=1.5,C=-1$$

7. 逆Z变换求解时间域响应

通过Z逆变换得到零状态响应：
$$y_{zs}[k]=A(0.25{)}^k+B(0.5{)}^k-C(1{)}^k$$

最终的零状态响应为：
$$y_{zs}[k]=(0.25{)}^k+1.5(0.5{)}^k-(1{)}^k$$

至此，我们完整地利用Z变换求解了题目的零状态响应。