巴特沃斯低通滤波器 matlab,基于MATLAB做巴特沃斯低通滤波器..doc
时间: 2023-09-11 10:05:02 浏览: 60
巴特沃斯低通滤波器是一种常用的数字滤波器,可用于信号处理中的低频信号滤波。MATLAB提供了butter函数来设计巴特沃斯低通滤波器。下面是一个基于MATLAB实现的巴特沃斯低通滤波器的例子。
首先,我们需要定义一些参数,包括采样频率、通带截止频率、阻带截止频率和通带最大衰减(dB):
```matlab
fs = 1000; % 采样频率
fpass = 100; % 通带截止频率
fstop = 150; % 阻带截止频率
Ap = 1; % 通带最大衰减(dB)
```
然后,我们可以使用butter函数来设计巴特沃斯低通滤波器:
```matlab
[n, Wn] = buttord(fpass/(fs/2), fstop/(fs/2), Ap, 60);
[b, a] = butter(n, Wn, 'low');
```
其中,`buttord`函数用于计算滤波器的阶数和截止频率,`butter`函数用于计算滤波器的系数。
最后,我们可以使用`filter`函数来应用滤波器:
```matlab
x = sin(2*pi*50*(0:1/fs:1));
y = filter(b, a, x);
```
其中,`x`是一个包含50Hz正弦波的信号,`y`是应用了巴特沃斯低通滤波器后的信号。
完整的MATLAB代码如下:
```matlab
% 定义参数
fs = 1000; % 采样频率
fpass = 100; % 通带截止频率
fstop = 150; % 阻带截止频率
Ap = 1; % 通带最大衰减(dB)
% 设计滤波器
[n, Wn] = buttord(fpass/(fs/2), fstop/(fs/2), Ap, 60);
[b, a] = butter(n, Wn, 'low');
% 应用滤波器
x = sin(2*pi*50*(0:1/fs:1));
y = filter(b, a, x);
% 绘制结果
subplot(2,1,1);
plot((0:length(x)-1)/fs, x);
title('原始信号');
xlabel('时间(秒)');
ylabel('幅值');
subplot(2,1,2);
plot((0:length(y)-1)/fs, y);
title('滤波后信号');
xlabel('时间(秒)');
ylabel('幅值');
```
运行此代码将生成一个包含原始信号和滤波后信号的图形。