数字信号处理第四章附加题

合集下载

数字信号处理-第4章习题

数字信号处理-第4章习题

习题四4.1 用窗函数法设计一线性相位FIR 高通滤波器,通带边界频率为0.6π,阻带边界频率为0.4π,要求阻带衰减不小于40dB ,窗函数从矩形窗、汉宁窗、海明窗和布莱克曼窗中选取,且要求滤波器的阶数最小。

求出该滤波器的单位脉冲响应h (n )的解析式。

解 求理想高通滤波器的边界频率n ω和过渡带宽ω∆:0.60.40.522c rn ωωππωπ++≈== 0.60.40.2c r ωωωπππ∆=-=-=选择窗函数和窗口长度。

阻带衰减不小于40dB ,因此选择汉宁窗。

根据过渡带宽求窗口长度N 和线性相位延迟常数α:6.20.231N Nππ=⇒= 1152N α-== 根据理想边界频率n ω和线性相位延迟常数α,求理想单位脉冲响应d ()h n :[][]()()d sin ()sin ()1()()d d 21n n n j n j n n n n n n h n e e n ωπωαωαπωαπαωαπαωωπωαπ----⎧---≠⎪⎪-⎡⎤=+=⎨⎢⎥⎣⎦⎪-=⎪⎩⎰⎰ 窗函数与理想单位脉冲响应相乘,即可得到线性相位FIR 高通滤波器的单位脉冲响应:[][]sin (15)sin (15)0.50.50.5cos ()15()15(15)0.515N n n n R n n h n n n ππππ⎧---⎡⎤⎛⎫-⋅⋅≠⎪ ⎪⎢⎥=-⎝⎭⎨⎣⎦⎪=⎩0,1,2,,30n =⋅⋅⋅6kHz ,阻带边界频率为2kHz 和8kHz ,采样频率为20kHz ,要求阻带衰减不小于50dB ,窗函数从矩形窗、汉宁窗、海明窗和布莱克曼窗中选取,且要求滤波器的阶数最小。

求出该滤波器的单位脉冲响应h (n )的解析式。

解 求理想带通滤波器的边界频率1ω、2ω和过渡带宽ω∆:12121212112212221120.4 , 20.620.2 , 20.80.3 , 0.7220.2c c c c s s r r r r s sc r c r r c c r f ff f f ff f ωππωππωππωππωωωωωπωπωωωωωπ========++≈=≈=∆=-=-= 选择窗函数和窗口长度。

数字信号处理第四章习题

数字信号处理第四章习题

