19、信号与系统：z变换及相关知识详解

信号与系统：z变换及相关知识详解

1. 复习问题

1.1 判断题与选择题

以下是一系列关于z变换的复习问题及答案：
| 题号 | 问题 | 选项 | 答案 |
| — | — | — | — |
| 7.1 | z变换可被视为伪装的拉普拉斯变换。 | (a) 正确，(b) 错误 | a |
| 7.2 | 图7.9所示序列的z变换是？ | (a) (2z + z^2)，(b) (2z^{-1} + z^{-1})，(c) (2z^{-1} + z^{-2})，(d) (z^{-1} + 2z^{-2}) | c |
| 7.3 | (\delta[n - 1])的z变换是？ | (a) 1，(b) z，(c) (z^{-1})，(d) (e^{-z})，(e) (e^{z}) | c |
| 7.4 | (u[n + 1])的z变换是？ | (a) (\frac{1}{z(1 - z)})，(b) (\frac{1}{1 - z})，(c) (\frac{z}{z - 1})，(d) (\frac{z^2}{z - 1}) | d |
| 7.5 | 设(X(z) = \frac{z + 1}{z^2 + 3z + 2})，初始值(x[0])是？ | (a) 0，(b) 1/2，(c) 1，(d) ∞ | a |
| 7.6 | 参考上一题的(X(z))，终值(x[\infty])是？ | (a) 0，(b) 1/2，(c) 1，(d) ∞ | a |
| 7.7 | 给定差分方程(y[n] + y[n - 1] + 2y[n - 2] = u[n])，(y[-1] = -1)，(y[-2] = -2)，(y[0])是？ | (a) -2，(b) -1，(c) 0，(d) 1，(e) 2 | a |
| 7.8 | 哪个MATLAB命令用于求z逆变换？ | (a) ztrans，(b) invztrans，(c) iztrans，(d) filter，(e) dstep | c |
| 7.9 | 哪个MATLAB命令用于确定系统的稳定性？ | (a) freqz，(b) residue，(c) residuez，(d) roots | d |
| 7.10 | 系统的极点图如图7.10所示，该系统是？ | (a) 稳定的，(b) 不稳定的，(c) 临界稳定的 | a |

1.2 问题分析

这些复习问题涵盖了z变换的基本概念、性质以及在MATLAB中的应用。例如，判断z变换与拉普拉斯变换的关系，求解序列的z变换、初始值和终值，利用差分方程求解系统响应，以及使用MATLAB命令进行z变换和逆变换等。

2. 问题求解

2.1 z变换的定义与收敛域

2.1.1 求序列的z变换

• 问题7.1 ：求图7.11所示序列的z变换。
• 问题7.2 ：确定图7.12所示信号的z变换。

2.1.2 确定求和收敛条件

• 问题7.3 ：
  • (a) 确定(\sum_{n = 1}^{\infty} (\frac{2}{3})^n z^n)收敛的条件。
  • (b) 确定(\sum_{n = 1}^{\infty} n (\frac{1}{5})^n z^{-n})收敛的条件。

2.1.3 求z变换及收敛域

• 问题7.4 ：求(x[n] = 2u[n] + (\frac{1}{2})^n u[n])的z变换及收敛域（ROC）。
• 问题7.5 ：求(x[n] = 2(\frac{2}{3})^n u[n] - 5(\frac{2}{5})^n u[n])的z变换及收敛域。

2.2 z变换的性质

2.2.1 求信号的z变换

• 问题7.6 ：
  • (a) 求(u[n - m])的z变换。
  • (b) 求(na^n u[n])的z变换。
  • (c) 求(a^n \cos(\pi n) u[n])的z变换。
• 问题7.7 ：
  • (a) 求(u[n] - u[n - 1])的z变换。
  • (b) 求(a^{n - 1} u[n - 1])的z变换。
  • (c) 求(n a^n u[n - 1])的z变换。
• 问题7.8 ：
  • (a) 求(x[n] = (\frac{1}{2})^n + (- \frac{1}{3})^n)，(n = 0, 1, 2, \cdots)的z变换。
  • (b) 求(x[n] = \begin{cases} (\frac{1}{4})^n, & n = 0, 1, 2, \cdots \ 0, & n < 0 \end{cases})的z变换。
  • (c) 求(x[n] = \begin{cases} 0.5^n, & 0 \leq n \leq 5 \ 0, & \text{otherwise} \end{cases})的z变换。
