matlab 帕多瓦数列 通项公式_分享matlab程序之——滤波器（带通，带阻）

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

1.带通滤波器

function y=bandp(x,f1,f3,fsl,fsh,rp,rs,Fs)

%带通滤波

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

%即，f1,f3,fs1,fsh,的值小于 Fs/2

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

% f 1：通带左边界

% f 3：通带右边界

% fs1：衰减截止左边界

% fsh：衰变截止右边界

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

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

%FS：序列x的采样频率

% f1=300;f3=500;%通带截止频率上下限

% fsl=200;fsh=600;%阻带截止频率上下限

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

% Fs=2000;%采样率

%

wp1=2*pi*f1/Fs;

wp3=2*pi*f3/Fs;

wsl=2*pi*fsl/Fs;

wsh=2*pi*fsh/Fs;

wp=[wp1 wp3];

ws=[wsl wsh];

%

% 设计切比雪夫滤波器；

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

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