Z变换实战：5个信号处理中的经典问题解析（附MATLAB代码）
Z变换实战：5个信号处理中的经典问题解析（附MATLAB代码）
在数字信号处理领域，Z变换就像一把瑞士军刀，它能将复杂的差分方程转化为简单的代数方程，让系统分析变得直观可控。无论是设计抗干扰的滤波器，还是评估控制系统的稳定性，Z变换都展现出强大的数学魅力。本文将聚焦工程实践中最常遇到的五大场景，通过MATLAB代码演示如何用Z变换这把利器解决实际问题。
1. 滤波器设计的Z变换实现
设计数字滤波器时，Z变换能将时域差分方程转换为直观的传递函数。以一个二阶低通滤波器为例，其差分方程为：
y[n] = 0.2x[n] + 0.3x[n-1] + 0.5y[n-1] - 0.2y[n-2]通过Z变换可得传递函数：
num = [0.2 0.3 0]; % 分子系数 den = [1 -0.5 0.2]; % 分母系数 zplane(num, den); % 绘制零极点图关键操作步骤：
- 用
tf2zpk转换传递函数形式：[z,p,k] = tf2zpk(num,den); - 通过零极点分布判断滤波器特性：
- 极点靠近单位圆时频响曲线出现峰值
- 零点在z=1处会抑制低频信号
注意：设计时应确保所有极点位于单位圆内以保证系统稳定
2. 系统稳定性判定技巧
Z变换的收敛域直接决定系统稳定性。某控制系统传递函数为： $$ H(z) = \frac{z}{(z-0.8)(z-1.2)} $$
稳定性分析流程：
- 计算极点位置：
poles = roots([1 -2 0.96]); % 解分母多项式 - 判定准则：
- 当所有极点满足|p|<1时系统稳定
- 本例中1.2>1，系统不稳定
改进方案： 通过反馈调节改变极点位置：
new_den = conv(den, [1 -0.5]); % 增加极点0.53. 差分方程快速求解
对于递归式数字系统： $$ y[n] - 0.6y[n-1] = x[n] $$
Z变换解法步骤：
- 转换到Z域：
syms z Y X eq = Y - 0.6/z*Y == X; - 求解传递函数：
H = solve(eq,Y)/X; - 逆变换得时域解：
h = iztrans(H);
对比传统迭代法，Z变换将计算复杂度从O(n)降到O(1)。
4. 信号重构中的逆变换实践
给定Z域表达式： $$ X(z) = \frac{z}{z-0.5} + \frac{z}{z-0.8} $$
三种逆变换实现方式：
| 方法 | MATLAB函数 | 适用场景 |
|---|---|---|
| 部分分式展开 | residue | 有理分式 |
| 幂级数展开 | iztrans | 符号运算 |
| 长除法 | deconv | 获取前N项 |
示例代码：
[num,den] = residue([1 1], [1 -1.3 0.4]); x = filter(num, den, [1 zeros(1,99)]); % 生成前100个样本5. 频响特性可视化分析
利用Z变换计算频率响应：
[H,w] = freqz(num,den); subplot(211); plot(w/pi, 20*log10(abs(H))); % 幅频特性 subplot(212); plot(w/pi, angle(H)); % 相频特性典型系统特性对比：
| 系统类型 | 极点分布特征 | 幅频曲线形状 |
|---|---|---|
| 低通 | 集中在z=1附近 | 高频衰减 |
| 高通 | 集中在z=-1附近 | 低频衰减 |
| 带通 | 共轭极点对 | 钟形曲线 |
在实际项目中，我常先用fvtool快速预览滤波器特性，再精细调整参数。例如设计心电图信号处理器时，通过移动极点位置成功抑制了50Hz工频干扰。