• 问题7.9 ：
  • (a) 求(x[n] = n a^n \cos(\Omega_0 n) u[n])的z变换。
  • (b) 求(x[n] = n^2 u[n])的z变换。
  • (c) 求(x[n] = e^n \sin(n) u[n])的z变换。
• 问题7.10 ：
  • (a) 求(x[n] = 3(\frac{2}{5})^n + 4(-\frac{1}{3})^n)的z变换。
  • (b) 求(y[n] = 1 - e^{-an})的z变换。

2.2.2 证明z变换对

• 问题7.11 ：
  • (a) 证明((n + 1) u[n] \leftrightarrow \frac{z^2}{(z - 1)^2})。
  • (b) 证明(n^2 a^n u[n] \leftrightarrow \frac{az(z + a)}{(z - a)^3})。
  • (c) 证明(a^n \sin(\Omega n) u[n] \leftrightarrow \frac{za \sin(\Omega)}{z^2 - 2za \cos(\Omega) + a^2})。

2.2.3 其他性质证明

• 问题7.12 ：证明(Z{n x[n]} = -z \frac{dX(z)}{dz})。
• 问题7.13 ：假设序列(x[n])是实值的，证明(X(z) = X^ (z^ ))。

2.2.4 计算参数

• 问题7.14 ：离散时间信号(x[n] = A \cos(\Omega n))的z变换为(X(z) = \frac{z^2 - 4z}{z^2 - 1.6z + 0.8})，计算(A)和(\Omega)的值。

2.2.5 求信号的z变换

• 问题7.15 ：已知(x[n] = 3^n u[n])，求(y[n] = n x[n] + x[n - 1] u[n - 1] + x[n + 1] u[n + 1])的z变换。

2.2.6 求初始值和终值

• 问题7.16 ：确定某些信号（问题7.24中给出z变换）的初始值。
• 问题7.17 ：求问题7.27中各情况的初始值。
• 问题7.18 ：
  • (a) 对于(X(z) = \frac{z}{(z - \frac{1}{6})(z - \frac{1}{4})(z - \frac{1}{5})})，求(x[n])的初始值和终值。
  • (b) 对于(X(z) = \frac{z^2 + 2}{z^2 + 3z + 2})，求(x[n])的初始值和终值。

2.2.7 计算特定值

• 问题7.19 ：离散时间信号(x[n])的z变换为(X(z) = \frac{z}{(z - 2)(z - 1)})，计算(x[0])，(x[1])和(x[105])。

2.2.8 求相关信号的z变换

• 问题7.20 ：已知离散时间信号的z变换为(X(z) = \frac{z^2 + 2}{z^2 + 3z + 1})，求以下信号的z变换：
  • (a) (y[n] = x[n - 1] u[n - 1])
  • (b) (y[n] = \sin(\frac{\pi n}{4}) x[n])
  • (c) (y[n] = n^n x[n])
  • (d) (y[n] = 2x[n] * x[n])

2.2.9 利用z变换求卷积

• 问题7.21 ：使用z变换确定序列(x[n] = [-1, 1, 1])和(h[n] = [-3, 2, 0, 2, 1])的卷积。

2.2.10 利用卷积性质求z变换

• 问题7.22 ：使用卷积性质求(x[n] = \sin(\frac{\pi}{4} n) u[n] * u[n])的z变换。

2.3 z逆变换

2.3.1 求逆z变换

• 问题7.23 ：
  • (a) 求(X_1(z) = \frac{z}{z^2 + \frac{5}{6}z + \frac{1}{6}})的逆z变换。
  • (b) 求(X_2(z) = \frac{z}{(z - 1)(z - 0.5)})的逆z变换。
  • (c) 求(X_3(z) = \frac{z}{(z - 2)(z - 8)})的逆z变换。
