指数高通滤波器代码_分享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);

你可能感兴趣的:(指数高通滤波器代码)