基于Matlab数字带通滤波设计

《数字信号处理课程设计报告》课程设计题目：数字带通滤波设计

学院：信息工程学院

专业：通信工程

班级：

学生姓名:

指导老师：

日期：2012年5月4日至17日

目录

1. 课程设计的目的和意义 (3)

1.1课程目的 (3)

1.2实验意义 (3)

2.课程设计题目描述及要求 (4)

3.实验内容 (5)

3.1数字滤波器的简介 (5)

3.2 椭圆数字带通滤波器设计原理 (5)

3.3实验流程框图 (7)

3.4实验步骤 (7)

3.5 程序代码 (8)

3.6实验仿真结果图 (9)

3.7实验结果分析 (10)

4.实验总结 (10)

参考文献 (11)

1.课程设计的目的和意义

1.1课程目的

(1) 学习椭圆模拟带通滤波器的设计。

(2) 使用椭圆模拟带通滤波器逼近数字带通滤波器。

(3) 双线性变换法的应用原理。

(4) 使用matlab对滤波器进行仿真和频谱分析。

1.2实验意义

数字滤波器是个离散时间系统（按预定的算法，将输入离散时间信号换为所

要求的输出离散时间信号的特定功能装置）。应用数字滤波器处理模拟信号时，

首先须对输入模拟信号进行限带、抽样和模数转换。数字滤波器输入信号的抽样

率应大于被处理信号带宽的两倍，其频率响应具有以抽样频率为间隔的周期重复特性，且以折叠频率即1/2抽样频率点呈镜像对称。为得到模拟信号，数字滤波器处理的输出数字信号须经数模转换、平滑。数字滤波器具有高精度、高可靠性、可程控改变特性或复用、便于集成等优点。数字滤波器在语言信号处理、图像信号处理、医学生物信号处理以及其他应用领域都得到了广泛应用。

利用Matlab设计滤波器,可以按照设计要求非常方便地调整设计参数,极大地减轻了设计的工作量,有利于滤波器设计的最优化。Matlab因其强大的数据处理功能被广泛应用于工程计算,其丰富的工具箱为工程计算提供了便利,利用Matlab信号处理工具箱可以快速有效地设计各种数字滤波器,设计简单方便。本文介绍了在Matlab环境下滤波器设计的方法和步骤。

2.课程设计题目描述及要求

用双线性变换法设计一个数字带通滤波器，使其指标接近接近如下技术指标

的椭圆模拟带通滤波器：ω

p1=100Hz, ω

p2

=200Hz, ω

s1

=50Hz, ω

s2

=250Hz,通带衰减

系数为Rp=0.5dB,组带衰减系数Rs=60dB,采样频率为20000Hz。

3.实验内容

3.1数字滤波器的简介

数字滤波器(digital filter)是由数字乘法器、加法器和延时单元组成的一种装置。其功能是对输入离散信号的数字代码进行运算处理，以达到改变信号频谱的目的。

数字滤波器可以理解为是一个计算程序或算法，将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列，并在转化过程中，使信号按预定的形式变化。数字滤波器有多种分类，根据数字滤波器冲激响应的时域特征，可将数字滤波器分为两种，即无限长冲激响应（IIR）滤波器和有限长冲激响应（FIR）滤波器。

IIR数字滤波器具有无限宽的冲激响应，与模拟滤波器相匹配。所以IIR滤波器的设计可以采取在模拟滤波器设计的基础上进一步变换的方法。FIR数字滤波器的单位脉冲响应是有限长序列。它的设计问题实质上是确定能满足所要求的转移序列或脉冲响应的常数问题，设计方法主要有窗函数法、频率采样法和等波纹最佳逼近法等。

3.2 椭圆数字带通滤波器设计原理

双线性变换法获得的数字滤波器频率响应特性中不会出现混叠现象，可以适用于高通、带通和带阻滤波器的设计。IIR数字滤波器的设计通常要借助于模拟

低通滤波器的设计，由原型低通滤波器到其他形式（高通、带通、带阻）IIR数字滤波器的频带变换有模拟频带变换法和数字频带变换法。本实验中用的是模拟频带变换法。

首先将给定的对数字滤波器(DF)的技术要求转换为一个低通模拟滤波器(AF)的技术要求，根据这种要求用某种逼近设计出原型的低通模拟滤波器(LP AF)，计算出模拟滤波器的阶数N、极点s i和传递函数)(s

H

a

，再按照双线性变换的变换

关系，将模拟滤波器的传递函数)(s

H

a 转换为数字滤波器的传递函数)(z

H。表1

给出了对数字滤波器(DF)的技术要求直接转换为对一个低通模拟滤波器(AF)的技术要求的频率预畸变校正关系和转换公式。

低通到带通的频率变换公式如下：

（1）

在p平面与s平面虚轴上的频率关系为：

(2)

式中，

，表示带通滤波器的带通宽度，分别为带通滤波器的带通下截止频率和通带截止频率，根据公式（2）的映射关系，频率λ=0映射为频率

，频率λ=。也就是说将低通滤波器G(p)的通带[]映射为带通带通滤波器的通带，同理，得出将Q（p）转换为带通滤波器的系统函数，即：

(3)

可以证明(4)

所以，带通滤波器的通带（阻带）边界频率关于中心频率几何对称，如果原指标给定的边界频率不满足公式（4），就要改变其中一个边界频率，以便满足公式（4），但要保证改变后的指标高于原指标，具体方法是，如果,则减小以使（4）得到满足，具体计算公式为：

或(5)

减小使带通宽度大于原指标要求的通带宽度，增大或减少都使左边的过渡带宽度小于原指标要求的过渡带宽度，反之，如果,，则减少使式（4）得到满足，而且在关于中心频率几何对称的两个正频率点上，带通滤波器的幅度值相等。

3.3实验流程框图

根据要求设计出椭

圆低通滤波器

将低通转化为

带通

使用matlab函数设计椭圆带通滤

波器

运用双线性法计算数

字带通滤波器参数

运用matlab画出数字带通滤波器的频

谱分析图和零极点分布图