Matlab和双线性变换的滤波器设计Word文件下载.docx

上传人:b****3 文档编号:6874329 上传时间:2023-05-07 格式:DOCX 页数:20 大小:210.81KB
下载 相关 举报
Matlab和双线性变换的滤波器设计Word文件下载.docx_第1页
第1页 / 共20页
Matlab和双线性变换的滤波器设计Word文件下载.docx_第2页
第2页 / 共20页
Matlab和双线性变换的滤波器设计Word文件下载.docx_第3页
第3页 / 共20页
Matlab和双线性变换的滤波器设计Word文件下载.docx_第4页
第4页 / 共20页
Matlab和双线性变换的滤波器设计Word文件下载.docx_第5页
第5页 / 共20页
Matlab和双线性变换的滤波器设计Word文件下载.docx_第6页
第6页 / 共20页
Matlab和双线性变换的滤波器设计Word文件下载.docx_第7页
第7页 / 共20页
Matlab和双线性变换的滤波器设计Word文件下载.docx_第8页
第8页 / 共20页
Matlab和双线性变换的滤波器设计Word文件下载.docx_第9页
第9页 / 共20页
Matlab和双线性变换的滤波器设计Word文件下载.docx_第10页
第10页 / 共20页
Matlab和双线性变换的滤波器设计Word文件下载.docx_第11页
第11页 / 共20页
Matlab和双线性变换的滤波器设计Word文件下载.docx_第12页
第12页 / 共20页
Matlab和双线性变换的滤波器设计Word文件下载.docx_第13页
第13页 / 共20页
Matlab和双线性变换的滤波器设计Word文件下载.docx_第14页
第14页 / 共20页
Matlab和双线性变换的滤波器设计Word文件下载.docx_第15页
第15页 / 共20页
Matlab和双线性变换的滤波器设计Word文件下载.docx_第16页
第16页 / 共20页
Matlab和双线性变换的滤波器设计Word文件下载.docx_第17页
第17页 / 共20页
Matlab和双线性变换的滤波器设计Word文件下载.docx_第18页
第18页 / 共20页
Matlab和双线性变换的滤波器设计Word文件下载.docx_第19页
第19页 / 共20页
Matlab和双线性变换的滤波器设计Word文件下载.docx_第20页
第20页 / 共20页
亲,该文档总共20页,全部预览完了,如果喜欢就下载吧!
下载资源
资源描述

Matlab和双线性变换的滤波器设计Word文件下载.docx

《Matlab和双线性变换的滤波器设计Word文件下载.docx》由会员分享,可在线阅读,更多相关《Matlab和双线性变换的滤波器设计Word文件下载.docx(20页珍藏版)》请在冰点文库上搜索。

Matlab和双线性变换的滤波器设计Word文件下载.docx

此报告重点介绍了用双线性不变法设计IIR数字滤波器的基本流程,比较了各种设计方法的优缺点,总结了模拟滤波器的性能特征。

最后以双线性不变法设计了—个高通巴特沃斯FIR数字滤波器,介绍了设计步骤,然后在Matlab环境下进行了仿真与调试,实现了设计目标。

关键字:

数字滤波器巴特沃斯MATLAB双线性变法

Abstract

Matlabfoilesearchandapplicationofdigitalfilteiprovidesanintuitive,efficientandconvenientweapon.Hematiixoperationsbasedonthecalculation,visualization,progiamdesignintegiatedintoanmteiactivewoikingenviiomnent・Iiipaniculai;

Matlabtoolboxtoenableleseaicheisinvanousfieldscanbeeasyandintuitiveapplicationofscientificxeseaichandengineering.

ThisreportintroducedwithemphasisofthebasicflowofdesigningtheIIRdigitfilterbythebilineaipoliticalrefbnn,comparedwitheachkindofdesignmethod'

sgoodandbadpoints,sunmianzedanalogfiltei^speiformancechaiactenstic.FuiallvdesignonebythebilinearpoliticallefbmitopassButtenvoithhighIIRdigitfiltei;

mtioducedthedesignpiocedure,thenhasearnedonthesmiulationandthedebuggingundertheMatlabenviiomnent,hasachievedtheprojectobjective・

Keyword:

digitalfilteiButtenvorthMATLAB

1数字滤波器

1.1数字滤波器的基本概念

数字滤波器是数字信号处理的重要基础,是对信号都是过滤检测与参数估计等处理过程中,它是使用最为广泛的一种线性系统。

数字滤波器处理的对象是经由采样期间将模拟信号转换而得到的数字信号。

数字滤波器是指完成信号滤波处理功能的,用有限精度算法实现的离散时间线性非时变系统。

