有限长单位脉冲响应滤波器设计
有限长单位脉冲响应(FIR)滤波器的设计方法-第一节
在通信系统中的应用
信号调制与解调
频分复用
FIR滤波器在通信系统的信号调制与解 调过程中起到关键作用，能够实现信 号的频谱搬移。
FIR滤波器可以实现频分复用，将多个 信号调制到不同的频段上，实现多路 信号同时传输。
信道均衡
在通信过程中，信号经过信道时会受 到各种干扰和失真，FIR滤波器可以用 于信道均衡，减小信号失真。
特点
稳定性好、易于实现、无递归结 构、相位线性等。
FIR滤波器的应用领域
01
02
03
信号处理
FIR滤波器广泛应用于信 号去噪、滤波、增强等处 理。
图像处理
在图像处理中，FIR滤波 器用于图像平滑、锐化等 操作。
通信系统
FIR滤波器用于通信系统 的调制、解调、信道均衡 等。
FIR滤波器与IIR滤波器的比较
群时延
群时延
群时延是描述滤波器对信号延迟影响 的参数。在通信和音频处理中，群时 延的稳定性非常重要，在设计FIR 滤波器时考虑使用特定的窗函数或优 化算法，以减小信号通过滤波器时的 延迟。
幅度响应
幅度响应
幅度响应描述了滤波器对不同频率信号的衰减程度。理想的幅度响应应该是对 所有频率成分具有相同的衰减，但在实际中很难实现。
VS
总结词：计算量较小，需要较高的存 储空间和通信开销，适用于大规模数 据和分布式系统。
05
FIR滤波器的应用实例
在音频处理中的应用
音频信号降噪
FIR滤波器能够有效地去除 音频信号中的噪声，提高 音频质量。
音频压缩
通过FIR滤波器对音频信号 进行压缩，可以减小音频 文件的大小，便于存储和 传输。
最优化方法
最优化方法是一种基于数学优化的FIR滤波器设计方法， 其基本思想是通过优化算法来求解FIR滤波器的系数，
课后习题及答案_第7章有限脉冲响应数字滤波器设计--习题（精品pdf）
第6章 有限脉冲响应(FIR)数字滤波器的设计习题1． 已知FIR 滤波器的单位脉冲响应为：（1） h (n )长度N =6h (0)=h (5)=1.5h (1)=h (4)=2h (2)=h (3)=3（2） h (n )长度N =7h (0)=－ h (6)=3h (1)=－ h (5)=－ 2h (2)=－h (4)=1h (3)=0试分别说明它们的幅度特性和相位特性各有什么特点。
2． 已知第一类线性相位FIR 滤波器的单位脉冲响应长度为16， 其16个频域幅度采样值中的前9个为：H g (0)=12， H g (1)=8.34， H g (2)=3.79， H g (3)~H g (8)=0根据第一类线性相位FIR 滤波器幅度特性H g (ω)的特点， 求其余7个频域幅度采样值。
3． 设FIR 滤波器的系统函数为求出该滤波器的单位脉冲响应h (n )， 判断是否具有线性相位， 求出其幅度特性函数和相位特性函数。
4． 用矩形窗设计线性相位低通FIR 滤波器， 要求过渡带宽度不超过π/8 rad 。
希望逼近的理想低通滤波器频率响应函数H d (e j ω)为（1） 求出理想低通滤波器的单位脉冲响应h d (n )；（2） 求出加矩形窗设计的低通FIR 滤波器的单位脉冲响应h (n )表达式， 确定)9.01.29.01(101)(4321−−−−++++=z z z z z Hα与N之间的关系；（3）简述N取奇数或偶数对滤波特性的影响。
5．用矩形窗设计一线性相位高通滤波器，要求过渡带宽度不超过π/10 rad。
希望逼近的理想高通滤波器频率响应函数H d(e jω)为（2）用h1(n)和h2(n)分别构成的低通滤波器是否具有线性相位？群延时为多少？题8图9．对下面的每一种滤波器指标，选择满足FIRDF设计要求的窗函数类型和长度。
（1）阻带衰减为20 dB，过渡带宽度为1 kHz，采样频率为12 kHz；（2）阻带衰减为50 dB，过渡带宽度为2 kHz，采样频率为20 kHz；（3）阻带衰减为50 dB，过渡带宽度为500 Hz，采样频率为5 kHz。
第3章 无限长单位脉冲响应(IIR)滤波器的设计方法
ˆ H a ( s)
[ha (t ) (t nT )]e
n
st
dt
n
h
ha (t ) (t nT )e st dt
nsT
n
a
(nT )e
n
h(n )e nsT
H ( z)
ut e e u nT e u n e u n
nT 3nT
T n 3T n
1 1 H ( z) 1 T 1 z e 1 z 1e 3T
z (e e ) 1 z 1 (e T e 3T ) e 4T z 2
数字滤波器的设计步骤： 1）按照实际需要确定滤波器的性能要求。 2）用一个因果稳定系统的 H(z) 或 h(n) 去逼近这个 性能要求，即求 h(n) 的表达式。 确定系数 a i 、 bi 或零极点 c i 、 d，以使滤波器 i 满足给定的性能要求——第三章、四章讨论 3）用一个有限精度的运算去实现这个系统函数。包 括 选择运算结构：如级联型、并联型、卷积型、频 率采样型以及快速卷积（FFT）型等; 选择合适的字长和有效数字的处理方法等(第五 章)。
这时数字滤波器的频响才能不失真地重现模拟滤 波器的频响（存在于折叠频率ΩS/2以内）
H (e
j
1 ) Ha ( j ) T T
但任何一个实际的模拟滤波器，其频响都不可能 是真正带限的，因此不可避免地存在频谱的交叠， 即混淆，如图，这时，数字滤波器的频响将不同于 原模拟滤波器的频响而带有一定的失真。模拟滤波 器频响在折叠频率以上衰减越大，失真则越小，这 时，采用脉冲响应不变法设计的数字滤波器才能得 到良好的效果。
实验四FIR数字滤波器设计与软件实现
实验四FIR数字滤波器设计与软件实现
实验目的：
掌握FIR数字滤波器的设计与软件实现方法，了解滤波器的概念与基
本原理。
实验原理：
FIR数字滤波器全称为有限脉冲响应数字滤波器，其特点是具有有限
长度的脉冲响应。
滤波器通过一系列加权系数乘以输入信号的延迟值，并
将这些值相加得到输出信号。
FIR滤波器的频率响应由滤波器系数所决定。
实验步骤：
1.确定所需的滤波器的设计规格，包括截止频率、通带波纹、阻带衰
减等。
2.选择适当的滤波器设计方法，如窗函数、最佳近似法、最小二乘法等。
3.根据所选方法，计算滤波器的系数。
4.在MATLAB环境下，使用滤波器的系数实现滤波器。
5.输入所需滤波的信号，经过滤波器进行滤波处理。
6.分析输出的滤波信号，观察滤波效果是否符合设计要求。
实验要求：
1.完成FIR数字滤波器的设计和软件实现。
2.对比不同设计方法得到的滤波器性能差异。
3.分析滤波结果，判断滤波器是否满足设计要求。
实验器材与软件：
1.个人电脑；
2.MATLAB软件。
实验结果：
根据滤波器设计规格和所选的设计方法，得到一组滤波器系数。
通过
将滤波器系数应用于输入信号，得到输出滤波信号。
根据输出信号的频率
响应、通带波纹、阻带衰减等指标，评估滤波器的性能。
实验注意事项：
1.在选择设计方法时，需要根据滤波器要求和实际情况进行合理选择。
2.在滤波器实现过程中，需要注意滤波器系数的计算和应用。
3.在实验过程中，注意信号的选择和滤波结果的评估方法。
第5章 有限长单位脉冲响应(FIR)滤波器的实现
17
4.
h(n)奇对称，h(n)=-h(N-1-n)，N为偶数
j N N 1 1 j 2 2 2
He
e
N 1 2hn sin n n 0 h(n) 2
3 0 1 2 4 5
h(n) 奇对称
图5.1.1 线性相位特性
8
5.1.2 线性相位FIR滤波器的幅度特性
对于线性相位FIR 滤波器h(n)有奇对 称和偶对称两种情 况，且每种情况包 含了N为奇数和偶 数，所以共分四种 情况:
h(n) h(n)
0
1
2
3 4
5
6ห้องสมุดไป่ตู้
n
0
1
2
3
4
5
n
(a) h(n) h(n)
(b)
4 1 2 3
16
N 1 2 H c(n) sin n n 1 所以 N 1 n c ( n) 2h 2
由于 对这些点也奇对称，有
点呈奇对称，所以 。
由于 时， 相当 于H(z)在 处有两个零点，不能用于 的滤波器设计，故不能用作低通、高通和带阻滤波器的设计。
( )
式中为常数，表示此时通过这一系统的各频率分量的时
延； 系统的群时延为:
d ( ) g d
5
线性相位FIR滤波器的DTFT为
H e j h n e jn H e j ( ) H e j
5
6
3
4
5
n
0
1
2
n
(c)
(d)
图5.1.2
线性相位FIR滤波器h(n)的四种对称形式
FIR
N 1 ( ) , 2 2 h(n) h( N 1 n), 0≤ n≤ N 1
（7.1.10）
由以上推导结论可知，如果要求单位脉冲响应为
h(n)、长度为N的FIR数字滤波器具有第二类线性相位 特性，则h(n)应当关于n=(N-1)/2点奇对称。N为奇数 和偶数时h(n)的对称情况分别如表7.1.1中情况3和情况 4
1． 线性相位FIR数字滤波器 对于长度为N的h(n)，频率响应函数为
j
H (e )
n 0
N 1
h(n)e jn
（7.1.1）
（7.1.2）
Hale Waihona Puke H (e j ) H g ()e j ( )
式中，Hg(ω)称为幅度特性; θ(ω)称为相位特性。注意，这
里Hg(ω)不同于|H(ejω)|，Hg(ω)为ω的实函数，可能取负值， 而|H(ejω)|总是正值。线性相位FIR滤波器是指θ(ω)是ω的线 性函数，即
低通滤波器的幅频特性技术指标示意图
ωp到ωs称为过渡带。
通带内和阻带内允许的衰减一般用分贝数表示，通带内允许
的最大衰减用p表示，阻带内允许的最小衰减用s表示。对
低通滤波器， 定义为：
max | H (e j ) | p 20lg dB , 0 p j min | H (e ) |
n 0
N 1
(7.1.6)
将（7.1.6）式中两式相除得到：
cos sin
h(n) cos n h(n)sin n
n0 n0 N 1
N 1
即
h(n) cos n sin h(n)sin n cos
有限脉冲响应滤波器
n
h1[n]
w[n]
h[n]
新的 n
因果脉冲响应
滤波器形状 - 0.0014 dB
10.11 kHz
- 75 dB
21.685 kHz
音频信号的带宽比较小，并且其频率主要集中在低频段，通常是 2 kHz 以下。 而噪声的带宽一般很大。为了提高音频信号的质量，往往需滤除高频的噪音。
Sample value x(n)
bk
FIR 滤波器可以保证稳定， 消除通带内的相位失真等， IIR 滤波器不具备的优点。
IIR滤波器，通常需要的系数较少。
9.2 滑动平均滤波器
滤波器项数越多，低通 滤波效果越强，滤除的 高频分量越多。
虚 部
截止频率约为：
实部
增益为0的最低数字频率：
（ ≈ 455 Hz ）
虚 部
实部
9.3 相位失真
矩形窗得到 的滤波器形状
如果通带增益为 1，21dB的衰减相当于 0.089125的增益。这个值似乎很小，但 一个信号与其相乘后仍然可以听到。 因此，阻带内对信号的衰减仅为 21dB 的滤波器是很差的，较好的滤波器应 具有更大的衰减。
汉宁窗脉冲响应
N=35 汉宁窗幅度响应
汉宁窗
汉宁窗得到 的滤波器形状
1
Amplitude |X(f)|
被 噪 声 污 染 的 音 频 信 号
0
-1 0 0.2 0.4 0.6 0.8 1 1.2 Number of samples 1.4 1.6 1.8 x 10 2
0.5
设计时 Ω1 的选 择 必须大于所要求 的通带边缘频率
截 断 的 脉 冲 响 应
截 断 并 移 位 的 脉 冲 响 应
汉宁窗
数字信号处理课程设计题目：用汉宁窗设计FIR数字低通滤波器对给定数据进行滤波院系：自动化与信息工程学院专业：通信工程班级: 通信 081学号: ************: *******:***职称: 副教授2011年6月26日－2011年6月30日一、设计原理：1.1 FIR 滤波器:有限长单位脉冲响应数字滤波器（Finite Impulse Response Digital Filter,缩写FIRDF ）：有限长单位冲激响应滤波器，是数字信号处理系统中最基本的元件，最大优点是可以实现线性相性滤波，它可以在保证任意幅频特性的同时具有严格的线性相频特性，同时其单位抽样响应是有限长的，因而滤波器是稳定的系统。
因此，FIR 滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。
FIR 滤波器的设计方法主要分为两类：第一类是基于逼近理想滤波器器特性的方法包括窗函数法、频率采样法、和等波纹最佳逼近法；第二类是最优设计法。
设FIRDF 的单位脉冲响应)(n h 的长度为N ，则其频率响应函数为∑-=-=10)()(N n n j j e n h e H ωω（1-1）一般将)(ωj e H 表示成如下形式：)()()(ωθωωj g j e H e H = （1-2）式中，)(ωg H 是ω的实函数（可以去负值）。
与前面的表示形式，即)()()(ωθωωj g j e H e H =相比， )(ωg H 与ω不同。
)(ωθ与 )(ωϕ不同。
为了区别于幅频响应函数)(ωj e H 和相频响应函数)(ωϕ，称)(ωg H 为幅频特性函数，称)(ωθ为相频特性函数。
第一类线性相位FIRDF 的相位特性函数是ω的严格线性函数：()ωτωθ-= （1-3）第二类线性相位FIRDF 的相位特性函数如下：()ωτθωθ-=0 （1-4）式中，τ是常数，0θ是起始相位。
2/0πθ-=在信号处理中很有实用价值（如希伯尔特变换器），这是FIRDF 除了线性相位滤波外，还具有真正交变换作用。
第6章无限脉冲响应滤波器的设计
Ha(
j)
2
1(
1
)2N
c
第六章 无限脉冲响应滤波器的设计
为了求出幅度平方函数，将|Ha(jΩ)|2写成s的函数：
Ha(s)Ha(s) 1(
1 s
)2N
此幅度平方函数有2N个极j点c ，极点
1
j(12k 1)
sk( 1)2N(j c) ce 2 2N
它们均匀的分布在半径为Ωc的圆上，间隔是π/N弧度。
3c
j2
j2
(sce 3 )(sc)(sce 3 )
第六章 无限脉冲响应滤波器的设计
由于各滤波器的幅频特性不同，为使设计统一，
将所有的频率归一化。这里采用对3dB截止频率Ωc归
一化，归一化后的Ha(s)表示为
Ha(s)
N 1
(
1 s
sk
)
k0 c c
(6.2.10)
式中，s/Ωc=jΩ/Ωc。
第六章 无限脉冲响应滤波器的设计
第六章 无限脉冲响应滤波器的设计
第六章 无限脉冲响应滤波器的设计
例6.2.1 已知通带截止频率fp=5kHz，通带最大衰减 αp=2dB ， 阻 带 截 止 频 率 fs=12kHz ， 阻 带 最 小 衰 减 αs=30dB。请按照以上技术指标设计巴特沃斯低通滤波 器。
H a (p ) p 5 3 .2
1 3 p 4 6 5 .213 p 3 6 5 .213 p 2 6 3 .213 p 1 6
系统函数的因式分解形式 1
H a (p ) (p 2 0 .6p 1 1 )8 p (2 1 .6p 1 1 )8 p ( 1 )
p 10 lg
有限脉冲响应数字滤波器设计实验报告
成绩：《数字信号处理》作业与上机实验（第二章）班级：学号：姓名：任课老师：完成时间：信息与通信工程学院2014—2015学年第1 学期第7章有限脉冲响应数字滤波器设计1、教材p238：19.设信号x(t) = s(t) + v(t),其中v(t)是干扰，s(t)与v(t)的频谱不混叠，其幅度谱如题19图所示。
要求设计数字滤波器,将干扰滤除，指标是允许|s(f)|在0≤f≤15 kHz频率范围中幅度失真为±2%（δ1 = 0.02）;f > 20 kHz,衰减大于40 dB（δ2=0.01）;希望分别设计性价比最高的FIR和IIR两种滤波器进行滤除干扰。
请选择合适的滤波器类型和设计方法进行设计，最后比较两种滤波器的幅频特性、相频特性和阶数。
题19图(1)matlab代码：%基于双线性变换法直接设计IIR数字滤波器Fs=80000;fp=15000;fs=20000;rs=40;wp=2*pi*fp/Fs;ws=2*pi*fs/Fs;Rp=-20*log10(1-0.02);As=40;[N1,wp1]=ellipord(wp/pi,ws/pi,Rp,As);[B,A]=ellip(N1,Rp,As,wp1);[Hk,wk1]=freqz(B,A,1000);mag=abs(Hk);pah=angle(Hk);%窗函数法设计FIR 数字滤波器 Bt=ws-wp;alph=0.5842*(rs-21)^0.4+0.07886*(rs-21); N=ceil((rs-8)/2.285/Bt); wc=(wp+ws)/2/pi;hn=fir1(N,wc,kaiser(N+1,alph)); M=1024;Hk=fft(hn,M); k=0:M/2-1;wk=(2*pi/M)*k;%画出各种比较结果图 figure(2);plot(wk/pi,20*log10(abs(Hk(k+1))),':','linewidth',2.5); hold onplot(wk1/pi,20*log10(mag),'linewidth',2); hold offlegend('FIR 滤波器','IIR 滤波器');axis([0,1,-80,5]);xlabel('w/\pi');ylabel('幅度/dB'); title('损耗函数'); figure(3)plot(wk/pi,angle(Hk(k+1))/pi,':','linewidth',2.5); hold onplot(wk1/pi,pah/pi,'linewidth',2); hold offlegend('FIR 滤波器','IIR 滤波器');xlabel('w/\pi');ylabel('相位/\pi'); title('相频特性曲线');(2)两种数字滤波器的损耗函数和相频特性的比较分别如图1、2所示：图1 损耗函数比较图 图2 相频特性比较图0.10.20.30.40.50.60.70.80.91-80-70-60-50-40-30-20-100w/π幅度/d B损耗函数FIR 滤波器IIR 滤波器0.10.20.30.40.50.60.70.80.91-1-0.8-0.6-0.4-0.200.20.40.60.81w/π相位/π相频特性曲线FIR 滤波器IIR 滤波器(3)IIR数字滤波器阶数：N=5FIR数字滤波器阶数：N=36(4)运行结果分析：由图2及阶数可见，IIR阶数低得多，但相位特性存在非线性失真，FIR具有线性相位特性。
第 9 章 有限脉冲响应滤波器
图 9.11a
n 1 sin(n h1[n]= nπ n 1
1)
=
1
π
sinc(n
1)
是滤波器的截止频率(cut-off 1是滤波器的截止频率
frequency)。 。
图9.12 理想低通滤波器脉冲响应
从实用来讲， 从实用来讲，这个脉冲响应要解决两个问题 1) 它的无限长度，两边响应值较小的部分截止 它的无限长度，
9.1 有限脉冲响应滤波器基础 9.1 FINITE IMPULSE RESPONSE FILTER BASICS 有限脉冲滤波器（ 有限脉冲滤波器（Finite impulse response filters) :(非递归滤波器 非递归滤波器nonrecursive filters)滤波 非递归滤波器 滤波 器的输出只与现在和过去的输入有关,跟过去的输出无 器的输出只与现在和过去的输入有关 跟过去的输出无 关。(它总是在有限长个非 0 采样值后，输出为 0 ， 它总是在有限长个非 采样值后， 脉冲响应有限项) 脉冲响应有限项) 差分方程: 差分方程 y[n]= b0+b1x[n-1]+…+bMx[n-M] 脉冲响应： 脉冲响应：h[n]= b0+b1 δ [n-1]+…+bM δ[n-M] 传输函数: 传输函数 H(z)= b0+b1z-1+…+bMz-M 频率响应: 频率响应 H( )= b0+b1e-j +…+bMe-jM
5 项滑动平均滤波器频率响应 H( )=0.2+0.2e-j +0.2e-j2 +0.2e-j3 +0.2e-j4 频率响应见图 9.2
滑动平均滤波器的特点： 滑动平均滤波器的特点： 1) 第一个零点出现在 2π/M 弧度 2) 截止频率（0.707对应的频率 约为第一个零 截止频率（ 对应的频率)约为第一个零 对应的频率 点增益频率的一半 π/M 3)项数越多，低通滤波器效果强，滤除高频分量 项数越多，低通滤波器效果强， 项数越多 多。
数字信号处理实验(吴镇扬)答案4
实验四 有限长单位脉冲响应滤波器设计朱方方 03 通信四班(1) 设计一个线性相位FIR 高通滤波器，通带边界频率为π，阻带边界频率为π，阻带衰减不小于40dB 。
要求给出h(n)的解析式，并用MATLAB 绘出时域波形和幅频特性。
解：(1)求数字边界频率:0.6 , 0.4c r ωπωπ== (2)求理想滤波器的边界频率:0.5n ωπ=(3)求理想单位脉冲响应:[]d sin ()sin[()]()()1n n n n n n h n n παωααπαωαπ⎧---≠⎪⎪-=⎨⎪-=⎪⎩(4) 选择窗函数。
阻带最小衰减为-40dB ，因此选择海明窗(其阻带最小衰减为-44dB)；滤波器的过渡带宽为ππ=π，因此6.210.231 , 152N N N ππα-=⇒=== (5) 求FIR 滤波器的单位脉冲响应h(n):[]31d sin (15)sin[0.5(15)]1cos ()15()()()15(15)115n n n R n n h n w n h n n n ππππ⎧---⎡⎤⎛⎫-⋅⋅≠⎪ ⎪⎢⎥==-⎝⎭⎨⎣⎦⎪=⎩程序：clear;N=31; n=0:N-1;hd=(sin(pi*(n-15))-sin*pi*(n-15)))./(pi*(n-15)); hd(16)=; win=hanning(N); h=win'.*hd;figure; stem(n,h);xlabel('n'); ylabel('h(n)'); grid;title('FIR 高通滤波单位脉冲响应h(n)'); [H,w]=freqz(h,1); H=20*log10(abs(H)); figure;3plot(w/pi,H);axis([0 1 -100 10]);xlabel('\omega/\pi'); ylabel('幅度/dB');grid;title('FIR 高通滤波器，hanning 窗，N=31');51015202530nh (n )FIR 高通滤波器的单位脉冲响应h(n)0.10.20.30.40.50.60.70.80.91-100-90-80-70-60-50-40-30-20-10010ω/π幅度/d BFIR 高通滤波器，hanning 窗，N=31分析：由图知阻带衰减最小值大于40，满足要求。
课后习题及答案_第7章有限脉冲响应数字滤波器设计--习题答案
∑ h ( n )e
n =0
N −1
− j ωm
1 [1 + 0.9e − jω + 2.1e − j2ω + 0.9e − j3ω + e − j4ω ] 10 1 j2ω (e + 0.9e jω + 2.1 + 0.9e − jω + e − j2ω )e − j2ω 10
1
=
1 ( 2.1 + 1.8 cos ω + 2 cos 2ω )e − j2ω 10
其中, a=(N－1)/2=10。 (2) 由 Hd(ejω)求得 hd(n)：
0
|ω | π
π 4
π <| ω | 4
π sin (n − 10) π − 4 1 4 hd (n) = e − jω10 e jωn dω = ∫ −π / 4 2π π(n − 10)
(3) 加窗得到 FIR 滤波器单位脉冲响应 h(n)： · 升余弦窗:
4π N
H(ejπ)=0，
不能实现高
π , 即 N≥40。取 N=41。 10
hd (n) =
1 π H d (e jω )e jω n dω ∫ −π 2π ωc + B 1 −ωc = e − jω a e jω m dω + ∫ω e − jω a e jω n dω ∫ ω − ( + B ) c 2π c
4
和［h2(n)］也可以得到同样的结论。
jθ ( ω ) jω 设 H1 (e ) = FT[h1 (n)] = H 1g (ω )e 1
H 2 (e jω ) = FT[h2 (n)] = H 2g (ω )e jθ 2 (ω )
dsp实验报告-有限脉冲响应滤波器(FIR)实验
实验四．数字信号处理算法实验实验4.1 ：有限脉冲响应滤波器（FIR ）算法实验一．实验目的1．掌握窗函数法设计FIR 滤波器的Matlab 实现，为CCS 提供滤波系数。
2．掌握采用C 语言在VC5509开发板上实现混频信号的FIR 滤波。
二．实验设备计算机，ICETEK-VC5509-A 实验箱及电源。
三．实验原理1. 窗函数法设计FIR 滤波器（详细理论请看《数字信号处理》原理书籍） 本实验要求：设计一个低通滤波器，通带截止频率fp=10kHz ，阻带截止频率fs1=22kHz ，阻带衰减ap=75dB ，采样频率fs=50kHz,计算出滤波系数fHn,并对混频信号（高频+低频正弦波）fIn 进行滤波,得输出波形fOut 。
解：过渡带宽度=fs1-fp=12kHz ；截止频率：f1=fp+(过渡带宽度)/2=16kHz f1对应的数字频率：Ω1=2πf1/fs=0.64π(rad) -理想低通滤波器单位脉冲响应：hd[n]=sin(0.64π(n-a))/(π(n-a)) 其中a=(N-1)/2 (n=0~N-1)-根据阻带衰减要求选择布莱克曼窗，窗函数长度N 为： N=5.98fs/过渡带宽度≈25则窗函数为：w[n]=0.42-0.5cos(2πn/24)+0.08cos(4πn/24) 滤波器脉冲响应为：h[n]=hd[n]w[n] (n=0~N-1) <1>-根据上面各式计算出h[n]。
2. FIR 滤波FIR 滤波器的差分方程为：1()()N i i y n h x n i -==-∑ <2>其中，h i ----滤波器系数；x(n)---滤波器的输入；y(n)--- 滤波输出。
根据公式<1><2>，得本例对应FIR 滤波器的差分方程为： y[n]=-0.001x[n-2]-0.002x[n-3]-0.002x[n-4]+0.01x[n-5]-0.009x[n-6]-0.018x[n-7]-0.049x[n-8]-0.02x[n-9] +0.11x[n-10]+0.28x[n-11]+0.64x[n-12] +0.28x[n-13]-0.11x[n-14]-0.02x[n-15]+0.049x[n-16]-0.018x[n-17]-0.009x[n-18]+0.01x[n-19] -0.002x[n-20]-0.002x[n-21]+0.001x[n-22] (n=0,1,2,...)采用线性缓冲区法（原理见备课笔记）解此差分方程，得FIR 滤波结果y(n)。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
n n
1 2
由此看出：
①由于
奇对称，所以 0,2 对
为偶对称。
也为奇对称，对
H ( ) 0,
②且由于
时，
处必有一零点。
③因此这种情况不能用于设计
时
的滤波器，如高通、带阻滤波器。
3. h(n)奇对称，N为奇数，h(n)=-h(N-1-n)
满足；
（5）无反馈运算，运算误差小。
（6）由于单位脉冲响应为有限长序列，因此可以 采用FFT快速算法实现信号过滤，可以提高效率。
缺点：（1）因为无极点，要获得好的过渡带特性， 需以较高的阶数为代价；
（2）无法利用模拟滤波器的设计结 果，一般无解 析设计 公式，要借助计算机辅助设 计程序完成。
三、FIR滤波器的设计方法
() (3)
利用类似的关系，可以得出新的解答为
hn Nh22N11 n (4)
（4）式是FIR滤波器具有（3）的线性相位的必要且 充分条件，它要求单位冲击响应的h（n）序列以
n (N 1)/ 2
为奇对称中心，此时，时
间延迟h(n)
(N 1)/ 2
等于
h(长n) 度N-1的一半，即
设计任务：选择有限长度的脉冲响 h(n)
应
H，(z得) 到系统函
数
，使幅频特性满足技术指标，
同时是相频特性达到线性相位。
• 设计方法：
• 窗函数法 • 频率采样法
• 切比雪夫等波纹逼近法
§4.1 线性相位FIR数字滤波器的特性
一、 线性相位的条件
1、线性相位的定义
线性相位意味着一个系统的相频特性是频率的 线性函数，即
2
2
2
N 1/ 2
则 H a(n) cos n n0
由于 ，因此 种滤波器的设计。
偶对称 对这些频率也呈偶对称。可以用于各
2．h(n)偶对称，N为偶数 h(n)=h(N-1-
n)
N 1
N 1
2
H e j
2
h n e jn
h N 1 n e j N 1n
n0
n0
N 1
2
分四种情况: 1、h(n)偶对称 、
N为偶数 2、h(n)偶对称
N为奇数
3、h(n)奇对称 N为偶数
4、h(n)奇对称 N为奇数
二、线性相位FIR滤波器的幅度特性 1．h(n) 偶对称，N为奇数
h(n)=h(N-1-n)
H e j H e j
N 1
h(n)e jn
n0
N 3
2
N 1
y(n) ai x(n i)
①
对应的系统函数
i0 N 1
H (z) ai zi
i0
因为它是一种线性时不变系统，可用卷积和形式表示
N 1
③
y(n) h(i)x(n i) i0
比较①、③得：
ai h(i)
N 1
H (z) h(i)zi
i0
二、FIR数字滤波器的特点(与IIR数字滤波器比较)
第四章 有限长单位脉冲响应（FIR） 滤波器的设计方法
序言 §4.1 线性相位FIR数字滤波器的特性 §4.2 窗口设计法（时间窗口法）设计FIR滤波 器§4.3 频率取样法设计FIR滤波器
§4.4 FIR数字滤波器的最优化设计 §4.5 IIR与FIR数字滤器的比较
序言
一、FIR滤波器的表示
FIR数字滤波器的差分方程描述
e
j
N 1 2
N 3
2 n0
2hnc
os
n
N 1 2
h
N 1 2
H
()
h
N 1 2
( N 3)/ 2 n0
2h(n)
cos
n
N21
() N 1
2
令
，则
H () h N 1 (N1)/2 2h( N 1 m) cosm
2 m1
2
令
a(0) h N 1, a(n) 2h N 1 n, n 1,2,, N 1
h
n
e jn
h
N
1
e
j
N 1 2
N 1
h n e jn
n0
2
n N 1
2
N 3
2
hn
e jn e jN 1n
h
N
1 e
j
N 1 2
n0
2
H (e
j
)
j
e
N 1 2
N 3 2
h
n0
n
j n N 1
(e 2
j n N 1
e 2 )
h
N 1 2
() (1)
式中为常数，此时通过这一系统的各频率分
量的时延为一相同的常数，系统的群时延为
g
d() d
2、线性相位的条件
FIR滤波器的DTFT为
N
H e j H e j hne jn
n
式中 H(ω)是正或负的实函数。等式中间和等式右
边的实部与虚部应当各自相等，同样实部与虚部的
比值应当相等:
个抽h(样N 周1)期 。h(在N 1)
满足 2
2
的这h种( N奇对1)称 0情况下，
2
，因
而
2 。这种线性相位情况和前一
种不同之处是，除了产生线性相位外，还有
的固定相移。0Biblioteka 2 022
(N 1)
(N 0.5)
h(n) 偶
对称
h(n) 奇
对称
图1 线性相位特性
二、 线性相位FIR滤波器的幅度特性
N
sin cos
hnsinn
n N
hncosn
n
将上式两边交叉相乘，再将等式右边各项移到左边， 应用三角函数的恒等关系
N
hnsin n
n
满足上式的条件（线性相位的条件）是
N 1
2
(2)
hn hN 1 n ,0 n N 1
（2）式是FIR滤波器具有（1）的线性相位的必要且 充分条件，它要求单位冲击响应的h（n）序列以
n (N 1)/ 2
延迟 h(n)
等于
为偶对称中心，此时，时间 (N 1)/ 2
长度N-1的一半，即
个抽样周期。N为偶数时，延时为整数，N为奇数时，
延时等于整h(n数) 加半个抽样周n期 (。N不管1)N/ 2为奇偶，此
时
都应满足对
轴呈偶
对称。
另外一种情况是，除了上述的线性相位外，还有一附 加的相位，即
h n e jn e j N 1n
n0
j
N 1
N 2
1
e 2 2h
n0
n
cos n
N 2 1
H
N / 21 n0
2h(n)
cos
n
N21
令
，则
H
N / 22h m1
N 2
1
m
cos
m
1 2
或写为：
H
N /2
b(n)
n1
b(n) 2h
cos
N 1 2
优点 ：（1）很容易获得严格的线性相位，避免被处理
的信号 产生相位失真，这一特点在 宽频带信
号处理、阵 列信号处理、数据传输等系统中
非常重要；
（2 ）可得到多带幅频特性；
性问题；
（3 ）极点全部在原点（永远稳定），无稳定
过一
（4 ）任何一个非因果的有限长序列，总可以通
因果性总是
定的延时，转变为因果序列， 所以