第四章习题4.1 (a) By expanding the equation()()[]()⎥⎦⎤⎢⎣⎡==⎰--∞→∞→2200021T T Ft j T xx T xx dt e t x T E lim F P E lim F 00πΓ taking the expected value, and finally taking the limit as ∞→0T ,show that the right-hand side converges to )(f xx Γ.(b) Prove that2102211)(1)(∑∑-=---+-==N n fn j fm j N N m xx en x N e m r ππ.4.2 For zero-mean, jointly Gaussian random variables, X 1, X 2, X 3, X 4, itis well known that)()()()()()()(3241423143214321X X E X X E X X E X X E X X E X X E X X X X E ++=. Use this result to derive the mean-square value of ()m r xx and the variance, given by()[][]()()()[]∑∞-∞=+-+-≈n xx xx xx xx m n m n n m N N m r γγγ*22varwhich is defined as[][][]22(()(var m r E m r E m r xx xx xx -=. 4.3 By use of the expression for the fourth joint moment for Gaussianrandom variables, show that(a)()()[]⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧⎥⎦⎤⎢⎣⎡--+⎥⎦⎤⎢⎣⎡+++=2212122121421)(sin )(sin )(sin )(sin 1f f N N f f f f N N f f f P f P E x xx xx ππππσ (b)[]⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧⎥⎦⎤⎢⎣⎡--+⎥⎦⎤⎢⎣⎡++=2212122121421)(sin )(sin )(sin )(sin )()(cov f f N N f f f f N N f f f P f P x xx xx ππππσ(c)[]⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧⎪⎪⎭⎫ ⎝⎛+=242sin 2sin 1)(var f N fN f P x xx ππσ under the condition that the sequence ()n x is a zero-mean white Gaussian noise sequence with variance 2x σ.4.4 Generalize the results in Problem 4.3 to a zero-mean Gaussian noiseprocess with power density spectrum )(f xx Γ, as given by()[]()⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫ ⎝⎛+Γ=222sin 2sin 1var f N fN f f P xx xx ππ (Hint: Assume that the colored Gaussian noise process is the output of a linear system excited by white Gaussian noise.)4.5 Show that the periodogram values at frequencies,1,1,0,/-==L k L k f k given by (4.1.35), can be computed by passing the sequence through a bank of L IIR filters, where each filter has an impulse response )()(/2n u e n h N nk j k π-= and then computing the magnitude-squared value of the filter outputs at n=N. Note that each filter has a pole on the unit circle at the frequency f k .4.6 The Bartlett method is used to estimate the power spectrum of asignal x(n). We know that the power spectrum consists of a single peak with a 3 dB bandwidth of 0.01 cycle per sample, but we do not know the location of the peak.(a) Assuming that N is large, determine the value of M=N/K so thatthe spectral window is narrower than the peak.(b) Explain why it is not advantageous to increase M beyond thevalue obtained in part (a).4.7 The N-point DFT of a random sequence x(n) is ∑-=-=10/2)()(N n N nk j e n x k X π.Assume that E[x(n)]=0 and E[x(n)x(n+m)]=)(2m w δσ (in other words,x(n) is a white noise process).(a) Determine the variance of X(k).(b) Determine the autocorrelation of X(k).4.8 An AR(2) process is described by the difference equation)()2(81.0)(n n x n x ω+-=, where w(n) is a white noise process withvariance 2ωσ.(a) Determine the parameters of the MA(2), MA(4), and MA(8)models that provide a minimum mean-sequare error fit to thedata x(n).(b) Plot the true spectrum and those of the MA (q), q=2,4,8spectra and compare the results. Comment on how well theMA(q) models approximate the AR (2) process.4.9 An MA (2) process is described by the difference equation )2(81.0)()(-+=n n n x ωω, where w(n) is a white noise process withvariance 2ωσ.(a) Determine the parameters of the AR(2), AR(4), and AR(8)models that provide a minimum mean-square error fit to the data x(n).(b) Plot the true spectrum and those of the AR(p), p=2,4,8, andcompare the results. Comment on how well the AR(p) modelsappoximate the MA (2) process.4.10 The autocorrelation sequence for an AR process x(n) ismxx m ⎪⎭⎫ ⎝⎛=41)(γ (a) Determine the difference equation for x(n)(b) Is your answer unique? If not, give any other possiblesolutions.4.11 Suppose that we represent an ARMA(p,q) process as a cascade ofan MA(q) followed by an AR(p) model. The input-output equation for the MA(q) model is ∑=-=qk k k n w b n v 0)()(, where w(n) is a whitenoise process. The input-output equation for the AR(p) model is∑==-+pk k n v k n x a n x 1)()()((a) By computing the autocorrelation of v(n), show thatq m d b b m mq k m w m k k w vv ≤≤==∑-=+0)(022σσγ(b) Show that 1)()(00=+=∑=a k m a m pk vx k vv γγ4.12 Suppose that the AR(2) process in Problem 4.8 is corrupted by anadditive white noise process v(n) with variance 2v σ. Thus, we havey(n)=x(n)+v(n)(a) Determine the difference equation for y(n) and thusdemonstrate that y(n) is an ARMA(2,2) process. Determinethe coefficients of the ARMA process.(b) Generalize the result in part (a) to an AR(p) process∑=+--=pk k n w k n x a n x 1)()()( and )()()(n v n x n y +=.4.13 The harmonic decomposition problem considered by Pisarenko maybe expressed as the solution to the equationa a a Γa H w yy H 2σ=The solution for a may be obtained by minimizing the quadratic form a Γa yy H subject to the constraint that a a H =1. The constraint can be incorporated into the performance index by means of a Lagrange multiplier. Thus the performance index becomes()a a a Γa H yy H 1-+=λζ.By minimizing ζ with respect to a , show that this formulation is equivalent to the Pisarenko eigenvalue problem given in (4.4.9), with the Lagrange multiplier playing the role of the eigenvalue. Thus,show that the minimum of ζ is the minimum eigenvalue 2w σ.4.14 The autocorrelation of a sequence consisting of a sinusoid withrandom phase in noise is)(2cos )(21m m f P m w xx δσπγ+=where 1f is the frequency of the sinusoidal, P its power, and 2w σthe variance of the noise. Suppose that we attempt to fit an AR(2) model to the data.(a) Determine the optimum coefficients of the AR(2) model as afunction of 2w σ and 1f .(b) Determine the reflection coefficients 1K and 2K correspondingto the AR(2) model parameters.(c) Determine the limiting values of the AR(2) parameters and (1K ,2K )as 02→w σ.4.15 This problem involves the use of cross-correlation to detect a signalin noise and estimate the time delay in the signal. A signal x(n) consists of a pulsed sinusoid corrupted by a stationary zero-mean white noise sequence. That is, 10),()()(0-≤≤+-=N n n w n n y n x ,where )(n w is the noise with variance 2w σ and the signal is⎩⎨⎧-≤≤=otherwise M n n A n y ,010,cos )(0ω. The frequency 0ω is known, but the delay 0n , which is a positiveinteger, is unknown, and is to be determined by cross-correlating x(n) with y(n). Assume that 0n M N +>. Let∑-=-=10)()()(N n xy n x m n y m rdenote the cross-correlation sequence between x(n) and y(n). In the absence of noise, this function exhibits a peak at delay 0n m =. Thus,0n is determined with no error. The presence of noise can lead toerrors in determining the unknown delay.(a) For 0n m =, determine ()[]0n r E xy . Also, determine thevariance ()[]0var n r xy , due to the presence of the noise. In bothcalculations, assume that the double-frequency term averages to zero. That is, 0/2ωπ>>M .(b) Determine the signal-to-noise ratio, defined as []{}[])(var )(020n r n r E SNR xy xy = (c) What is the effect of the pulse duration M on the SNR?。

数字信号处理习题及答案4

数字信号处理习题及答案4

习题及答案4一、填空题(每空1分, 共10分)1.序列()sin(3/5)x n n π=的周期为 .2.线性时不变系统的性质有 律、 律、 律. 3.对4()()x n R n =的Z 变换为 ,其收敛域为 。

4.抽样序列的Z 变换与离散傅里叶变换DFT 的关系为 。

5.序列x(n )=(1,-2,0,3;n=0,1,2,3), 圆周左移2位得到的序列为 。

6.设LTI 系统输入为x (n ) ,系统单位序列响应为h (n ),则系统零状态输出y (n )= . 7.因果序列x (n ),在Z →∞时,X(Z )= 。

二、单项选择题(每题2分, 共20分)1.δ(n)的Z 变换是 ( )A.1 B 。

δ(ω) C 。

2πδ(ω) D 。

2π2.序列x 1(n)的长度为4,序列x 2(n )的长度为3,则它们线性卷积的长度是 ( )A. 3 B. 4 C 。

6 D 。

73.LTI 系统,输入x (n)时,输出y (n );输入为3x (n —2),输出为 ( ) A. y (n-2) B.3y (n —2) C 。

3y (n ) D.y(n )4.下面描述中最适合离散傅立叶变换DFT 的是 ( )A.时域为离散序列,频域为连续信号B 。

时域为离散周期序列,频域也为离散周期序列C 。

时域为离散无限长序列,频域为连续周期信号D 。

时域为离散有限长序列,频域也为离散有限长序列5.若一模拟信号为带限,且对其抽样满足奈奎斯特条件,理想条件下将抽样信号通过 即可完全不失真恢复原信号 ( )A.理想低通滤波器 B 。

理想高通滤波器 C 。

理想带通滤波器 D 。

理想带阻滤波器6.下列哪一个系统是因果系统 ( )A.y(n)=x (n+2) B 。

y (n )= cos (n+1)x (n ) C 。

y (n )=x (2n ) D.y (n)=x (— n )7.一个线性时不变离散系统稳定的充要条件是其系统函数的收敛域包括 ( ) A 。

数字信号处理 第四章04

数字信号处理 第四章04

同理:∀n1 , 0 ≤ n ≤ L − 1
详见(4-38) P.142
L
M
例:N=12=4×3, M=4 , L=3 算法流图:图4-20,P.144
这是一个蝶形 三点蝶形
这仍然是一个蝶形 四点蝶形
x(n)={x(0), x(1), x(2)…x(11)}
x(0) x(1) x(5) x(9)
= ∑ x( n)WNkn = = =
=
n=0 M −1 L −1
n0 = 0 n1 = 0 M −1 L −1
N −1
x(Mn1+n0)
一列一列 求DFT
= DFTn1 [ x( n1 , n0 )] 0 ≤ k0 ≤ L − 1, ∀n0

n1 = 0
( Mn1 + n0 )( Lk1 + k0 ) x ( n , n ) W ∑∑ 1 0 N Mn1k0 Lk1n0 k0 n0 MLk1n1 x ( n , n ) W W W W ∑∑ 1 0 N N N N
n1 = 0 L −1
n1 = 0,1,..., L − 1
行号
k0 = 0,1,..., L − 1
′ k 0 n0 (3) X 1 (k0 , n0 ) = X 1 (k0 , n0 )WN
0 ≤ k0 ≤ L − 1
0 ≤ n0 ≤ M − 1 (4) ∀k0 , 0 ≤ k0 ≤ L − 1 (针对每一行) M −1 ′ kn ′ X 2 (k 0 , k1 ) = DFTn0 [ X 1 (k 0 , n0 )] = ∑ X 1 (k0 , n0 )WM , k 0 = 0,1,..., M − 1
0 ≤ n0 ≤ M − 1, ∀k

数字信号处理课后习题答案(全)1-7章

数字信号处理课后习题答案(全)1-7章
因此系统是非线性系统。
第 1 章 时域离散信号和时域离散系统
(6) y(n)=x(n2)
令输入为
输出为
x(n-n0)
y′(n)=x((n-n0)2) y(n-n0)=x((n-n0)2)=y′(n) 故系统是非时变系统。 由于
T[ax1(n)+bx2(n)]=ax1(n2)+bx2(n2) =aT[x1(n)]+bT[x2(n)]
第 1 章 时域离散信号和时域离散系统
题4解图(一)
第 1 章 时域离散信号和时域离散系统
题4解图(二)
第 1 章 时域离散信号和时域离散系统
题4解图(三)
第 1 章 时域离散信号和时域离散系统
(4) 很容易证明: x(n)=x1(n)=xe(n)+xo(n)
上面等式说明实序列可以分解成偶对称序列和奇对称序列。 偶对称序列可 以用题中(2)的公式计算, 奇对称序列可以用题中(3)的公式计算。
(2) y(n)=x(n)+x(nN+1)k 0
(3) y(n)= x(k)
(4) y(n)=x(n-nn0)n0
(5) y(n)=ex(n)
k nn0
第 1 章 时域离散信号和时域离散系统
解:(1)只要N≥1, 该系统就是因果系统, 因为输出 只与n时刻的和n时刻以前的输入有关。
如果|x(n)|≤M, 则|y(n)|≤M, (2) 该系统是非因果系统, 因为n时间的输出还和n时间以 后((n+1)时间)的输入有关。如果|x(n)|≤M, 则 |y(n)|≤|x(n)|+|x(n+1)|≤2M,
=2x(n)+x(n-1)+ x(n-2)
将x(n)的表示式代入上式, 得到 1 y(n)=-2δ(n+2)-δ(n+1)-0.5δ(2n)+2δ(n-1)+δ(n-2)

北京邮电大学《数字信号处理》门爱东-数字信号处理第四章附加习题.docx

北京邮电大学《数字信号处理》门爱东-数字信号处理第四章附加习题.docx

5. 一个线性吋不变因果系统由下列差分方程描述y(n) = x(n)-x (AZ-l)-0.5y(n-l)(1) 系统函数H(Z),判断系统屈于FIR 和11R 屮的哪一类以及它的滤波特性。

第四章附加题 1.请推导出三阶巴特沃思低通滤波器的系统函数,设Sl c =\rad!s o 2.设计一个满足下列指标的模拟Butterworth 低通滤波器,要求通带的截止频率 f p = 6kHz,,通带最大衰减Ap=3dB,,阻带截止频率X = 12kHz,,阻带的最 小衰减4 = 25dB ,求出滤波器的系统函数。

3. 设计一个模拟切比雪夫低通滤波器,要求通带的截止频率f P =3kHz,通带衰 减要不大于0.2dB,阻带截止频率f s = 12kHz,阻带衰减不小于50dBo 4. 数字滤波器经常以下图描述的方式来处理带限模拟信号。

在理想情况下,模 —数变换器把模拟信号取样,产生序列x(n) = x a (nT),而数一模变换器乂将 取样)心)变成限带波形 D sin (龙/T)(—M) (”)加)(,") 整个系统等效于一个线性吋不变模拟系统。

如果系统/?(〃)的截止频率是Tr/Srad/s, i/T = iOkHz,等效模拟滤波器 的截止频率是多少? "=一8 模-数变换器 心)、 2 (周期T) 2 1/T = 20kHz ,重复(1)。

h (n) 数-模变换器 儿⑴、 (周期T) 2y 何〉⑵ 若输入x©) = 2cos(0.5勿)+ 5(7?>0),求系统稳态输岀的最大幅值。

