巴特沃斯切比雪夫滤波器的仿真与实现要点Word文件下载.docx

巴特沃斯切比雪夫滤波器的仿真与实现要点Word文件下载.docx

《巴特沃斯切比雪夫滤波器的仿真与实现要点Word文件下载.docx》由会员分享，可在线阅读，更多相关《巴特沃斯切比雪夫滤波器的仿真与实现要点Word文件下载.docx（26页珍藏版）》请在冰豆网上搜索。

1）完成巴特沃斯、切比雪夫滤波器的设计与仿真；

2）撰写课程设计说明书。

进

度

安

排

起止日期

工作内容

16周周一

任务分配和讨论。

16周周一至17周周一

掌握滤波器的组成与工作原理原理，并进行性能分析

17周周一至17周周三

制定滤波器的设计方案

17周周四至18周周三

完成滤波器的实现与性能仿真

18周周四

调试和分析

18周周五

完成课程设计说明书

参考资料

1、《移动通信基础》杨家玮著电子工业出版社

2、《移动通信》郭梯云著西安电子科技大学出版社

3、《基于MATLAB的移动通信信道建模与仿真》刘岚，

4、《MATLAB通信仿真及应用实例详解》邓华等，人民邮电出版社.2003-9

摘要：

利用MATLAB设计滤波器,可以按照设计要求非常方便地调整设计参数,极大地减轻了设计的工作量,有利于滤波器设计的最优化。

MATLAB因其强大的数据处理功能被广泛应用于工程计算,其丰富的工具箱为工程计算提供了便利,利用MATLAB信号处理工具箱可以快速有效地设计各种数字滤波器,设计简单方便。

本文介绍了在MATLABR2009a环境下滤波器设计的方法和步骤。

关键词：

滤波器，MATLAB

1．前言

1.1MATLAB

MATLAB是美国MathWorks公司开发的一种功能极其强大的高技术计算语言和内容极其丰富的软件库，集数值计算、矩阵运算和信号处理与显示于一身。

该软件最初是由美国教授CleveMoler创立的。

1980年前后，他在教线性代数课程时，发现用其他高级语言编程时极不方便，便构思开发了MATLAB，即矩阵实验室（MatrixLaboratory）。

该软件利用了当时代表数值线性代数领域最高水平的EISPACK和LINPACK两大软件包，并且利用Fortran语言编写了最初的一套交互式软件系统，MATLAB的最初版本便由此产生了。

最初的MATLAB由于语言单一，只能进行矩阵的运算，绘图也只能用原始的描点法，内部函数只有几十个，因此功能十分简单。

1984年该公司推出了第一个MATLAB的商业版，并用C语言作出了全部改写。

现在的MATLAB程序是MathWorks公司用C语言开发的，第一版由steveBangert主持开发编译解释程序，SteveKleiman完成图形功能的设计，JohnLittle和CleveMoler主持开发了各类数学分分析的子模块，撰写用户指南和大部分的M文件。

接着又添加了丰富的图形图像处理、多媒体功能、符号运算和与其它流行软件的接口功能，使MATLAB的功能越来越强大。

MTALAB系统主要由以下五个部分组成：

（1）MATALB语言体系。

MATLAB是高层次的矩阵／数组语言，具有条件控制、函数调用、数据结构、输入输出、面向对象等程序语言特性。

利用它既可以进行小规模端程，完成算法设计和算法实验的基本任务，也可以进行大规模编程，开发复杂的应用程序。

（2）MATLAB工作环境。

这是对MATLAB提供给用户使用的管理功能的总称。

包括管理工作空间中的变量据输入输出的方式和方法，以及开发、调试、管理M文件的各种工具。

（3）图形句相系统。

这是MATLAB图形系统的基础，包括完成2D和3D数据图示、图像处理、动画生成、图形显示等功能的高层MATLAB命令，也包括用户对图形图像等对象进行特性控制的低层MATLAB命令，以及开发GUI应用程序的各种工具。

（4）MATLAB数学函数库。

这是对MATLAB使用的各种数学算法的总称。

包括各种初等函数的算法，也包括矩阵运算、矩阵分析等高层次数学算法。

（5）MATLAB应用程序接口（API）。

