第7章 有限长单位抽样响应数字滤波器的设计

合集下载

《数字信号处理教程》程佩青(第三版)清华大学出版社课后答案

《数字信号处理教程》程佩青(第三版)清华大学出版社课后答案

结果 y (n ) 中变量是 n ,


∑ ∑ y (n ) =
x ( m )h (n − m ) =
h(m)x(n − m) ;
m = −∞
m = −∞
②分为四步 (1)翻褶( -m ),(2)移位( n ),(3)相乘,
(4)相加,求得一个 n 的 y(n) 值 ,如此可求得所有 n 值的 y(n) ;
10
T [ax1(n)+ bx2 (n)] =
n

[ax1
(n
)
+
bx2
(n
)]
m = −∞
T[ax1(n) + bx2(n)] = ay1(n) + by2(n)
∴ 系统是线性系统
解:(2) y(n) =
[x(n )] 2
y1(n)
= T [x1(n)] = [x1(n)] 2
y2 (n) = T [x2 (n)] = [x2 (n)] 2
(3) y(n) = δ (n − 2) * 0.5n R3(n) = 0.5n−2 R3(n − 2) (4) x(n) = 2n u(−n −1) h(n) = 0.5n u(n)
当n ≥ 0 当n ≤ −1
∑ y(n) = −1 0.5n−m 2m = 1 ⋅ 2−n
m = −∞
3
y(n) = ∑n 0.5n−m 2m = 4 ⋅ 2n
+ 1)

x1 (n
+ 1)]
=
−a n
综上 i) , ii) 可知: y1 (n) = −a nu(−n − 1)
(b) 设 x(n) = δ (n − 1)
i)向 n > 0 处递推 ,

课后习题及答案_第7章有限脉冲响应数字滤波器设计--习题(精品pdf)

课后习题及答案_第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。

第七章_有限长单位冲激响应(FIR)数字滤波器的设计

第七章_有限长单位冲激响应(FIR)数字滤波器的设计
j e H d ( e j ) 0
| | c
c | |
表示其群时延
2.冲激响应序列
1 j jn hd ( n) H ( e ) e d d 2 1 c j jn e e d c 2 s in[(n ) c ] (n ) c s in[(n ) c ] c (n )
• H (0) 0 ,
0,2 奇对称,关于 偶对称. • H ( ) 关于
可用于设计:
•高通滤波器 •带通滤波器 4种不同的幅度特性中,以第一种幅度特性最好,因而在FIR滤 波器的设计中,通常都采用第一种幅度特性。称第一类FIR滤波 器。
h( n) 偶对称,N为奇数
7.1.3 线性相位FIR数字滤波器的零点分布特点
且 h( n) 关于 N 1 偶对称或奇对称
2 偶对称 h(n) h( N n 1)
奇对称 h(n) h( N n 1)
1. h( n)为偶对称
H ( z ) h( n) z
n 0 N 1 n N 1 n 0
h( N n 1) z n
WR ( ) 为矩形窗频率响应幅度函数 8 主瓣宽度 N 第一旁瓣比主瓣低 31dB
4 N
0
WHan (Biblioteka )4 N4.海明窗 5.布莱克曼窗

n h ( N n 1 ) z n 0
N 1

m N n 1
( N m 1) ( N 1) m ( N 1) 1 h ( m ) z z h ( m ) z z H ( z ) m 0 N 1 N 1
H ( z)
m 0

7有限长单位冲激晌应(FIR)数字滤波器的设计方法讲解

