Arduino RP2040录制3kHz音频去噪失败：报错及无声音问题求助

核心原因 #

1. 截止频率违反奈奎斯特采样定理：你的音频采样率为3kHz，根据奈奎斯特定理，有效频率上限为采样率的一半（即1500Hz）。但代码中设置的highcut=4000Hz和lowcut=2000Hz都远高于这个上限，导致滤波器无法生成有效系数，滤波后的数据出现大量无效值（如NaN）。
2. 类型转换时的数值异常：无效值转换为int16类型时触发RuntimeWarning，最终生成的音频文件全是无效数据，自然没有声音。

解决措施 #

• 修正截止频率范围：将lowcut和highcut调整到1500Hz以内，比如设置为lowcut=300.0、highcut=1200.0，可根据你要保留的实际音频频段灵活调整。
• 添加数值校验与处理：滤波后去除NaN或无穷大值，将浮点数据归一化到int16的合法范围（-32768至32767），避免转换时的异常。
• 优化滤波器参数：降低滤波器阶数（比如从6调整到4），避免高阶滤波器带来的数值不稳定；或使用filtfilt代替lfilter，实现零相位滤波，同时提升数值稳定性。
• 支持立体声处理：如果原音频是立体声，需要分别处理左右声道，避免维度不匹配导致的问题。

修正后的示例代码 #

from scipy.signal import lfilter, butter
from scipy.io import wavfile
import numpy as np

# 调整为符合3kHz采样率的有效频段
lowcut = 300.0
highcut = 1200.0


def butter_bandpass(lowcut, highcut, fs, order=5):
    nyq = 0.5 * fs
    # 确保归一化后的频率在有效区间内（0~1）
    low = max(lowcut / nyq, 0.01)
    high = min(highcut / nyq, 0.99)
    b, a = butter(order, [low, high], btype='band', fs=fs)
    return b, a

def butter_bandpass_filter(data, lowcut, highcut, fs, order=5):
    b, a = butter_bandpass(lowcut, highcut, fs, order=order)
    y = lfilter(b, a, data)
    # 替换无效值为0
    y = np.nan_to_num(y)
    return y

def bandpass_filter(buffer):
    return butter_bandpass_filter(buffer, lowcut, highcut, sr, order=4)


sr, data = wavfile.read('3000.wav')
# 处理立体声/单声道情况
if len(data.shape) > 1:
    filtered = np.column_stack([bandpass_filter(data[:, i]) for i in range(data.shape[1])])
else:
    filtered = bandpass_filter(data)

# 归一化到int16的合法范围并转换类型
filtered = np.clip(filtered, -32768, 32767)
filtered = filtered.astype(np.int16)

wavfile.write('3000_filter.wav', sr, filtered)

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