巴特沃斯低通滤波器matlab代码

function [y,t] = butterworth_lowpass_filter(x, fs, cutoff_freq, order)

% Butterworth lowpass filter implementation in MATLAB
% x - input signal
% fs - sampling frequency
% cutoff_freq - cutoff frequency
% order - filter order

% Design the filter
Wn = cutoff_freq / (fs/2);
[b,a] = butter(order, Wn, 'low');

% Filter the signal
y = filter(b,a,x);

% Time vector
t = 0:1/fs:(length(x)-1)/fs;

% Plot the input and filtered signals
figure;
subplot(2,1,1);
plot(t,x);
xlabel('Time (s)');
ylabel('Amplitude');
title('Input Signal');
subplot(2,1,2);
plot(t,y);
xlabel('Time (s)');
ylabel('Amplitude');
title('Filtered Signal');

end

相关问题

巴特沃斯低通滤波器 matlab,基于MATLAB做巴特沃斯低通滤波器..doc

巴特沃斯低通滤波器是一种常用的数字滤波器，可用于信号处理中的低频信号滤波。MATLAB提供了butter函数来设计巴特沃斯低通滤波器。下面是一个基于MATLAB实现的巴特沃斯低通滤波器的例子。

首先，我们需要定义一些参数，包括采样频率、通带截止频率、阻带截止频率和通带最大衰减（dB）：

fs = 1000;   % 采样频率
fpass = 100; % 通带截止频率
fstop = 150; % 阻带截止频率
Ap = 1;      % 通带最大衰减（dB）

然后，我们可以使用butter函数来设计巴特沃斯低通滤波器：

[n, Wn] = buttord(fpass/(fs/2), fstop/(fs/2), Ap, 60);
[b, a] = butter(n, Wn, 'low');

其中，`buttord`函数用于计算滤波器的阶数和截止频率，`butter`函数用于计算滤波器的系数。

最后，我们可以使用`filter`函数来应用滤波器：

x = sin(2*pi*50*(0:1/fs:1));
y = filter(b, a, x);

其中，`x`是一个包含50Hz正弦波的信号，`y`是应用了巴特沃斯低通滤波器后的信号。

完整的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('幅值');

运行此代码将生成一个包含原始信号和滤波后信号的图形。

巴特沃斯低通滤波器matlab

要在MATLAB中实现巴特沃斯低通滤波器，可以使用butter函数和filter函数。巴特沃斯滤波器是一种IIR滤波器，用于滤除信号中的高频成分，保留低频成分。

首先，使用butter函数设计巴特沃斯滤波器。这个函数有两个参数，第一个参数是滤波器的阶数，第二个参数是截止频率。阶数越高，滤波器的陡峭度越高。截止频率是一个0到1之间的值，表示截止频率占取样频率的比例。

例如，使用以下代码可以设计一个阶数为4，截止频率为0.5的巴特沃斯低通滤波器：
[b, a] = butter(4, 0.5, 'low');

然后，可以使用filter函数应用设计好的滤波器到信号上。这个函数有三个参数，第一个参数是滤波器的系数，第二个参数是输入信号，第三个参数是输出信号的初始状态。

例如，使用以下代码可以将滤波器应用到信号上：
filtered_signal = filter(b, a, input_signal);

这样就可以得到经过巴特沃斯低通滤波器滤波后的信号filtered_signal。