7有限长单位冲激晌应(FIR)数字滤波器的设计方法讲解
(7-8) 式和 (7-12 〉式 , 即h(n)= 士 h(N-1-n)
因而系统函数可表示为
在这一公式中,方括号内有“±”号。当取“十”号时,h(n) 满 足 h(n)=h(N-1-n) 偶 对 称 ; 当 取 “ 一 ” 号 时 ,h(n) 满 足 h(n)=-h(N-1-n)奇对称。下面对应这两种情况分别讨论它们的 频率响应。
的这种奇对称情况下,满足
h
N 1 2
h
N 1种线性相位情况和前一种不同之处是
,

了产生线性相位外 , 还有士π /2 的固定相移。
由于 h(n) 有上述奇对称和偶对称两种 , 而 h(n)的点数 N 又有奇数、偶数两种情况 . 因而 h(n)可以有 4 种类型 , 如图 7-1和图 7-2所示 , 分别对应于4种线性相位FIR数字滤波器。
2. h(n) 奇对称 此时 , 由 (7-15) 式可知 , 频率响应为
图 7-3 h(n)偶对称时的 线性相位特性
图7-4 h(n)奇对称时的90° 相移线性相位特性
将此式与 (7-13) 式比较 , 可得幅度函数为
相位函数为 幅度函数可正可负,相位函数既是线性相位的,又包 括π/2相移, 如图7-4所示,看出不仅有(N-1)/2个抽 样间隔的延时,而且还产生一个90°的相移。这种使 所有频率的相移皆为90°的网络,称为90°移相器, 或称正交变换网络,它和理想低通滤波器、理想微分 器一样,有着极重要的理论和实际意义。
两式相除,可得
(7-8 〉式是 FIR 滤波器具有 (7-3a) 式的线性相 位的必要且充分条件 , 它要求单位冲激 响应的 h(n) 序列以 n=(N-1)/2 为偶对称中心 , 此时时间延时τ等 于以 h(n)长度 N-1 的 一半 , 即为τ =(N-1)/2 个抽样 周期。 N 为偶数时 , 延时为整数 ;N 为奇数时 , 延时 为整 数加半个抽样周期。不管 N 为奇偶 , 此时 h( n ) 都应满足对 n=(N-1)/2 轴呈偶对称。

信号与系统课件--第七章有限长单位脉冲响应FIR数字滤波器的设计方法

信号与系统课件--第七章有限长单位脉冲响应FIR数字滤波器的设计方法

z FIR数字滤波器的 H(z)为 1 的多项式 ,设计方法
有: 1、直接近似法-----窗函数法 2、频率抽样法、等波纹逼近法
§7.1 线性相位FIR滤波器的特点
一,实现线性相位的条件
N 1
H (e j ) h(n)e jn H g ( )e j ( ) n0
第一类线性相位: ()=- 第二类线性相位: ()=0-
N为奇数 h(n)
0
C(n) 1
0
N- 1 n
N 1
n
2
( N 1) / 2
H () c(n)sin(n)
n1
H( )
2
o
N为偶数 h(n)
0
d(n) 1
0
N- 1 n
N
n2H来自()N /2
d
n1
(n)
s i n
n
1 2
H( )
o
2
四,线性相位FIR的零点
1,一般 : 四个零点为一组
∵ h(n) 实序列 有:
二,线性相位的特点
系统函数 1、h(n)偶对称
H (z)
( N 1)
z2
N 1
h(n)
z
N 1n 2
z
N 1n 2
n0 2
线性相位
H (e
j
)
j ( N 1)
e2
N 1
h(n) cos[(n
N
1) ]
n0
2
H g
( )
N 1 n0
h(n)
c
os
n
N 1
2
() 1 (N 1)
第七章 有限长单位脉冲响应(FIR)
数字滤波器的设计方法

第七章.有限冲激响应数字滤波器设计

第七章.有限冲激响应数字滤波器设计

FIR数字滤波器设计的窗函数法
B
A
D
取N=10,矩形窗频谱:
N
π
2

N
π
2
2、FIR数字滤波器设计的窗函数法2、FIR数字滤波器设计的窗函数法
Bartlett窗
时域波形
频谱(N=21)
2、FIR数字滤波器设计的窗函数法2、FIR数字滤波器设计的窗函数法
FIR数字滤波器设计的窗函数法数字滤波器设计的窗函数法
汉宁窗的幅度特性
FIR数字滤波器设计的窗函数法
4. 汉明(Hamming)窗——改进的升余弦窗
2、FIR数字滤波器设计的窗函数法
2、FIR数字滤波器设计的窗函数法
其幅度函数由五部分组成,它们都是移位不同,且幅度也不同的W R(ω)函数,使旁瓣再进一步抵消。

阻带衰减进一步增加,这样可以得到更低的旁瓣,但主瓣宽度却不得不加宽到矩形窗的3倍,即B=12π/N,第一旁瓣的值比主瓣小58dB,即A=-58dB,旁瓣谱峰渐近衰减速度为D=-18dB/oct。

常用的窗函数2、FIR数字滤波器设计的窗函数法
常用窗函数的幅
度特性 
(a)矩形窗;
(b)巴特利特窗
(三角形窗);
(c)汉宁窗; 
(d)哈明窗;
(e)布莱克曼窗
2、FIR数字滤波器设计的窗函数法
理想低通加窗后的
幅度特性
(N=51,ωc=0.5π)
(a)矩形窗;
(b)巴特利特窗(三
角形窗);
(c)汉宁窗; 
(d)哈明窗;
(e)布莱克曼窗
低通幅度特性。

有限长单位冲激响应(FIR)数字滤波器的设计方法.

有限长单位冲激响应(FIR)数字滤波器的设计方法.
(2) 零点zi在单位圆上,但不在实轴上,即 zi rie ji , ri 1,i 0或
系统函数的基本因子为
主要参数:N=3,=(N-1)/2=1;
z2 1 0
z2
Re[ z] 1
7.2 线性相位FIR数字滤波器的特点
四、 零点位置
第一类和第二类线性相位的系统函数满足等式:H (z) z(N1)H (z1)
h(n)sin[( N 1 n)]
n0
2
因此
H g ( )
N 1 n0
h(n)sin[( N 1 n)]
2
Q( ) ( N 1)
2
2
, , 为常数 N 1
2
2
这种使所有频率的相移皆为90的网络,称为90移相器,或称正交变换网络,它和理
想低通滤波器、理想微分器一样,有着极重要的理论和实际意义。
N 1/ 2
H () c(n)sin n
其中
c(n) 2h N 1 n, n 1,2, , N
n0
2
2
特点: H(0)=0, H()=0, H(2)=0; 对= 0,, 2奇对称;
7.2 线性相位FIR数字滤波器的特点
三、 幅度函数的特点
4. h(n)奇对称,N为偶数
N 1
设希望设计的滤波器传输函数为Hd(ejω), hd(n)是与其对应的单位脉冲响应,因此
设计可在时域进行,使设计的FIR滤波器的h(n)逼近hd(n),以低通为例:
Hd
(e
j
)
e ja
0,c
, c
0 c hd(n) Nhomakorabea1
2
c c
e jae jnd sin(c (n a)) (n a)

有限脉冲响应数字滤波器设计实验报告

有限脉冲响应数字滤波器设计实验报告

成绩:《数字信号处理》作业与上机实验(第二章)班级:学号:姓名:任课老师:完成时间:信息与通信工程学院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具有线性相位特性。

第7章 FIR-有限长单位脉冲响应 及其滤波器设计

第7章 FIR-有限长单位脉冲响应 及其滤波器设计

3 0 1 2
4
5
h(n) 奇 对称, N为偶 数
n
图7-2 h(n) 奇对称
第7章有限长单位脉冲响应(FIR)数字滤波器的设计方法 7.2.2 线性相位FIR滤波器频率响应特点 先看h(n)偶对称的情况: h(n)=h(N-1-n) 其系统函数为
N −1 n=0
0≤n≤N-1
N −1 n=0
(7-13)
h(n) 偶 对称, N为奇 数
n
h(n) 偶 对称, N为偶 数
(N-1)/2 偶对称 中心
h(n)
0
1
图7-1 h(n) 偶对称
2
3
4
5
n
第7章有限长单位脉冲响应(FIR)数字滤波器的设计方法
h(n)
(N-1)/2 奇对称 中心
4 0 1 2 3
5
6
h(n) 奇 对称, N为奇 数
n
h(n)
(N-1)/2 奇对称 中心

(7-20)
式中,grd(group
delay)为群延迟函数。由式(7-20)可
知,当h(n)满足偶对称时,FIR数字滤波器具有(N-1)/2个 采样的延时, 它等于单位脉冲响应h(n)长度的一半。也 就是说,FIR数字滤波器的输出响应整体相对于输入延时 了(N-1)/2个采样周期。 再看h(n)奇对称的情况: h(n)=-h(N-1-n) 0≤n≤N-1 (7-21)
第7章有限长单位脉冲响应(FIR)数字滤波器的设计方法 •再者,只要经过一定延时,任何非因果有限长序列都可以 变成因果有限长序列,因而总能用因果系统来实现。 •还有,FIR滤波器由于单位冲激响应是有限长的,因而可 利用FFT实现信号过滤,从而使运算效率得到大幅提升。 •但是,要取得很好的衰减特性,FIR滤波器H(z) 的阶次比 IIR滤波器的要高。 •IIR滤波器设计中的各种变换法对FIR滤波器设计是不适用 的,这是因为那里是利用有理分式的系统函数,而FIR滤 波器的 系统函数只是z-1的多项式。 •从以上讨论可见,我们最感兴趣的是具有线性相位的FIR 滤波器。对于非线性相位的FIR滤波器,一般可用IIR滤波 器来代替,IIR滤波器所需阶数比FIR滤波器的要少很多。

第7章有限长单位抽样响应数字滤波器的设计

第7章有限长单位抽样响应数字滤波器的设计

-1
0
1
2
3
n1
/
不适合低通滤波器、高通滤波器、带阻滤波器
第7章 有限长单位抽样响应数字滤波器的设计
27 /178
4、h(n)为奇对称,N为偶数
H ( )
N 1
h(n) sin[ (
N
1
n)]
n0
2
N 1
2
2h(n)sin[( N 1 n)]
n0
2
N
2 2h( N m)sin[(m 1)]
3 /178
FIR数字滤波器的优点: 1)可以做到严格线性相位 2)可以具有任意的幅度特性
3)对一个信号滤波,相当于 x n hn
可以用FFT计算 4)总可以用一个因果系统来实现 5) 无反馈运算,运算误差小。
第7章 有限长单位抽样响应数字滤波器的设计
4 /178
FIR数字滤波器的缺点
1)不能用模拟滤波器的设计方法,无解 析设计公式,要借助计算机辅助设计;
第7章 有限长单位抽样响应数字滤波器的设计
18 /178
Ⅰ、第一类线性相位
H ( )
N 1 n0
h(
n
)
c
o
s
n
N
2
1
( ) ( N 1 )
2
Ⅱ、第二类线性相位
H ( )
N 1
h(n ) sin
n0
n
N
2
1
( ) ( N 1 )
22
第7章 有限长单位抽样响应数字滤波器的设计
第7章 有限长单位抽样响应数字滤波器的设计
1 /178
第7章 有限脉冲响应数字滤波器的设计
第7章 有限长单位抽样响应数字滤波器的设计