数字滤波器的输入是一组数字量。

它本身既可以是用数字硬件装配而成的一台用于完成给定运算的专用数字计算机,也可以是将所需的运算编写的程序通过计算机来执行。

数字滤波器具有稳定性高、精度高、灵活性大等突出优点。

随着数字技术的发展,用数字技术实现滤波器的功能愈来愈受到人们的重视,并得到了广泛的应用。

数字滤波器的数学运算通常有两种实现方式。

一种是频域法,即利用FFT快速运算方法对输入信号进行离散傅里叶变换,分析其频谱,然后再根据所希望的频率特性进行滤波,再利用傅里叶反变换得到时域信号。

这种方法具有较好的频域选择特性和灵活性,并且由于信号频率与所希望的频率特性是简单的相乘关系,所以它比计算等价的时域卷积要快得多。

另一种方法是时域法,这种方法是通过离散的抽样数据做差分数学运算来达到滤波目的的。

1.2数字滤波器的分类

数字滤波器按照不同的分类方法,可分为许多种,但总体来讲可以分成两大类。

—类称为经典滤波器.即一般滤波器,特点是输入信号中有用的频率成分和希望滤除的频率成分各占有不同的频带,通过一个合适的选频滤波器达到滤波的目的。

但对于一般滤波器,如果信号和干扰的频带互不重叠,则不能完成对干扰的有效滤除,这是需要采用另一类所谓的现代滤波器•例如维纳滤波器、卡尔曼滤波器、自适应滤波器等最佳滤波器,这些滤波器可按照随机信号内部的一些统计分布规律,从干扰中最佳地提取信号。

从功能上,一般数字滤波器可以分为低通、高通、带通、带阻和全通等,此种分类方法是和模拟滤波器一样的。

1.3数字滤波器的MATLAB实现

MATLAB工具箱未滤波器的设计应用提供了丰富而简便的方法,如函数方法和图形工具方法等、使原来非常繁琐复杂的程序设计变成了简单的函数调用,为滤波器的设计和实现开辟了广阔的天地。

数字滤波器的一般设计过程为:

按照实际需要,确定滤波器的性能要求;

用一个因果的、稳定的离散线性时不变系统,去逼近这一性能指标;

用有限精度的运算实现所设计的系统;

通过模拟,验证所设计的系统是否符合给定性能要求。

2双线性变换法

2.1双线性变换法知识简介

脉冲响应不变法的主要缺点是产生频率响应的混叠失真。

这是因为从S平面到Z平面是多值的映射关系所造成的。

为了克服这一缺点,可以采用非线性频率压缩方法,将整个频率轴上的频率范围压缩到-兀/T〜MT之间,再用z=esT转换到Z平面上。

也就是说,第一步先将整个S平面压缩映射到S1平面的-tt/T〜MT—条横带里;

第二步再通过标准变换关系z=esIT将此横带变换到整个Z平面上去。

这样就使S平面与Z平面建立了一一对应的单值关系,消除了多值变换性,也就消除了频谱混叠现象。

由图2-1看出,在零频率附近,模拟角频率Q与数字频率co之间的变换关系接近于线性关系;

但当Q进一步增加时,co增长得越来越慢,最后当Qtoc时,co终止在折叠频率3=兀处,因而双线性变换就不会出现由于高频部分超过折叠频率而混淆到低频部分去的现象,从而消除了频率混叠现象。

图2-1

2.2双线性变换法设计数字滤波器原理

将S平面jQ轴压缩变换到si平面jG轴上的-MT到k/T一段,可以采用以下变换关系]

I

J

2tailIT这样±

cc=G变换到T

0=Q

变换到0

=2可将上式写成

221111Tj

Tj

e

ee

jQ-GdQ+-=

s

llsj=G,解析延拓到整个s平面和si平面,可得

221111TsTsTsT

再将Is平面通过以下变换关系映射到z平面•即

Tse

Z1=

从而得到S平面和Z平面的单值映射关系为

111-+-=

Z

ZSS

SZ-+=

11

—般来说,为了使模拟滤波器的某一频率与数字滤波器的任一频率有对应关系,可引入待定常数C,

T

sTse

ec

s1111-+-=

将T

se

z

匸代入到上式,可得

111一十-=z

zc

cscz-+=

在MATLAB中,双线性Z变换可以通过bilinear函数实现淇调用格式为:

