第 9–10 章 · 拉普拉斯与 Z 变换:极点零点决定一切
傅里叶变换要求信号绝对可积,$\sin(\omega t)$、$e^{at}$ 这类"不安分"的信号进不了门。拉普拉斯(连续)与 Z 变换(离散)给指数因子 $e^{-\sigma t}/z^{-n}$ 加权后再变换——收敛域(ROC)与极点零点把系统分析推到全能模式。
拉普拉斯变换:带衰减的傅里叶
存在条件变成 ROC($\sigma > $ 某值);系统函数 $H(s) = Y(s)/X(s)$(LCCDE 两边变换即得)。核心语言:
拉氏反变换用部分分式(先求极点再展开);一阶系统 $\frac{1}{s + a}$、二阶系统 $\frac{\omega_n^2}{s^2 + 2\zeta\omega_ns + \omega_n^2}$(欠阻尼振铃、超调)是控制工程的第一语言。
- 极点(分母零点):$h$ 的自然模式 $e^{p_it}$——极点位置定波形(左半平面衰减、右半平面发散、虚轴等幅振荡);
- 零点:幅度谱的谷;
- 稳定性 ⟺ 所有极点在左半平面(因果系统);因果 ⟺ ROC 右于最右极点。
Z 变换:离散的对应
单位圆 $z = e^{j\Omega}$ 上就是离散傅里叶——Z 是拉氏的离散表亲。稳定性 ⟺ 极点在单位圆内(因果);ROC 环形($r_- < |z| < r_+$)。离散滤波器谱型速记:极点靠近单位圆 → 该频率凸峰(谐振/陷波的构造术)、零点在圆上 → 该频率挖零(陷波器 $H(z) = \frac{1 - \cos\Omega_0 z^{-1}}{\dots}$ 类)。
系统结构与实现
- 有理系统函数 ↔ 差分方程 ↔ 抽头结构(直接型、级联型、并联型)——同一 $H(z)$ 的多种实现,数值特性不同(定点敏感性);
- FIR(纯零点,恒稳定、可线性相位——滤波器设计的主流)与 IIR(有极点,效率高但需判稳);
- 双线性变换 $s \leftrightarrow \frac{2}{T}\frac{1 - z^{-1}}{1 + z^{-1}}$ 把连续设计搬到离散(频率轴扭曲需预矫)——模拟滤波器遗产的继承通道。
老师的提醒
- 极点-零点草图读频响的三步:标极零 → 看单位圆(或虚轴)上各点与极零的距离 → 距离近极点处增益升、近零点处增益降——十秒读出滤波器形状。
- ROC 与因果/稳定的配对表要烂熟:右序列、左序列、双边序列的 ROC 各在哪儿。
自检
- 由 LCCDE 求系统函数并判断稳定性。
- 极点在单位圆附近对频响的影响?
- 双线性变换的频率扭曲与预矫公式?
更深一层:极点零点是"系统的性格签名"
拉普拉斯与 Z 变换的统一语言是极点-零点:极点(分母零点)决定自然模式(左半平面衰减、右半平面发散、虚轴等幅振荡)、零点决定幅度谱的谷。稳定性 ⟺ 所有极点在左半平面(连续)或单位圆内(离散)——一个代数检查替代全部稳定性分析。极点-零点图的"读图术":极点靠近单位圆 → 该频率凸峰(谐振/陷波的构造术)、零点在圆上 → 该频率挖零(陷波器)。极零图是滤波器与控制器的性格签名——看极零分布就知道系统脾气。
收敛域(ROC)是拉氏/Z 变换的一部分:同一代数式配不同 ROC 是不同信号(因果/反因果之分)——ROC 不是附录,是变换的一部分。因果 ⟺ ROC 右于最右极点(连续)/最外极点之外(离散)。
推导细读:从 LCCDE 到系统函数的一步代数
线性常系数微分/差分方程两边变换即得系统函数 $H(s) = Y(s)/X(s)$——微分方程变代数方程是变换域方法的全部价值。一阶系统 $\frac{1}{s+a}$(惯性环节)、二阶系统 $\frac{\omega_n^2}{s^2 + 2\zeta\omega_n s + \omega_n^2}$(欠阻尼振铃、超调)是控制工程的第一语言——传递函数的分母阶数与系数,直接读出系统的动态性格。这个"一步代数"的威力在滤波器设计(模拟原型的传递函数 → 双线性变换到数字)与控制综合(期望极点 → 反馈增益)里反复兑现。
在真实世界里:双线性变换的"模拟遗产继承"
双线性变换 $s \leftrightarrow \frac{2}{T}\frac{1-z^{-1}}{1+z^{-1}}$ 把连续设计搬到离散——模拟滤波器的遗产(巴特沃斯/切比雪夫/椭圆)通过它进入数字世界。频率轴的扭曲($\omega = \frac{2}{T}\tan\frac{\Omega}{2}$)需要预矫(先把模拟设计的临界频率按扭曲公式预弯)——这个细节的遗漏是数字滤波器"设计对了但频率错了"的头号事故。双线性变换 = 模拟滤波器设计手册与数字信号处理器之间的桥。
专题深入:拉普拉斯变换与"收敛域的灵魂"完整版
第九章的拉普拉斯要抓住收敛域(ROC)这个灵魂:$X(s) = \int x(t)e^{-st}dt$ 的积分不是总收敛——ROC 是使积分收敛的 s 平面条带,"X(s) 相同、ROC 不同 = 不同信号"($(1/2)^n u[n]$ 与 $-(1/2)^n u[-n-1]$ 共享同一 X(z) 但 ROC 不同——ROC 是变换的第二张身份证)。ROC 的结构定理(有理 X(s) 的 ROC 是平行于虚轴的条带、不含极点、左边信号的 ROC 在最左极点之左、右边信号在其之右、双边信号是中间条带);因果性与稳定性的 s 域判据(因果 ⟺ ROC 是右半平面区域且延伸到无穷;稳定 ⟺ ROC 包含虚轴——两个判据合取:因果稳定 ⟺ 全部极点在左半平面——控制系统稳定性判据的信号论根源);零极点图的读图术(极点位置 ↔ 时域模式:左实极点衰减指数、虚轴极点等幅振荡、右半平面发散、共轭对 ↔ 振荡频率——"极点在撒时域的种子")。
推导细读:Z 变换与"系统函数统一论"完整版
Z 变换的对应结构($X(z) = \sum x[n]z^{-n}$、ROC 圆环、因果稳定 ⟺ 全部极点在单位圆内)与拉氏完全平行——"z = e^{sT} 把两个平面焊在一起"(s 左半平面 ↔ 单位圆内、jω 轴 ↔ 单位圆——数字滤波器设计的映射基础)。系统函数的统一论:$H(s) = Y(s)/X(s)$ = 传输算子 = 零极点 + 增益——三种实现的等价结构(直接型(对极点零点敏感、系数量化误差大)、级联型(二阶节串接——"分组配对极点零点、定点实现的首选")、并联型——系数量化敏感性是数字滤波器结构选择的工程判据);一阶/二阶系统的 s 域参数化(ω0 与 Q 值:Q 大 ⟺ 极点靠近虚轴 ⟺ 谐振尖——滤波器与谐振电路共享参数语言)。频率响应的几何读法($H(j\omega) = $ 零点矢量积/极点矢量积——"ω 扫过时测到各零极点的距离":极点附近 H 大(峰)、零点附近 H 小(陷波)——"图形化 Bode"是设计滤波器的直觉引擎)。
工程案例细节:变换域设计的产业出口工程版
九-十章的工程出口:模拟滤波器的原型设计(巴特沃斯/切比雪夫的极点布局("极点均匀分布在 s 左半平面的圆上"——巴特沃斯的极点几何)+ 去归一化——有源滤波器(Sallen-Key 拓扑)的 RC 网络实现);数字滤波器设计的两条路线(IIR:模拟原型 → 双线性变换(频率预畸变补偿非线性映射)——阶数省、相位非线性;FIR:窗函数/等波纹(Parks-McClellan)——线性相位、阶数多——"相位要求高用 FIR、阶数敏感用 IIR"的选型律);控制系统的先行课(传递函数 + 零极点语言直接对接根轨迹/频域设计——"信号与系统是控制论的前置"的精确含义);开关电容与 Sigma-Delta(z 域分析开关电容滤波器、噪声整形的 z 域推导——ΔΣ 调制器的稳定性 = z 平面极点配置问题);通信与音频(IIR 均衡器、梳状滤波器(极点等间隔分布)做 NTSC 消串色与声学回声抵消)。"零极点语言是模拟与数字的普通话"——两个世界的设计师靠它沟通。
例题
离散系统 $H(z) = \frac{1}{1 - 0.9z^{-1}}$。(1)极点在哪、因果稳定吗?(2)粗画频率响应(极点角 0)。(3)$z = e^{j0}$ 处增益多大?
解答
(1) 极点 $z = 0.9$(单位圆内)⟹ 因果系统稳定 ✓(ROC $|z| > 0.9$ 含单位圆)。
(2) 极点在 $\Omega = 0$ 附近靠单位圆 ⟹ 低频增益高(低通性格);$\Omega = \pi$(最远点)增益低:
直流放大 10 倍、奈奎斯特频率只剩 0.53 倍——一阶低通滤波器,极点离圆越近峰越尖。