山东大学数字信号处理课程试题答案(A卷)

2
,α =
N −1 ，故该 FIR 滤波器具有线性相位特性。 2
6
证明： 冲激响应不变法： 使数字滤波器 h(n)等于原模拟滤波器冲激响应 ha (t ) 的等间隔抽样 ha (nT ) H(z)与H(s)的关系：H(z)=Z[h(n)]=Z[£-1[Ha(s)t=nT]] 分解： H ( s ) =
∑ h( n) z − n =
∞
∞
n = −∞
∑ h (nT ) z
a N
∞
−n
＝
∑ ∑ Ak e sk nT z −n = ∑ Ak ] [∑ (e skT z −1 ) n ] = ∑
n = 0 k =1 k =1 n =0
∞
N
N
Ak sk T −1 z k =1 1 − e
（2）
7
即得证。 答：IIR 滤波器设计的基本方法是利用复值映射将大家熟知的模拟滤波器变换为数字滤 波器。 基本步骤： （1）设计模拟低通滤波器 （2）研究并实行滤波器变换以得到数字低通滤波器 （3）研究并实行频带变换以便从数字低通滤波器得到其他数字滤波器 以巴特滤波器为例， 由参数 Ω p , R p , Ω s As 给出， 实质是为了求其阶次 N 和截止频率 Ω c （1） Ω = Ω p ⇒ R p = −10 lg[
− wc )
] = −(e − jwc + e jwc )
（1）
两边同乘以 e
得：
⇒
′ wc + wc cos( ) 2 α =− ′ w − wc cos( c ) 2
（2）
即：由低通到高通的公式：
′ wc + wc cos( ) z −1 + α 2 z −1 = − ， α = − ′ 1 + αz −1 wc − wc cos( ) 2
n =0 − ( N −1− n )
( N −1) / 2
=
∑ h( n ) z − n –
n =0
( N − 3) / 2
由于 h(n)是奇对称的，故 h (
( N −3) / 2
N −1 )=0 2
− ( N −1− n )
H(z)=
∑ h( n ) z
n =0 ( N −3) / 2
−n
( N − 3) / 2
′
z −1 − α ′ 得：G(1)=–1,G(–1)=1, wc 等效于高通滤波器的 π + wc ，而 −1 1 − αz
– wc 等效为高通滤波器的 wc ，则：
′
e − j ( − wc ) = −
′
e − jwc + α 1 + αe − jwc
c '
⇒
α [1 + e j ( w
j − ( wc ' − wc ) 2
答：
4
IIR:用较少的阶数获得很高的选择性，存储单元少，运算次数少，效率高，代价 为相位非线性 FIR：严格线性相位，较多的阶数，效率低 （2 分） （1） 若按相同的选择性和相同的相位线性要求， IIR 加全通网络进行相位校正， 要大大增加滤波器的节数和复杂性。 （2） 若要求严格线性相位，FIR 性能上和经济上优于 IIR。 （3） 从设计上看， IIR 借助于模拟滤波器的结果可利用现成的设计公式来设计， 成熟方便。 （4） FIR 无闭式解，用窗口函数设计，大多数设计都用迭代法，实现时对硬件 要求高。 (8 分) 证 明 ： FFT 算 法 ： 将 一 个 N 点 的 DFT 运 算 分 解 ， 不 妨 设
∑s−s
k =1
N
Ak
，单极点多项式和
k
N
对 H a ( s ) 进行拉式反变换，得： ha (t ) =
∑A e
k =1 k N k =1
sk t
u (t )
h(n) = ha (nT ) = ∑ Ak e sk nT u (nT )
对 h(n) 求 Z 变换： H ( z ) =
（1）
n = −∞
N −1
N −1
n 为偶数
N / 2 −1
n 为奇数
N / 2 −1 r =0
＝
∑ g (r )WN2rk ＋
r =0
∑ h(r )W
k N / 2 −1 r =0
( 2 r +1) k N
N / 2 −1
＝
∑ g (r )W
r =0
−j 2π ⋅2 r N
2 rk N
＋ WN
2π ⋅r N /2
∑ h(r )W
∑ h( n ) z − n +
n = ( N +1) / 2
∑ h( n) z
− ( N −1− n )
N −1
−n
令第二项：n′=N-1-n ，得：
( N −1) / 2
H(z)
=
∑ h( n ) z − n +
n =0
( N −3) / 2 n =0
∑ h( N − 1 − n) z
∑ h( n) z
N −1 2
N −1 2
（2）
X(
N ) = ∑ x( n )WN 2 n =0 =
N −1 2 r =0
N −1
j
N n 2
= ∑ x( n )( −1) n =
n=0 N −1 2 r =0
N −1
∑ x(2r ) − ∑ x(2r + 1)
r =0 r =0
N −1 2
N −1 2
∑ x( N − 1 − 2r ) − ∑ x(2r + 1) ∑ x(2k + 1) − ∑ x(2r + 1) =0
–jw(N-1) / 2
− jwn
− e − jw( N −1− n ) ]
jw ( N −1 −n) 2
N / 2 −1 n =0
=e
∑ h(n)[e
−e
− jw (
N −1 −n ) 2
]
=e –jw(N-1) / 2
π N −1 j ( −w ) 2 2
N / 2 −1 n =0
∑ h(n) ⋅ 2 j sin[w(
–
∑ h( n) z
n =0
H(e
jw
)=
∑ h(n)[e
n =0
− jwn
− e − jw( N −1− n ) ]
jw ( N −1 −n ) 2 − jw ( N −1 −n ) 2
( N −3) / 2
= e –jw(N-1) / 2
∑
n =0
h(n)[e
−e
]
( N −3) / 2
= e –jw(N-1) / 2
∑1−W
H(
N −1
zk 1− r N z −N )= r N
∑ 1 − rW
k =0
N −1
H (k )
−k N
z −1
+
N / 2 −1 H (k ) = ∑ ∑ − k −1 k =0 k = 0 1 − rW N z
H (k ) 1 − re
2πk j N
z −1
k=N / 2
∑
N −1
H (k ) 1 − re
r ( N / 2+ k )
N / 2 −1 N / 2 −1 r (k + ) N rk ∴ G ( k + ) = ∑ g ( r )W N / 2 2 = ∑ g ( r )W N / 2 = G (k ) 2 r =0 r =0
N
同理： H (k +
( N +k )
N ) = H (k ) 2
N / 2 −1
− ( N −1− n )
N / 2 −1
=
∑ h( n) z
n =0
−n
–
∑ h( n ) z
n =0
−( N −1− n )
N / 2 −1
=
∑ h(n)[ z
n =0
−n
− z −( N −1− n ) ]
将z=e jw代入，得： H(e jw)=
N / 2 −1 n =0
∑ h(n)[e
1+ (
1 Ωp Ωc
] )
2N
（2） Ω = Ω s ⇒ As = −10 lg[
1 ] Ωs 2N 1+ ( ) Ωc
−1
⇒ N , Ω c 的值即可。
8
解：由已知，低通到高通滤波器的变换： z
=−
z −1 + α ， 1 + αz −1
设 wc 为原型低通滤波器的截止频率， wc 为高通滤波器的截至频率 由 G( z) = −
N = 2 M , x(2r ) = g (r ), x(2r + 1) = h(r ), r=0,1,…
则：x(n)的 DFT 可分解为：
nk X (k ) = DFT [ x(n)] = ∑ x(n)WN n =0 N −1
N −1 2
＝
∑ x(n)WNnk ＋ ∑ x(n)WNnk
n =0 n =0
1 + 1 / 3z −1 1 • −1 1 − 1/ 4z 1 − 1 / 2 z −1
z –1 1/3 z -1 1/2
1/4
d)
并联型：H(z)=
−7/3 10 / 3 + −1 1 − 1/ 4z 1 − 1 / 2 z −1
-7/3 z
–1
x(n)
1/4 10/3 z 1/2
–1
y(n)
3
k =0 r =0 N −1 2 N −1 2
N－1－2 r = 2 k +1 ⎯令 ⎯ ⎯⎯⎯ ⎯→ =