巴特沃斯低通滤波程序matlab

以下是一个简单的巴特沃斯低通滤波器的 MATLAB 程序：

% 设计巴特沃斯低通滤波器
fc = 1000; % 截止频率
fs = 8000; % 采样频率
order = 4; % 阶数

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

% 读取音频文件
[x, fs] = audioread('audio.wav');

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

% 播放原始音频和滤波后的音频
sound(x, fs);
pause(length(x)/fs);
sound(y, fs);

% 绘制频率响应曲线
freqz(b, a);

其中，`fc` 是滤波器的截止频率，`fs` 是采样频率，`order` 是滤波器的阶数。`b` 和 `a` 分别是滤波器的分子和分母系数。