有限脉冲响应数字滤波器的设计

有限脉冲响应数字滤波器的设计

切比雪夫最佳一致逼近准则
设希望设计的滤波器幅度特性为Hd(ω),实际设计的滤波器幅度特性为Hg(ω),其加权误差E(ω)用下式表示:
E(ω)=W(ω)[Hd(ω)-Hg(ω)] (7.4.3)
为设计具有线性相位的FIR滤波器,其单位脉冲响应h(n)或幅度特性必须满足一定条件。假设设计的是h(n)=h(n-N-1),N=奇数情况,
凯塞—贝塞尔窗(Kaiser-Basel Window)
01
式中
02
I0(x)是零阶第一类修正贝塞尔函数,可用下面级数计算:
03
一般I0(x)取15~25项,便可以满足精度要求。α参数可以控制窗的形状。一般α加大,主瓣加宽,旁瓣幅度减小,典型数据为4<α<9。当α=5.44时,窗函数接近哈明窗。α=7.865时,窗函数接近布莱克曼窗。凯塞窗的幅度函数为
图7.2.4 常用的窗函数
图7.2.5 常用窗函数的幅度特性 矩形窗;(b)巴特利特窗(三角形窗);(c)汉宁窗; 哈明窗;(e)布莱克曼窗
图7.2.6 理想低通加窗后的幅度特性(N=51,ωc=0.5π) 矩形窗;(b)巴特利特窗(三角形窗);(c)汉宁窗; 哈明窗;(e)布莱克曼窗
汉宁(Hanning)窗——升余弦窗
图7.2.3 汉宁窗的幅度特性
4.哈明(Hamming)窗——改进的升余弦窗
(7.2.11)
其频域函数WHm (e jω)为
其幅度函数WHm(ω)为
当N>>1时,可近似表示为
(7.2.13)
1
其频域函数为
2
其幅度函数为
3
(7.2.14)
4
5.布莱克曼(Blackman)窗
(7.1.16)