• 问题7.24 ：
  • (a) 求(X(z) = \frac{z}{z^2 - z + 2})的逆z变换。
  • (b) 求(Y(z) = \frac{z^2 + 2}{(z + 1)(z - 1)})的逆z变换。
  • (c) 求(H(z) = \frac{z^2 + 4z + 1}{(z + 1)(z + \frac{1}{2})})的逆z变换。
• 问题7.25 ：求(X(z) = \frac{z^2 + 2}{(z - 1)(z + 2)(z + 3)})的逆z变换。
• 问题7.26 ：已知(x[n])的z变换为(X(z) = \frac{z^2 + 2}{(z + 1)(z - 1)})，
  • (a) 求终值(x[\infty])。
  • (b) 确定(X(z))的逆z变换并验证(a)中的结果。
• 问题7.27 ：
  • (a) 求(X_1(z) = \frac{z}{z^2 - z - 0.75})的逆z变换。
  • (b) 求(X_2(z) = \frac{z}{z^2 + z - 0.8 + 0.64})的逆z变换。
• 问题7.28 ：使用长除法求(X(z) = \frac{z}{z^2 - z + 2})的逆z变换。

2.4 应用

2.4.1 求解系统响应

• 问题7.29 ：离散时间系统由(y[n] + 2y[n - 1] + 0.36y[n - 2] = 0)，(y[0] = 0)，(y[1] = 1)描述，求响应(y[n])。
• 问题7.30 ：系统的差分方程为(y[n] - 6y[n - 1] + 15y[n - 2] = 0)，(y[-2] = 0)，(y[-1] = 1)，求(y[n])。
• 问题7.31 ：使用z变换求解差分方程(y[n] + y[n - 1] = (1.5)^n)，(y[0] = 1)。

2.4.2 求系统的传递函数

• 问题7.32 ：求由(y[n] + 2y[n - 1] + 0.5y[n - 2] = x[n] + x[n - 1] + 2x[n - 2])表示的离散时间系统的传递函数。

2.4.3 求系统的完整响应

• 问题7.33 ：系统由(y[n] - y[n - 1] + 3y[n - 2] = u[n])，(y[-2] = 2)，(y[-1] = 0)，(y[0] = 1)表示，求完整响应(y[n])。

2.4.4 求数字滤波器的阶跃响应

• 问题7.34 ：数字滤波器的单位脉冲响应为(h[n] = (\frac{1}{2})^{n - 1} u[n - 1])，求该滤波器的阶跃响应。

2.4.5 求离散时间线性系统的传递函数

• 问题7.35 ：离散时间线性系统的单位脉冲响应为(h[n] = 2 + 3(0.6)^{-n} - 5(0.8)^{-n})，确定其传递函数。

2.4.6 求系统对单位阶跃输入的响应

• 问题7.36 ：离散时间系统的传递函数为(H(z) = \frac{z + 1}{z^2 - z + \frac{1}{2}})，当输入为单位阶跃函数(u[n])时，求系统响应(y[n])。

2.4.7 求系统输出的z变换

• 问题7.37 ：求图7.13所示系统输出的z变换。

2.4.8 求反馈系统的传递函数

• 问题7.38 ：确定图7.14所示反馈系统的传递函数。

2.4.9 求系统的总体传递函数

• 问题7.39 ：确定图7.15所示系统的总体传递函数，设(h_1[n] = \delta[n] + \delta[n - 1])，(h_2[n] = n u[n])，(h_3[n] = (\frac{1}{3})^n u[n])。

2.4.10 求离散时间线性系统的传递函数和差分方程

• 问题7.40 ：图7.16所示离散时间线性系统，
  • (a) 求系统的传递函数(H(z))。
  • (b) 求描述该系统的差分方程。

2.4.11 求系统对不同输入的响应

• 问题7.41 ：系统的传递函数为(H(z) = \frac{0.6z}{(z + 0.2)(z - 0.8)})，当输入为：
  • (a) (x[n] = u[n])
  • (b) (x[n] = 2^n u[n])时，求系统响应。

2.4.12 求系统的传递函数和单位脉冲响应

• 问题7.42 ：离散时间系统的输入(x[n] = \delta[n] + \delta[n - 1])，输出(y[n] = \delta[n] + \delta[n - 1] + \delta[n - 3])，求(H(z))和(h[n])。

2.4.13 求离散时间系统的脉冲响应和阶跃响应

