四阶巴特沃斯滤波器 matlab

在Matlab中，你可以使用`butter`函数来设计四阶巴特沃斯滤波器。该函数具有以下语法：

[b, a] = butter(n, Wn, 'ftype')

其中：
- `n`表示滤波器的阶数，对于四阶滤波器，`n`应设置为4。
- `Wn`是归一化截止频率，其范围为0到1。截止频率是指在此频率处滤波器的增益下降3dB。对于四阶滤波器，你需要提供两个截止频率，例如 `[0.2, 0.5]`。
- `'ftype'`是滤波器类型，对于巴特沃斯滤波器，你可以使用 `'bandpass'` 表示带通滤波器。

该函数的输出是滤波器的分子和分母系数，你可以将这些系数传递给`filter`函数来应用滤波器。

以下是一个使用四阶巴特沃斯滤波器的示例：

% 设计四阶巴特沃斯滤波器
n = 4;
Wn = [0.2, 0.5];
[b, a] = butter(n, Wn, 'bandpass');

% 生成输入信号
t = 0:0.01:10;
x = sin(2*pi*1*t) + sin(2*pi*5*t) + sin(2*pi*10*t);

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

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

在这个示例中，我们设计了一个四阶巴特沃斯带通滤波器，截止频率为0.2和0.5。然后，我们生成了一个包含三个频率分量的输入信号，并使用`filter`函数应用滤波器。最后，我们绘制了原始信号和滤波后的信号。你可以根据自己的需求修改截止频率和输入信号的频率分量。