定点 DSP 的 Q15 与饱和运算：音频算法为什么离不开它

做音频、电机控制、滤波器时，DSP 工程师经常在定点和浮点之间做选择。很多人知道"定点便宜、浮点好写"，但少有人讲透：定点 DSP 真正的核心竞争力，不是"便宜"，而是一整套**硬件级保护机制**——Q 格式与饱和运算。

## Q15 是什么

Q15 是定点小数最常见的格式：16 位二进制数表示 [-1, 1) 范围内的小数。最高位是符号位，剩下 15 位是小数位。比如十进制 0.5 写成 Q15 就是 0x4000（16384）。

为什么选这个范围？因为滤波器系数、音频采样值、PID 参数天然都落在 [-1, 1) 附近。用 Q15，一个 16 位乘法器的单次乘加就能完成 FIR 滤波的核心运算：

```c
// Q15 乘法：两个 Q15 相乘，结果是 Q30，需要右移 15 位回到 Q15
acc += (int32_t)a * b >> 15;
```

这行代码在 DSP 上是一条 MAC（乘加）指令，单周期完成。浮点 DSP 同样做得到，但**定点 DSP 的 MAC 单元面积更小、功耗更低、主频利用率更高**，在大批量低成本场景（家电、电机驱动）里是硬优势。

## 饱和运算：定点的"安全气囊"

浮点数溢出会变成 inf，还能继续算；定点数溢出会**回绕**——比如 Q15 里 0.9 + 0.9 = 1.8，超出 [-1, 1) 范围，16 位回绕后变成 -0.2。音频里这意味着刺耳的爆音，电机里意味着转矩突变。

饱和运算就是硬件在溢出时**钳位**：1.8 被钳成 0.99997（最大值），-1.8 钳成 -1。输出失真，但方向正确、不会突变极性。DSP 的 ALU 通常有一位饱和标志位（SAT），配套的加减乘指令（如 `SSUB16`）在硬件层面做钳位，**零额外周期**。

## 工程实践要点

1. **标定（Scaling）是定点算法的第一道工序**：每个变量先仿真出动态范围，再定 Q 格式。Q15 不够用就上 Q31（32 位），或者用块浮点（block floating point）。
2. **累加器位宽要留余量**：FIR 滤波几十阶累加，中间结果很容易爆 32 位。DSP 通常配 40 位累加器，多出的 8 位就是"保护位"（guard bits），最后一步再饱和输出。
3. **除法尽量避开**：定点除法慢且精度损失大，能改成乘以倒数的都改掉；倒数用查表 + 牛顿迭代精修。
4. **从浮点原型迁移时**：先在 PC 上用浮点跑通算法，再逐变量转定点、用 MATLAB/Python 的定点工具箱做量化噪声分析，别直接上手写 Q15。

定点不是"穷人的浮点"，它是一套为信号处理优化过的完整运算体系。理解了 Q 格式和饱和，你才算真正入门了定点 DSP。

补充：定点算法里还需要注意乘法结果的 Q 格式变化（Q15×Q15=Q30，需要移位回 Q15）、累加器保护位（如 40 位累加器）以及舍入方式（截断 vs 四舍五入）对噪底的影响。你们在把浮点算法转成 Q15 时，是用 MATLAB 定点工具箱做仿真对比，还是直接在目标板上测？