• 问题7.43 ：离散时间系统的传递函数为(H(z) = \frac{0.8z}{(z - 0.6)(z - 2)})，求其脉冲响应和阶跃响应。

2.5 MATLAB计算

2.5.1 求信号的z变换

• 问题7.44 ：使用MATLAB求以下信号的z变换：
  • (a) (x[n] = (0.8)^n u[n])
  • (b) (x[n] = n(0.6)^n u[n])

2.5.2 求逆z变换

• 问题7.45 ：使用MATLAB求(X(z) = \frac{z}{z - 0.6})的逆z变换。

2.5.3 绘制系统的阶跃响应

• 问题7.46 ：线性离散时间系统的传递函数为(H(z) = \frac{z^2 + 2z + 3}{z^3 + 3z^2 + 2z})，使用MATLAB绘制该系统的阶跃响应。

2.5.4 确定传递函数的极点和零点

• 问题7.47 ：确定传递函数(H(z) = \frac{z^3 + 2z^2 + 4z + 3}{z^3 + 6z^2 + 3z - 10})的极点和零点。

2.5.5 绘制z变换的极点和零点位置

• 问题7.48 ：系统的z变换为(X(z) = \frac{3z^4 + 2z^3 + 4z^2 + 3z + 2}{0.6z^4 + 0.2z^3 + z^2 + 0.5z + 0.5})，绘制(X(z))的极点和零点位置。

2.5.6 确定系统的稳定性

• 问题7.49 ：确定传递函数(H(z) = \frac{3z^4 + 2z^3 + 4z^2 + 3z + 2}{3z^4 + 1.25z^3 + 0.5z^2 - 0.375z - 0.2})表示的系统的稳定性。
• 问题7.50 ：检查传递函数(H(z) = \frac{-3z^5 + 2z^4 + 5z^3 - 4z^2 + 2z - 6}{z^5 + 2z^4 + 5z^3 - 4z^2 + z - 1})描述的系统的稳定性。

2.6 数学公式

2.6.1 二次公式

二次方程(ax^2 + bx + c = 0)的根为(x_{1,2} = \frac{-b \pm \sqrt{b^2 - 4ac}}{2a})。

2.6.2 三角函数恒等式

• (\sin(-x) = -\sin(x))，(\cos(-x) = \cos(x))
• (\sec x = \frac{1}{\cos x})，(\csc x = \frac{1}{\sin x})，(\tan x = \frac{\sin x}{\cos x})，(\cot x = \frac{1}{\tan x})
• (\sin(x \pm 90^{\circ}) = \pm \cos x)，(\cos(x \pm 90^{\circ}) = \mp \sin x)，(\sin(x \pm 180^{\circ}) = -\sin x)，(\cos(x \pm 180^{\circ}) = -\cos x)
• (\sin^2 x + \cos^2 x = 1)
• 正弦定理：(\frac{a}{\sin A} = \frac{b}{\sin B} = \frac{c}{\sin C})
• 余弦定理：(a^2 = b^2 + c^2 - 2bc \cos A)
• 正切定理：(\frac{\tan(\frac{A - B}{2})}{\tan(\frac{A + B}{2})} = \frac{a - b}{a + b})
• (\sin(x \pm y) = \sin x \cos y \pm \cos x \sin y)，(\cos(x \pm y) = \cos x \cos y \mp \sin x \sin y)，(\tan(x \pm y) = \frac{\tan x \pm \tan y}{1 \mp \tan x \tan y})
• (2\sin x \sin y = \cos(x - y) - \cos(x + y))，(2\cos x \cos y = \cos(x - y) + \cos(x + y))，(2\sin x \cos y = \sin(x + y) + \sin(x - y))
• (\sin 2x = 2\sin x \cos x)，(\cos 2x = \cos^2 x - \sin^2 x = 2\cos^2 x - 1 = 1 - 2\sin^2 x)，(\tan 2x = \frac{2\tan x}{1 - \tan^2 x})

2.6.3 双曲函数

