切比雪夫1型滤波器概要.docx
《切比雪夫1型滤波器概要.docx》由会员分享,可在线阅读,更多相关《切比雪夫1型滤波器概要.docx(15页珍藏版)》请在冰豆网上搜索。
切比雪夫1型滤波器概要
1课题描述,,,,,,,,,,,,,,,,,,,,,,,,,,1
2设计原理,,,,,,,,,,,,,,,,,,,,,,,,,,1
2.1滤波器的分类,,,,,,,,,,,,,,,,,,,,,,1
2.2模拟滤波器的设计指标,,,,,,,,,,,,,,,,,,1
2.3切比雪夫1型滤波器,,,,,,,,,,,,,,,,,,,2
2.3.1切比雪夫1型滤波器的设计原理,,,,,,,,,,,3
2.3.2切比雪夫1型滤波器的设计步骤,,,,,,,,,,,3
3脉冲响应不变法
3.1脉冲响应不变法原理,,,,,,,,,,,,,,6
4设计内容,,,,,,,,,,,,,,,,,,,,,,,,,,,6
14.1设计步骤6
4.2用MATLA编程实现6
4.3设计结果分析,,,,,,,,,,,,,,,,,,,,,,10
5总结,,,,,,,,,,,,,,,,,,,,,,,,,,,,,10
6
10
参考文献
1课题描述
数字滤波器是数字信号处理的重要工具之一,它通过数值运算处理改变输入信号所含频率成分的相对比例或者滤出某些频率成分的数字器件或程序,二数字滤波器处理精度较高,体积小,稳定,重量轻,灵活,不存在阻抗匹配问题,可以实现模拟滤波器无法实现的特殊功能。
故本书课题使用MATLAB信号处理箱和运算用切比雪夫法设计数字低通滤波器。
利用脉冲响应不变法设计切比雪夫I数字低通滤波器,通带截止频率
100hz,阻带截止频率150Hz,采样频率1000hz,通带最大衰减为0.5dB,阻带最小衰减为10dB,画出幅频、相频响应曲线,并设计信号验证滤波器设计的正确性。
设计原理
2.1.滤波器的分类
(1)从功能上分;低、带、高、带阻。
(2)从实现方法上分:
FIR、IIR
(3)从设计方法上来分:
Chebyshev(切比雪夫),Butterworth(巴特沃斯)
(4)从处理信号分:
经典滤波器、现代滤波器
2.2模拟滤波器的设计指标
设ha(j?
)是一个模拟滤波器的频率响应,则基于平方幅度响应J(?
)=Ha(j?
)的低通滤波器技术指标为:
0<1Ha(j?
)I<1/AA2,?
sI
其中&为通带波动系数,P?
和s?
是通带和阻带边缘频率。
A为阻带衰减系数
从图知必须满足i'
其中参数&和A是数字滤波器指标
2.3切比雪夫1型滤波器
2.3.1切比雪夫滤波器介绍
在巴特沃兹滤波器中,幅度响应在通带和阻带内都是单调的。
因此,若滤波器的技术要求是用最大通带和阻带的逼近误差来给出的话,那么,在靠近通带低频端和阻带截止频率以上的部分都会超出技术指标。
一种比较有效的途
径是使逼近误差均匀地分布于通带或阻带内,或同时在通带和阻带内都均匀
*其中1€|
分布,这样往往可以降低所要求的滤波器阶次。
通过选择一种具有等波纹特性而不是单调特性的逼近方法可以实现这一点。
切比雪夫型滤波器就具有这种性质:
其频率响应的幅度既可以在通带中是等波纹的,而在阻带中是单调的(称为I型切比雪夫滤波器),也可以在通带中是单调的,而在阻带中是等波纹的(称为II型切比雪夫滤波器)。
丨型切比雪夫滤波器的幅度平方函数是
率(注意:
常用的以幅度下降3分贝的频率点作为截止频率的定义不适用于切比雪夫滤波器!
)
从定义切比雪夫多项式可以直接得出由Tn(x)和Tn-1(x)求Tn+1(x)
的递推公式。
将三角恒等式得Tn+1(x)=2(x)T(x)-Tn-1(x)
当01时,cos-是虚数,
所以Tn(x)像双曲余弦一样单调地增加;IHn(w)I对于0Ww/wCK1呈现出在1和1/(2£A2)之间的波动;而对于w/w0〉1单调地减小。
需要用三个参量来确定该滤波器:
&,w0和在典型的设计中,用容许的通带波纹来确定&,而用希望的通带截止频率来确定w(X然后选择合适的阶次n,以便阻带的技术要求得到满足。
2.3.2切比雪夫1型滤波器的设计原理
切比雪夫滤波器(又译车比雪夫滤波器)是在通带或阻带上频率响应幅度等波纹波动的滤波器。
在通带波动的为“I型切比雪夫滤波器”,在阻带波动的为“II型切比雪夫滤波器”。
切比雪夫滤波器在过渡带比巴特沃斯滤波器的衰减快,但
频率响应的幅频特性不如后者平坦。
切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小,但是在通频带内存在幅度波动。
切比雪夫I型为例介绍其设计原理幅度平方函数用丨H(j门)丨表示
£为小于1的正数,表示通带内幅度波动的程度,£越大,波动幅度也越大。
?
P称为通带截止频率。
令入=?
/?
P,称为对?
P的归一化频率。
定义允许的通带内最大衰减ap用下式表示
定义允许的通带内最大衰减pa用下式表示ap=10lgmax|Ha(j?
)|八2/min
IHa(j?
)|A2|?
|
p
式中max|Ha(j?
)|八2=1min|Ha(j?
)|A2=1/(1+£a?
)
因此ap=10lg(1+£A2)£A2=10A(0.1ap)-1
这样,可以根据通带内最大衰减ap,可求出参数£。
阶数N影响过渡带的宽度,同时也影响通带内波动的疏密,因为N等于通带内的最大值和最小值的总个数。
设阻带的起点频率为?
S,则有
令久,=9由入>I,有
Cv(&)=M(NarchA^)=—
可以解出
arch
列囤g)『
——archN
3dB截止频率用Q,.表示,
忆(血)卜%
按照(2-3-1)式,有
斗丄]
Ine)
经过一系列推论得归一化系统函数为
G4/?
)=
山归一化的系统函数为
Rd」")
才1!
!
"叫)
233切比雪夫低通滤波器的设计步骤
⑴确定低通滤波器的技术指标:
边带频率11p,它们满足
dp—)Olg
丹"(丿吋
匕-lOlg'r
Ha<(jny
(2)求滤波器阶数N和参数;
archki1
先求出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和阻带频率wsc。
wp.ws是实际模拟角频率。
(4)[B,A]=cheby1(N,Rs,wso,'ftype')
该格式用于计算N阶切比雪夫1型数字滤波器系统函数的分子和分母的多项式系数向量B和Ao调用参数N和wso分别为切比雪夫1型数字滤波器的阶数和阻带截止频率的归一化值。
(5)[B,A]=cheby1(N,Rp,wso,'ftype','s')
该格式用于计算N阶切比雪夫1型模拟滤波器系统函数的分子和分母多项式系数向量B和Ao调用参数N和wso分别为切比雪夫1型模拟滤波器的阶数和阻带截止频率的归一化值。
Impinvar
功能:
用“脉冲响应不变法”将模拟滤波器离散化。
[BZ,AZ]=impinvar(B,A,Fs)
Plot(X,Y)功能:
绘制线性二维图形
gridon
功能:
对当前坐标系添加主要网格线。
3脉冲响应不变法脉冲响应不变法原理
脉冲响应不变法是实现模拟滤波器数字化的一种直观而常用的方法。
它特别适合于对滤波器的时域特性有一定要求的场合。
脉冲响应不变法的核心原理是通过对连续函数ha(t)等间隔采样得到离散序
列ha(nT)。
令h(n)=ha(nT),T为采样间隔。
它是一种时域上的转换方法。
一个模拟滤波器的传递函数可以用有理分式表达式表示为:
H')b。
bs…bws"
(3-1)
a。
Qs■…as"
通过反拉普拉斯变换我们就可以得到它的冲激相应:
-X
ha(t)二L[Ha(s)]
(3-2)
脉冲响应不变法就是要保证脉冲响应不变,即:
h(n)小(nT)
(3-3)
对上式的冲激相应序列h(n)作Z变换,就可以得到数字滤波器的传递函数:
H⑵二ZT[ha(nT)]
(3-4)
一般来说,Ha(S)的分母多项式阶次总是大于分子多项式的阶次。
假定Ha(S)
的没有多重极点,则式
(1)就可分解为:
二A
Ha(S)=、
j-s-Si
(3-5)式中Si,A均为复数,Si是Ha(S)的极点。
其拉氏变换为
二Si.Tn
ha(t)二為Aeu(n)
i±
一般来说,Ha(s)的分母多项式阶次总是大于分子多项式的阶次。
数字滤波
器的传递函数H(z)经过合并简化,成为一般形式的有理分式传递函数
H⑵二
c0-C|Z‘-...-qz丄
dod才...-(3-6)
在讨论采样序列z变换与模拟信号拉氏变换之间关系的有关章节中,我们已
经知道
1孑2兀
H(z)|z』Ha(sjK)
TkT
(3-7)
按照z=eST的关系,每一个s平面上宽度为2二/T的水平条带将重迭映射到
z平面上。
因此脉冲响应不变法将S平面映射到z平面,不是一个简单的对
应的关系。
对于高采样频率(T小)的情况,数字滤波器在频域可能有极高的增益。
为此我们采用
H⑵八FA—r
y1—ez
SiT
(3-8)
N
h(n)=Tha(nT)T^eSJnu(n)(3-9)
i4
在脉冲响应不变法设计中,模拟频率与数字频率之间的转换关系是线性的
(.」;*T)。
同时,它可以保持脉冲响应不变,h(t)=ha(nT)。
因此,这一方法往往用于低通时域数字滤波器设计及相应的模拟系统数字仿真设计。
4实验内容
4.1设计步骤
(1)按一定规则将数字滤波器的技术指标转换为模拟低通滤波器的技术指标;
(2)根据转换后的技术指标使用滤波器阶数函数,确定滤波器的最小阶数N和
截止频率Wc
(3)利用最小阶数N产生模拟切比雪夫低通滤波原型
(4)利用冲激响应不变法或双线性不变法把模拟切比雪夫滤波器转换成数字切
比雪夫低通滤波器。
4.2试验程序
低通滤波器程序
fp=100;fs=150;
Rp=0.5;Rs=10;
计算阶数N和截止频率wc
Fs=1000;Wp=2*pi*fp;Ws=2*pi*fs;[N,wc]=cheb1ord(Wp,Ws,Rp,Rs,'s');%N
[z,p,k]=cheb1ap(N,Rp);%[Bap,Aap]=zp2tf(z,p,k);%的系数
[b,a]=lp2lp(Bap,Aap,wc);%[bz,az]=impinvar(b,a,Fs);%换为数字低通滤波器
wc
计算系统函数零极点和增益
将系统函数的零极点转化为系统函数一般形式
将模拟低通滤波器原型,转换为模拟低通滤波器采用脉冲响应不变法,将模拟低通滤波器,转
figure
(1)
计算频率H和频率响应W
[H,W]=freqz(bz,az,Fs);%
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');ylabel('相位/Ao');gridon;
欢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
N
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‘);gridon;
subplot⑵lj2)plot(W/pi,180/pi*unwrap(angle(H)))
xlabeK'相频/Hz,);ylabel('相位/o’);gridon;
wc=
62&3185
»
Figure1
DesktopWindowHelp
FileEditViewInsertTools
□Q爭In|直Q①⑥丨遲I□因|■回
0
-50
2-100
畀-150
200
0010.20.30.4050.&0.708091
相频心n
-250
■J
验证滤波器的正确性程序
figure
(2);
f仁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);title('f2输入信号');
subplot(2,3,3)
x=sin(2*pi*f1*t)+sin(2*pi*f2*t);
plot(t,x);title('叠加后输入信号');
y=filtfilt(bz,az,x);
subplot(2,3,4)
plot(t,y);title('滤波之后的信号')
subplot(2,3,5)
plot(t,sin(2*pi*f1*t));title('希望输出的信号');
入MATLAB
:
ileEditDebugDesktopWindow
Help
21Q
CurrentDirectory:
e:
\MATLAB7\work
Z1』⑥
shortcuts
臣]HowtoAdd回WhafsNew
>>fp=100;fs=150;
^p=0.5;Rs=10;
:
s=1000;Wp=2*pi*fp;Ws=2*pi*fs;
[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);%采用脉冲响应不变法,捋複拟低通慮皮器,转换为埶宁低通痢皮器
figured)
[H,W]=freqz(bz,az,Fs);%计算频率H和频率响应W
mbplot(2,1,1),plot(W/pi,20*logl0(abs(H)))dabelCfl频/Hz,);ylabelC度/dB、);gridon;
subplot(2,1,2),plot(W/pi,180/pi*unwrap(angle(H)))
dabelC相频/Hz,);ylabel(?
相位/o')fgridon;
:
igure
(2);
:
l=50;
i=150;m=0:
(n-1);t=m/Fs;%采样点数
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|
战扎ET号輪AH咎
鸯59彎出的re号
El
|d何辽1
FihtdittfiurInartTookLwEklnpVrindowH»lp
□门口凸h險直杠毎崔□□■□
4.3实验结果分析
求得阶数N=3,通带边界频率wc=628.3125,频率在0-100hz(0-0.2二)完全通过,在100-150hz(0.2二-0.3二)通过幅度递减,150hz(0.3二-后被截止。
验证信号f仁50hz,f2=200hz,观察验证后的波形,可以看出fl通过,而f2被截止,数字低通滤波器的阻带截止频率为150hz,所以设计的数字低通滤波器
符合技术指标
5实验总结与体会
通过对数字信号处理的课程设计,熟悉了MATLA的运行环境,了解更多有关于Matlab软件的知识,初步掌握了MATLA语言在数字信号处理中一些基本库函数的调用和编写基本程序等应用;通过对数字低通滤波器的设计让我熟悉了滤波器设计的一般原理,对滤波器有了一个感性的认识;学会了数字低通滤波器设计的一般步骤;加深了对脉冲响应不变法的理解和认识。
在设计的过程中,我也认识到了自己所学知识的不足。
以前上课都是学一些基本的东西,自以为知识掌握的很熟练,实际只是了解了最表面的东西,好多理论只是不能很好的应用于实践,现在运用学到得的东西做出有实际应用价值的东西,对所学知识点进一步的理解,并进行系统化。
这也让我再次认识到知识是无尽的,只有不断的充实自己、完善自己的知识理论体系,才能够更好的胜任自己以后的工作。
设计过程中知识的不足也让我更加坚定了终身学习的决心。
6参考文献数字信号处理(第三版)高西全丁玉美
MATLA基础及应用(第三版)于润伟朱晓慧
MATLAB辅助现代工程数字信号处理(第二版)李益华