6. 设% (/)表示一模拟滤波器的单位冲激响应,用冲激响应不变法,将此模拟滤波器转化成数字滤波器(/2(町表示单位取 样响应,即h(n) = Th a (nT)\确定系统函数H(z),并把T 作为参数,T 为任何 值时,数字滤波器是稳定的,并说明数字滤波器近似为低通滤波器还是高通滤波 器。

《现代数字信号处理》第4章习题答案

《现代数字信号处理》第4章习题答案

(a)试求
AR(2)模型的系数 a2
=
⎡⎣1, a2 (1), a2 (2)T
⎤⎦
(表示为 w0 ,
σ
2 w

P
的函数形式。)
(b)求AR(2)模型对应的反射系数Γ1和Γ2。
(c)当 σ
2 w

0
时,AR(2)参数和反射系数的极限值是多少?
解:(a)
rx (0) =
P
+
σ
2 w
,
rx (1) =
P cosω0,
{ } E
ei− (n) x∗ (n − k )
=
E
⎧⎪⎡ ⎨⎢
x
(
n
− i) +
i
∑ ai∗
(
j)
x(n
−i
+
j)⎤⎥ ⋅
x∗
(n

k )⎫⎪⎬
⎪⎩⎣
j =1

⎪⎭
i
= rx (k − i) + ∑ ai∗ ( j) rx (k − i + j) j =1
=
⎡ ⎢rx
(i

k
)
+
i
∑ ai
(
j)
rx
1 6
2 3
⎤ ⎥ ⎦

且:
b
(0)
=
rx
(
0)
+
a
(1)
rx
(1)
+
a
(
2)
rx
(
2)
=
1

1 6
×
1 2

2 3

《数字信号处理》课后答案

《数字信号处理》课后答案

数字信号处理课后答案 1.2 教材第一章习题解答1. 用单位脉冲序列()n δ及其加权和表示题1图所示的序列。

解:()(4)2(2)(1)2()(1)2(2)4(3) 0.5(4)2(6)x n n n n n n n n n n δδδδδδδδδ=+++-+++-+-+-+-+-2. 给定信号:25,41()6,040,n n x n n +-≤≤-⎧⎪=≤≤⎨⎪⎩其它(1)画出()x n 序列的波形,标上各序列的值; (2)试用延迟单位脉冲序列及其加权和表示()x n 序列; (3)令1()2(2)x n x n =-,试画出1()x n 波形; (4)令2()2(2)x n x n =+,试画出2()x n 波形; (5)令3()2(2)x n x n =-,试画出3()x n 波形。

解:(1)x(n)的波形如题2解图(一)所示。

(2)()3(4)(3)(2)3(1)6() 6(1)6(2)6(3)6(4)x n n n n n n n n n n δδδδδδδδδ=-+-+++++++-+-+-+-(3)1()x n 的波形是x(n)的波形右移2位,在乘以2,画出图形如题2解图(二)所示。

(4)2()x n 的波形是x(n)的波形左移2位,在乘以2,画出图形如题2解图(三)所示。

(5)画3()x n 时,先画x(-n)的波形,然后再右移2位,3()x n 波形如题2解图(四)所示。

3. 判断下面的序列是否是周期的,若是周期的,确定其周期。

(1)3()cos()78x n A n ππ=-,A 是常数;(2)1()8()j n x n e π-=。

解:(1)3214,73w w ππ==,这是有理数,因此是周期序列,周期是T=14; (2)12,168w wππ==,这是无理数,因此是非周期序列。

5. 设系统分别用下面的差分方程描述,()x n 与()y n 分别表示系统输入和输出,判断系统是否是线性非时变的。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

第四章附加题
1. 由三阶巴特沃思低通滤波器的幅度平方函数推到其系统函数,设
1/c rad s Ω=。

2. 设计一个满足下列指标的模拟Butterworth 低通滤波器,要求通带的截止频率
6,p f kHz =,通带最大衰减3,p A dB =,阻带截止频率12,s f kHz =,阻带的最
小衰减25s A dB =,求出滤波器的系统函数。

3. 设计一个模拟切比雪夫低通滤波器,要求通带的截止频率 f p =3kHz ,通带衰
减要不大于0.2dB ,阻带截止频率 f s = 12kHz ,阻带衰减不小于 50dB 。

4. 数字滤波器经常以下图描述的方式来处理限带模拟信号。

(1)
如果系统()h n 的截止频率是8rad s π,110T kHz =,等效模拟滤波器的截止频率是多少?
(2)
设120kHz =,重复(1)。

()
()
()
()
()
()
()
T T a x
t x n y n y t a
h n −−−→
−−−→
−−−→
−−−→模-数变换器
数-模变换器
采样周期采样周期
5. 一个线性时不变因果系统由下列差分方程描述
()()()()10.51y n x n x n y n =----
(1) 系统函数()H Z ,判断系统属于FIR 和IIR 中的哪一类以及它的滤波特性(低通、高通等)。

(2) 若输入()()2cos 0.55x n n π=+ ()0n ≥,求系统输出信号达到稳态后的最大幅度値。

6. 设()a h t 表示一模拟滤波器的单位冲激响应,
()0.9,0
0,0
t a e t h t t -⎧≥=⎨<⎩
用冲激响应不变法,将此模拟滤波器转化成数字滤波器(()h n 表示单位取样响应,即()()a h n Th nT =)。

确定系统函数()H z ,并把T 作为参数,T 为任何值时,数字滤波器是稳定的,并说明数字滤波器近似为低通滤波器还是高通滤波器。

7. 用冲激响应不变法将以下
变换为
,抽样周期为T 。

(1) 。

(2)
,n 为任意正整数。

8. 用冲激响应不变法设计一个离散时间低通滤波器,连续时间巴特沃斯滤波器
的幅度平方函数为
()()
2
21
1/a N
c H j j j Ω=
-ΩΩ
滤波器的技术指标为
()11
0j p p H e ωδωω-≤≤≤≤
()j s
s H e ωδωωπ≤≤≤
假设没有混叠,问冲激响应不变法中所用的采样周期值对设计结果是否有影响,并说明理由。

9. 图示是由RC 组成的模拟滤波器,写出其传输函数()a H s ,并选用一种合适
的转换方法,将()a H s 转化成数字滤波器()H z ,最后画出网络结构图。

10. 试用双线性变换法设计一低通数字滤波器,并满足技术指标如下: (1) 通带和阻带都是频率的单调下降函数,而且没有起伏; (2) 频率在0.5π处的衰减为 3.01dB -;
(3) 频率在0.75π处的衰减至少为15dB 。

11. 设计一个数字滤波器,要求3dB 截止频率为π2.0,频率在π5.0到π之间的
阻带衰减至少为15dB ,系统取样频率为500Hz ,用双线性变换法设计满足指标的最低阶巴特沃思滤波器的传递函数,并用正准型结构实现。

12. 设计一个数字低通滤波器,要求其截止频率 f c = 50Hz ,在 fc 处的衰减为α
c = 3dB 。

系统的取样频率为
500Hz ,要求从二阶 Butterworth 模拟低通滤波
器中用双线性变换法导出数字滤波器的 H(z)。

13. 冲激响应不变法和双线性变换法是两种滤波器设计方法。

这两种方法将s 左
半平面的极点映射到z 平面单位圆内部,从而保持模拟滤波器的稳定性。

如果一个模拟滤波器的所有极点和零点都在s 左半平面内,那么这个滤波器具有最小相位。

(1)确定冲激响应不变法是否可以将最小相位模拟滤波器映射为最小相位
离散时间系统。

(2)对于双线性变换法重复(1)。

14. 阶数2N ≥时,连续时间滤波器()a H s 的系统函数可以用两个低阶系统的级
联来表示
()()()12a a a H s H s H s =
所以,一个数字滤波器可以通过将变换直接应用到()a H s 来设计,也可以通过分别将()1a H s 、()2a H s 变换为()1H z 、()2H z 来设计,然后以级联方式实现
()()()12H z H z H z =
(1) 如果()1H z 、()2H z 用冲激响应不变法由()1a H s 、()2a H s 设计,与用冲激响应不变法由()a H s 直接设计的滤波器()H z ,式()()()12H z H z H z =是否成立。

(2) 对双线性变换法重复(1)。

15. 假设某模拟滤波器()a H s 是一个低通滤波器,又知()()1
1
a z s z H z H s +=
-=,数字
滤波器()H z 的通带中心位于下面哪种情况?并说明原因。

(1) 0ω=(低通)。

(2) ωπ=(高通)。

(3) 除0或π以外的某一频率(带通)。

16. 设计一个数字巴特沃斯高通滤波器,要求其通带截止频率0.8p ωπ=rad ,通
带衰减不大于3dB ,阻带截止频率0.5s ωπ=rad ,阻带衰减不小于18dB ,采样间隔为2s 。

17. 证明u z =-(旋转变换)是一个低通→高通的稳定变换。

18. 把模拟低通滤波器传递函数中的s 用1/s 代替,就得到模拟高通滤波器。


若)(s G a 是低通滤波器的传递函数,)(s H a 是高通滤波器的传递函数,则
)/1()(s G s H a a =。

另外,数字滤波器还可以借助双线性变换1
1
+-=
z z s 从模拟滤波器映射得到(为方便起见,设2=T )。

在这种映射下,虽然频率刻度有
了畸变,但保留了幅度特性的特征。

下图的网格表示一个截止频率为
2/
πω=L
的低通滤波器。

常数A 、B 、C 、D 都是实数。

试问为了得到截止频
率为2/πω=H 的高通滤波器,应如何修改这些系数?
19. 利用双线性变换法设计一个满足下面指标要求的数字带通巴特沃思滤波器:
通带上下边缘频率各为 200Hz 和300Hz ,通带波动3dB ,阻带上下边缘频率分别为 50Hz 和450Hz ,阻带衰减20dB ,取样频率为1kHz 。

请用模拟滤
ω
波器频率变换实现本设计。

20. 设滤波器差分方程为
)2(6
1
)1(65)1(21)()(---+-+=n y n y n x n x n y
试用正准型及一阶节的级联型、并联型结构实现此差分方程。

相关文档
最新文档