matlab 滤波器代码

以下是一个简单的MATLAB滤波器代码示例：

% 设计一个50Hz的低通滤波器
fs = 1000; % 采样率
fc = 50; % 截止频率
[b, a] = butter(4, fc/(fs/2), 'low'); % 4阶Butterworth滤波器设计
% b为分子系数，a为分母系数

% 产生一个随机信号
t = 0:1/fs:1; % 时间向量
x = sin(2*pi*100*t) + sin(2*pi*200*t) + sin(2*pi*300*t) + randn(size(t));

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

% 绘图
subplot(2,1,1);
plot(t, x);
title('原始信号');
subplot(2,1,2);
plot(t, y);
title('滤波后信号');

在上面的示例中，我们使用Butterworth滤波器设计函数butter来创建一个低通滤波器。我们使用randn函数创建一个随机信号，并使用filter函数将其通过滤波器进行滤波。最后，我们使用subplot和plot函数绘制原始信号和滤波后信号的图形。