• (\sinh x = \frac{e^x - e^{-x}}{2})，(\cosh x = \frac{e^x + e^{-x}}{2})，(\tanh x = \frac{\sinh x}{\cosh x})，(\coth x = \frac{1}{\tanh x})，(\text{csch} x = \frac{1}{\sinh x})，(\text{sech} x = \frac{1}{\cosh x})
• (\sinh(x \pm y) = \sinh x \cosh y \pm \cosh x \sinh y)，(\cosh(x \pm y) = \cosh x \cosh y \pm \sinh x \sinh y)，(\tanh(x \pm y) = \frac{\tanh x \pm \tanh y}{1 \pm \tanh x \tanh y})

2.6.4 导数公式

如果(U = U(x))，(V = V(x))，(a)为常数，则：
- (\frac{d}{dx}(aU) = a \frac{dU}{dx})
- (\frac{d}{dx}(UV) = U \frac{dV}{dx} + V \frac{dU}{dx})
- (\frac{d}{dx}(\frac{U}{V}) = \frac{V \frac{dU}{dx} - U \frac{dV}{dx}}{V^2})
- (\frac{d}{dx}(aU^n) = naU^{n - 1} \frac{dU}{dx})
- (\frac{d}{dx}(a^U) = a^U \ln a \frac{dU}{dx})
- (\frac{d}{dx}(e^U) = e^U \frac{dU}{dx})
- (\frac{d}{dx}(\sin U) = \cos U \frac{dU}{dx})
- (\frac{d}{dx}(\cos U) = -\sin U \frac{dU}{dx})
- (\frac{d}{dx}(\tan U) = \sec^2 U \frac{dU}{dx})

2.6.5 不定积分公式

• (\int adx = ax + C)
• (\int UdV = UV - \int VdU)（分部积分法）
• (\int U^n dU = \frac{U^{n + 1}}{n + 1} + C)（(n \neq -1)）
• (\int \frac{dU}{U} = \ln |U| + C)
• (\int a^U dU = \frac{a^U}{\ln a} + C)（(a > 0)，(a \neq 1)）
• (\int e^{ax} dx = \frac{1}{a} e^{ax} + C)
• (\int xe^{ax} dx = \frac{e^{ax}}{a^2}(ax - 1) + C)
• (\int \ln x dx = x \ln x - x + C)
• (\int \sin ax dx = -\frac{1}{a} \cos ax + C)
• (\int \cos ax dx = \frac{1}{a} \sin ax + C)
• (\int \sin^2 ax dx = \frac{x}{2} - \frac{\sin 2ax}{4a} + C)
• (\int \cos^2 ax dx = \frac{x}{2} + \frac{\sin 2ax}{4a} + C)
• (\int x \sin ax dx = \frac{1}{a^2}(\sin ax - ax \cos ax) + C)
• (\int x \cos ax dx = \frac{1}{a^2}(\cos ax + ax \sin ax) + C)
• (\int x^2 \sin ax dx = \frac{1}{a^3}(2ax \sin ax + (2 - a^2 x^2) \cos ax) + C)
• (\int x^2 \cos ax dx = \frac{1}{a^3}(2ax \cos ax - (2 - a^2 x^2) \sin ax) + C)
• (\int e^{ax} \sin bx dx = \frac{e^{ax}}{a^2 + b^2}(a \sin bx - b \cos bx) + C)
• (\int e^{ax} \cos bx dx = \frac{e^{ax}}{a^2 + b^2}(a \cos bx + b \sin bx) + C)
• (\int \sin ax \sin bx dx = \frac{\sin((a - b)x)}{2(a - b)} - \frac{\sin((a + b)x)}{2(a + b)} + C)（(a \neq \pm b)）
• (\int \sin ax \cos bx dx = -\frac{\cos((a + b)x)}{2(a + b)} - \frac{\cos((a - b)x)}{2(a - b)} + C)（(a \neq \pm b)）
• (\int \cos ax \cos bx dx = \frac{\sin((a + b)x)}{2(a + b)} + \frac{\sin((a - b)x)}{2(a - b)} + C)（(a \neq \pm b)）
• (\int \frac{dx}{a^2 + x^2} = \frac{1}{a} \tan^{-1}(\frac{x}{a}) + C)
• (\int \frac{x dx}{a^2 + x^2} = \frac{1}{2} \ln(a^2 + x^2) + C)
• (\int \frac{dx}{(a^2 + x^2)^2} = \frac{x}{2a^2(a^2 + x^2)} + \frac{1}{2a^3} \tan^{-1}(\frac{x}{a}) + C)

