巴特沃斯低通滤波的matlab代码

以下是巴特沃斯低通滤波器的MATLAB代码示例：

% 设置滤波器参数
fc = 1000; % 截止频率
fs = 8000; % 采样频率
order = 4; % 滤波器阶数

% 计算归一化频率
Wn = fc / (fs/2);

% 计算滤波器系数
[b, a] = butter(order, Wn, 'low');

% 读取音频文件
[y, Fs] = audioread('example.wav');

% 应用滤波器
y_filtered = filter(b, a, y);

% 播放原始音频
sound(y, Fs);

% 等待音频播放完毕
pause(length(y)/Fs);

% 播放滤波后的音频
sound(y_filtered, Fs);

其中，`fc` 是滤波器的截止频率，`fs` 是采样频率，`order` 是滤波器的阶数。`Wn` 是归一化频率，`b` 和 `a` 是滤波器的系数。通过 `filter` 函数将滤波器应用到原始音频文件上，得到滤波后的音频。最后使用 `sound` 函数播放原始音频和滤波后的音频。