数字butterworth滤波器的设计与实现
- 格式：doc
- 大小：192.39 KB
- 文档页数：18
,.
课程设计任务书
学生姓名： 专业班级：
指导教师： 工作单位：
题 目: 数字Butterworth滤波器的设计
初始条件：
1. Matlab软件；
2. 课程设计辅导资料：“Matlab语言基础及使用入门”、“信号与系统”、“数字信号处理原理与实现”、“Matlab及在电子信息课程中的应用”等；
3. 先修课程：信号与系统、数字信号处理、Matlab应用实践及信号处理类课程等。
要求完成的主要任务：（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）
1. 课程设计时间：1周；
2. 课程设计内容：数字Butterworth滤波器的设计，具体包括：基本Butterworth滤波器的设计，数字高通、带通滤波器的设计，冲激响应不变法和双线性变换法的应用等；
3. 本课程设计统一技术要求：研读辅导资料对应章节，对选定的设计题目进行理论分析，针对具体设计部分的原理分析、建模、必要的推导和可行性分析，画出程序设计框图，编写程序代码（含注释），上机调试运行程序，记录实验结果（含计算结果和图表），并对实验结果进行分析和总结，按要求进行实验演示和答辩等；
4. 课程设计说明书按学校“课程设计工作规范”中的“统一书写格式”撰写，具体包括：
① 目录；
② 与设计题目相关的理论分析、归纳和总结；
③ 与设计内容相关的原理分析、建模、推导、可行性分析；
④ 程序设计框图、程序代码（含注释）、程序运行结果和图表、实验结果分析和总结； ,.
⑤ 课程设计的心得体会（至少500字）；
⑥ 参考文献（近5年的、不少于5篇）；
⑦ 其它必要内容等。
时间安排：1周
附——具体设计内容：
1. 低通巴特沃斯模拟滤波器设计。设计一个低通巴特沃斯模拟滤波器：指标如下：通带截止频率：pf=3400HZ, 通带最大衰减：pR=3dB 阻带截至频率：sf=4000HZ，阻带最小衰减：SA=40dB
2. 模拟低通转换为数字低通滤波器已知一模拟滤波器的系数函数为
1000()1000asSH 分别用冲激响应不变法和双线性变换法将()asH转换成数字滤波器系统函数()Hz，并图示()asH和 ()Hz的幅度相应曲线。分别取采样频率sF＝1000Hz和sF＝500Hz，分析冲激响应法中存在的频率混叠失真和双线性变换法存在的非线性频率失真等。
3.设计一个10阶的Butterworth滤波器，通带为100～250Hz，采样频率为1000Hz，绘出滤波器的单位脉冲响应。
指导教师签名： 年 月 日
系主任（或责任教师）签名： 年 月 日 ,.
目 录
1 设计总体思路 .................................................................................................. 1
2 低通巴特沃斯模拟滤波器设计 .................................................................... 2
2.1 MATLAB提供的设计函数................................................................. 2
2.1.1巴特沃斯滤波器阶数选择函数 .............................................. 2
2.1.2 完整的巴特沃斯滤波器滤波器设计函数 ............................ 3
2.2 设计步骤 .............................................................................................. 3
2.3 MATLAB程序及运行结果................................................................. 4
2.4 结果分析 .............................................................................................. 5
3 模拟低通转换为数字低通滤波器设计 ........................................................ 7
3.1冲激响应不变法 ................................................................................... 7
3.2 双线性变换法 ...................................................................................... 7
3.3 MATLAB提供的设计函数................................................................. 8
3.4 MATLAB程序及运行结果................................................................. 8
3.5 结果分析 ............................................................................................ 10
4 设计10阶巴特沃斯滤波器 ........................................................................ 11
4.1设计原理分析 ..................................................................................... 11
4.2 matlab 设计程序及运行结果 ........................................................ 12
4.3 结果分析 ............................................................................................ 13
5 心得体会 ........................................................................................................ 14
6 参考文献 ........................................................................................................ 15 ,.
1 设计总体思路
IIR数字滤波器的设计，主要采用间接法，即：首先设计出低通模拟滤波器H(S)；进行频率变换，将其转换为高通、带通、带阻滤波器；再用脉冲响应不变法或双线性变换法从模拟滤波器转换为数字滤波器。
对单极点的N阶H(S)用部分展开式：
）－（k1)(SSASHKNK
冲激响应不变法取H(S)的单阶极点Sk的指数函数sFKSe作为H（Z）的极点
）－（111)(zeAzHsFksKNK
双线性变换法是用11112ZZTS代换H（S）中的S得到H(Z),双线性变换法可完全消除频率混叠失真但存在非线性频率失真，而冲激响应不变法存在混叠失真。
在不同的设计阶段MATLAB的信号处理工具箱都给出了相应的滤波器设计函数，这些函数代表了不同类型的逼近函数的滤波器，常用的有巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器类型。
模拟域 冲激响应不变法
频率变换 双线性变换法
图1.1 设计整体思路
归一化模拟低通原型 模拟高通带通带阻滤波器 数字高通带通带阻滤波器 ,.
2 低通巴特沃斯模拟滤波器设计
2.1 MATLAB提供的设计函数
2.1.1巴特沃斯滤波器阶数选择函数
（1）求最小阶数N的函数buttord
调用格式1：[N, Wn] = BUTTORD(Wp, Ws, Rp, Rs, 's')
功能：求出巴特沃斯模拟滤波器的最小阶数N和频率参数Wn
说明：Wp, Ws,为通带截止频率和阻带截止频率如式所示；
sstopsffW2，为时域采样频率其中sspasspfffW,2
Rp, Rs为通带最大衰减 Apass(dB)和阻带最小衰减Astop(dB)如式
)(lg20)d()()(lg20220ppasspassfjfjjeHBeHeHA)(lg20)d()()(lg20220sstopstopfjfjjeHBeHeHA
's'对应模拟滤波器。
（2）模拟低通滤波器设计函数buttap
调用格式：[Z,P,K] = BUTTAP(N)
功能：按给定的阶数N设计出巴特沃斯模拟低通滤波器
说明：Z,P,K为返回的N阶模拟滤波器的零、极点和增益系数。如： ,.
NiiNjjpZZZH11)()z(k)(＝
其中k为增益常数，zj和pi为系统函数H(Z)的N个零点和N个极点。
2.1.2 完整的巴特沃斯滤波器滤波器设计函数
（1）合为一步的数字滤波器设计函数butter
调用格式1：[B,A] = BUTTER(N,Wn)
功能：设计N阶低通数字滤波器，系数矢量B，A按Z的降幂排列，当Wn为二元矢量时，Wn = [W1 W2]，返回一个2N阶带通滤波器，通带范围：W1 < W < W2。
调用格式2：[B,A] = BUTTER(N, Wn, 'high') 设计高通数字滤波器.
调用格式3：[B,A] = BUTTER(N, Wn, 'stop') 设计带阻滤波器 Wn = [W1,W2].
调用格式4：[Z,P,K] = BUTTER(N,Wn) 返回低通数字滤波器的N阶零、极点矢量.
（2）幅频响应FREQS
调用格式： [H,W] = FREQS(B,A,M)
功能：此函数可以求出系统频率响应的数值解，并列出此系统的幅
频及相频响应曲线。
2.2 设计步骤
（1） 确定滤波器设计指标通带内截止频率wp、阻带内截止频率 ws、通带内最大衰减 rp、阻带内最大衰减 rs。
