切比雪夫1型滤波器讲解文档格式.docx
- 文档编号:20669502
- 上传时间:2023-01-25
- 格式:DOCX
- 页数:11
- 大小:320.46KB
切比雪夫1型滤波器讲解文档格式.docx
《切比雪夫1型滤波器讲解文档格式.docx》由会员分享,可在线阅读,更多相关《切比雪夫1型滤波器讲解文档格式.docx(11页珍藏版)》请在冰豆网上搜索。
)=Ha(j?
)的低通滤波器技术指标为:
0≤∣Ha(j?
)∣≤1/A^2,?
s≤∣?
∣
其中ε为通带波动系数,p?
和s?
是通带和阻带边缘频率。
A为阻带衰减系数
从图知必须满足
其中参数ε和A是数字滤波器指标
2.3切比雪夫1型滤波器
2.3.1
切比雪夫滤波器介绍在巴特沃兹滤波器中,幅度响应在通带和阻带内都是单调的。
因此,若滤波器的技术要求是用最大通带和阻带的逼近误差来给出的话,那么,在靠近通带低频端和阻带截止频率以上的部分都会超出技术指标。
一种比较有效的途径是使逼近误差均匀地分布于通带或阻带内,或同时在通带和阻带内都均匀分布,这样往往可以降低所要求的滤波器阶次。
通过选择一种具有等波纹特性而不是单调特性的逼近方法可以实现这一点。
切比雪夫型滤波器就具有这种性质:
其频率响应的幅度既可以在通带中是等波纹的,而在阻带中是单调的(称为I型切比雪夫滤波器),也可以在通带中是单调的,而在阻带中是等波纹的(称为II型切比雪夫滤波器)。
I型切比雪夫滤波器的幅度平方函数是
率(注意:
常用的以幅度下降3分贝的频率点作为截止频率的定义不适用于切比雪夫滤波器!
)
从定义切比雪夫多项式可以直接得出由Tn(x)和Tn-1(x)求Tn+1(x)的递推公式。
将三角恒等式得Tn+1(x)=2(x)T(x)-Tn-1(x)当0<
x<
1时,Tn^2(x)在0和1之间变化;
当x>
1时,cos-是虚数,所以Tn(x)像双曲余弦一样单调地增加;
∣Hn(w)∣对于0≤w/w0≤1呈现出在1和1/(2ε^2)之间的波动;
而对于w/w0〉1单调地减小。
需要用三个参量来确定该滤波器:
ε,w0和N。
在典型的设计中,用容许的通带波纹来确定ε,而用希望的通带截止频率来确定w0。
然后选择合适的阶次n,以便阻带的技术要求得到满足。
2.3.2切比雪夫1型滤波器的设计原理
切比雪夫滤波器(又译车比雪夫滤波器)是在通带或阻带上频率响应幅度等波纹波动的滤波器。
在通带波动的为“I型切比雪夫滤波器”,在阻带波动的为“II型切比雪夫滤波器”。
切比雪夫滤波器在过渡带比巴特沃斯滤波器的衰减快,但频率响应的幅频特性不如后者平坦。
切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小,但是在通频带内存在幅度波动。
切比雪夫I型为例介绍其设计原理幅度平方函数用︱H(j)︱表示
ε为小于1的正数,表示通带内幅度波动的程度,ε越大,波动幅度也越大。
?
p称为通带截止频率。
令λ=?
/?
p,称为对?
p的归一化频率。
定义允许的通带内最大衰减αp用下式表示
定义允许的通带内最大衰减pα用下式表示αp=10lgmax︱Ha(j?
)︱^2/min︱Ha(j?
)︱^2︱?
︱≤?
p
式中max︱Ha(j?
)︱^2=1min︱Ha(j?
)︱^2=1/(1+ε^2)
因此αp=10lg(1+ε^2)ε^2=10^(0.1αp)-1
这样,可以根据通带内最大衰减αp,可求出参数ε。
阶数N影响过渡带的宽度,同时也影响通带内波动的疏密,因为N等于通带内的最大值和最小值的总个数。
设阻带的起点频率为?
s,则有
令久,=9由入>I,有
Cv(&
)=M(NarchA^)=—
可以解出
arch
列囤g)『
——archN
3dB截止频率用Q,.表示,
忆(血)卜%
按照(2-3-1)式,有
斗丄]
Ine)
经过一系列推论得归一化系统函数为
2.3.3切比雪夫低通滤波器的设计步骤
(1)确定低通滤波器的技术指标:
边带频率p,它们满足
(2)求滤波器阶数N和参数
先求出1/k1,再求出阶数N,取符合条件的N的最小整数。
(3)求归一化系统函数Ga(p)
4)将去归一化,得到实际的Ha(s)
函数说明:
用MATLAB设计切比雪夫1型低通滤波器
(1)[z,p,k]=cheb1ap(N,Rs)
该格式用于计算N阶切比雪夫1型归一化模拟滤波器的零极点和增益因子,返回长度为N的列向量Z和P,分别给出N个零点和极点的位置。
Rs是阻带最小衰减。
(2)[N,wso]=cheb1ord(wp,ws,Rp,As)该格式用于计算切比雪夫1型数字滤波器的阶数N和阻带截止频率wso,调用参数分别为数字滤波器的通带频率和阻带边界频率的归一化值。
(3)[N,wso]=cheb1ord(wp,ws,Rp,As)
该格式用于计算切比雪夫1型模拟滤波器的阶数N和阻带频率wso。
wp.ws是实际模拟角频率。
(4)[B,A]=cheby1(N,Rs,wso,'
ftype'
)该格式用于计算N阶切比雪夫1型数字滤波器系统函数的分子和分母的多项式系数向量B和A。
调用参数N和wso分别为切比雪夫1型数字滤波器的阶数和阻带截止频率的归一化值。
(5)[B,A]=cheby1(N,Rp,wso,'
'
s'
)该格式用于计算N阶切比雪夫1型模拟滤波器系统函数的分子和分母多项式系数向量B和A。
调用参数N和wso分别为切比雪夫1型模拟滤波器的阶数和阻带截止频率的归一化值。
Impinvar
功能:
用“脉冲响应不变法”将模拟滤波器离散化。
[BZ,AZ]=impinvar(B,A,Fs)Plot(X,Y)功能:
绘制线性二维图形
gridon功能:
对当前坐标系添加主要网格线。
3脉冲响应不变法脉冲响应不变法原理脉冲响应不变法是实现模拟滤波器数字化的一种直观而常用的方法。
它特别
适合于对滤波器的时域特性有一定要求的场合。
脉冲响应不变法的核心原理是通过对连续函数ha(t)等间隔采样得到离散序列ha(nT)。
令h(n)=ha(nT),T为采样间隔。
它是一种时域上的转换方法。
Hb0b1s...bNs
(3-1)
一个模拟滤波器的传递函数可以用有理分式表达式表示为:
a0a1s...aMs
通过反拉普拉斯变换我们就可以得到它的冲激相应:
1
ha(t)L1[Ha(s)]
(3-2)
脉冲响应不变法就是要保证脉冲响应不变,即:
h(n)ha(nT)
(3-3)
对上式的冲激相应序列h(n)作Z变换,就可以得到数字滤波器的传递函数:
H(z)ZT[ha(nT)]
(3-4)
一般来说,Ha(s)的分母多项式阶次总是大于分子多项式的阶次。
假定Ha(s)的没有多重极点,则式
(1)就可分解为:
NAi
Ha(s)A
i1ssi
(3-5)式中si,Ai均为复数,si是Ha(s)的极点。
其拉氏变换为
Nsi.Tn
ha(t)Aieu(n)
i1
一般来说,Ha(s)的分母多项式阶次总是大于分子多项式的阶次。
数字滤波
器的传递函数H(z)经过合并简化,成为一般形式的有理分式传递函数
H(z)
c0c1z1...cLzL
d0d1z1...dNzN
(3-6)
在讨论采样序列z变换与模拟信号拉氏变换之间关系的有关章节中,我们已经知道
12
H(z)|zesTHa(sjK)
TKT
(3-7)
按照z=eST的关系,每一个s平面上宽度为2/T的水平条带将重迭映射到z平面上。
因此脉冲响应不变法将s平面映射到z平面,不是一个简单的一一对应的关系。
对于高采样频率(T小)的情况,数字滤波器在频域可能有极高的增益。
为此我们采用
N
h(n)Tha(nT)TAieSITnu(n)(3-9)
i1在脉冲响应不变法设计中,模拟频率与数字频率之间的转换关系是线性的
(*T)。
同时,它可以保持脉冲响应不变,h(t)=ha(nT)。
因此,这一方法往
往用于低通时域数字滤波器设计及相应的模拟系统数字仿真设计。
subplot(2,1,1),plot(W/pi,20*log10(abs(H)))xlabel('
幅频/Hz'
);
ylabel('
幅度/dB'
gridon;
subplot(2,1,2),plot(W/pi,180/pi*unwrap(angle(H)))xlabel('
相频/Hz'
相位/^o'
欢MATLAB
FileEditDebugDesktopWindow
Help
DW
Wof
CurrentDirectory:
&
WIATLAB7\work二]』也
Shortctts
上]HowtoAdd上]What'
sNew
»
fp=100;
fs=150;
Rp=O.5;
Rs=10;
Fs=1030;
Wp=2*pi*fp;
Ws=2*pi*fs;
[N,wc]二cheblord(Wp,Ws,Rp,Rs,'
s'
乐计茸阶数N和截止频至wc
wc
[z,p,k]=cheb1ap(N,Rp);
%计尊系统函数零极点和増益
[Bap,Aap]=zp2tf(z,p,k);
%将系统函数的寥极点转化为系统因数一般形式的系数
[b,a]=lp21p(Bap,Aap,wc);
%将複拟低通廳波器原型,转换为模拟低通1虑波器
[bz,az]=impinvar(b,a,Fs);
%采用脉冲响应不娈法,将模拟.低鱼虑波器,转换为新宇低通慮皮器
figured)
[H,W]=freqz(bz,az,Fs);
%计算频奉H和频率响应甲
subplot(2j1,1),plot(I/pij20*log10(abs(H)))
xlabelC幅频/Hz'
);
ylabel(*幅度/dB‘);
subplot⑵lj2)plot(W/pi,180/pi*unwrap(angle(H)))
xlabeK'
相频/Hz,);
相位/o’);
wc=
62&
3185
验证滤波器的正确性程序
figure
(2);
f1=50;
n=150;
m=0:
(n-1);
t=m/Fs;
%采样点数subplot(2,3,1)
x=sin(2*pi*f1*t);
plot(t,x);
title('
f1输入信号'
f2=200;
x=sin(2*pi*f2*t);
subplot(2,3,2)plot(t,x);
f2输入信号'
subplot(2,3,3)x=sin(2*pi*f1*t)+sin(2*pi*f2*t);
叠加后输入信号'
y=filtfilt(bz,az,x);
subplot(2,3,4)
plot(t,y);
滤波之后的信号'
)subplot(2,3,5)
plot(t,sin(2*pi*f1*t));
希望输出的信号'
入MATLAB
:
ileEditDebugDesktopWindow
21Q
e:
\MATLAB7\work
Z1』⑥
shortcuts
臣]HowtoAdd回WhafsNew
>
^p=0.5;
Rs=10;
s=1000;
[N,wc]=cheb1ord(Wp,Ws^Rp,Rs,s'
);
%计算阶数N和截止频率初c
VC
[z,p,k]=cheblap(N,Rp);
%计算系统函数霎扱点和增缢
[Bap,Aap]=zp2tf(z,p,k);
%将系纟充函数的零极点转化为杀绒函数一般形式的系数
[b,a]=lp21p(bap,Aap,wc);
%将模拟低通戏波器原型>
牧换为模拟低通滤披畚
[bz,az]=impinvar(b,a,Fs);
%采用脉冲响应不变法,捋複拟低通慮皮器,转换为埶宁低通痢皮器
[H,W]=freqz(bz,az,Fs);
%计算频率H和频率响应W
mbplot(2,1,1),plot(W/pi,20*logl0(abs(H)))dabelCfl频/Hz,);
ylabelC度/dB、);
subplot(2,1,2),plot(W/pi,180/pi*unwrap(angle(H)))
dabelC相频/Hz,);
ylabel(?
相位/o'
)fgridon;
igure
(2);
l=50;
i=150;
(n-1);
%采样点数
mbplot(2,3,1)
c=sin(2*pi*f1*t);
^lot(t,x);
titleCf1输入信号’);
2=200;
c=sin(2*pi*f2*t);
subplot(2,3,2)
)lot(t,x);
titleCf2输入信号’);
subplot(2,3,3)
c=sin(2*pi*f1*t)+sin(2*pi*f2*t);
lot(t,x)iJitleC输入信号’);
^=filtfilt(bz$az^x);
subplot(2,3,4)
Jlot(t,y);
titleC滤波之后的信号’)
subplot(2,3,5)
)1ot(t,sin(2*pi*fl*t)):
title('
幕望输出的信号’):
肌Start|
4.3实验结果分析
求得阶数N=3,通带边界频率wc=628.3125,频率在0-100hz(0-0.2)完全通过,在100-150hz(0.2-0.3)通过幅度递减,150hz(0.3-∞)后被截止。
验证信号f1=50hz,f2=200hz,观察验证后的波形,可以看出f1通过,而f2被截止,数字低通滤波器的阻带截止频率为150hz,所以设计的数字低通滤波器符合技术指标
5实验总结与体会
通过对数字信号处理的课程设计,熟悉了MATLAB的运行环境,了解更多有关于Matlab软件的知识,初步掌握了MATLAB语言在数字信号处理中一些基本库函数的调用和编写基本程序等应用;
通过对数字低通滤波器的设计让我熟悉了滤波器设计的一般原理,对滤波器有了一个感性的认识;
学会了数字低通滤波器设计的一般步骤;
加深了对脉冲响应不变法的理解和认识。
在设计的过程中,我也认识到了自己所学知识的不足。
以前上课都是学一些基本的东西,自以为知识掌握的很熟练,实际只是了解了最表面的东西,好多理论只是不能很好的应用于实践,现在运用学到得的东西做出有实际应用价值的东西,对所学知识点进一步的理解,并进行系统化。
这也让我再次认识到知识是无尽的,只有不断的充实自己、完善自己的知识理论体系,才能够更好的胜任自己以后的工作。
设计过程中知识的不足也让我更加坚定了终身学习的决心。
6参考文献数字信号处理(第三版)高西全丁玉美MATLAB基础及应用(第三版)于润伟朱晓慧
MATLAB辅助现代工程数字信号处理(第二版)李益华
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 滤波器 讲解
![提示](https://static.bdocx.com/images/bang_tan.gif)