数字信号处理实战：从零极点图到系统特性分析

1. 零极点图数字信号处理的X光片第一次接触零极点图时我完全不明白这些散落在复平面上的小圆圈和叉叉有什么用。直到有次调试音频滤波器当我把一个极点的位置向单位圆外移动了0.1喇叭里立刻传出刺耳的啸叫声——那一刻我才真正理解零极点图就像是系统的X光片能直观揭示系统的内在特性。零极点图的横轴代表实部纵轴代表虚部。其中零点用○表示使系统传递函数值为零的点极点用×表示使系统传递函数值趋向无穷大的点举个例子假设有个系统的传递函数是H(z)(z-0.5)/(z-0.8)那么零点在z0.5处极点在z0.8处用Python可以快速绘制这个零极点图import matplotlib.pyplot as plt from scipy import signal # 定义系统分子(z-0.5)分母(z-0.8) sys signal.TransferFunction([1, -0.5], [1, -0.8], dt1) plt.figure() signal.zplane(sys.zeros, sys.poles) plt.title(零极点图示例) plt.show()2. 系统稳定性别让极点越狱在工程实践中判断系统稳定性是首要任务。有次我设计了一个IIR滤波器仿真一切正常但实际运行时输出却不断增大最终溢出——后来发现是一个极点的模略大于1。稳定性判据很简单当且仅当所有极点都位于单位圆内即模小于1时系统才是稳定的。这就像给极点画了个警戒线单位圆内极点贡献衰减的分量单位圆上极点产生等幅振荡单位圆外极点导致输出发散通过零极点图可以直观判断画出单位圆|z|1检查所有极点是否都在圆内例如下面这个系统就存在稳定性风险H(z) 1 / (z - 1.2)因为极点z1.2在单位圆外。用向量法分析更直观当z在单位圆上旋转时极点向量长度会周期性变化但永远不会为零。3. 频率响应零极点的引力场零极点分布决定了系统的频率响应特性这就像不同天体引力场影响周围空间。我曾用这个原理快速设计了一个带阻滤波器效果比传统方法更好。幅频响应的计算公式|H(e^jw)| K * ∏(零点到e^jw的距离) / ∏(极点到e^jw的距离)具体分析步骤在单位圆上取测试点e^jw测量到各零极点的距离按公式计算增益当测试点接近零点时分子项趋近零形成陷波接近极点时分母趋近零产生峰值。例如# 带通滤波器设计示例 zeros [1, -1] # 单位圆上零点阻隔低频和高频 poles [0.9*exp(1j*pi/4), 0.9*exp(-1j*pi/4)] # 增强中频实测发现极点越靠近单位圆峰值越尖锐零点越靠近单位圆陷波越深。4. 特殊系统全通与最小相位在实际项目中我遇到过需要校正相位失真的情况这时全通系统就派上用场了。而有次设计均衡器时最小相位系统帮我避免了不必要的相位畸变。4.1 全通系统特点是幅频响应恒为1只改变相位。其零极点呈共轭倒数对称H(z) (z^-1 - a*)/(1 - a z^-1), |a|1MATLAB实现示例[b,a] tfdata(allpass(0.8),v); zplane(b,a)4.2 最小相位系统所有零极点都在单位圆内具有最小群延时最小相位滞后唯一性给定幅频响应时设计技巧将单位圆外零点用其倒数替换例如把零点z1.2改为z0.833。5. 实战案例噪声抑制滤波器设计去年处理ECG信号时我设计了一个基于零极点分析的50Hz工频陷波器。关键步骤零点配置在50Hz对应频率处ω2π*50/fs放置零点theta 2 * np.pi * 50 / 1000 # 假设fs1kHz zeros [np.exp(1j*theta), np.exp(-1j*theta)]极点配置在相同角度但半径0.95处放置极点确保窄带抑制poles [0.95 * np.exp(1j*theta), 0.95 * np.exp(-1j*theta)]稳定性验证检查所有极点模是否小于1print(np.abs(poles)) # 应输出[0.95, 0.95]实测效果显示该滤波器在保持信号形态的同时有效去除了工频干扰。通过调整极点半径还能控制阻带宽度——半径越大阻带越窄。