2.6.6 定积分公式

如果(m)和(n)为整数，则：
- (\int_{0}^{2\pi} \sin ax dx = 0)
- (\int_{0}^{2\pi} \cos ax dx = 0)
- (\int_{0}^{\pi} \sin^2 ax dx = \int_{0}^{\pi} \cos^2 ax dx = \frac{\pi}{2})
- (\int_{0}^{\pi} \sin mx \sin nx dx = \int_{0}^{\pi} \cos mx \cos nx dx = 0)（(m \neq n)）
- (\int_{0}^{\pi} \sin mx \cos nx dx = \begin{cases} \frac{\pi}{2}, & m = n \text{ 且 } m \text{ 为偶数} \ 0, & m \neq n \ \frac{\pi}{2}, & m = n \text{ 且 } m \text{ 为奇数} \end{cases})
- (\int_{0}^{2\pi} \sin mx \sin nx dx = \int_{0}^{2\pi} \cos mx \cos nx dx = \begin{cases} 0, & m \neq n \ 2\pi, & m = n \end{cases})
- (\int_{0}^{\infty} \frac{\sin ax}{x} dx = \begin{cases} \frac{\pi}{2}, & a > 0 \ 0, & a = 0 \ -\frac{\pi}{2}, & a < 0 \end{cases})
- (\int_{0}^{\infty} \frac{\sin^2 x}{x^2} dx = \frac{\pi}{2})
- (\int_{0}^{\infty} \frac{\cos bx}{x^2 + a^2} dx = \frac{\pi}{2a} e^{-ab})（(a > 0)，(b > 0)）
- (\int_{0}^{\infty} \frac{\sin bx}{x(x^2 + a^2)} dx = \frac{\pi}{2a^2}(1 - e^{-ab}))（(a > 0)，(b > 0)）
- (\int_{0}^{\infty} \sin cx dx = \int_{0}^{\infty} \cos cx dx = \frac{1}{2})
- (\int_{0}^{\pi} \sin^2 nx dx = \int_{0}^{\pi} \cos^2 nx dx = \frac{\pi}{2})（(n)为整数）

2.6.7 洛必达法则

如果(f(0) = 0 = h(0))，则(\lim_{x \to 0} \frac{f(x)}{h(x)} = \lim_{x \to 0} \frac{f’(x)}{h’(x)})，其中撇号表示求导。

2.6.8 泰勒和麦克劳林级数

• (f(x) = f(a) + (x - a)f’(a) + \frac{(x - a)^2}{2!} f’‘(a) + \cdots)
• (f(x) = f(0) + xf’(0) + \frac{x^2}{2!} f’‘(0) + \cdots)

2.6.9 幂级数

• (e^x = 1 + x + \frac{x^2}{2!} + \frac{x^3}{3!} + \cdots)
• (\sin x = x - \frac{x^3}{3!} + \frac{x^5}{5!} - \cdots)
• (\cos x = 1 - \frac{x^2}{2!} + \frac{x^4}{4!} - \cdots)
• (\tan x = x + \frac{x^3}{3} + \frac{2x^5}{15} + \cdots)
• ((1 + x)^n = 1 + nx + \frac{n(n - 1)}{2!} x^2 + \cdots)（(|x| < 1)）
• (\frac{1}{1 - x} = 1 + x + x^2 + \cdots)（(|x| < 1)）

2.6.10 求和公式

• (\sum_{k = 1}^{N} k = \frac{N(N + 1)}{2})
• (\sum_{k = 1}^{N} k^2 = \frac{N(N + 1)(2N + 1)}{6})
• (\sum_{k = 1}^{N} k^3 = (\frac{N(N + 1)}{2})^2)
• (\sum_{k = 0}^{N} a^k = \frac{1 - a^{N + 1}}{1 - a})（(a \neq 1)）
• (\sum_{k = M}^{N} a^k = \frac{a^M - a^{N + 1}}{1 - a})（(a \neq 1)）
• ((a + b)^N = \sum_{k = 0}^{N} \binom{N}{k} a^{N - k} b^k)，其中(\binom{N}{k} = \frac{N!}{k!(N - k)!})

