数字信号处理第七章
- 格式：doc
- 大小：2.28 MB
- 文档页数：27
下载文档原格式
合集下载
数字信号处理第七章
数字信号处理第七章第七章数字滤波器设计7.1：无限脉冲响应滤波器的阶数估计q7.1用mattab确定一个数字无限冲激响应低通滤波器所有四种类型的最低阶数。
指标如下:40khz的抽样率，,4khz的通带边界频率，8khz的阻带边界频率，0.5db的通带波纹，40db的最小阻带衰减。
评论你的结果。
答：标准通带边缘角频率wp是：标准阻带边缘角频率WS为：理想通带波纹rp是0.5db理想阻带波纹rs是40db1.使用这些值，巴特沃斯低通滤波器的最低阶数为n=8，相应的标准通带边缘频率wn 为0.24692.使用这些值得到切比雪夫1型低通滤波器最低阶数n=5，相应的标准通带边缘频率wn是0.2000.3/使用这些值，切比雪夫2型低通滤波器n=5的最低阶数和相应的标准通带边缘频率wn为0.40004.使用这些值得到椭圆低通滤波器最低阶数n=8，相应的标准通带边缘频率wn是0.2000.从以上结果中观察到椭圆滤波器的阶数最低，并且符合要求。
问题7。
2.用MATLAB确定四种数字无限冲激响应高通滤波器的最低阶数。
指标如下：3500hz采样率、1050hz通带边界频率、600Hz阻带边界频率、1dB通带纹波和50dB最小阻带衰减。
对结果的评论a：标准通带边缘角频率WP为：标准阻带边缘角频率ws是：理想通带纹波RP为1dB，理想阻带纹波RS为50dB1.使用这些值得到巴特沃斯高通滤波器最低阶数n=8，相应的标准通带边缘频率wn是0.5646.2.使用这些值，切比雪夫1高通滤波器的最低阶数为n=5，相应的标准通带边缘频率wn为0.60003.使用这些值得到切比雪夫2型高通滤波器最低阶数n=5，相应的标准通带边缘频率wn是0.3429.4.使用这些值，椭圆低通滤波器的最低阶数n=4，相应的标准通带边缘频率wn为0.6000。
从上述结果可以看出，椭圆滤波器的阶数最低，满足要求。
q7.3用matlab确定一个数字无限冲激响应带通滤波器所有四种类型的最低阶数。
数字信号处理第七章有限单位冲激响应FIR数字滤波器的设计方法(共95张PPT)
线性相位分析
H (z)z (N 2 1 )N n 0 1h (n ) 1 2Z (n (N 2 1 )) 1 2Z (n (N 2 1 ))
H (ej)e e j j(( N )2 1) N n 0 1 h( n) c o s(n (N 2 1 ) ) (1) H ()
m 0
即 H (z) z (N 1 )H (z 1 )
H (z) z (N 1 )H (z 1 )
所以有： h (z) 1H (z) z (N 1 )H (z 1 ) 2
1N 1h (n )z nz (N 1 )zn 2n 0
z (N 2 1 )N n 0 1 h (n ) 1 2Z (n (N 2 1 )) 1 2Z (n (N 2 1 ))
m1
(N 1)/2a(n)con s)(
n0
其中: a ( 0 ) h (N 1 ),a ( n ) 2 h ( n N 1 ),( n 1 )
2
2
由于con s对 0,,2
是偶对称的。
因此，H()对0,,2
为偶对称。
线性相位滤波器的幅度特点
2、h(n)偶对称，N为偶数
对（1）式与如上合并项，注意到由于N为偶数， h(N 1) 项即为0，则
四种线性相位滤波器
偶对称单位冲激响应
h (n ) ＝h (N－ 1－n )
相位响应
( ) N 1 2
情
况
( )
1
o
－ N( － 1)
N为 奇 数 h (n )
0 a (n )
N－ 1 n
0
N 1
n
2
( N 1) / 2
H ( ) a (n) cos n
n0
精品课件-数字信号处理(第三版) 刘顺兰-第7章
第７章数字信号处理中的有限字长效应
7.1.2 定点制误差分析 1. 数的定点表示 定点制下，一旦确定了小数点在整个数码中的位置，在整个
运算过程中即保持不变。因此，根据系统设计要求、 数值范围来 确定小数点处于什么位置很重要，这就是数的定标。 数的定标有Q表示法和S表示法两种。Q表示法形如Qn，字母Q后的 数值n表示包含n位小数。如Q0表示小数点在第0位的后面，数为整 数；Q15 表示小数点在第15位的后面，0～14位都是小数位。S表 示法则形如Sm.n，m表示整数位，n表示小数位。以16位DSP为例， 通过设定小数点在16位数中的不同位置，可以表示不同大小和不 同精度的小数。表7.1列出了一个16位数的16种Q表示、 S表示及 它们所能表示的十进制数值范围。
小的正数： (01.000..0)2×2－127=1×2－127≈5.9×10－39
(4) 当S=1，E=－127，F的23位均为1时，表示的浮点数为绝 对值最小的负数：
(10.111..1)2×2－127=(－1－2－23)×2－127≈－5.9×10－39 双精度浮点数占用8个字节(64位)存储空间，包括1位符号位、 11位阶码、 52位尾数，数值范围为1.7E-308～1.7E+308。
第７章数字信号处理中的有限字长效应
乘除运算时，假设进行运算的两个数分别为x和y，它们的Q 值分别为Qx和Qy，则两者进行乘法运算的结果为xy，Q值为Qx+Qy， 除法运算的结果为x/y，Q值为Qx－Qy。
在程序或硬件实现中，上述定标值的调整可以直接通过寄存 器的左移或右移完成。若b>0，实现x×2b需将存储x的寄存器左 移b位；若b<0，实现x×2b则需将存储x的寄存器右移|b|位即可。
称为小数点位置。
第七章 模拟滤波器的设计(数字信号处理)
1 (
s
c
)
2N
10
a s / 10
(7.2.15)
由(7.2.14)和(7.2.15)式得到：
(
p
s
)
N
10 10
a p / 10 a s / 10
1 1
令
sp s / p , k sp
10 10
a p 10 as 10
1 1
,则N由下式表示：
N
1
1
1
1
0
fC a ) 低通
f
0
fC b ) 高通
f
0
fC1 c) 带通
fC2
f
0
fC1 d ) 带阻
fC2 f
7.1 理想滤波器
无过渡带且在通频带内满 足不失真测试条件的滤波 器称为理想滤波器。理想 滤波器的频率响应函数为：
|H(f)| A0
-fc
A e j 2 p ft 0 0 H(f) 0 f fc 其它
lg k sp lg sp
(7.2.16)
用上式求出的N可能有小数部分，应取大于等于N
的最小整数。关于3dB截止频率Ωc，如果技术指标中没 有 给 出 ， 可 以 按 照 (7.2.14) 式 或 (7.2.15) 式 求 出 ， 由
图7.2.2 低通滤波器的幅度特性
滤波器的技术指标给定后，需要设计一个传输函
数Ha(s)，希望其幅度平方函数满足给定的指标αp和αs， 一般滤波器的单位冲激响应为实数，因此
H a ( j )
2
H a ( s )G ( s )
s j
H a ( j ) H a ( j )
s
c
)
2N
10
a s / 10
(7.2.15)
由(7.2.14)和(7.2.15)式得到：
(
p
s
)
N
10 10
a p / 10 a s / 10
1 1
令
sp s / p , k sp
10 10
a p 10 as 10
1 1
,则N由下式表示：
N
1
1
1
1
0
fC a ) 低通
f
0
fC b ) 高通
f
0
fC1 c) 带通
fC2
f
0
fC1 d ) 带阻
fC2 f
7.1 理想滤波器
无过渡带且在通频带内满 足不失真测试条件的滤波 器称为理想滤波器。理想 滤波器的频率响应函数为：
|H(f)| A0
-fc
A e j 2 p ft 0 0 H(f) 0 f fc 其它
lg k sp lg sp
(7.2.16)
用上式求出的N可能有小数部分，应取大于等于N
的最小整数。关于3dB截止频率Ωc，如果技术指标中没 有 给 出 ， 可 以 按 照 (7.2.14) 式 或 (7.2.15) 式 求 出 ， 由
图7.2.2 低通滤波器的幅度特性
滤波器的技术指标给定后，需要设计一个传输函
数Ha(s)，希望其幅度平方函数满足给定的指标αp和αs， 一般滤波器的单位冲激响应为实数，因此
H a ( j )
2
H a ( s )G ( s )
s j
H a ( j ) H a ( j )
精品课件-数字信号处理-第7章
xA (n) x(n) jxˆ(n)
(7-7)
式中 xˆ(n) 是时间离散信号x(n)
xˆ(n) x(n) h(n)
（7-8）
解析信号对实信号来说就是有一阶导数的连续信号。由此意 义来说，任何序列都不是解析信号，因为它是一个以整数为变量 的函数，但xA(n)是xA(t)的采样，如果xA(t)是解析的， 我们仍认 为xA(n)也是解析的， 这是对解析信号的修正。
第七章 离散希尔伯特变换
７
h(n) |H(j )|
－
n) π
2
－ π 2
(a)
－7 －6 －5 －4 －3 －2－1
1 2 3 4 5 6 7 n
(b)
图7.1 (a) 频域特性； (b) 时域特性
第七章 离散希尔伯特变换
７
7.4 因果序列傅里叶变换下的希尔伯特变换
当xA(n)是解析序列时，其实部和虚部成希尔伯特变换关系。 它对应的频谱则是单边的。如果把频谱看成解析的，即其实部与 虚部成希尔伯特变换关系，则对应的时域序列应是单边的， 即 因果的。本节主要讨论因果序列傅里叶变换的希尔伯特变换。
第七章 离散希尔伯特变换
７
7.2 时间连续信号的希尔伯特变换
给定一时间连续信号x(t)，其希尔伯特变换 xˆ(t)定义为
xˆ(t) 1 x( )d 1 x(t )d x(t) 1
π t
π
πt
（7-1）
xˆ(t) 可以看成是x(t)通过单位冲激响应 h(t) 1 滤波器
的输出。
第七章 离散希尔伯特变换
７
在时间连续信号处理中解析信号是一个重要的概念，本章我 们将其推广到时间离散信号。从形式上说不能把复时间离散信号 或复序列看成是解析函数，因为它是一个以整数为变量的函数， 但是也可以按照类似的处理方式，将复序列之实部和虚部联系起 来使复序列的频谱在单位圆上的-π≤ω<0范围内为零。用类似 的方法也可以将周期性（或有限时宽）序列的傅里叶变换之实部 和虚部联系起来，在这种情况下，“因果性”条件是，该周期序 列在各周期的后半部为零。根据对偶关系，对于时间序列呈单边 特性的因果序列，在频域(其实部与虚部)也应存在某种变换关系。 最小相位序列是一类很重要的信号， 其傅里叶变换幅度和相位 之间存在希尔伯特变换关系。
第七章 数字信号处理中的有限字长效应
设系数采用b位量化长度和舍入方式进行量化，系数量化误
差为e(n)，其变化范围 ( / 2, / 2) ，均值为0，方差为 2 /12
则实际系数为：
ˆ h(n) h(n) e(n)
0 n ( N 1) / 2
ˆ 且量化后 h(n) 也一定满足偶对称，即
ˆ ˆ h(n) h( N 1 n)
2.有限字长效应对信号量化的影响;
3.有限字长效应对系统参数表示的影响
4.有限字长效应在运算过程中的影响
7.1
数字信号处理中的有限长效应
有限字长效应:
在实际的处理过程中，数字信号和系统都不是无限精度的，而是有 限精度，精度的大小则有字长的大小决定，正是由于有限精度，从而给 原有的数字信号处理系统带来了影响，这种影响称为数字信号处理中的 有限字长效应。
z1 0.85 j 0.15
求得a2对z1和z2的影响
z2 0.85 j 0.15
z1 1 j 900 3.3333e a2 z1 z2
z2 1 j 900 3.3333e a2 z2 z1
可见， a2对z1和z2的影响是相同的。因而
z2 z2 a2 a2
i 1 i 1
b
b1
i b 1
b1
ai 2 i
故截尾误差满足：
0 ET (2b 2b1 ), x 0
即
0 ET , x 0
②对于反码负数
b
x 1 2 b1 ai 2 i
i 1
b1
ET QT [ x] x 1 2 ai 2 (1 2
若采用截尾处理，试分别求出原码负数1.1001、反码负数1.1100
数字信号处理讲义第7章滤波器的设计方法
第７章滤波器的设计方法教学目的1．掌握由连续时间滤波器设计离散时间IIR滤波器的方法，包括冲激响应不变法，双线性变换法等；2．了解常用的窗函数，掌握低通IIR滤波器的频率变换法、用窗函数法设计FIR滤波器的方法；3．掌握FIR滤波器的逼近原理与设计方法。
教学重点与难点重点：本章是本课程的重中之重，滤波器的设计是核心内容之一。
1．连续时间滤波器设计离散时间IIR滤波器的方法，包括冲激响应不变法，双线性变换法等；2．常用的窗函数，掌握低通IIR滤波器的频率变换法、用窗函数法设计FIR滤波器的方法；3．掌握FIR滤波器的逼近原理与设计方法。
难点：1.冲激响应不变法，双线性变换法2.用窗函数法设计FIR滤波器FIR滤波器的逼近原理与设计方法基本概念7.0.1 选频滤波器的分类数字滤波器是数字信号处理的重要基础。
在对信号的过滤、检测与参数的估计等处理中, 数字滤波器是使用最广泛的线性系统。
数字滤波器是对数字信号实现滤波的线性时不变系统。
它将输入的数字序列通过特定运算转变为输出的数字序列。
因此，数字滤波器本质上是一台完成特定运算的数字计算机。
我们已经知道，一个输入序列x(n)，通过一个单位脉冲响应为h(n)的线性时不变系统后，其输出响应y(n)为∑∞-)(y))()()(n(nn=m*=xmhnhx将上式两边经过傅里叶变换，可得式中，Y (e j ω)、X (e j ω)分别为输出序列和输入序列的频谱函数， H (ejω)是系统的频率响应函数。
可以看出，输入序列的频谱X (e j ω)经过滤波后，变为X (e j ω)H (e j ω)。
如果|H (e j ω)|的值在某些频率上是比较小的，则输入信号中的这些频率分量在输出信号中将被抑制掉。
因此，只要按照输入信号频谱的特点和处理信号的目的，适当选择H (ej ω)，使得滤波后的X (e j ω)H (e j ω)符合人们的要求，这就是数字滤波器的滤波原理。
和模拟滤波器一样，线性数字滤波器按照频率响应的通带特性可划分为低通、高通、带通和带阻几种形式。
数字信号处理第三版第七章
h
N 1 2
0
，得到:
Hg(ω)关于 ω=0, π, 2π三 点奇对称
情况3只能实 现带通滤波器
当ω=0，π, 2π时，
sin［ω(n－τ)］=0, 且sin ［ω(n－τ)］
关于过零点奇对称
M1
Hg() 2h(n)sin[(n)] n0
情况4： h(n)=－h(N－n－1), N为偶数。
吉布斯效应是由于将hd(n)直接截断引起的，因此，也 称为截断效应。
图7.2.2 吉普斯效应
Hd(ejω)是一个以2π为周期的函数，可以展为傅里叶级数，
即
Hd(ej) hd(n)ejn
n
hd(n)，当然就是Hd(ejω)对应的
单位脉冲响应。设计FIR滤波器就是根据要求找到N个傅
2. 线性相位FIR的时域约束条件 线性相位FIR滤波器的时域约束条件是指满足线性
相位时，对h(n) 1) 第一类线性相位对h(n)
相位函数θ(ω)=－ωτ，由式（7.1.1）和（7.1.2）得 到:
N1
H(ej) h(n)ejnHg()ej n0 N 1
h (n )(c o s n jsin n ) H g ( )(c o s jsin )
加矩窗形后窗的幅H滤度(e波特j器性) 的We2R1 幅πgj(度ωπ)π特的2H 1π性卷d等积g(ππ于。H )e理djg(想W )低W R 通Rg(滤g (波器))e的dj幅(度特)d性Hdg(ω)与
将H(ejω)写成H(ejω)=Hg(ω)e－jω ，则
2) 第二类线性相位对h(n)的约束条件
相位函数θ(ω)=－π/2－ωτ，由式（7.1.1）和（7.1.2），
数字信号处理第七章
H d (e j ) h d (n ) h d (n )w (n )
H(ej)h(n)
Hd (e j)为理想低通
滤波器的传输函数。
数字信号处理第七章
h (n )h d(n )R N (n )
如果对截取后的信号进行傅里叶变换，假设采用矩形窗截
取，对截取后信号进行傅里叶变换得：
频域卷积定理
H(ej) 1
Hd
(e
j
)
1 e 0
j
c
c
：低通滤波器的延时
hd(n)
1
2
Hd(ej)ejnd
1
2
c ej
c
ejnd
1
2
c ej(n)d
c
1
2
1
j(n)
ej(n)
|c c
s
in( c(n)) (n)
数字信号处理第七章
理想特性的hd(n)和Hd(ω)
hd
(n)
sin(c(n ) (n )
hd(n)的最大 值是多少？
ej 2 1 H d()W R ()d
H(ej)H()ej 数字信号处理第七章
则实际FIR滤波器的幅度函数H (ω) 为
H ()2 1 H d()W R()d
取样函数
矩形窗
正好是理想滤波器幅度函数与窗函数幅度函数的卷积。
数字信号处理第七章
H(0) 0.5H(0) H(ω)max H(ω)min
③N增加,过渡带宽减小,肩峰值不变。 因主瓣附近
（a）
（b）
hd(n)是一个以(N-1)/2为中心的偶对称的无限长非因果序列， 如果截取一段n=0～N-1的hd(n)作为h(n)，则为保证所得到的
H(ej)h(n)
Hd (e j)为理想低通
滤波器的传输函数。
数字信号处理第七章
h (n )h d(n )R N (n )
如果对截取后的信号进行傅里叶变换，假设采用矩形窗截
取，对截取后信号进行傅里叶变换得：
频域卷积定理
H(ej) 1
Hd
(e
j
)
1 e 0
j
c
c
：低通滤波器的延时
hd(n)
1
2
Hd(ej)ejnd
1
2
c ej
c
ejnd
1
2
c ej(n)d
c
1
2
1
j(n)
ej(n)
|c c
s
in( c(n)) (n)
数字信号处理第七章
理想特性的hd(n)和Hd(ω)
hd
(n)
sin(c(n ) (n )
hd(n)的最大 值是多少？
ej 2 1 H d()W R ()d
H(ej)H()ej 数字信号处理第七章
则实际FIR滤波器的幅度函数H (ω) 为
H ()2 1 H d()W R()d
取样函数
矩形窗
正好是理想滤波器幅度函数与窗函数幅度函数的卷积。
数字信号处理第七章
H(0) 0.5H(0) H(ω)max H(ω)min
③N增加,过渡带宽减小,肩峰值不变。 因主瓣附近
（a）
（b）
hd(n)是一个以(N-1)/2为中心的偶对称的无限长非因果序列， 如果截取一段n=0～N-1的hd(n)作为h(n)，则为保证所得到的
精品文档-数字信号处理(吴瑛)-第7章
第7章 IIR数字滤波器
因此，复幅度平方函数Q(s)的极点分布是以4个(对应复数 极点)或2个(对应实数极点)一组出现的，如图7.2.1所示。当 极点不在坐标轴上时，这时复幅度平方函数的极点是呈4个一 组成对出现的，如图7.2.1(a)所示; 当极点在横坐标轴上时， 复幅度平方函数的极点是呈2个一组成对出现的，如图7.2.1(b) 所示。为保证系统稳定，此时将左半平面的极点归于Ha(s)，零 点可以任意选一半的零点。对于稳定系统，系统函数在虚轴上
系统频率响应如图7.2.3所示。在λ=1处，通带衰减为 3 dB，在λ=4处阻带衰减大于40 dB
第7章 IIR数字滤波器 图7.2.3 例7.2.3巴特沃思型归一化模拟低通滤波器的幅频响应
第7章 IIR数字滤波器
巴特沃思型归一化模拟低通滤波器的幅频响应在通带和阻 带都是随频率单调变化的，因而如果在通带边缘满足指标，则 在通带内肯定有富裕量，也即会超出指标要求。如果将指标均 匀分布在通带或阻带内，在相同的指标要求下，可以设计出阶
b1s N a1s N
1 1
bN 1s aN 1s
bN aN
第7章 IIR数字滤波器
例7.2.3的Matlab程序如下:
Wp=1; Ws=4; Rp=3; Rs=40；
％设置滤波器技术
［N，Wc］=buttord(Wp，Ws，Rp，Rs，′s′)；
％计算滤波器阶数和3dB
［z，p，k］=buttap(N)；
Ha
(
p)
(
p
e j35π
)(
p
e j45π
)(
p
1 ejπ )(
p
e j65π
)(
p
e j75π
数字信号处理程佩青第七章
3）h(n)奇对称，N为奇数
幅度函数：
H
(
)
N 1
h(n) sin
n0
N 1 2
n
sin
N
2
1
(
N
1
n)
sin
n
N 2
1
sin
N 1 2
n
sin
N
2
1
n
对
N 1 2
呈奇对称
h(n)奇对称且N为奇数
h
N 1 2
0
N -3
令
N
1
H ( )
nm
2 n0
2h(n) sin
1 2r cosi z1 z2 N 3 N 1 1
2
3） zi rie ji ri 1 i 0或 ，即零点在实轴上
零点： ri
1 ri
Hi (z)
1 ri z1
1
1 ri
z 1
1
ri
1 ri
z 1
z 2
" " i 负实轴上 " " i 0 正实轴上 N 3 N 1 1
呈偶对称
N -3
H
(
)
h
N
2
1
2 n0
2h(n)
cos
N
2
1
n
令 N 1 n m
2
N 1
h
N 2
1
2 m1
2h
N 2
1
m
cos(m
)
N 1 2
H () a(n)cos(n) n0
其中：
a(0)
h
N 1 2
a(n)
2h
数字信号处理第七章离散希尔伯特变换
离散希尔伯特变换与连续希尔伯特变换的关系
连续希尔伯特变换是离散希尔伯特变 换在时间上的连续形式，两者在数学 表达和性质上有很多相似之处。
离散希尔伯特变换是连续希尔伯特变 换的离散化，因此在应用上也有很多 相似之处，如信号分析、滤波器设计 等。
02
离散希尔伯特变换的基本原理
离散时间信号的表示
离散时间信号
快速算法的优缺点
优点是计算量较小，效率较高；缺点是需要对算法进行深入理解，实 现难度较大。
离散希尔伯特变换的软件实现
软件实现的基本思想
利用现有的软件库或编程语言，编写程序实现离散希尔伯特变换 的功能。
软件实现的步骤
首先选择合适的软件库或编程语言，然后编写程序实现离散希尔伯 特变换的功能，最后进行测试和验证。
02
通过傅里叶变换、快速傅里叶变换等算法对离散时间信号进行
频谱分析。
离散时间信号的频域特性
03
包括频率范围、频率分辨率、频率分辨率与采样频率的关系等。
离散希尔伯特变换的定义与性质
离散希尔伯特变换的定义
将一个实数序列通过一定的数学运算转换为复 数序列的过程。
离散希尔伯特变换的性质
包括线性性、时移性、频移性、共轭对称性等。
在时间上离散取值的信号，可以用序列的形式表示。
离散时间信号的数学描述
可以用离散时间实数序列或离散时间复数序列来表示。
离散时间信号的分类
根据取值是否连续，可以分为离散时间连续信号和离散时间离散信 号。
离散时间信号的频谱分析
离频率分量的大小和分布情况。
离散时间信号的频谱分析方法
离散希尔伯特变换的优缺点及未来发展方向 离散希尔伯特
变换的优缺点
01
数字信号处理(程佩青)
第七章 二维信号处理的一般方法§1 引言实践中不少信号是二维的，图象信号是一个典型的例子。
早期的图象处理技术采用的是信号处理的方法，其中不少技术至今仍广泛地应用着。
但是后来人们发现信号处理所得到的一个好的图片，并不一定能让人看着舒服。
这是因为人的视觉对图象的感受和信号处理中所采用的质量指标并不协调。
20世纪80年代，图象处理技术从采用信号处理的方法转向了采用人工智能、模式识别的方法，形成了一个新的技术领域——计算机视觉。
然而本章仍只讨论二维信号的信号处理方法。
这固然是因为由此而建立起来的许多技术还被广泛地应用着，另外也因为它是计算机视觉的研究基础。
二维信号可以通过扫描变成一维信号——电视信号就是一个典型例子——这种信号的处理，本质上仍是一维的，这里不再作讨论。
我们只讨论直接对二维信号进行处理的方法，它们是从一维的方法中推广过来的，但并不是所有的一维处理技术都能推广到二维中来。
这一点将在以后的讨论中予以说明。
本章讨论的内容是把一维信号处理中的时域和频域技术推广到二维中来。
本节则先把各种术语和变换推广过来。
无论是二维信号还是二维线性定常的系统，在时域里都表示成为一个二维序列f(1n ,2n )。
因此我们从介绍基本的二维序列开始我们的讨论。
1. 单位样本序列定义为：⎩⎨⎧===δ0n n 1 0)n ,n (2121其余 (7.1) 即仅在（0，0）点取1值而在其它点均为零的序列。
此序列作用到线性定常系统后的输出，即称为该系统的脉冲响应h(1n ,2n )。
这里的定常性是指无论冲激作用到哪一点（比如(1m ,2m )点），所得到的输出都是同形状的，只不过中心点的位置不同（由(1m ,2m )给定）罢了，即当输入为)m n ,m n (2211--δ时，输出为：)m n ,m n (h 2211--的系统称为定常的。
（附带说明一句：本教材一直采用“定常”这一术语，读者应明确：对一维连续系统它指的是time invariant （时不变），对一维离散系统它指的是Shift invariant （移不变），对二维系统指的是Space invariant （空间不变）——这因为对图象来说(1m ,2m )表示了空间点的位置）。
数字信号处理第7章
没复共轭；只有倒数。
例如，
Zi 1/ 2,1/ Zi 2
Zi
1/ Zi
01
1
2
2
（4） Zi 既在实轴上也在单位圆上。此时，
只有一个零点，且有两种可能，或位于Z=1， 或位于Z=-1。
Zi 1
N为偶数时的偶对称
H() 0, Z 1
为其零点；N为偶数奇对称 H（0）=0，有Z=1零点；
N为奇数奇对称H(0) H() 0, 有零点Z=1，和Z= -1。
Zi 1
§7-3 窗函数设计法
一、设计方法
1、设计思想
先给定理想filter的频响 Hd (ej)，所要求设计一个 FIR的filter的频响为 H(ej) ，使 H(ej) 逼近 Hd (ej)
2、设计过程
设计是在时域进行的，先用傅氏反变换求出理
想filter的单位抽样响应 hd (n) ，然后加时间窗 w(n)
c
2
N
时，随
增加，WR ( )
左边 旁瓣的起伏部分扫过通带，卷积 H()
也随着 WR ( ) 的旁瓣在通带内的面积
变化而变化，故 H()将围绕着零值而波动。
(6)当 c 2N时，WR ( )的右边旁瓣将进入
Hd ()的通带，右边旁瓣的起伏造成H() 值围绕 H(0) 值而波动。
H() / H(0)
理想低通filter的频响
为
1
0
为群延时
0
hd (n)
F1[Hd (e j)]
1 2
c eje jnd
c
1
2
e j(n)d 1
c
1 j2(n
)
e j(n)
c c
c sin[( n )c ] (n )c
例如，
Zi 1/ 2,1/ Zi 2
Zi
1/ Zi
01
1
2
2
（4） Zi 既在实轴上也在单位圆上。此时，
只有一个零点，且有两种可能，或位于Z=1， 或位于Z=-1。
Zi 1
N为偶数时的偶对称
H() 0, Z 1
为其零点；N为偶数奇对称 H（0）=0，有Z=1零点；
N为奇数奇对称H(0) H() 0, 有零点Z=1，和Z= -1。
Zi 1
§7-3 窗函数设计法
一、设计方法
1、设计思想
先给定理想filter的频响 Hd (ej)，所要求设计一个 FIR的filter的频响为 H(ej) ，使 H(ej) 逼近 Hd (ej)
2、设计过程
设计是在时域进行的，先用傅氏反变换求出理
想filter的单位抽样响应 hd (n) ，然后加时间窗 w(n)
c
2
N
时，随
增加，WR ( )
左边 旁瓣的起伏部分扫过通带，卷积 H()
也随着 WR ( ) 的旁瓣在通带内的面积
变化而变化，故 H()将围绕着零值而波动。
(6)当 c 2N时，WR ( )的右边旁瓣将进入
Hd ()的通带，右边旁瓣的起伏造成H() 值围绕 H(0) 值而波动。
H() / H(0)
理想低通filter的频响
为
1
0
为群延时
0
hd (n)
F1[Hd (e j)]
1 2
c eje jnd
c
1
2
e j(n)d 1
c
1 j2(n
)
e j(n)
c c
c sin[( n )c ] (n )c
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第七章数字滤波器设计7.1：无限冲激响应滤波器的阶数的估计Q7.1用MATTAB确定一个数字无限冲激响应低通滤波器所有四种类型的最低阶数。
指标如下:40 kHz的抽样率，,4 kHz的通带边界频率，8 kHz的阻带边界频率，0.5 dB的通带波纹，40 dB的最小阻带衰减。
评论你的结果。
答：标准通带边缘角频率Wp是：标准阻带边缘角频率Ws是：理想通带波纹Rp是0.5dB理想阻带波纹Rs是40dB1.使用这些值得到巴特沃斯低通滤波器最低阶数N=8，相应的标准通带边缘频率Wn是0.2469.2.使用这些值得到切比雪夫1型低通滤波器最低阶数N=5，相应的标准通带边缘频率Wn是0.2000.3/使用这些值得到切比雪夫2型低通滤波器最低阶数N=5，相应的标准通带边缘频率Wn是0.4000.4.使用这些值得到椭圆低通滤波器最低阶数N=8，相应的标准通带边缘频率Wn是0.2000. 从以上结果中观察到椭圆滤波器的阶数最低，并且符合要求。
Q7.2用MATLAB确定一个数字无限冲激响应高通滤波器所有四种类型的最低阶数。
指标如下:3500Hz的抽样率，1050 Hz的通带边界频率,600 Hz的阻带边界频率，1 dB 的通带波纹，50 dB的最小阻带衰减。
评论你的结果答：标准通带边缘角频率Wp是：标准阻带边缘角频率Ws是：理想通带波纹Rp是1dB理想阻带波纹Rs是50dB1.使用这些值得到巴特沃斯高通滤波器最低阶数N=8，相应的标准通带边缘频率Wn是0.5646.2.使用这些值得到切比雪夫1型高通滤波器最低阶数N=5，相应的标准通带边缘频率Wn是0.6000.3.使用这些值得到切比雪夫2型高通滤波器最低阶数N=5，相应的标准通带边缘频率Wn是0.3429.4.使用这些值得到椭圆低通滤波器最低阶数N=4，相应的标准通带边缘频率Wn是0.6000.从以上结果中观察到椭圆滤波器的阶数最低，并且符合要求。
Q7.3用MATLAB确定一个数字无限冲激响应带通滤波器所有四种类型的最低阶数。
指标如下:7 kHz的抽样率,1.4 kHz和2.1 kHz的通带边界频率，1.05 kHz和2.45 kHz的阻带边界频率，,0 .4 dB的通带波纹,50 dB的最小阻带衰减。
评论你的结果。
答：标准通带边缘角频率Wp是：标准阻带边缘角频率Ws是：理想通带波纹Rp是0.4dB理想阻带波纹Rs是50dB1.使用这些值得到巴特沃斯带通滤波器最低阶数2N=18，相应的标准通带边缘频率Wn是[0.3835 0.6165].2.使用这些值得到切比雪夫1型带通滤波器最低阶数2N=12，相应的标准通带边缘频率Wn 是[0.4000 0.6000].3.使用这些值得到切比雪夫2型带通滤波器最低阶数2N=12，相应的标准通带边缘频率Wn 是[0.3000 0.7000].4.使用这些值得到椭圆带通滤波器最低阶数2N=8，相应的标准通带边缘频率Wn是[0.4000 0.6000].从以上结果中观察到椭圆滤波器的阶数最低，并且符合要求。
Q7.4用MATLAB确定一个数字无限冲激响应带阻滤波器所有四种类型的最低阶数。
指标如下:12 kHz的抽样率，2.1 kHz和4.5 kHz的通带边界频率,2.7 kHz和3.9 kHz的阻带边界频率，0.6 dB的通带波纹,45 dB的最小阻带衰减。
评论你的结果。
答：标准通带边缘角频率Wp是：标准阻带边缘角频率Ws是：理想通带波纹Rp是0.6dB理想阻带波纹Rs是45dB1.使用这些值得到巴特沃斯带阻滤波器最低阶数2N=18，相应的标准通带边缘频率Wn是[0.3873 0.7123].2.使用这些值得到切比雪夫1型带阻滤波器最低阶数2N=10，相应的标准通带边缘频率Wn 是[0.3500 0.7500].3.使用这些值得到切比雪夫2型带阻滤波器最低阶数2N=10，相应的标准通带边缘频率Wn 是[0.4500 0.6500].4.使用这些值得到椭圆带阻滤波器最低阶数2N=8，相应的标准通带边缘频率Wn是[0.3500 0.7500].从以上结果中观察到椭圆滤波器的阶数最低，并且符合要求。
7.2：无限冲激响应滤波器设计程序P7.1说明巴特沃斯带阻滤波器的设计。
% 巴特沃斯带阻滤波器的设计Ws = [0.4 0.6]; Wp = [0.2 0.8]; Rp = 0.4; Rs = 50;% 估计滤波器的阶数[N1, Wn1] = buttord(Wp, Ws, Rp, Rs);% 设计滤波器[num,den] = butter(N1,Wn1,'stop');% 显示传输函数disp('分子系数是 ');disp(num);disp('分母系数是 ');disp(den);% 计算增益响应[g, w] = gain(num,den);% 绘制增益响应plot(w/pi,g);gridaxis([0 1 -60 5]);xlabel('\omega /\pi'); ylabel('增益， dB');title('巴特沃斯带阻滤波器的增益响应');Q7.5通过运行程序P7. 1来设计巴特沃兹带阻滤波器。
写出所产生的传输函数的准确表达式。
滤波器的指标是什么，你的设计符合指标吗，使用MATLAB，计算并绘制滤波器的未畸变的相位响应及群延迟响应。
答：表达式是：滤波器参数是：Wp1=0.2π，Ws1=0.4π，Ws2=0.6π，Wp2=0.8π，Rp=0.4dB,Rs=50dB.设计的滤波器增益响应如下：从图中可以总结出设计符合指标。
滤波器的未畸变的相位响应及群延迟响应如下：Q7.6修改程序P7.1来设计符合习题Q7.1所给指标的切比雪夫1型低通滤波器。
写出所产生的传输函数的准确表达式。
你的设计符合指标吗?使用MATLAB,计算并绘制滤波器的未畸变的相位响应及群延迟响应。
答：表达式如下：设计的滤波器增益响应如下：从图中可以总结出设计符合指标。
滤波器的未畸变的相位响应及群延迟响应如下：Q7.7修改程序P7.1来设计符合习题Q7.2所给指标的切比雪夫2型高通滤波器。
写出所产生的传输函数的准确表达式。
你的设计符合指标吗?使用MATLAB,计算并绘制滤波器的未畸变的相位响应及群延迟响应。
答：表达式如下：设计的滤波器增益响应如下：从图中可以总结出设计符合指标。
滤波器的未畸变的相位响应及群延迟响应如下：Q7.8修改程序P7.1来设计符合习题Q7.3所给指标的椭圆带通滤波器。
写出所产生的传输函数的准确表达式。
你的设计符合指标吗，使用MATLAB，计算井绘制滤波器的未畸变的相位响应及群延迟响应。
答：表达式如下：设计的滤波器增益响应如下：从图中可以总结出设计符合指标。
滤波器的未畸变的相位响应及群延迟响应如下：7.3：吉布斯现象Q7.9使用函数sinc编写一个MATLAB程序，以产生截止频率在Wc= 0.4π处、长度分别为81,61,41和21的四个零相位低通滤波器的冲激响应系数，然后计算并画出它们的幅度响应。
使用冒号“:”运算符从长度为81的滤波器的冲激响应系数中抽出较短长度滤波器的冲激响应系数。
在每一个滤波器的截止频率两边研究频率响应的摆动行为。
波纹的数量与滤波器的长度之间有什么关系?最大波纹的高度与滤波器的长度之间有什么关系?你将怎样修改上述程序以产生一个偶数长度的零相位低通滤波器的冲激响应系数?答：长度为81时幅度响应如下：长度分别为61,41和21的幅度响应如下：从中可以观察到由于吉布斯现象产生的幅度响应的摆动行为。
波纹的数量与滤波器的长度之间的关系——波纹的数量减少与长度成正比。
最大波纹的高度与滤波器的长度之间的关系——最大波纹的高度与长度无关。
Q7.10使用函数sinc编写一个MATLAB程序，以产生一个截止频率在Wc= 0.4π处、长度为45的零相位高通滤波器的冲激响应系数，计算并画出其幅度响应。
在每一个滤波器的截止频率两边研究频率响应的摆动行为。
你将怎样修改上述程序以产生一个偶数长度的零相位高通滤波器的冲激响应系数?答：长度为45时幅度响应如下：从中可以观察到由于吉布斯现象产生的幅度响应摆动行为。
在这种情况下你不能改变长度。
原因：这是一个零相位滤波器，这意味着它也是一个线性相位滤波器，因为零相是一种特殊的线性相位的子集。
现在，理想的有限脉冲响应长度甚至有对称的中点h[n]。
使其成了一个线性相位FIR滤波器。
二型滤波器不可能是高通滤波器，因为必须在z=-1处有零点，意味着w=+-π。
Q7.11编写一个MATLAB程序，以产生长度分别为81,61,41和21的四个零相位微分器的冲激响应系数，计算并画出它们的幅度响应。
下面的代码段显示了怎样产生一个长度为2M+1的微分器。
n=1:M;b=cos(pi*n)./n;num=[-fliplr(b) 0 b];对于每种情况,研究微分器的频率响应的摆动行为。
波纹的数量与微分器的长度之间有什么关系，最大波纹的高度与滤波器的长度之间有什么关系?答：幅度响应分别如下：从中可以观察到由于吉布斯现象产生的幅度响应的摆动行为。
波纹的数量与微分器的长度之间的关系——两者成正比。
最大波纹的高度与滤波器的长度之间的关系——两者间没有关系。
Q7.12编写一个MA11AB程序，以产生长度分别为81,61.41和21的四个离散时间希尔伯特变换器的冲激响应系数，计算并画出它们的幅度响应。
下面的代码段显示了怎样产生一个长度为2M十1的希尔伯特变换器。
n=1:M;c=sin(pi*n)./2;b=2*(c.*c)./(pi*n);num=[-fliplr(b) 0 b];对于每种情况，研究希尔伯特变换器的频率响应的摆动行为。
波纹的数量与希尔伯特变换器的长度之间有什么关系?最大波纹的高度与滤波器的长度之间有什么关系?答：幅度响应如下：从中可以观察到由于吉布斯现像产生的幅度响应的摆动行为。
波纹的数量与希尔伯特变换器的长度之间的关系——两者成正比。
最大波纹的高度与滤波器的长度之间的关系——两者无关系。
7.4：有限冲激响应滤波器的阶数估计Q7.13 线性相位低通FIR滤波器的阶数估算，参数如下:p =2 kHz, s =2.5 kHz, p = 0.005, s = 0.005, F T = 10kHz使用 kaiord 的结果为N = 46使用 ceil 命令的目的是朝正方向最接近整数方向取整。
使用nargin命令的目的是表明函数M文件体内变量的数目。
Q7.14 (a)线性相位FIR滤波器的阶数估算，其中采样频率改为F T = 20 kHz ，则结果为 N=91。
(b) 线性相位FIR滤波器阶数的估计，其中通带波纹改成p= 0.002和s= 0.002 结果为 N=57。