避坑指南：MATLAB采样函数resample/decimate的5个常见错误用法（附正确代码示例）

MATLAB信号重采样实战：避开resample、decimate与downsample的五大陷阱

信号处理工程师的日常工具箱里，MATLAB的采样率转换函数是高频使用的利器。无论是处理生物医学信号、音频数据还是通信波形，resample、decimate和downsample这几个函数都扮演着关键角色。然而，正是这些看似简单的函数，在实际工程应用中埋下了不少“坑”。新手容易混淆参数，老手也可能因惯性思维而犯错，最终导致信号失真、分析结果偏离，甚至引发难以察觉的系统性误差。这篇文章，我想结合自己处理EEG（脑电图）和振动信号时踩过的“雷”，聊聊这几个函数最常见的五种错误用法，并提供经过验证的正确代码示例。我们的目标不是复述手册，而是让你在下一个项目中，能写出更稳健、更可靠的代码。

1. 参数顺序之殇：resample(x, p, q) 到底谁是新谁是旧？

这可能是最经典、也最令人困惑的错误。resample(x, p, q)的官方描述是“通过有理因子p/q改变采样率”。但问题来了：p代表目标采样率吗？还是q代表目标采样率？一个顺序搞反，采样率可能就变成了原来的平方关系，结果完全错误。

错误重现： 假设我们有一段采样频率为200 Hz的EEG信号x_eeg，我们希望将其重采样到100 Hz。一个直觉性的错误写法是：

Fs_original = 200; % Hz
Fs_target = 100;   % Hz
% 错误写法：误以为p是目标频率，q是原始频率
y_wrong = resample(x_eeg, Fs_target, Fs_original);

这段代码看似合理，实则大错特错。在MATLAB的语法中，y = resample(x, p, q)的含义是：将信号x的采样率乘以p/q。也就是说，新的采样率 = 原始采样率 * (p / q)。

因此，正确的对应关系是：

• 你希望新采样率是原来的 Fs_target / Fs_original 倍。
• 所以，p/q = Fs_target / Fs_original。
• 对于从200 Hz降到100 Hz，Fs_target / Fs_original = 100/200 = 1/2。
• 我们需要找到一对整数p和q，使得 p/q = 1/2。最简单的就是 p=1, q=2。

正确代码示例：

Fs_original = 200; % Hz
Fs_target = 100;   % Hz
% 计算有理因子 p/q = Fs_target / Fs_original
% 使用 rat 函数获取最接近的有理近似，tol是容忍度
[p, q] = rat(Fs_target / Fs_original);
fprintf('重采样因子: p = %d, q = %d\n', p, q); % 输出 p=1, q=2
y_correct = resample(x_eeg, p, q);
% 验证：新采样率 = Fs_original * (p/q) = 200 * (1/2) = 100 Hz

注意：rat函数并非总是返回最简分数（如1/2），它返回的是在指定容忍度内的近似。对于非整数倍采样率转换（如从22050 Hz转到16000 Hz），rat函数非常有用，但务必检查返回的p和q值，避免它们过大导致计算量激增或数值不稳定。

为了更清晰，我们对比一下直觉错误与正确理解下的参数：