2.6.11 对数恒等式

• (\log(xy) = \log x + \log y)
• (\log(\frac{x}{y}) = \log x - \log y)
• (\log(x^n) = n \log x)
• (\log_{10} x)为常用对数，(\ln x)为自然对数，(\log_e x = \ln x)

2.6.12 指数恒等式

• (e^x = 1 + x + \frac{x^2}{2!} + \frac{x^3}{3!} + \cdots)（(e = 2.7182)）
• (e^x e^y = e^{x + y})
• ((e^x)^n = e^{nx})
• (e^{\ln x} = x)

2.6.13 近似公式

如果(|x| \ll 1)，则：
- ((1 \pm x)^n \approx 1 \pm nx)
- (e^x \approx 1 + x)
- (\ln(1 + x) \approx x)
- (\sin x \approx x)（(x \to 0)）
- (\cos x \approx 1)
- (\tan x \approx x)

3. 总结

通过上述复习问题和求解过程，我们对z变换的基本概念、性质、逆变换以及在系统分析中的应用有了更深入的理解。同时，掌握了使用MATLAB进行z变换和逆变换的方法，以及如何利用数学公式解决相关问题。在实际应用中，这些知识可以帮助我们分析和设计离散时间系统，处理各种信号处理问题。

以下是一个简单的mermaid流程图，展示求解系统响应的一般步骤：

graph TD;
    A[给定差分方程] --> B[进行z变换];
    B --> C[求解Y(z)];
    C --> D[求逆z变换];
    D --> E[得到系统响应y[n]];

这些内容涵盖了z变换的多个方面，从理论知识到实际应用，再到使用MATLAB进行计算，为我们学习和研究信号与系统提供了全面的指导。在后续的学习和实践中，我们可以进一步深入探讨这些知识，解决更复杂的问题。

4. 深入分析与拓展

4.1 z变换性质的应用拓展

z变换的性质在信号处理和系统分析中具有广泛的应用。例如，线性性质使得我们可以对多个信号的线性组合进行z变换，从而简化计算。在问题7.15中，已知(x[n] = 3^n u[n])，求(y[n] = n x[n] + x[n - 1] u[n - 1] + x[n + 1] u[n + 1])的z变换。我们可以利用z变换的线性性质，分别求出(n x[n])、(x[n - 1] u[n - 1])和(x[n + 1] u[n + 1])的z变换，然后将它们相加。

再如，时域移位性质在处理具有时间延迟或提前的信号时非常有用。在问题7.20(a)中，求(y[n] = x[n - 1] u[n - 1])的z变换，就可以直接应用时域移位性质。

4.2 系统稳定性分析

系统的稳定性是系统分析中的一个重要问题。在复习问题7.10中，通过系统的极点图判断系统的稳定性。一般来说，对于离散时间系统，若其所有极点都位于单位圆内，则系统是稳定的；若有极点位于单位圆外，则系统是不稳定的；若有极点位于单位圆上且为单极点，则系统是临界稳定的。

在MATLAB中，可以使用 roots 命令来确定传递函数的极点，进而判断系统的稳定性。例如，对于传递函数(H(z) = \frac{z^3 + 2z^2 + 4z + 3}{z^3 + 6z^2 + 3z - 10})，可以使用以下代码来确定其极点：

num = [1, 2, 4, 3];
den = [1, 6, 3, -10];
poles = roots(den);
disp(poles);

通过分析这些极点的位置，就可以判断系统的稳定性。

4.3 逆z变换的方法总结

求逆z变换有多种方法，如部分分式展开法、长除法等。在问题7.23 - 7.28中，分别使用了这些方法来求解逆z变换。
- 部分分式展开法 ：将复杂的z变换函数分解为简单的分式之和，然后根据已知的z变换对求出逆z变换。例如，对于(X(z) = \frac{z}{z^2 + \frac{5}{6}z + \frac{1}{6}})，先将其进行部分分式展开，然后再求逆z变换。
- 长除法 ：通过长除法将z变换函数展开为幂级数，从而得到逆z变换的序列。在问题7.28中，使用长除法求(X(z) = \frac{z}{z^2 - z + 2})的逆z变换。