数字信号处理导论

数字信号处理导论
成果如右图
接上例:N=10
分别用矩形窗
和Hamming 窗
使用Hamming 窗后,阻带衰 减变好,但过 渡带变宽。
高通:
令:
H
d
(e
j
)
e
j
N 1 2
0
c 0 c
hd
(n)
sin[(n
N 1) ] sin[(n
2
(n N 1)
N 2
1)c
]
2
相当于用一种截止频率在 处旳低通滤波器
理想微分器 x(t) 旳频率特征:
H (s)
y(t)
y(t) dx(t) dt
H (s) s H ( j) j
令: x(n) x(t) tnTs , y(n) y(t) tnTs
x(n) H (z)
y(n)
理想差分器
旳频率特征: H (e j ) j,
Hd (e j ) j
奇对称, 纯虚函数
主瓣宽度最宽:12
N
旁瓣幅度最小
汉宁窗-布拉克曼窗比较
矩形窗-汉宁窗-布拉克曼窗比较
矩形窗-三角形窗比较
矩形窗-海明窗-凯泽窗比较
六种窗函数基本参数比较
窗函数
矩形窗 三角形窗
汉宁窗 海明窗 布拉克曼窗 凯泽窗
窗谱性能指标
加窗后滤波器性能指标
旁瓣峰值 主瓣宽度 过渡带宽 阻带最小衰减
/dB

