带通滤波器的MATLAB代码

资 源 简 介

带通滤波器的MATLAB代码

详 情 说 明

以下是使用MATLAB编写的带通滤波器代码的示例：

```matlab

% 设计带通滤波器

fs = 1000; % 采样频率

f1 = 50; % 通带低频边缘

f2 = 200; % 通带高频边缘

f3 = 250; % 阻带低频边缘

f4 = 300; % 阻带高频边缘

rp = 1; % 通带最大衰减

rs = 60; % 阻带最小衰减

% 设计滤波器

[n, wn] = buttord([f1 f2]/(fs/2), [f3 f4]/(fs/2), rp, rs);

[b, a] = butter(n, wn, 'bandpass');

% 生成输入信号

t = 0:1/fs:1; % 时间向量

x = sin(2*pi*100*t) + sin(2*pi*300*t) + randn(size(t));

% 应用滤波器

y = filter(b, a, x);

% 绘制结果

subplot(2,1,1);

plot(t, x);

title('原始信号');

xlabel('时间');

ylabel('幅值');

subplot(2,1,2);

plot(t, y);

title('滤波后的信号');

xlabel('时间');

ylabel('幅值');

```

以上代码展示了如何使用MATLAB设计和应用带通滤波器。带通滤波器用于保留指定频率范围内的信号，并丢弃其他频率的信号。在此示例中，我们使用了巴特沃斯滤波器设计方法来创建带通滤波器，并将其应用于输入信号。通过绘制原始信号和滤波后的信号，我们可以观察到滤波器对信号的影响。

请注意，您可以根据需要调整代码中的参数和频率范围以适应您的应用场景。

```