4.4 卷积性质的应用

卷积性质在信号处理中非常重要。在问题7.21中，使用z变换确定序列(x[n] = [-1, 1, 1])和(h[n] = [-3, 2, 0, 2, 1])的卷积。根据卷积性质，两个序列的卷积的z变换等于它们各自z变换的乘积。因此，可以先求出(x[n])和(h[n])的z变换(X(z))和(H(z))，然后计算(Y(z) = X(z)H(z))，最后求(Y(z))的逆z变换得到卷积结果。

5. 案例分析

5.1 离散时间系统的设计与分析

假设我们要设计一个离散时间系统，其输入为(x[n])，输出为(y[n])，满足差分方程(y[n] + 2y[n - 1] + 0.36y[n - 2] = x[n])。我们可以按照以下步骤进行分析和设计：
1. 求系统的传递函数 ：对差分方程进行z变换，得到(Y(z)(1 + 2z^{-1} + 0.36z^{-2}) = X(z))，则传递函数(H(z) = \frac{Y(z)}{X(z)} = \frac{1}{1 + 2z^{-1} + 0.36z^{-2}} = \frac{z^2}{z^2 + 2z + 0.36})。
2. 分析系统的稳定性 ：使用 roots 命令求传递函数的极点：

num = [1, 0, 0];
den = [1, 2, 0.36];
poles = roots(den);
disp(poles);

判断极点是否都位于单位圆内，以确定系统的稳定性。
3. 求系统的单位脉冲响应和阶跃响应 ：
- 单位脉冲响应 ：令输入(x[n] = \delta[n])，则(X(z) = 1)，(Y(z) = H(z))，求(Y(z))的逆z变换得到单位脉冲响应(h[n])。
- 阶跃响应 ：令输入(x[n] = u[n])，则(X(z) = \frac{z}{z - 1})，(Y(z) = H(z)X(z))，求(Y(z))的逆z变换得到阶跃响应。

5.2 信号处理中的应用

在信号处理中，z变换可以用于滤波、信号恢复等。例如，设计一个数字滤波器，其单位脉冲响应为(h[n] = (\frac{1}{2})^{n - 1} u[n - 1])。我们可以通过以下步骤分析该滤波器：
1. 求滤波器的传递函数 ：对(h[n])进行z变换，得到(H(z) = \frac{z^{-1}}{1 - \frac{1}{2}z^{-1}} = \frac{1}{z - \frac{1}{2}})。
2. 分析滤波器的频率响应 ：使用 freqz 命令绘制滤波器的频率响应：

num = [1];
den = [1, -0.5];
freqz(num, den);

通过频率响应图，我们可以了解滤波器对不同频率信号的滤波特性。

6. 总结与展望

6.1 总结

本文详细介绍了z变换的基本概念、性质、逆变换以及在系统分析中的应用。通过复习问题和求解过程，我们掌握了z变换的计算方法，包括求序列的z变换、确定收敛域、利用性质求z变换等。同时，学会了使用MATLAB进行z变换和逆变换的计算，以及判断系统的稳定性。此外，还介绍了z变换在系统设计、信号处理等方面的应用案例。

6.2 展望

z变换作为信号与系统中的重要工具，在未来的研究和应用中具有广阔的前景。随着科技的发展，离散时间系统在通信、控制、图像处理等领域的应用越来越广泛。我们可以进一步研究z变换在这些领域的应用，例如设计更高效的数字滤波器、优化离散时间控制系统等。同时，结合其他信号处理方法，如傅里叶变换、小波变换等，可以更好地解决复杂的信号处理问题。

以下是一个mermaid流程图，展示设计离散时间系统的一般步骤：

graph TD;
    A[确定系统要求] --> B[建立差分方程];
    B --> C[求传递函数];
    C --> D[分析系统稳定性];
    D --> E[求系统响应];
    E --> F[评估系统性能];
    F --> G{是否满足要求};
    G -- 是 --> H[完成设计];
    G -- 否 --> B;

通过不断学习和实践，我们可以更好地掌握z变换的知识，将其应用到实际问题中，为解决各种信号处理和系统分析问题提供有力的支持。