指数高通滤波器代码_分享matlab程序之——滤波器篇（高通，低通）

快毕业了，把自己写的现成的matlab函数分享给有需要的人，由于个人水平有限，写的不好请见谅，愿意拍砖的尽管拍好了。目前还不考虑读博，所以写的程序仍了可惜，所以就拿出来分享。好了不废话了，开始正题。

以下两个滤波器都是切比雪夫I型数字滤波器，不是巴特沃尔滤波器，请使用者注意！

1.低通滤波器

使用说明：将下列代码幅值然后以m文件保存，文件名要与函数名相同，这里函数名：lowp。

function y=lowp(x,f1,f3,rp,rs,Fs)

%低通滤波

%使用注意事项：通带或阻带的截止频率的选取范围是不能超过采样率的一半

%即，f1,f3的值都要小于 Fs/2

%x:需要带通滤波的序列

% f 1：通带截止频率

% f 3：阻带截止频率

%rp：边带区衰减DB数设置

%rs：截止区衰减DB数设置

%FS：序列x的采样频率

% rp=0.1;rs=30;%通带边衰减DB值和阻带边衰减DB值

% Fs=2000;%采样率

%

wp=2*pi*f1/Fs;

ws=2*pi*f3/Fs;

% 设计切比雪夫滤波器；

[n,wn]=cheb1ord(wp/pi,ws/pi,rp,rs);

[bz1,az1]=cheby1(n,rp,wp/pi);

%查看设计滤波器的曲线

[h,w]=freqz(bz1,az1,256,Fs);