这是MATLAB为用户提供的一个函数库，使得用户能够在MATLAB环境中使用C程序或FORTRAN程序，包括从MATLAB中调用于程序（动态链接），读写MAT文件的功能。

除此之外，MATLAB系统还具有如下特点：

（1）具有易学易用的语言体系；

（2）具有交互式的工作环境；

（3）具有多层面的图像处理系统；

（4）具有丰富高效的MATLAB工具箱；

（5）具有便利的程序接口（API）;

（6）应用领域广泛；

（7）嵌入了面向对象编程语言。

1.2滤波器的概念 

滤波器是一种用来消除干扰杂讯的器件，将输入或输出经过过滤而得到纯净的交流电。

您可以通过基本的滤波器积木块----二阶通用滤波器传递函数，推导出最通用的滤波器类型：

低通、带通、高通、帯阻和椭圆型滤波器。

传递函数的参数——f0、d、hHP、hBP和hLP，可用来构造所有类型的滤波器。

转降频率f0为s项开始占支配作用时的频率。

设计者将低于此值的频率看作是低频，而将高于此值的频率看作是高频，并将在此值附近的频率看作是带内频率。

阻尼d用于测量滤波器如何从低频率转变至高频率，它是滤波器趋向振荡的一个指标。

实际阻尼值从0至2变化。

高通系数hHP是对那些高于转降频率的频率起支配作用的分子的系数。

带通系数hBP是对那些在转降频率附近的频率起支配作用的分子的系数。

低通系数hLP是对那些低于转降频率的频率起支配作用的分子的系数。

设计者只需这5个参数即可定义一个滤波器。

1.2.1滤波器的原理

凡是可以使信号中特定的频率成分通过，而极大地衰减或抑制其他频率成分的装置或系统都称之为滤波器，相当于频率“筛子”。

滤波器的功能就是允许某一部分频率的信号顺利的通过，而另外一部分频率的信号则受到较大的抑制，它实质上是一个选频电路。

滤波器中，把信号能够通过的频率范围，称为通频带或通带；

反之，信号受到很大衰减或完全被抑制的频率范围称为阻带；

通带和阻带之间的分界频率称为截止频率；

理想滤波器在通带内的电压增益为常数，在阻带内的电压增益为零；

实际滤波器的通带和阻带之间存在一定频率范围的过渡带。

1.2.2理想滤波器与实际滤波器

理想滤波器使通带内信号的幅值和相位都不失真,阻喧内的频率成分都衰减为零的滤波器，其通带和阻带之间有明显的分界线。

如理想低通滤波器的频率响应函数或理想滤波器实际上并不存在。

实际滤波器的特性需要以下参数描述：

1）恒部平均值A0：

描述通带内的幅频特性；

波纹幅度：

d。

2）上、下截止频率：

以幅频特性值为A0/2时的相应频率值WC1，WC2作为带通滤波器的上、下截止频率。

带宽。

因为所以也称“-3dB”带宽

3）选择性：

实际滤波器过渡带幅频曲线的倾斜程度表达了滤波器对通带外频率成分的衰减能力，用信频程选择性和滤波器因素描述。

信频程选择性：

与上、下截止频率处相比，频率变化一倍频程时幅频特性的衰减量，即倍频程选择性：

=或=信频程选择性总是小于等于零，显然，计算信量的衰减量越大，选择性越好。

滤波器因素：

－60dB处的带宽与－3dB处的带宽之比值，即越小,选择性越好.分辨力：

即分离信号中相邻频率成分的能力，用品质因素Q描述。

Q越大，分辨率越高。

实际带通滤波器的形式

1）恒定带宽带通滤波器:

B=常量，与中心频率f0无关。

2）恒定百分比带通滤波器：

在高频区恒定百分比带通滤波器的分辨率比恒定带宽带通滤波器差。

1.2.3滤波器的分类

从大的方面分，滤波器分为模拟滤波器和数字滤波器。

模拟滤波器有电阻，电容，电感，及由原器件构成；

实际中数字滤波器应用的比较广泛。

从实现方法上分，数字滤波器分为IIR和FIR，即无限冲激响应滤波器和有限冲激响应滤波器；

