如何消除FFTW实现NFM解调生成WAV文件中的咔哒声？

1. 定位咔哒声来源的步骤 #

• 排查WAV格式与波形问题：用Audacity这类音频工具打开生成的WAV文件，查看波形图——如果存在突然的电平跳变（比如从0突变到高电平、或电平骤升骤降），基本可以锁定是数据不连续或写入错误导致。同时核对WAV的采样率、位深是否与你输出的音频参数完全匹配，参数不匹配也会引发播放异常。
• 验证IQ数据连续性：在解调前，把原始IQ数据导出为文件，或实时打印相邻数据块的首尾样本差值。如果相邻块的样本值出现突变，说明SDR设备接收时存在丢包，或是SoapySDR的readStream调用未处理好欠载/过载情况。
• 隔离处理环节：跳过解调步骤，直接将IQ数据的实部/虚部生成单声道WAV文件。如果此时仍有咔哒声，说明问题出在SDR数据读取或WAV写入环节；若没有，则再聚焦解调算法排查。
• 检查时钟同步状态：确认SDR设备时钟是否稳定，外接时钟需检查连接，内部时钟可通过SoapySDR API查询偏移量——时钟不同步会导致采样率不准，累积后引发数据错位。

2. 消除或减轻咔哒声的步骤 #

• 修复数据读取连续性：调用SoapySDR的readStream时，必须校验返回的样本数是否等于请求数量。若返回值小于请求数，说明发生欠载，此时不能直接将不完整数据块送入处理，可选择丢弃或做平滑补零；同时调整SDR硬件缓冲区参数（比如增大SOAPY_SDR_BUFFLEN）。
• 优化WAV写入逻辑：避免每次处理小块数据就调用写入函数，尽量累积足够大的数据块后再写入；如果是边处理边写入，要在最后修正WAV文件头的总样本数，避免头信息与实际数据长度不匹配。
• 添加块间平滑过渡：在相邻数据块的边界处，应用短时长的淡入淡出（比如首尾各10个样本做线性增益过渡），即使存在小幅电平突变，也能掩盖咔哒声。
• 校正时钟偏移：若为时钟同步问题，可启用SDR的自动频率控制（AFC），或手动校准采样率偏移——比如用已知信号的频率计算实际采样率与设置值的偏差，再调整SoapySDR的采样率参数。
• 防止样本溢出截断：确保解调后的音频样本落在WAV格式的合法范围内（比如16位WAV是-32768到32767），超出范围时不能直接截断，需添加限幅器将数值限制在合法区间，避免电平突变。

3. FFT或FM解调是否会导致咔哒声？ #

从你描述的“多调制模式都有咔哒声、仅WFM不明显”来看，大概率不是FFT或FM解调算法本身的问题——不同调制的解调逻辑差异极大，共用问题的可能性极低，但仍可排查以下几点：

• 如果FFT是解调环节的一部分（比如部分数字解调方案），需确认是否使用了合适的窗口函数（如汉明窗），若直接对数据块做FFT后截断，边界的频谱泄漏可能引入突变，但这种情况一般不会全程出现咔哒声。
• FM解调常用的鉴频器（如希尔伯特变换法、差分法），若处理数据块时未保留块间相位连续性，会导致每个块起始相位突变，进而产生咔哒声——比如用希尔伯特变换时，需保留上一块的末尾相位信息作为下一块的初始参考。
• 检查解调后的音频样本是否存在NaN或无穷大值，这类异常值会直接引发播放咔哒声，可能是解调计算中出现除以零等错误导致。

内容的提问来源于stack exchange，提问作者StefenMA