带通滤波器matlab

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

% 设计带通滤波器
Fs = 1000;       % 采样率
f1 = 50;         % 通带截止频率1
f2 = 200;        % 通带截止频率2
Wp = [f1 f2]/(Fs/2);  % 截止频率归一化
Rp = 1;          % 通带最大衰减量
Rs = 60;         % 阻带最小衰减量
[n, Ws] = buttord(Wp(1), Wp(2), Rp, Rs); % 计算阶数和截止频率
[b, a] = butter(n, Ws, 'bandpass'); % 滤波器系数

% 使用滤波器过滤信号
t = 0:1/Fs:1;    % 时间序列
x = sin(2*pi*100*t) + sin(2*pi*300*t) + sin(2*pi*500*t); % 带噪声的信号
y = filter(b, a, x); % 过滤后的信号

% 绘制原始信号和滤波后的信号
figure;
subplot(2,1,1);
plot(t,x);
title('原始信号');
xlabel('时间 (s)');
ylabel('幅度');
subplot(2,1,2);
plot(t,y);
title('滤波后的信号');
xlabel('时间 (s)');
ylabel('幅度');

该示例使用了 Butterworth 滤波器设计方法，并将其应用于一个包含三个频率分量的信号中。可以通过调整截止频率、通带最大衰减量和阻带最小衰减量来优化滤波器效果。