2. 降采样基础

在进入具体的降采样算法之前，先把几个基础概念理清楚，否则后面很容易出现“数据点少了但波形变样了”“滤波之后幅值不对了”这类问题。本节围绕四个关键词展开：采样定理、抗混叠滤波、抽取与插值、降采样因子选择。它们共同决定了降采样是否“安全”、结果是否可信。

一、采样定理：降采样的合法性边界

奈奎斯特—香农采样定理指出：对一个最高频率为 f_max 的带限信号，只要采样频率 f_s 大于 2f_max，就可以由采样序列无失真地恢复原信号。这里的 2f_max 称为奈奎斯特率，f_s/2 称为奈奎斯特频率。

降采样本质上是降低采样率。设原始采样率为 f_s，降采样因子为 M，则降采样后的采样率为 f_s' = f_s/M，新的奈奎斯特频率变为 f_s/(2M)。如果原信号中仍含有高于 f_s/(2M) 的频率成分，这些成分在降采样后会“折叠”到低频区，形成混叠。混叠一旦发生，就无法通过后续处理还原，这是降采样中最需要防范的错误。

举个教学示例：某振动信号以 f_s = 10 kHz 采集，主要关注 0～1 kHz 的轴承特征频率，但传感器同时拾取到 3.6 kHz 的电机干扰。若直接取 M = 4，降采样后 f_s' = 2.5 kHz，奈奎斯特频率为 1.25 kHz，3.6 kHz 干扰会折叠到 |3.6 − 2.5| = 1.1 kHz 处，恰好落在关注频带内，造成误判。这说明降采样前必须先处理带外成分。

二、抗混叠滤波：先滤后抽

抗混叠滤波器的任务是在抽取之前把高于新奈奎斯特频率的成分衰减掉。工程上通常使用低通滤波器，截止频率设在 f_s/(2M) 附近，并留出过渡带。滤波器的阻带衰减、过渡带宽度和阶数需要权衡：阶数越高，过渡带越窄、抑制越好，但计算量和相位延迟也越大。

实际操作顺序非常关键：先滤波，再抽取。如果先抽取再滤波，混叠已经发生，滤波只能去掉折叠后的假频，无法恢复原信号。对于实时系统，常采用多级抽取结构：第一级用较高采样率配合较宽松的滤波器，逐级降低采样率，既降低单级滤波器阶数，又减少总运算量。

三、抽取与插值：降采样与升采样的基本运算

抽取（Decimation）指每隔 M 个点保留一个点，数学上记为 y[n] = x[nM]。它直接降低数据量，但单独使用会引入混叠，因此“抽取”在工程语境中常指“抗混叠滤波 + 下采样”的完整过程。插值（Interpolation）则相反，是在样本之间插入新点以提高采样率，典型做法是先补零再低通滤波，或采用线性、样条等插值核。

抽取和插值常配合使用，用于采样率转换。例如把 48 kHz 音频转为 16 kHz，可先设计抗混叠低通，再按 M = 3 抽取；若要把 16 kHz 转为 48 kHz，则按 L = 3 插值后再滤波。对于非整数倍转换，可采用 L/M 的有理数比例结构，先插值 L 倍再抽取 M 倍，中间滤波器截止频率取两者奈奎斯特频率的较小值。

四、降采样因子怎么选

降采样因子 M 不是越大越好，需要同时满足几个约束。第一，关注频带必须完整保留：若实验真正关心的最高频率为 f_interest，则必须满足 f_s/(2M) ≥ f_interest，即 M ≤ f_s/(2f_interest)。第二，抗混叠滤波器可实现：过渡带越窄，滤波器阶数越高，需评估实时性与相位延迟。第三，数据量与后续处理匹配：降采样后仍需保证每个关注周期内有足够点数用于峰值检测或包络分析。

以教学示例说明：原始 f_s = 10 kHz，关注频带 0～1 kHz，则 M ≤ 10/(2×1) = 5。若取 M = 4，新奈奎斯特频率 1.25 kHz，过渡带为 1～1.25 kHz，相对较窄，需要较高阶滤波器；若取 M = 2，新奈奎斯特频率 2.5 kHz，过渡带 1～2.5 kHz 较宽，滤波器容易实现，但数据量只降到一半。工程上常先按关注频带确定 M 的上限，再根据滤波器资源和延迟在 2～5 之间折中。

还需要注意：降采样因子最好取整数，便于定点实现和索引对齐；若必须非整数倍，采用 L/M 结构并保证中间滤波器截止频率取 f_s·min(1/(2M), 1/(2L)) 的保守值。对于多级抽取，总因子等于各级因子之积，每级都应按本级输出采样率重新核算抗混叠截止频率。

把本节四个概念串起来：采样定理给出降采样的合法性边界，抗混叠滤波保证边界内不被污染，抽取与插值提供采样率转换的基本运算，降采样因子则是在保留信息、滤波器成本和数据量之间的折中决策。后续的平滑处理，正是在这个已经“安全降采样”的数据基础上进行的。