目录
1. 数字滤波器的设计任务及要求 (2)
2. 数字滤波器的设计及仿真 (2)
2.1数字滤波器的设计 (3)
2.2数字滤波器的性能分析 (3)
3. 数字滤波器的实现结构对其性能影响的分析 (8)
3.1数字滤波器的实现结构一及其幅频响应 (10)
3.2数字滤波器的实现结构二及其幅频响应 (12)
3.3 数字滤波器的实现结构对其性能影响的小结 (12)
4. 数字滤波器的参数字长对其性能影响的分析 (13)
4.1数字滤波器的实现结构一参数字长及幅频响应特性变化
4.2数字滤波器的实现结构二参数字长及幅频响应特性变化
4.3 数字滤波器的参数字长对其性能影响的小结 (16)
5. 结论及体会 (16)
5.1 滤波器设计、分析结论 (16)
5.2 我的体会 (16)
5.3 展望 (16)
1.数字滤波器的设计任务及要求
1. 设计说明
每位同学抽签得到一个四位数，由该四位数索引下表确定待设计数字滤波器的类型及其设计方法，然后用指定的设计方法完成滤波器设计。
要求：滤波器的设计指标： 低通：
(1)通带截止频率πrad (id)
pc 32
ln =
ω (2)过渡带宽度πrad )
(i d 160
10log tz ≤∆ω
(3)滚降dB αroll 60=
其中，i d — 抽签得到那个四位数（学号的最末四位数），本设计中i d =0201。 2. 滤波器的初始设计通过手工计算完成；
3. 在计算机辅助计算基础上分析滤波器结构对其性能指标的影响（至少选择两种以上合适的滤波器 结构进行分析）；
4. 在计算机辅助计算基础上分析滤波器参数的字长对其性能指标的影响；
5. 以上各项要有理论分析和推导、原程序以及表示计算结果的图表；
6. 课程设计结束时提交设计说明书。
2.数字滤波器的设计及仿真
2.1数字滤波器（编号0201）的设计
数字滤波器是数字信号处理的重要工具之一，它通过数值运算处理改变输入信号所含频率成分的相对比例或者滤出某些频率成分的数字器件或程序，而数字滤波器处理精度高、体积小、稳定、重量轻、灵活、不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊功能。
本次课程设计使用MATLAB 信号处理箱和运用切比雪夫法设计数字滤波器，将手工计算一个切比雪夫I 型的IIR 的低通模拟滤波器的系统函数，并在MATLAB 的FDATool 设计工具分析其性能指标。
其设计流程如图所示：
图1.1
切比雪夫滤波器的幅频特性具有等波纹特性。切比雪夫I 型滤波器振幅特性在通带内是等波纹的，在阻带内是单调递减的。
图1.2
2.2数字滤波器的性能分析
手工计算：
数字低通滤波器设计指标： (1)通带截止频率πrad πrad (id)
pc 1657.032
ln ==ω (2)过渡带宽度πrad πrad )
(i d 014.016010log tz =≤∆ω
（3）阻带截止频率πrad
st 1801.0=ω
(4)滚降dB αroll 60= 2.2.1滤波器计算步骤
（1）数字低通滤波器指标预畸为样本模拟低通滤波器指标 转换关系为：
计算模拟低通滤波器的阶数N ，计算滤波器截止频率
将指标转换为模拟低通滤波器
模拟低通滤波器指标转换成数字低通滤波器
计算数字滤波器指标，并读入指标参数
2
tan 2ω
T Ω=
，其中，令T=2s ,计算得
①通带截止频率573002
tan 2.ωT Ωp
p ==
②阻带截止频率5922.02
tan
2==st
st ωT Ω ③通带允许最大衰减dB R p 1= ④dB As 61=
⑤0335.15730.05922
.0===
p st ΩΩλ
（2）求样本模拟低通切比雪夫滤波器的阶次N
16
1465.15)
(]
1101
10[arcch 1.010≈=--=
p
st
Rp
A .ΩΩarch N s （3）求归一化原型样本模拟低通滤波器的系统函数 ）
（k k N N an p p Πε（s）H -=
=-1
1
21
5088.0110
1.0=-=p
R ε
由公式求出归一化极点，k=1,2，…33.
2*33)12(cos 33)12(sin
π
ξπξ-+--=k jch k ch p k
0893.01
161==
εξarsh
（s）
H an
16
1
1
16)
(2
*5088.01
)(=--=
i i
a p p p G
（4）求样本模拟低通滤波器 （将 去归一化）
0003.00039.00291.01139.01
*
4582.00105.17498.29454.31120.88952.71
*
0013.134512.85428.116028.43332.50020.12716
.0*10|)(234
5678910
111213141516
3++++++++++++++=Ω=-=s s s s s s s s s s s s s s s p G s H p s
p s α
（5）求所需数字低通滤波器的系统函数（双线性变换法）
1
111--+-=
=z z s s （s）
H z H ）（
由于阶数较高，为了简化计算，便用matlab 工具辅助计算，将模拟低通滤波器系统函数系数转换为数字低通滤波器的系统函数。
程序如下：
B1 =1.0e-003 *[ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.2716]
A1 =[ 1.0000 1.0020 5.3332 4.6028 11.5428 8.4512 13.0013 7.8952 8.1120 3.9454 2.7498 1.0105 0.4582 0.1139 0.0291 0.0037 0.0003] [Bz,Az]=bilinear(B1,A1,0.5) [H,f]=freqz(Bz,Az,512,0.5); plot(f,20*log10(abs(H)))
xlabel('频率：HZ)');ylabel('幅度(db)') title('图2.1') axis([0,0.26,-350,5]) grid on
程序运行结果：
（s）H lp （s）H an