c
2
N
处出现肩峰值,两侧形成起伏振
荡,振荡旳幅度和多少取决于旁瓣旳幅度和多少
变化N只能变化窗谱旳主瓣宽度,但不能变化主
瓣与旁瓣旳相对百分比。其相对百分比由窗函数 形状决定,称为Gibbs效应
例1.设计低通 FIR

7-有限长数字滤波器的设计

7-有限长数字滤波器的设计

的通带内,这时应出现正的肩峰。
( 4)
时,主瓣全部在通带外, 出现负的肩峰。
c 2 / N
(5)当
时,随
增加,
左边 旁瓣的起伏部分扫过通带,卷积
也随着 变化而变化,故 的旁瓣在通带内的面积 将围绕着零值而波动。
(6)当
为幅度函数,
(3)h(n)的频响
其中, 为相位函数。
为幅度函数,
4、窗函数频响产生的影响 从几个特殊频率点的卷积过程就可看出其影响:
(1 )
时,
也就 在 到 全部面积的积分。 因此,H(0)/H(0)=1(用H(0)归一化)。
0

0
(2 )
时,
正好与 。
的一半相重叠。这时有
(3)
时,
的主瓣全部在
4、FIR的系统函数是Z-1的多项式,故IIR的方法不适用。 5、FIR的相位特性可以是线性的,因此,它有更广泛的 应用,非线性的FIR一般不作研究。
7-2
线性相位FIR DF的特点
一、线性相位的条件 如果FIR DF的单位抽样响应h(n)为实数,而且 满足偶对称h(n)=h(N-1-n),或满足奇对称 h(n)=-h(N-1-n),其对称中心在 n N 1 处,可证 2 明filter就具有准确的线性相位。
§ 7-1
引言
一、IIR DF的特点 1、DF的设计依托AF的设计,有图表可查,方便简单。 2、相位的非线性 H(Z)的频响: H (e j ) H ( Z ) j H (e j ) e j ( ) ,
Z e
j H ( e ) 是幅度函数, 其中, ( ) 是相位函数。

