切比雪夫Ⅰ型滤波器和切比雪夫Ⅱ型滤波器的对比

设计流程图如下：
设计思想：
首先设计一个源信号和一个混合信号，通过其频谱对比得出最大和最小通带,最大和最小阻带；然后再根据得到的参数来设计切比雪夫滤波器,最后通过切比雪夫Ⅰ型滤波器和切比雪夫Ⅱ型滤波器的对比来得出那种效果好。

切比雪夫滤波器设计原理：
切比雪夫滤波器的振幅具有等波纹特性,它有两种形式：1)振幅特性在通带内是等波纹的、在阻带内是单调的切比雪夫I 型滤波器;2)振幅特性在通带内是单调的、在阻带内是等波纹的切比雪夫I Ｉ型滤波器,采用何种形式的切比雪夫滤波器取决于实际用途.
切比雪夫滤波器的设计方法就是将逼近精确度均匀分布在整个通带内，或者均匀分布在整个阻带内,或者均匀分布在两者之内,这样就可以使滤波器阶数大大降低。

切比雪夫I 型滤波器平方幅度响应函数表示为:
2
)(Ωj G =[1+2εC 2
N (Ω)]2/1- 其中 ε<1(正数)，它与通带波纹有关，ε越大,波纹也越大；C N (Ω) 是切比雪夫多项式,它被定义为:
C N （Ω)=cos （N ａr ｃｃo ｓ(Ω))，Ω≤1, C N (Ω)＝cosh （Ｎarcosh(Ω)),Ω
＞1．
而切比雪夫II 型滤波器平方幅度响应函数表示为：
)(Ωj G 2={１+2ε{ Ｃ2
N （Ω)/[2N (Ω／c Ω)］2｝}
1- 其中 ε<１(正数)，表示波纹变化情况；c Ω为截止频率;N 为滤波器的阶次,也是 C N （N ΩΩ/) 的阶次。

源信号编码及其图形：
t=-１：0．０1:1
ｙ＝(ｃos(2＊pi*10*t)+ｃos(2*pi*4０*t ））；
Ｎ＝leng ｔh(y);
fx=fft(y);
df=10０/Ｎ；
n=0:N ／2;
f=n*d ｆ;
ｓu ｂplot(2,1,１)；
p ｌot(f,ab ｓ(f ｘ(ｎ+1))＊2/Ｎ）; g ｒid;
titｌe('源波形频谱')
图（一)
混合信号编码及其图形：
ｔ=-1:0.01：１;
X＝(ｃos(2*pｉ*１０*ｔ）＋cｏs（2*ｐi＊25*t)+coｓ(2*ｐi*40*t)); N＝lｅngth(X);
ｆx=fft(X);
dｆ=100/N;
ｎ＝0:N／2；
f=n*df;
subｐlｏｔ（２,1,2);
ploｔ(f，abｓ(fx(ｎ＋1）)*２／Ｎ）; ｇrｉd；
ｔitle('混合波形频谱')
图(二)
从图（一)和图（二）对比可以得出：为了能达到和满足我们的要求,我们取以下的参数,最大通带wp２:0．5,最小通带ｗp1：0．0５,最大阻带ｗs2:０.3,最小阻带ws１:0．1。

切比雪夫Ⅰ型滤波器设计如下：
ws1＝0.１*pｉ；ws2=0.3*pi; %滤波器的阻带截止频率
wp１=0.０5*pi;wｐ２=0.5*pi；％滤波器的通带截止频率
Ｒp=1;As=２0；%滤波器的通阻带衰减指标
%转换为模拟滤波器的技术指标
T=0.01;Fs=1／T;
Ｏｍgp1=(2／Ｔ）*tａn(wｐ１/2)；Omｇp２=(2/T）＊taｎ(ｗp2/2)；
Omgｐ=[Omgp1,Omｇp2];
Oｍgｓ１＝(2／T)*tan(ws1／2）;Ｏmgs2＝(2／Ｔ）*taｎ（ws2／2)；
Oｍgs=[Omｇs１,Oｍｇs２];
bw＝Ｏｍgp2－Omgp１;w0=sqｒt(Ｏmgｐ1*Ｏmｇp2); ％模拟通带带宽和中心频率
ripｐｌｅ=10^(-Ｒp／２0); %滤波器的通带衰减对应的幅度值
Attｎ=１0^(-As／20)；%滤波器的阻带衰减对应的幅度值
%模拟原型滤波器计算
[n,Oｍgn］=cｈeb１ord（Omgp,Omｇs,Rp，Aｓ,'ｓ') %计算阶数ｎ和截止频率
[z０,p0,k０]=cheb1ap（ｎ,Rp）; %设计归一化的模拟滤波器原型
ｂａ1=ｋ0*rｅal(ｐoly(ｚ0）)；％求原型滤波器的系数ｂ
aa1=reａl(pｏｌｙ(ｐ０)); %求原型滤波器的系数ａ
[ba,aａ］=lｐ2bs（ｂa1,aa１,ｗ0，bw)；
%用双线性变换法计算数字滤波器系数
[bd,ad]=ｂiｌinｅar(bａ,aa，Fs)
%求数字系统的频率特性
[Ｈ,w]=frｅqｚ（bd,ａｄ);
dbH=２0*log10((ａbs(H)+eps)／maｘ(abs(H)））;
subplot(２,2,1);pｌot(ｗ/ｐi,aｂs(Ｈ)）；
ylabeｌ('｜H|');ｘlabel('频率(\pi)')；tiｔlｅ(＇幅度响应＇);aｘｉs(［0,1,0，1.１]); sｅt(gca,'ＸＴｉckMｏde＇,'maｎual'，'XＴｉｃｋ＇,[０,0.２,0．3，0.7,0.8］）; sｅt(gca，'ＹTickMｏde','ｍanｕal＇,'YＴick',[0，Attn,ripple，1]）;griｄ
图(三）
ｎ＝３
Oｍgn=16.１402 200.00００
bd =０.１698-０.8703 1.9961－2.５870 1.9９６1-０.87０３0．169８
ad =1.000０－2.5４50 2.533２-１．73５6 0．９60５-0.04６９-0.1６19
分析：由图（三）运行结果可知,最大通带0.5，最小通带０.０5,最大阻带0．3,最小阻带0.１;切比雪夫Ⅰ型滤波器的设计的个性技术指标精确度是均匀分布的。

而其幅度特性在通带内是等波纹的，在阻带内是单调下降的。

虽然达到了所设计滤波器的要求，滤除了我们不需要的波形,但是存在一些问题,如误差等。

切比雪夫Ⅱ型滤波器设计如下：
ｗｓ１=0.1＊pi;wｓ2=0.3*pi; %滤波器的阻带截止频率
ｗp１＝0.０5＊pi;ｗp2=０.５*pｉ; ％滤波器的通带截止频率
Rp=1;Aｓ=20;%滤波器的通阻带衰减指标
%转换为模拟滤波器的技术指标
T＝0.0１;Fs＝1/T;
Ｏｍgp1＝(2/T）＊taｎ(wｐ１/2);Omgp2=(2/T)＊ｔaｎ（wp2/2)；
Ｏmgp=[Oｍｇp1,Oｍｇp2]；
Omｇs1=(２/T)*ｔaｎ（ｗs１/2);Omgs2=(２/T)*tan(ws2/２)；
Ｏmgs=[Omgs1,Omgs2];
bw=Omgp２－Ｏｍgp１；w０=sqｒｔ(Omgp1＊Oｍgp2); %模拟通带带宽和中心频率
rｉｐｐlｅ＝10＾(-Rp／2０)；%滤波器的通带衰减对应的幅度值
Attｎ=10^(-As/２0）; ％滤波器的阻带衰减对应的幅度值
%模拟原型滤波器计算
[n,Omｇn]＝cｈeb2oｒd(Oｍｇｐ,Omgｓ,Ｒp,As，'s＇)
[z0,p0,k0］=chｅb2ap（ｎ，Ａs); ％设计归一化的ｃhｅb2型模拟滤波器原型ba1=ｋ0*real(poly（z0))；%求原型滤波器的系数b
ａa1=reａl(ｐoly(ｐ０))；％求原型滤波器的系数a
［bａ,aa]＝lp2ｂs(bａ1，aa１,w0,bw）；%变换为模拟带通滤波器
％用双线性变换法计算数字滤波器系数
［bd,ａd］=biｌｉneａr(ba,ａa,Fs）
%求数字系统的频率特性
[H,ｗ]=fｒeqz（bd,ad);
ｄbＨ=20*loｇ10（(abｓ(H)+eps)／max(abs（H））);
ｓubplot(2,2,1);plot(w/pi,ａbs(Ｈ));
ｙlabel('|H|＇）;xｌabｅl(＇频率（\pi)')；tiｔlｅ('幅度响应');axｉｓ([０,1,0，１．１]); set（gca,'XＴickＭｏde','mａnual','XＴick',［０,０.２，0.3,0．7,0.8]);
ｓet(gca，＇ＹTickMode'，'maｎｕａl',＇ＹTiｃk'，[０,Attｎ,rｉｐｐｌe，1]）；g ｒiｄ
图(四)
n =3
Omgn ＝25.７3４7 １2５．4348
bd ＝0．21７9 -0．8５32 1.５78１－1.88１9 1.5781 -0.8５3２0.2179
ad=1．０00０-2.2979 １.7４5６-0.9１83０.８465 -0.37２1 -0.0002
分析:由图（四)的运行结果可知，最大通带0.５,最小通带0.05,最大阻带0.3,最小阻带0.１;所设计的切比雪夫Ⅱ型滤波器的各项技术指标精确度是均匀分布的。

而其幅度特性在通带内是单调的,在阻带内是等波纹的。

虽然也达到了所设计滤波器的要求，滤除了我们不需要的并行,但是也还是存在误差,这是可以理解的。

由图(三)和图（四）的运行结果对比来看,切比雪夫Ⅰ型滤波器在阻带跟通带过渡时要比切比雪夫Ⅱ型滤波器过渡时平稳许多,所以切比雪夫Ⅰ型滤波器要比切比雪夫Ⅱ型滤波器所做出来的效果要好。

在实际的生活用途中,我们要根据我们的不同需求和要工作的不同场合，来选择切比雪夫Ⅰ型滤波器或者切比雪夫Ⅱ型滤波器,这样才能达到我们的要求，来满足我们的各项需要。