完成后应当能做到
- 掌握相关、能量谱、功率谱与匹配滤波入口
- 能识别 FIR、IIR 和多采样率题的第一步
- 能在时域、频域、s 域和 z 域之间选择最短路径
第 12 章 系统综合题
12.1 综合题常见组合
清华 830 的专业课题可能不是单一公式题,而是把多个模块串起来:
- 先由系统性质判断 LTI、因果、稳定。
- 再由冲激响应写卷积。
- 再转到傅里叶、拉普拉斯或 Z 域。
- 最后根据极点和 ROC 判断系统特性。
12.2 解题总流程
- 判断时间类型:连续还是离散。
- 判断系统类型:是否 LTI。
- 若 LTI,找 、 或方程。
- 输入输出关系优先选最短路径:
- 时域短:用卷积。
- 频域短:用乘法。
- 方程短:用系统函数。
- 最后检查因果、稳定、初值和单位。
12.3 典型例题框架
例 1:由微分方程求频率响应
若:
零状态拉普拉斯变换:
所以:
若系统因果且 ,则:
频率响应:
例 2:由差分方程判断稳定
若:
系统函数:
极点 。若因果,ROC 为 ,包含单位圆,所以稳定。
12.4 最后复盘清单
考前你必须能独立完成:
- 判断连续/离散信号周期。
- 求能量和平均功率。
- 用阶跃和冲激表示分段信号。
- 判断系统线性、时不变、因果、稳定。
- 做矩形、阶跃、指数、冲激卷积。
- 求傅里叶级数系数并画离散谱。
- 熟练使用傅里叶变换性质。
- 解释采样频谱复制和混叠。
- 用调制性质画频谱。
- 根据 、ROC 判断连续系统因果稳定。
- 根据 、ROC 判断离散系统因果稳定。
- 从微分方程和差分方程求系统函数。
24. 典型真题式例题
下面这些例题不是为了刷题量,而是把前面的二级结论串起来。每题都给出“落笔第一步”。
24.1 平移矩形卷积
求:
令:
已知:
又因为第二个信号是 ,卷积结果整体平移:
所以:
考场触发器:卷积中某个输入平移,输出跟着同向平移。
24.2 门控载波频谱
求:
先写矩形频谱:
乘余弦等于左右搬移各半:
考场触发器:看到“窗函数乘载波”,立刻画两个 sinc,中心在 。
24.3 系统函数和 ROC 反变换
已知:
先部分分式:
若 ROC 为:
两项都取右边信号:
若 ROC 为:
两项都取左边信号:
若 ROC 为:
第一项取左边,第二项取右边:
考场触发器:同一个代数式,ROC 不同,时域信号不同。
24.4 采样混叠判断
设 带限:
以采样周期 采样:
采样角频率:
采样后频谱:
相邻两份频谱中心间隔为 ,每份宽度为 。无混叠条件:
若题目给 Hz:
考场触发器:不要只背奈奎斯特,先画复制谱,复制谱不重叠才可理想重建。
24.5 DFT 做线性卷积
若 长度为 , 长度为 ,线性卷积长度为:
用 点 DFT 做卷积时,实际得到的是 点圆周卷积。要让圆周卷积等于线性卷积,必须补零到:
证明思路:线性卷积的非零结果只占 到 。若 ,圆周卷积会把 的尾部折回前面相加,产生时域混叠;若 ,没有尾部折回,所以两者相同。
考场触发器:看到“用 DFT 快速卷积”,第一步先算 。
25. 更多真题式例题:调制、响应与初值
这一节继续补考场常见题型。重点不是题目本身,而是每类题的“第一步写什么”。
25.1 DSB-SC 同步解调
已知抑制载波双边带信号:
接收端乘本振 :
利用:
得到:
若 是低通信号,第二项搬到 附近,低通滤波后:
考场触发器:同步解调本振若写成 ,低通后刚好恢复原幅度。
25.2 同步解调有相位误差
若本振有相位误差:
则:
用积化和差:
得到:
低通后:
二级结论:相干解调相位误差会造成幅度衰减;当 时,低通输出为零。
25.3 普通 AM 的包络不失真条件
普通调幅:
高频载波的包络为:
若希望包络检波后直接得到 的形状,必须保证括号内不变号:
等价于:
若定义调幅指数:
则无过调制条件是:
二级结论:普通 AM 能用包络检波的前提是“不翻包络”,也就是 始终为正。
25.4 零输入响应与零状态响应
考虑一阶连续系统:
初始条件:
输入:
用单边 Laplace:
代入微分方程:
所以:
第一项只由初始条件产生,叫零输入响应:
第二项只由输入产生,叫零状态响应:
总响应:
考场触发器:带初始条件的题,用单边变换;系统函数 只描述零状态响应,不包含初始条件。
25.5 用系统函数求零状态响应
还是系统:
零初始条件下:
所以系统函数:
若输入为单位阶跃:
零状态输出:
部分分式:
所以:
二级结论:看到“系统函数 + 输入”,默认求的是零状态响应。
25.6 单边 Z 解一阶差分方程
已知:
初始条件:
输入:
单边 Z 变换中:
代入:
整理:
所以:
第一项是零状态响应,第二项是零输入响应。
当 :
而:
因此:
考场触发器:单边 Z 的移位会带出初始项,别把它当成双边 Z 的纯 。
25.7 单边 Z 初值定理快速验算
单边 Z 定义:
令 ,则 ,所有 的项消失:
以上一题结果为例:
取极限:
也可由差分方程在 直接验证:
所以:
二级结论:做完单边 Z 题,用初值定理验算 ,很容易抓出移位漏项。
25.8 零输入响应的特征根写法
若离散系统齐次方程为:
设:
代入:
两边除以 :
这就是特征方程。若两个根为 且不相等,则:
若重根为 ,则:
常数由初始条件决定。
二级结论:零输入响应只看系统自然模态,也就是特征根;输入只影响零状态响应和总响应中的强迫部分。
29.4 相关函数、能量谱和功率谱
相关不是卷积的替代品,它描述的是两个信号在不同平移下的相似程度。
连续互相关:
连续自相关:
离散互相关:
常见性质:
- ,即零延时自相关等于能量。
- 。
- 自相关在 处通常取最大值。
能量信号的能量谱密度:
维纳-辛钦思想:
功率信号对应功率谱密度,周期信号的功率谱是离散谱线。
考场二级结论:相关题先写“滑动内积”,谱密度题先看是能量信号还是功率信号,不要把能量谱和功率谱混用。
继续加深:相关和卷积最容易混的地方在“翻转对象”和“物理含义”。
卷积是系统响应:
它的关键词是“输入经过系统”。相关是相似度:
它的关键词是“平移后有多像”。如果题目出现“检测延时、同步、相似程度、最大相关峰、白噪声中检测已知信号”,优先想到相关和匹配滤波。
能量谱和功率谱的区别:
| 信号类型 | 看什么量 | 频域名称 | 常见对象 |
|---|---|---|---|
| 能量信号 | 总能量有限 | 能量谱密度 $ | X(j\omega) |
| 功率信号 | 平均功率有限 | 功率谱密度 | 周期信号、平稳随机信号 |
小例题模板:若 是实偶信号,则 也是实偶函数,且 。求自相关时,不要一上来做傅里叶变换;如果波形简单,先用“重叠长度或重叠面积”算滑动内积,通常更快。
29.5 匹配滤波器的最小推导模板
匹配滤波器用于在白噪声背景下最大化某个采样时刻的输出信噪比。若已知待检测信号为 ,希望在 处采样,则匹配滤波器冲激响应为:
频域形式:
输出:
在 处,输出等于信号和自身的相关峰,因此最容易和噪声区分。
考场落笔:
- 先确定采样时刻 。
- 写 。
- 若题目在频域给 ,写 。
- 检查是否需要归一化增益 。
二级结论:匹配滤波器不是普通低通滤波器,它的核心是“时间反转 + 共轭 + 延时”。
继续加深:匹配滤波器的“匹配”不是匹配频带,而是匹配波形。若输入为
其中 是白噪声,则滤波输出在 采样:
代入 后,信号项变成 与自身的内积,输出峰值正比于信号能量。白噪声通过滤波器后的噪声功率与滤波器能量有关,柯西不等式说明当滤波器形状取 时,采样信噪比最大。
考场不要展开完整证明,只要会写三句话:
- 目标是在 处最大化输出信噪比。
- 最优冲激响应是待检测信号的共轭时间反转并延时。
- 输出峰值等于相关峰,峰值位置给出同步时刻。
易错点:如果信号是实信号,共轭可以省略;如果题目给的是复基带信号,共轭不能省略。
29.8 IIR 数字滤波器设计:两条主路
IIR 数字滤波器常从模拟原型滤波器变换而来。
第一条:脉冲响应不变法。
核心思想:
优点:时域冲激响应对应清楚。
缺点:模拟频谱会按采样频率周期复制,可能混叠,所以更适合低通或带限较好的原型。
第二条:双线性变换法。
映射:
优点:左半平面映到单位圆内,不产生频谱混叠。
缺点:频率轴非线性映射,有频率畸变:
预畸变:
二级结论:脉冲响应不变法怕混叠;双线性变换法不混叠但要预畸变。
继续加深:IIR 设计题最重要的是先判断“模拟指标怎么变成数字指标”。若题目给数字角频率 ,用双线性变换法时必须先预畸变:
然后按模拟滤波器设计流程求 ,最后用
替换得到 。
两种方法对比:
| 方法 | 映射对象 | 优点 | 风险 |
|---|---|---|---|
| 脉冲响应不变法 | 时域直观 | 频谱混叠 | |
| 双线性变换法 | 稳定性保持,不混叠 | 频率畸变 |
易错点:脉冲响应不变法不是把 直接替换成某个 表达式;双线性变换法不是线性频率映射,必须记住 。
29.9 FIR 数字滤波器设计:线性相位和窗函数
FIR 系统:
若冲激响应满足对称:
或反对称:
则 FIR 可能具有广义线性相位。
窗函数法步骤:
- 根据理想频响写理想冲激响应 。
- 由于 无限长,用窗函数截断:
- 根据通带/阻带指标选择窗类型和长度。
常见取舍:
- 矩形窗主瓣窄、旁瓣高。
- 汉宁/汉明窗旁瓣低一些,但主瓣更宽。
- 窗越长,过渡带越窄。
频率采样法则是直接在频域取样,再通过 IDFT 得到 。
二级结论:FIR 设计题先看是否要求线性相位;若要求线性相位,第一反应就是 的对称性。
继续加深:FIR 题里的“阶数、长度、对称性、群时延”要统一起来。若 FIR 长度为 ,阶数为 。线性相位系统的相位通常含有
所以群时延为
四类线性相位 FIR:
| 类型 | 长度 | 对称性 | 常见限制 |
|---|---|---|---|
| I 型 | 奇数 | 偶对称 | 最通用 |
| II 型 | 偶数 | 偶对称 | 处必为 0 |
| III 型 | 奇数 | 奇对称 | 处为 0 |
| IV 型 | 偶数 | 奇对称 | 处为 0 |
窗函数法的考场句子:
- 理想频响突变,对应无限长 sinc 型冲激响应。
- 实际 FIR 必须有限长,所以乘窗截断。
- 乘窗会带来过渡带和波纹,窗越长过渡带越窄,旁瓣由窗类型决定。
小例题模板:若题目要求线性相位低通且 ,第一反应是 I 型 FIR,群时延为 个采样点。若写出的 不满足 ,就说明构造错了。
29.10 多采样率:抽取、插值和有理倍采样
抽取 倍:
抽取前必须低通抗混叠。
插值 倍先插零:
插值后必须低通抗镜像。
有理倍采样率变换:
标准顺序是:
- 先插值 倍。
- 低通滤波。
- 再抽取 倍。
中间低通滤波器要同时完成抗镜像和抗混叠。
二级结论:变采样率题不要只改时间下标;必须画频谱,检查镜像和混叠。
继续加深:抽取和插值不是简单的“点变少、点变多”,它们都会改变频谱排列。
抽取 倍的频域风险:频谱被压缩后相加,若原信号不先限带,就会混叠。抗混叠低通的截止频率通常按新 Nyquist 频率来定。
插值 倍的频域风险:插零后频谱会出现镜像,需要低通滤掉镜像。抗镜像低通保留原基带,去掉周期副本。
有理倍 的标准结构:
中间滤波器的通带和截止要同时满足两件事:
- 去掉插值产生的镜像。
- 保证后续抽取不会混叠。
考场检查:若题目问从 变到 ,不要只写最终采样率;还要说明先升采样到 ,滤波后再降采样到 。
29.12 本章综合题型:从题干到工具的选择
第 29 章的专题跨度大,最后要形成一个工具选择表。看到题干时先抓关键词:
| 题干关键词 | 优先工具 | 第一行 |
|---|---|---|
| 初始电容电压、初始电感电流 | 单边 Laplace | 写元件单边变换式 |
| 极点、单位圆、频响峰 | 零极点图 | 先画极点和零点 |
| 相似、延时检测、同步 | 相关函数 | 写滑动内积 |
| 白噪声、已知波形检测 | 匹配滤波 | 写 |
| 有限长采样数据频谱 | DFT 谱分析 | 写 |
| 快速算法 | FFT | 说明 FFT 是 DFT 算法 |
| 模拟滤波器转数字 | IIR 设计 | 判断是否预畸变 |
| 线性相位数字滤波器 | FIR 设计 | 写 |
| 抽取、插值、采样率转换 | 多采样率 | 画频谱并加低通 |
| 复杂反馈图 | Mason 公式 | 数前向通路和回路 |
综合题的稳妥流程:
- 先定域: 域、 域、 域、 域、频域还是 DFT 栅格。
- 再定对象:信号、系统、滤波器、状态、谱估计还是采样率变换。
- 写第一行公式,第一行必须能体现题型。
- 算完后做条件检查:ROC、稳定性、采样率、窗函数、线性相位、初始条件。
- 最后用单位和极限情况验算。
31.4 滤波器题先看指标,再决定工具
滤波器最常考的不是“名字会不会背”,而是你能不能把指标翻译成设计动作。
先记四类最常见滤波器:
- 低通:保留低频,抑制高频。
- 高通:保留高频,抑制低频。
- 带通:保留某个频带。
- 带阻或陷波:压低某个频带。
再记六个常用指标:
- 通带
- 阻带
- 过渡带
- 截止频率
- 通带纹波
- 阻带衰减
若题目从模拟原型变数字滤波器,常见两条路:
- 脉冲响应不变法:保留模拟冲激响应的采样对应关系,但可能混叠。
- 双线性变换法:把左半 平面映到单位圆内,不发生混叠,但会引入频率非线性压缩,因此要预畸变。
双线性变换最该先写出的式子是
对应的数字频率预畸变关系是
FIR 与 IIR 的课堂判别也要清楚:
- IIR 阶数通常更省,但稳定性和相位要额外盯紧。
- FIR 更容易做到线性相位,若题目强调“波形不失真”“群时延固定”,优先联想到线性相位 FIR。
例题 3:数字低通设计题给出通带边缘数字角频率 ,并说明采用双线性变换法。第一行应写什么?
参考入口:
第一行先写预畸变,而不是直接把 当作模拟频率:
因为双线性变换下频率映射不是线性的。
31.7 模拟原型滤波器三兄弟:Butterworth、Chebyshev、Bessel
前面我们已经讲了数字滤波器的两条设计主路,但很多学生一到“先选模拟原型”这一步就空掉。这里把最常见的三个原型补成一句话能分清的程度。
先记总原则:
- 先看题目更在意幅度平坦、过渡带陡峭,还是相位保真。
- 原型滤波器的选择,本质上是在“幅频指标”和“相频指标”之间取舍。
三个最常见原型如下。
Butterworth:
- 通带单调平滑,没有纹波。
- 过渡带下降速度中等。
- 适合题目强调“通带尽量平坦”,又没有特别苛刻的过渡带要求。
Chebyshev:
- 用通带或阻带纹波换更陡的过渡带。
- 同阶数下,选择性通常比 Butterworth 更强。
- 适合题目强调“阶数尽量低、过渡带尽量窄”,且允许纹波。
Bessel:
- 相位更接近线性,群时延更平坦。
- 波形失真较小。
- 幅频滚降不如前两者陡。
所以考场可以直接这么判:
- 要“幅度平坦”,先想 Butterworth。
- 要“过渡带更陡、可接受纹波”,先想 Chebyshev。
- 要“波形不失真、相位更稳”,先想 Bessel。
例题 5:若题目强调“脉冲波形尽量不失真传输”,三种模拟原型里该优先考虑哪一个?
参考入口:
优先考虑 Bessel,因为它更接近线性相位,群时延变化更小,时域波形更不容易被拉扯变形。
31.8 数字滤波器设计时,先把指标翻译成动作
很多题其实不是不会设计,而是第一眼没把文字指标翻成操作步骤。这里给一个最小流程。
若题目给的是低通指标,例如通带边缘、阻带边缘、通带最大衰减、阻带最小衰减,则第一步要做四件事:
- 判定是低通、高通、带通还是带阻。
- 看是模拟指标还是数字指标。
- 若是双线性变换法,先做预畸变。
- 再决定选 Butterworth、Chebyshev 还是其他原型。
一句考场版口诀:
先判类型,再看域,再做预畸变,最后选原型。
例题 6:题目给出数字低通指标 ,并要求“通带尽量平坦”。你第一段思路应如何写?
参考入口:
先说明这是数字低通题;若采用双线性变换法,先把 预畸变为模拟频率;再因为题目强调“通带尽量平坦”,优先选 Butterworth 原型。