N /2
呈奇对称,而对
呈偶对称。

有限冲激响应滤波器的设计

有限冲激响应滤波器的设计
2
(7.11) (7.12)
I型线性相位滤波器的幅度函数和相位函数的特点:
幅度函数对 N 1 偶对称,同时对 0, , 2 也呈偶对称;
2 相位函数为准确的线性相位。
2021/4/22
9
2.Ⅱ型线性相位滤波器
一个Ⅱ型线性相位滤波器,由于N是偶数,所以, h(n) 的对称
中心在半整数点 N 1 。其频率响应可以表示为:
(7.9)
2021/4/22
6
把满足式(7.7)、(7.8)和式(7.9)的奇对称条件的FIR 滤波器分别称为Ⅲ型线性相位滤波器和Ⅳ型线性相位滤波 器。
2021/4/22
7
7.1.2 线性相位滤波器频率响应的特点
1.I型线性相位滤波器
由于偶对称性,一个I型线性相位滤波器的频率响应可表示为
其中
( N 1) / 2
(7.24)
式中“ ”号相当于i ,零点在负实轴上,“ ”相当i 于0
零点在正实轴上。
在此情况下, N 3
N 1 1 2
2021/4/22
19
(4)零点 z zi 既在单位圆上,但在实轴上,即ri 1,i 0或,
此时零点只有两种情况,即 z 1或z 1 ,这时零点既是自己
的复共轭,又是倒数,其基本因子为
1 ) ]
(7.19)
n1
2
其中 d(k) 2h( N k) k 1, 2,..., N
2
2
2021/4/22
14
幅度函数为 相位函数为
H () N /2 d (k) sin[(k 1)]
n1
2
() (N 1)
2
2
(7.20) (7.21)
Ⅳ型线性相位滤波器的幅度函数和相位函数的特点:

第7章 有限长单位抽样响应数字滤波器的设计

第7章 有限长单位抽样响应数字滤波器的设计

12 /94
N 1 2 h( n) h( N 1 n) 0 n N 1
第7章 有限长单位抽样响应数字滤波器的设计
13 /94
(2)满足第二类线性相位的条件是: h(n)是实序列 且对 (N-1)/2奇对称,即 (7.1.10) N 1 N 1 证明: n n
H ( ) cos() h(n) cos(n)
n 0 N 1
H ( ) sin( ) h(n) sin(n)
n 0
N 1
第7章 有限长单位抽样响应数字滤波器的设计
11 /94
将上述两式相除,得
sin( ) cos( )
n 0 N 1 n 0
h(n) sin(n)
N 1
h(n) cos(n)
化简,得
N 1 n 0
h(n) sin[( n)] 0
n ( N 1) / 2 观察上式中的正弦函数,它是关于 奇 h( n) 对称的,为了使上式成立,则 必是偶对称的, 即
第7章 有限长单位抽样响应数字滤波器的设计
第7章 有限长单位抽样响应数字滤波器的设计
14 /94
同样可以表示为
N 1 1 1 H ( z ) [ H ( z ) z ( N 1) H ( z 1 )] h(n)[z n z ( N 1) z n ] 2 2 n 0
z
j

N 1 N 1 2
N 1 n 1 n N21 2 h ( n ) [ z z ] 2 n 0
第7章 有限长单位抽样响应数字滤波器的设计
9 /94
将z=e jω代入上式,得到:
H (e ) e
j j( N 1 N 1 ) 2
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