其中IIR网络中有反馈回路，FIR网络中没有反馈回路。

从小的方面分：

1）按所处理的信号分为模拟滤波器和数字滤波器两种。

2）按所通过信号的频段分为低通、高通、带通和带阻滤波器四种。

低通滤波器：

它允许信号中的低频或直流分量通过，抑制高频分量或干扰和噪声。

高通滤波器：

它允许信号中的高频分量通过，抑制低频或直流分量。

带通滤波器：

它允许一定频段的信号通过，抑制低于或高于该频段的信号、干扰和噪声。

带阻滤波器：

它抑制一定频段内的信号，允许该频段以外的信号通过。

1）巴特沃斯滤波器

巴特沃斯滤波器是电子滤波器的一种，巴特滤波器的特点是通频带的频率响应曲线最平滑，没有起伏，而础在阻带则逐渐下降为零。

在振幅的对数对角频率的波特图上，从某一边界角频率开始，振幅随着角频率的增加而逐步减少，趋向负无穷大。

2）切比雪夫滤波器

切比雪夫滤波器是在通带或阻带上频率响应幅度等波纹波动的滤波器。

切比雪夫滤波器在过渡带比巴特沃斯滤波器的衰减快，但频率响应的幅度特性不如后者平坦。

切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小，但是在通频带内存在幅度波动。

I型切比雪夫滤波器在通带上频率响应幅度等波纹的滤波器称为“I型切比雪夫滤波器”，II型切比雪夫滤波器在阻带上频率响应幅度等波纹的滤波器称为“II型切比雪夫滤波器”.

2.设计目的

1）掌握滤波器的组成与工作原理原理，学会MATLAB编程仿真，并进行性能分析。

2）巩固移动通信所学知识。

3.设计原理

3.1.模拟低通滤波器的设计指标及逼近方法

模拟低通滤波器的设计指标有αp,Ωp,αs和Ωs。

其中Ωp和Ωs分别称为通带截止频率和阻带截止频率，αp是通带Ω（=0~Ωp）中的最大衰减系数，αs是阻带Ω≥Ωs的最小衰减系数，αp和αs一般用dB数表示。

对于单调下降的幅度特性，可表示成：

（式3-1）

（式3-2）

如果Ω=0处幅度已归一化到1，即|Ha（j0）|=1,αp和αs表示为

（式3-3）

（式3-4）

以上技术指标用图3-1表示。

图中Ωc称为3dB截止频率，因

（式3-5）

图3-1低通滤波器的幅度特性

滤波器的技术指标给定后，需要设计一个传输函数Ha（s），希望其幅度平方函数满足给定的指标αp和αs，一般滤波器的单位冲激响应为实数，因此

Ha（s）必须是因果稳定系统，因此其所有极点必须落在s平面的左半平面。

3.2.巴特沃斯低通滤波器的设计方法

巴特沃斯低通滤波器的幅度平方函数|Ha（jΩ）|2用下式表示：

（式3-6）

图3-2巴特沃斯幅度特性和N的关系

用s替换j,将幅度平方函数|Ha（jΩs）|2写成s的函数：

（式3-7）

此式表明幅度平方函数有2N个极点，极点sk用下式表示：

（式3-8）

为形成稳定的滤波器，2N个极点中只取s平面左半平面的N个极点构成Ha（s），而右半平面的N个极点构成Ha（-s）。

Ha（s）的表示式为

（式3-9）

图3-3三阶巴特沃斯滤波器极点分布

设N=3，极点有6个，它们分别为

取s平面左半平面的极点s0,s1,s2组成Ha（s）：

由于各滤波器的幅频特性不同，为使设计统一，将所有的频率归一化。

这里采用对3dB截止频率Ωc归一化，归一化后的Ha（s）表示为

（式3-10）

令，p＝s/Ωc=+jλ。

令λ=Ω/Ωc，λ称为归一化频率；

p称为归一化复变量，这样归一化巴特沃斯的传输函数为

（式3-11）

式中，pk为归一化极点，用下式表示：

（式3-12）

将极点表示式（3-12）代入（3-11）式，得到的Ha（p）的分母是p的N阶多项式，用下式表示：

（