基于切比雪夫滤波器设计及matlab仿真Word下载.docx
- 文档编号:21746776
- 上传时间:2023-02-01
- 格式:DOCX
- 页数:10
- 大小:61.38KB
基于切比雪夫滤波器设计及matlab仿真Word下载.docx
《基于切比雪夫滤波器设计及matlab仿真Word下载.docx》由会员分享,可在线阅读,更多相关《基于切比雪夫滤波器设计及matlab仿真Word下载.docx(10页珍藏版)》请在冰豆网上搜索。
计并按要求编写课程设计报告书。
1.3
进一步学习低通滤波器的原理,在通信系统仿真软件MATLAB平台上,设计出IIR模拟低通滤波器,对设计项目进行调试,对程序进行仿真,对结果结合
理论进行分析。
2IIR2.1
IIR滤波器具有无限长脉冲响应,因此能够与模拟滤波器相匹敌;
一般来说,
所有的模拟滤波器都有无限长脉冲响应。
因此,IIR滤波器设计的基本方法是利用复值映射将大家熟知的模拟滤波器变换为数字滤波器。
这一方法的优势在于各
种模拟滤波器设计(AFD)表格和映射在文献中普遍能够获得。
这个基本方法也
称为A/D(模拟-数字)滤波器变换。
然而,AFD表格仅对低通滤波器适用,而
同时也想要设计其他频率选择滤性波器。
为此,需要对低通滤波器实行频带变换,
这些变换也是复值映射,在各种文献中也能得到。
这种IIR滤波器设计的基本方法存在两种主要途径。
途径一:
1.设计模拟低通滤波器,2.实行频带变换,3.实
行滤波器变换。
途径二:
1.设计模拟低通滤波器,2.实行滤波器变换,3.实行频带变换。
模拟滤波器的设计指标以及在模拟滤波器的技术要求中所用到的幅度平方
响应,这将导致三种广泛采用的模拟滤波器特性,即巴特沃兹(Butterworth),切比雪夫(Chebyshev)和椭圆(Elliptic)滤波器。
然后再将这些原型模拟滤波
器转换到不同的频率选择性数字滤波器。
2.2
H()j,设是某个模拟滤波器的频率响应,那么低通滤波器在幅度平方响应a
上的技术指标给出为:
12,,,,|H()|1j,||,,pa2,1,
120|H()|,,,j,,,||,sa2A
,式中是通带波纹参数,是通带截止频率一rad/s(弧度/秒)计,A是阻带衰p
减参数,以及是阻带截止频率以rad/s计。
这些参数如图2.1所示。
s
图2.1模拟低通滤波器技术指标要求
由图可见,
2|H()|j,必须满足a
12,,,|H()|j,,,pap21,,
12|H()|j,,,,,,sap2A
参数,RA和A是分别与以dB计的参数和有关的,这些关系是ps
Rp110R,,,,,,10log101(2-1)p102,,1
As120AA,,,,10log10(2-2)10s2A
,,波纹和的绝对标尺是通过下式与和A有关的:
12
2,1,,111,,,,211,,,,1,,21
,1,121,,,A1,A,,12
3I3.1
切比雪夫I型滤波器的振幅平方函数为:
122A()=H(),,,ja,,,221,v,,,N,,,C
式中,,,,为有效带通截止频率,是与带通波纹有关的参量,大,波纹大,0<
<
1;
C
为N阶切比雪夫多项式。
(),N
有关参数的确定:
(1)预先给定带通截止频率,;
(2)确定,。
通带波纹表示成:
22HjHj()(),,1aamaxmax,,,10lg20lg20lg,221HjHj()(),,aa2minmin1,,
所以
220.1,,,,,,,,10lg
(1),101
2给定通带波纹值,,()dB分贝数后,可求得。
(3)由阻带的边界条件确定阶数N.
2设,,,为事先给定的边界条件,即在阻带中的频率点处,要求滤波器频响Ass
1衰减达到以上。
2A
21即,,,Hj(),,,2saA
11由此得,2A,,,221,,v,,N,C,,
2A,1因此,,,,,v,,N,,C,,
X,,1,()1,,cosh(arccos)Nhx,=N
2A,1arccos()h,N,,sarccos()h,c
因此,要求阻带边界频率出衰减越大,,也越大。
参数,,给定后,查,NNc
阅有关模拟滤波器手册,就可以求得系统函数。
Hs()a
MATLAB提供了[Z,P,K]=cheblap(N,Rp)函数来设计N阶通带波纹为Rp的归一化切比雪夫I型模拟原型滤波器,并得到数组Z和P(即零点和极点)以及增益值K。
我们需要一个具有任意,的非归一化的切比雪夫I型滤波c
器。
这可以通过将这个归一化滤波器的数组P倍乘以,来完成。
和巴特沃兹原c
型滤波器相类似,这个滤波器没有零点。
新的增益K利用式子:
Hj(0)=1,Na
1为奇数或(0)Hj,,N为偶数确定这可通过将原来的K倍乘以非归一化a21,,
分母多项式对归一化分母多项式的比在s=0的求值来完成。
下面的函数称为
U-chblap(N,Rp,Omegac),用于设计一个非归一化的切比雪夫I模拟原型滤波器,得到直接的Hs()。
a
function[b,a]=uchblap(N,Rp,Omegac);
%[b,a]=uchblap(N,Rp,Omegac),
%b=numeratorpolynomialcoefficients
%a=denominatorpolynomialcoefficients
%N=orderoftheellipticfilter
%Rp=passbandroppleindB;
Rp>
0
%Omegac=cutofffrequencyinradians/sec
%[z,p,k]=cheblap(N,Rp);
%a=real(poly(p));
%aNn=a(N+1);
%p=p*Omegac;
%aNu=a(N+1);
%k=k*aNu/aNn;
%b0=k;
%B=real(poly(z));
%b=k*B;
3.2
R已知,,和,有三个参数就可以确定一个切比雪夫I型滤波器,A,spps
根据相对线性标尺中的2-1与2-2式我们可得到:
RAps1020,,,101和A,10根据上面以前讨论的性质有:
s,,r,=和,pc,p阶N给出为
2A,1g,2,
2,,,,log1gg,,10,,,,N,2,,,,10log1,,,,rr,,,,,,
44.1
设计一个低通的切比雪夫I型低通滤波器满足:
通带截止频率:
R=0.2π;
通带波纹:
=1dBpp
0.1*1阻带截止频率:
A=0.3π;
阻带波纹:
=16dB,,,101s我们首先计算必要的参数:
160.1*12010=0.5088,A==6.3096,,,101
0.3,,,c==0.2π,==1.5,rp0.2,
2A,1=12.2429,N=4g,2,
现在能确定Hs(),a
111,,=+=0.417022,,
N1a=0.5(N,,)=0.3646,
N1Nb=0.5(,,)=1.0644,
有四个极点:
,,,,,,,=(a)(b)=-0.0877j0.6179Pcos,,sin,,c,c,0.3,,,,2828,,,,
,3,,3,,,,cos,sin,P=(a)(b)=-0.2117j0.2559,c,c,,1.2,,,,2828,,,,
因此,
K0.891250.11030.3895,,Hs(),=a322(0.17540.3895)(0.42340.1103)ssss,,,,sp(),k,k0,
分子是要使得有:
1(0)Hj,=0.89125a21,,
4.2MATLAB
已知滤波器的技术要求,利用U-chblap函数,我们提供了一个称为afd-chb1
的函数用于设计一个模拟切比雪夫I型低通滤波器。
这个函数给出如后面的目录
所示。
为了展示模拟滤波器的频域图,我们提供的一个称为freqs_m的函数,它是
有MATLAB提供的函数freqs的修正形式。
这个函数以绝对值和相对值dB尺度
计算幅度响应以及相位响应,在freqs_m函数中响应是一直要计算到最大频率,。
freqs_m函数见附录。
max
这个模拟滤波器的脉冲响应h(t)是用MATLAB的impulde函数计算出来的。
通过MATLAB实现的系统函数为
0.891250.11030.3895,,Hs(),a22(0.17540.3895)(0.42340.1103)ssss,,,,
这个H,c(s)与计算的结果稍有不同,这是因为在计算中用了=0.5,而在asfd_butta
函数中,,是按满足在的要求选定的,有关这个滤波器所画出的图如下所示。
pc
图4-1幅度响应
图4-2幅度
图4-3相位响应
图4-4脉冲响应
5
在这次课程设计中,通过查找相关书籍和相关资料,使得自己增长了不少相关知识。
了解到了滤波器分为两大类,确定信号滤波器和随机信号滤波器。
通常
我们用滤波器是来降低噪声,估计信号的,不过滤波器也可以用来处理一些其它
的问题,如图象压缩等。
滤波器本质上是一个输入输出系统,可以是线性的,也
可以是非线性的。
确定信号滤波器是指在滤波器设计时没有用到统计理论,不考
虑信号与噪声的统计特性。
象常见的低通,高通,带通滤波器等都属于确定信号
滤波器,它们针对确定的信号来进行处理。
这三种频域滤波器设计的原理是根据
噪声与有用信号的频率分布是基本分开的,通过确定一个频率开关值,将噪声的
频率去掉,再进行傅立叶反变换就得到有用信号。
上面所提到的低通,高通,带
通滤波器仅适用于有用信号与干扰信号在频谱上具有分离特性的场合,实际中有
另外一些情况,比如干扰是由于信号本身的回声所引起的,这时回声干扰除了振
幅和相位的差别外,基本上跟原信号是一样的,因此频率滤波器无法消除这种干
扰。
但是我们可以遵循广泛适用的最小平方准则设计最小平方滤波器,来达到消
除干扰的目的。
通过这次课程设计,加强了自己掌握和理解书本知识的能力,培养了自己的实际动手能力与综合设计能力,并提高了自己的技术素质。
基本达到了了解信号
处理课程设计的任务,明确了滤波器设计的基本原则,掌握了滤波器设计的基本
方法与任务。
加深了自己对数字信号处理这门课程的理解。
希望自己在以后的学
习生活中不断加强自我学习的能力,努力完善自己。
[1](美)恩格尔,(美)普罗克斯著,刘树棠译.数字信号处理.陕西:
西安交通大学出版社,2002
[2]韩利竹,王华.MATLAB电子仿真与应用.北京:
国防工业出版社,2003[3]罗军辉.MATLAB7.0在数字信号处理中的应用.北京:
机械工业出版社,2005[4]孙祥,徐流美,吴清.MTLAB7.0基础教程.北京:
清华大学出版社,2006[5]聂祥飞,王宝海,谭泽富.MATLAB程序设计及其在信号处理中的应用.成都:
西南交通大学出版社,2005
附录:
源程序清单
Wp=0.2*pi;
Ws=0.3*pi;
Rp=1;
As=16;
Ripple=10^(-Rp/20);
Attn=10^(-As/20);
%模拟滤波器设计:
[b,a]=afd_chb1(Wp,Ws,Rp,As);
%****切比雪夫-I型滤波器阶次=4
%计算二阶环节部分:
[C,B,A]=sdir2cas(b,a)
%%C=0.038
%%B=001
%%A=1.00000.42330.1103
%%1.00000.17530.3895
%计算频率响应:
[db,mag,pha,w]=freqs_m(b,a,0.5*pi);
%计算脉冲响应:
[ha,t]=impz(b,a);
%Plots
subplot(2,2,1);
plot(w/pi,mag);
title('
幅度响应'
)
xlabel('
模拟频率(单位:
pi)'
);
ylabel('
|H|'
axis([0,0.5,0,1.1])set(gca,'
XTickMode'
'
manual'
XTick'
[0,0.2,0.3,0.5]);
set(gca,'
YTickmode'
YTick'
[0,Attn,Ripple,1]);
grid
subplot(2,2,2);
plot(w/pi,db);
幅度(db)'
)xlabel('
分贝数'
axis([0,0.5,-30,5])
[-30,-As,-Rp,0]);
YTickLabelMode'
YTickLabels'
[30,'
16'
1'
0'
])
subplot(2,2,3);
plot(w/pi,pha/pi);
相位响应'
弧度'
axis([0,0.5,-1,1])
[-1,-0.5,0,0.5,1]);
subplot(2,2,4);
plot(t,ha,[0,max(t)],[0,0]);
脉冲响应'
时间秒'
ha(t)'
axis([0,max(t),min(ha),max(ha)])
AFD_CHEB1.M
Function[b,a]=afd_cheb1(Wp,Ws,Rp,Rs)
%切比雪夫I型模拟低通滤波器的设计
%b为分子多项式的系数Hs()a
%a为分母多项式的系数Hs()a
%Wp为以弧度/秒为单位的通带边缘频率;
Wp>
0%Ws为以弧度/秒为单位的阻带边缘频率;
Ws>
0%Rp为通带中的振幅波动的dB数
%As为阻带衰减的dB数
%ifWp<
=0
Error(‘通带必须大于0’)
end
ifWs<
=Wp
Error(‘阻带边缘必须大于通带边缘’)end
if(Rp<
=0)|(As<
0)
Error(‘通带波纹或阻带衰减必须大于0’)end
%ep=sqrt(10^(Rp/10)-1);
A=10^(As/20);
OmegaC=Wp;
OmegaR=Ws/Wp;
g=sqrt(A*A-1)/ep;
N=ceil(log10(s+sqrt(g*g-1))/log10(OmegaR+sqrt(OmegaR*OmegaR-1)));
Fprintf(‘\n***切比雪夫I型滤波器的阶次=%4.0f\n’,N)
[b,a]=cheblap_o(N,Rp,Rp,Rs)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 切比雪夫滤波器 设计 matlab 仿真