[B乙Az]=bilmeai(B5A,Fs;

其中B,A为模拟滤波器传递函数G(s的分子分母多项式的系数向量,而Bz,Az为数字滤波器的传递函数H(z的分子分母多项式的系数向量。

3设计任务及方案选择

3.1设计任务与要求

利用MATLAB仿真软件系统结合双线性变换法设计一个数字巴特沃斯高通、低通和带通FIR滤波器。

假设一个信号:

tf

tftx2

12cos5.02sin(兀兀+=

其中fl=100Hz.f2=550Hzo试将该信号分别通过低通、高通、带通滤波器,对输出信号进行分析比较。

3.2方案设计论证与选择

本次滤波器设计的主要方法是先设计低通模拟滤波器,然后转换为高通、带通或带阻数字滤波器。

对于其他如高通、带通,则通过频率变换转换为设计相应的高通、带通等。

在设计的全过程的各个步骤,matlab都提供相应的工具箱函数,使得IIR数字滤波器设计变得非常简单。

总的来说,我的设计思路主要有以下两种:

方案一「从归一化模拟低通原型出发,先在模拟域内经频率变换成为所需类型的模拟滤波器;

然后进行双线性变换,由S域变换到Z域,而得到所需类型的数字滤波器。

图3-1先频率变换再离散

方案二:

先进行双线性变换,将模拟低通原型滤波器变换成数字低通滤波器;

然后在Z域内经数字频率变换为所需类型的数字滤波器。

图3-2先离散再频率变换

以上两种方案都可以,我最后选择了第一种方案进行设计,即先在模拟域内经频率变换成为所需类型的模拟滤波器;

然后进行双线性变换,由S域变换到Z域,而得到所需类型的数字滤波器。

因为高通、带通滤波器的设计方法是先将要设计的滤波器的技术特性指标通过频率转换成模拟低通滤波器的技术指标,再根据这些性能指标设计出低通滤波器的传递函数,为了设计程序的简洁,故选择了方案一。

4程序设计与调试

4.1设计要求

若实验设计要求得巴特沃斯高通滤波器的性能指标为:

通带边界频率fjp:

fp=400Hz;

通带波纹Rp:

Rp=ldB;

阻带边界频率fs:

fs=300Hz;

阻带衰减Rs:

Rs=30dB;

采样频率Fs:

Fs=2000Hz;

要求巴特沃斯低通滤波器通带边界频率与阻带边界频率与高通的正好相反。

模拟带通滤波器中心频率为«

)=200Hz,带宽为Bw=400Hz,

4.2程序设计与说明

4.2.1模拟低通滤波

首先应用双线性变换作为数字低通滤波器到模拟高通滤波器的频率预畸变处理:

2tanT

ppco=Q2cot22ss

pwT

fT

Jico2=

利用上式确定的技术指标和已知的参数Rp、Rs,就可以进行模拟原型低通滤波器设计了。

该部分的程序设计与注释如下:

Fs=2000;

Ts=l/Fs;

%采样频率和采样周期

wp=400*2*pi;

ws=300*2*pi;

%低通、高通数字滤波器截止频率

Rp=l;

Rs=30;

%通带波纹和阻带我减

wpl=Fs*2*tan(ws*Ts/2;

%按频率转换公式进行转换wsl=Fs*2*tan(wp*Ts/2;

[N,Wii]=buttofd(wpl,wsLRp,RsJst;

%计算模拟滤波器的最小阶数

[z,p,k]=buttap(N;

%设计模拟原型低通滤波器

[b,a]=zp2tf(z,p,k;

%变零极点增益形式为系统传递函数形式[H,w]=fieqs(b,a;

%求模拟低通滤波器的频率响应figine;

plot(w,abs(H;

grid;

%绘制低通频率特性皿凤模拟低通原型'

xlabel(W;

ylabelC幅度'

;

422模拟到数字的转换

(1从低通向低通的转换[b,a]=lp21p(Bap,Aap,Wn;

*能是把模拟滤波器原型转换成截至频率为Wn的低通滤波器。

再由得到的模拟低通滤波器通过函数[Bz,Az]=bilineai(B,A,Fs得到数字低通滤波器。

该部分的程序设计如下:

[el,fl]=lp21p(b,a,300*2*pi;

%原型低通转换为低通模拟滤波器[hl,wl]=freqs(el,fl;

%求模拟低通滤波器的频率响应figure;

plot(wl/(2*pi,abs(hl;

%绘制高通频率特性gild;

mleC模拟低通滤波器:

xlabelC频率Hz:

ylabe«

幅度:

[rl,ql]=bilineai(el,fl,1000;

%模拟低通滤波器转换为数字滤波器[sl,w2]=fieqz(il,ql;

%求数字低通滤波器的频率响应figuie;

plot(w2/pi,20*logl0(abs(sl;

giid;

%绘制数字低通滤波器的频率响应

mleC数字低通滤波器:

xlabel('

w/pr;

ylabel('

幅度(dB:

(2从低通向高通的转换[b,a]=lp2hp(Bap、Aap,Wn;

功能是把模拟滤波器原型转换成截至频率为Wn的模拟高通滤波器,再进一步转换成数字高通滤波器。

程序设计如下:

[e2,f2]=lp2hp(b,a,300*2*pi;

%原型低通模拟滤波器转换为高通模拟滤波器[Ii2,w3]=freqs(e2,f2;

%求模拟高通滤波器的频率响应

figure;

plot(w3/(2*pi,abs(h2;

MleC模拟高通滤波器•;

ylabelC幅度

[i2,q2]=bilineai(e2,f2,1000;

%模拟高通滤波器转换为数字滤波器[s2,w4]=fieqz(i2,q2;

%求数字高通滤波器的频率响应figuie;

plot(w4/pi,20*logl0(abs(s2;

%绘制数字高通滤波器的频率响应titleC数字高通滤波器:

xlabel(\v/pr;

ylabel(,幅度(dB:

(3从低通向带通的转换[b,a]=lp2bp(Bap,Aap,W0.Bw;

功能是把模拟滤波器原型转换成中心频率为W0、带宽为Bw的模拟带通滤波器,在进一步转换成数字带通滤波器。

程序设计如下:

W0=200*2*pi;

Bw=200;

[e3,f3]=lp2bp(b,a,W0.Bw;

%原型低通转换为带通模拟滤波器[Ii3,w5]=fieqs(e3,f3;

%求模拟带通滤波器的频率响应figuie;

plot(w5/(2*pi,abs(h3;

%绘制带通频率特性gild;

titleC模拟带通滤波器•;

xlabel('

频率Hz'

幅度*;

[r3,q3]=bilineai(e3,f3,1000;

%模拟滤波器转换为数字滤波器[s3,w6]=fieqz(r3,q3;

%求数字带通滤波器的频率响应figure;

plot(w6/pi,20*logl0(abs(s3;

%绘制数字带通滤波器的频率响应title。

数字带通滤波器'

幅度(dB'

4丄3对输入信号进行滤波

读取输入信号程序]

n=0J80;

dt=l/Fs;

t=n*dt;

%时间序列fl=100;

f2=550;

%输入信号频率x=sin(2*pi*fl*t+cos(2*pi*f2*t+sin(2*pi*1*t;

%输入信号figtue;

plot(t,x;

%绘制输入信号titleC输入信号

xlabelC时间t'

ylabel(*幅度:

将输入信号分别通过得到的低通、高通、带通滤波器,调用函数y=filtei(bz,az,x对输入信号进行滤波。

v=filtei(il,ql,x;

%对输入信号进行低通滤波figure;

plot(t,y;

%绘制输出信号“啦低通滤波输出信号:

ylabel(*幅度'

再进行高通、带通滤波只需将调用函数兔⑹中参量改为12、q2和13、q3即可。

4.3仿真结果与分析

将设计好的程序在matlab软件下进行仿真,得到仿真波形。

4.3.1仿真结果记录

图4-1模拟低通原型幅频响应

(a低通模拟滤波器的频率特性曲线(b巴特沃斯低通数字滤波器的频率特性曲线图4-2

JFigure1

摸拟低通原型

14

06

04

02

Fil«

IdttVie*InsertIdolsIindo*Help

«

■——AM"

*■—]OB■

Dc^lH昌ft題戟口⑥遲□Sh□

0.8

10

 

(a高通模拟滤波器的频率特性曲线(b巴特沃斯高通数字滤波器的频率特性曲线图4-3

(a带通模拟滤波器的频率特性曲线(b巴特沃斯带通数字滤波器的频率特性曲线图4・4

图4-5输入信号波形

JFigure8

FileEditV*iewInsertToolsDesktop¥

indo皆Kelp

G!

寻fe愆®

◎⑥退口匡1盲叵

图4・6经过巴特沃斯低通滤波器滤波后的输出波形

FileEditVie-«

InsertToolsDesktopWindowHelp□rasn致竄裁⑥空□sVq

图4・7经过巴特沃斯高通滤波器滤波后的输出波形

□回冈

Figure11

图4-8经过巴特沃斯高通滤波器滤波后的输出波形

4.3.2仿真结果分析

由仿真结果可知,设计出的数字低通,高通,带通的幅频响应均符合各自的特性,从模拟滤波器特性曲线可以看出,所设计的高通滤波器在大于400Hz为通带,其衰减均小于ldB;

小于300Hz为阻带,其衰减大于30dB,而设计的低通滤波器情况正好相反,符合本次设计的要求。

当滤波器输入100Hz和550Hz两种信号后,低通滤波器可以滤除频率为550Hz的高频信号,而低频信号100Hz的信号可以通过滤波器,验证了该低通滤波器的功能。

而高通滤波器可以滤除频率为100Hz的信号,而高频信号550Hz的信号可以通过滤波器。

4.3.3程序调试中出现的问题

带通滤波器的仿真结果虽然符合其带通特性,但要求得到的性能指标存在较大误差,由于时间紧迫,将在课后进一步寻找原因及解决方案。

—个线性是不变系统的典型表达式有状态空间型、传递函数型、零极增益型、极点留数型共四种,在开始的设计过程中我发现在把零极点增益型向传递函数型的转换过程中多用了一步,即先把零极点增益型转换为状态空间型,进行相应的频率变换和双线性变换后,再由状态空间型转换为传递函数型,进而绘制出其频率响应图,但是经仿真分析后发现,性能并不理想,所以并没有采用。

5心得体会

本次课程设计是基于数字滤波器的MATLAB软件应用,在本学期中,我在数字信号处理课程中学习了关于本次课程设计的一些内容,对于MATLAB软件在数字滤波器的设计中的一些函数也有所认识,通过着一周的时间更是让自己在这个方面的学习更进一步、使自己在这方面的学习深度有所加深。

本次用双线性变换法对巴特沃斯低通,高通及带通滤波器的设计,要求我能够用所学理论知识对巴特沃斯高通滤波器进行设计,并能够对设计的滤波器的频率特性进行分析。

在设计过程中•首先,我通过查找资料等途径,对关于所要设计的题目进行更加全面深入的学习,认真分析其中的原理,弄明白设计的思路;

其次,根据自己选择的性能指标要求,结合自己通过学习掌握的知识和资料进行程序设计,根据自己提出的方案利用软件实现;

最后,运行程序,分析所得结果是否符合自己的设计要求,对结果进行分析,实现由理论到实践,再由实践结果验证理论的过程。

总的来说,这次课程设计让我对MATLAB有了更深刻的了解,对数字滤波器的设计流程有了大致的了解•掌握了一些设计滤波器的基本方法,提高了理论用于实践的能力.掌握了更多专业相关的使用知识与技能。

同时,也暴露了我很多的不足,在以后的学习中,将进一步发扬有点,克服缺点。

参考文献

[1]袁润章,自蔓延高温合成技术研究进展.武汉:

武汉工业大学出版社,1994

[2]SchachtE.Iiidustnalpolysaccharides.Ainsteidam:

ElsevierScience,1987[3]郭景坤•陶瓷材料的强化与增韧新途径的探索•无机材料学报,1998.13910;

23-26

[4]SpnggsGE.Ahistoiyoffinegiamedhaidmetal.LitJofRefiactoiyMetalandHard

Material,1995J3:

241-255

[5]王为民.二硼化钛陶瓷的自蔓延高温合成与加工:

[博士论文].武汉工业大学,1988[6]吴中庆,张基田.S1C颗粒与液体铝的界石层•广州「92秋季中国材料科学研讨会,C-MRS,1992.603-609.

附录

本次课设总程序如下:

%采样频率和采样周期wp=400*2*pi;

%低通、高通数字滤波器截止频率Rp=l;

%通带波纹和阻带衰减wpl=Fs*2*tan(ws*Ts/2;

%按频率转换公式进行转换wsl=Fs*2*tan(wp*Ts/2;

[N,Wn]=buttofd(wpl,wsLRp,RsJst;

%计算模拟滤波器的最小阶数[z,p,k]=buttap(N;

%设计模拟原型低通滤波器[b,a]=zp2tf(z,p,k;

%变零极点增益形式为系统传递函数形式[H,w]=fieqs(b,a;

%绘制低通频率特性titleC模拟低通原型'

xlabel(濒率Hz:

ylabel(,幅度

%求数字低通滤波器的频率响应figure;

%原型低通模拟滤波器转换为高通模拟滤波器[Ii2,w3]=fieqs(e2,f2;

%求模拟高通滤波器的频率响应figiue;

皿凤模拟高通滤波器

[r2,q2]=bilineai(e2,f2,1000;

%求数字高通滤波

展开阅读全文
相关资源
猜你喜欢
相关搜索
资源标签

当前位置:首页 > PPT模板 > 商务科技

copyright@ 2008-2023 冰点文库 网站版权所有

经营许可证编号:鄂ICP备19020893号-2