fir带通滤波器去噪matlab

以下是使用MATLAB实现FIR带通滤波器去噪的示例代码：

% 生成含有3种频率和白噪声的信号
fs = 1000; % 采样率
t = 0:1/fs:1; % 时间序列
f1 = 50; % 第一个频率成分
f2 = 150; % 第二个频率成分
f3 = 300; % 第三个频率成分
x = sin(2*pi*f1*t) + sin(2*pi*f2*t) + sin(2*pi*f3*t) + randn(size(t));

% 设计FIR带通滤波器
N = 100; % 滤波器阶数
fcuts = [100 200]; % 截止频率
mags = [0 1 0]; % 通带增益和阻带增益
devs = [0.01 0.05 0.01]; % 通带和阻带的最大允许波动
fsamp = fs; % 采样率
b = firpm(N, fcuts, mags, devs, fsamp);

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

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

这段代码首先生成了一个包含三种频率和白噪声的信号。然后使用`firpm`函数设计了一个FIR带通滤波器，指定了滤波器的阶数、截止频率、通带增益和阻带增益等参数。最后，使用`filter`函数将滤波器应用到原始信号上，得到滤波后的信号。绘制了原始信号和滤波后的信号的波形图。