第7章 有限长单位抽样响应数字滤波器的设计
40 20
26 /178
H ( )
( N 1) / 2
Hr()
0 -20 -40 -1

c ( n ) sin( n )
n1
0
1 /
2
3
不适合低通滤波器、高通滤波器、带阻滤波器
第7章 有限长单位抽样响应数字滤波器的设计
27 /178
4、h(n)为奇对称,N为偶数
( N 1)/ 2 N 1 N 1 h 2h m cos( m ) 2 2 m 1
N 1 a(0) h( ) 2 a( n) 2h( N 1 n), n 1, 2, 3, , N 1 2 2
H (z )
1
N 1 j N 1 2 2
N 1 h( n) sin n 2 n0
第7章 有限长单位抽样响应数字滤波器的设计
17 /178
幅度函数 相位函数
N 1 H ( ) h( n ) sin n 2 n0
N 1 H ( ) h( n)sin[ ( n)] 2 n 0
N 1
N 1 2h( n)sin[ ( n)] 2 n0
N 1 2h( m )sin[ (m )] 2 2 m 1
N 2
N 1 2
第7章 有限长单位抽样响应数字滤波器的设计
40 20
奇对称
h(n)sin[( n)] 0
( N 1) h(n)关于 n 2 偶对称
第7章 有限长单位抽样响应数字滤波器的设计
16 /178
通过类似的推导,可以得到满足第二类线性相 位条件的系统函数
H (z) z
频率响应函数函数
H (e ) e
j j
( N 1)
N 1
N 1 2h( n)cos n 2 n 0
N 1 2
第7章 有限长单位抽样响应数字滤波器的设计
40 20
23 /178
1 H ( ) b( n)cos n 2 n 1
N /2
13 /178
证明(第一类线性相位条件):
N 1 幅度函数 H ( ) h( n ) cos n 2 n0
N 1
相位函数
N 1 ( ) ( ) 2
0
( N 1 ) 2
( )

2

严格的线性相位
( N 1)

n 1
N 1 N 1 2 h( n ) sin n h( ) 2 2 =0
H ( )
( N 1) / 2

c ( n ) sin( n )
n1
N 1 N 1 c ( n ) 2 h( n ), n 1, 2, , 2 2
N 1
( )

2 N ( 2) 2
0
N 1 ( ) ( ) 2 2
广义的线性相位

2

3 ( N ) 2
第7章 有限长单位抽样响应数字滤波器的设计
18 /178
结论
Ⅰ、第一类线性相位
N 1 N 1 H ( ) h ( n ) cos n ( ) ( ) 2 2 n0
第7章 有限长单位抽样响应数字滤波器的设计
1 /178
第7章 有限脉冲响应数字滤波器的设计
第7章 有限长单位抽样响应数字滤波器的设计
2 /178
IIR数字滤波器的优点:

可以利用模拟滤波器的设计结果
IIR数字滤波器的缺点:

非线性相位 若需线性相位,要采用全通网络进行相 位校正
第7章 有限长单位抽样响应数字滤波器的设计
h(m)z
m0
N 1
( N m 1)
h(m)z
m
第7章 有限长单位抽样响应数字滤波器的设计
11 /178
证明(第一类线性相位条件):
H (z) z H (z ) 1 ( N 1) 1 H (z) H ( z ) z H ( z ) 2 1 N 1 n ( N 1) n h( n) z z z 2 n0
第一类线性相位成立的充分必要条件
N 1 2 h( n) h( N 1 n) 0 n N 1
第二类线性相位成立的充分必要条件 N 1
2 2 h( n) h( N 1 n) 0 n N 1
第7章 有限长单位抽样响应数字滤波器的设计
10 /178
证明(第一类线性相位条件):
充分性 FIR滤波器的系统函数
H ( z ) h(n)z
n 0
N 1 n 0
N 1
n
h( N n 1)z n
m N 1 n n N 1 m
z
( N 1) N 1 m0
偶对称
N 1 cos 偶对称 ( N 1 n) 2
以(N-1)/2为中心,把两两相等的项进行合并
第7章 有限长单位抽样响应数字滤波器的设计
21 /178
( N 3)/ 2 N 1 N 1 H ( ) h 2h n cos( n ) 2 2 n0
20 /178
7.1.3 线性相位FIR滤波器幅度特性
1、h(n)为偶对称,N为奇数
N 1 H ( ) h ( n ) cos n 2 n0
N 1
N 1 (N 1) N 1 ( N 1) cos n cos n n n 关于 关于 2 2 2 2
e e
j(
N 1 N 1 ) 2
N 1 h( n )cos n 2 n0 N 1 h( n )cos n n0 2
j(
N 1 N 1 ) 2
第7章 有限长单位抽样响应数字滤波器的设计
N 1 n 0
h(n)sin( n)
第7章 有限长单位抽样响应数字滤波器的设计
15 /178
证明(第一类线性相位条件):
N 1
sin( ) cos( )
N 1 n 0
n 0 N 1 n 0
h(n)sin( n)
h(n)cos( n)
( N 1) 关于 n 2
第7章 有限长单位抽样响应数字滤波器的设计
7 /178
7.1.1. 线性相位的定义
FIR滤波器的频率响应函数
H (e j ) h(n)e j n H (e j ) e j ( )
n 0
N 1
H (e j ) e j ( ) H ( )e j ( )
第7章 有限长单位抽样响应数字滤波器的设计
14 /178
证明(第一类线性相位条件):
必要性
H (e j ) H ( )e j h(n)e jn
n 0 N 1
实部、虚部分别相等
H ( )cos( )
两式相除
N 1 n 0
h(n)cos( n)
H ( )sin( )
28 /178
Hr()
0 -20 -40 -1
0
1
2
3
/
不适合于低通滤波
器、带阻滤波器 h(n)为奇对称时,有90o 相移,适用于微分器和 90o移相器,而选频滤 波器采用h(n)为偶对称
第7章 有限长单位抽样响应数字滤波器的设计
0 -20 -1 0
n0
Hr()
H ( )
( N 1) / 2

40
a ( n ) cos( n )
20
1 /
2
3
第7章 有限长单位抽样响应数字滤波器的设计
22 /178
2、h(n)为偶对称,N为偶数
N 1 H ( ) h( n)cos n 2 n 0
2)由于没有极点,同样幅频特性,要获 得较好的过渡带特性,需较高的阶数做为代 价。
第7章 有限长单位抽样响应数字滤波器的设计
5 /178
FIR滤波器的设计方法 设计任务:选择有限长度的脉冲响应h(n),得 到系统函数H(z),使幅频特性满足技术指标, 同时使相频特性达到线性相位。
设计方法:
窗函数法Βιβλιοθήκη Hr()0 -20 -40 -1
0
1
2
3
N b( n ) 2h n n=1,2, 3, …, N/2 不适合设计高通滤 2
/
波器、带阻滤波器 1 cos n 0 , H(ω)对ω=π呈奇对称 ω=π时, 2
1 cos n 1, H(ω)对ω=0呈偶对称 ω=0时, 2 H(ω)对ω=2π呈偶对称
频率采样法
切比雪夫等波纹逼近法
第7章 有限长单位抽样响应数字滤波器的设计
6 /178
7.1 线性相位FIR数字滤波器的特性
本节主要介绍FIR滤波器具有线性相位的条 件及幅度特性以及零点、网络结构的特点。 FIR滤波器的单位冲激响应: 系统函数:
在 z 平面有N –1 个零点
在 z = 0 处是N –1 阶极点
第7章 有限长单位抽样响应数字滤波器的设计
24 /178
3、h(n)为奇对称,N为奇数
N 1 H ( ) h ( n ) sin n 2 n0
相关文档
最新文档