Question

MATLAB中如何实现巴特沃斯滤波器处理信号，在线等。。。

一段正弦电信号，中间有电压跌落，想把该信号通过2阶，100Hz截止频率的巴特沃斯滤波器处理下，对这个东西一点也不懂，请问应怎么实现。。。网上其他的看不懂。。。

t=0:0.0002:0.2;
u=sin(314*t);
u1=0.5*(314*t+pi/6);
u(300:700)=u1(300:700);

现在就是想对u处理一下，画出个处理后的图出来。。。请高手指教。。。

Answer 1

说实话，网上一大堆现成的，随便找本matlab信号处理的书也是现成的代码
给你贴一个我自己写的吧，不过我用的是切比雪夫，你自己改成巴...就行了
x是输入
function y=bandp(x,fp1,fp2,rp,rs,Fs)
%带通滤波
%x：需要滤波的序列
%fp1：通带左边界
%fp2：通带右边界
%fs1：衰减截止左边界
%fs2：衰减截止右边界
%rp：通带区衰减DB数设置
%rs：截止区衰减DB数设置
%Fs：序列x的采样频率
fp1=2;fp2=35;
fs1=1;fs2=40;
rp=3;rs=15; %%%%通带区衰减DB数设置 截止区衰减DB数设置
Fs=500;
wp1=0.08*pi;
wp2=0.14*pi;
ws1=0.01*pi;
ws2=0.16*pi;
wp=[wp1 wp2];
ws=[ws1 ws2];
%设计切比雪夫滤波器
[n,wn]=cheb1ord(ws/pi,wp/pi,rp,rs);
[bz1,az1]=cheby1(n,rp,wp/pi);
%查看设计滤波器的曲线
[h,w]=freqz(bz1,az1,256,Fs);
h=20*log10(abs(h));
figure;
ploy(w,h);
title('所设计滤波器的通带曲线');
grid on;
y=filter(bz1,az1,x);
end