数字信号处理习题与答案_数字信号处理习题答案
数字信号处理习题与答案由刀豆文库小编整理,希望给你工作、学习、生活带来方便,猜你可能喜欢“数字信号处理习题答案”。
3.已知
单位抽样响应为,通过直接计算卷积和的办法,试确定的线性移不变系统的阶跃响应。
9.列出下图系统的差分方程,并按初始条件
求输入为
时的输出序列,并画图表示。
解:系统的等效信号流图为:
解:根据奈奎斯特定理可知:
6.有一信号,它与另两个信号
和的关系是:
其中,已知,解:根据题目所给条件可得:
而
所以
8.若是因果稳定序列,求证:
证明:
∴
9.求的傅里叶变换。
解:根据傅里叶变换的概念可得:
13.研究一个输入为
和输出为的时域线性离散移不变系
统,已知它满足
并已知系统是稳定的。试求其单位抽样响应。解:
对给定的差分方程两边作Z变换,得:,为了使它是稳定的,收敛区域必须包括
即可求得
16.下图是一个因果稳定系统的结构,试列出系统差分方程,求系统函数。当
时,求系统单位冲激响应 , 画出系统零极点图和频率响应曲线。
由方框图可看出:差分方程应该是一阶的
则有
因为此系统是一个因果稳定系统;所以其收敛
17.设是一离散时间信号,其z变换为
求它们的z变换:,对下列信
号利用(a),这里△记作一次差分算子,定义为:
(b)(c)解:(a){
(b),(c)
由此可设
1.序列x(n)是周期为6的周期性序列,试求其傅立叶级数的系数。
~解: X(k)n05~x(n)W6nkn05j2nk~x(n)e6 j2k1412e6j22k10e6j23k8e6j24k6e6j25k10e6
计算求得:
~2.设x(n)R4(n),x(n)x((n))6.~~ 试求X(k)并作图表示~x(n),X(k)。~~~X(0)60;X(1)9j33;X(2)3j3;~~~X(3)0;X(4)3j3;X(5)9j33。
~解: X(k)n0x(n)W6nk~5n0j~x(n)e52nk6
~~~计算求得:X(0)4;X(1)j3;X(2)1;~~~ X(3)0;X(4)1;X(5)j3。jk1e3j2ke3ejk
n1,0n43.设x(n),h(n)R4(n2),0,其它n~令~x(n)x((n))6,h(n)h((n))4,~试求~x(n)与h(n)的周期卷积并作图。解:在一个周期内的计算
~~~y(n)~x(n)*h(n)h(nm)~~~y(n)~x(n)*h(n)h(nm)7x(n), 0n5设有两序列 x(n)0, 其他ny(n), 0n14 y(n)0, 其他n各作15点的DFT,然后将两个DFT相乘,再求乘积的IDFT,设所得结果为f(n),问f(n)的哪些点对应于x(n)y(n)应该得到的点。
解:序列x(n)的点数为N16,y(n)的点数为N215故又x(n)*y(n)的点数应为:NN1N2120f(n)为x(n)与y(n)的15点的圆周卷积,即L15所以,混叠点数为NL20155。用线性卷积结果 以15 为周期而延拓形成圆周卷积序列 f(n)时,一个周期 内在n0到n4(NL1)这5点处发生混叠,即f(n)中只有n5到n14的点对应于x(n)*y(n)应该得到的点。
8.已知x(n)是N点有限长序列,X(k)DFT[x(n)]。现将长度变成rN点的有限长序列y(n)x(n), 0nN-1y(n)0, NnrN-1试求DFT[y(n)](rN点DFT)与X(k)的关系。解: X(k)DFTxn Y(k)DFTy(n)
x(n)n0rN1N1j2nkeNN1n00kN1n0nky(n)WrNx(n)WnkrNn0N1j2πnkx(n)eNrkX()rklr(l0,1,N1)在一个周期内,Y(k)的抽样点数是X(k)的r倍(Y(k)的周期为Nr),相当于在X(k)的每两个值之间插入(r1)个其他的数值k(不一定为零),而当k为r的整数l倍时,Y(k)与X()相等。r 9已知x(n)是长为N点的有限长序列,X(k)DFT[x(n)]现将x(n)的每两点之间补进r1个零值点,得到一个长为rN点的有限长度x(n/r), nir, 0iN序列y(n), y(n)0, 其他n试求rN点DFT[y(n)]与X(k)的关系。解: X(k)DFTxn Y(k)DFTy(n)
N1n0n0nkx(n)WN,0kN1rN1nky(n)WrNN1i0x(ir/i0N1irkr)WrNx(i)WikN,0krN1Y(k)X((k))NRrN(k)Y(k)是将X(k)(周期为N)延拓r次形成的,即Y(k)周期为rN。
10.频谱分析的模拟信号以8kHz被抽样,计算了512个抽样的DFT,试确定频谱抽样之间的频率间隔,并证明你的回答。
证明 : s2fssF00fsF002其中s是以角频率为变量 的 频谱的周期,0是频谱抽样之间的频谱间隔。fssNF00F0对于本题:fsNfs8KHzN512 8000F015.625Hz51211.设有一谱分析用的信号处理器,抽样点数必须为2的整数幂,假定没有采用任何殊数据处理措施,要求频率分辨力10Hz,如果采用的抽样时间间隔为0.1ms,试确定(1)最小记录长度;(2)所允许处理的信号的最高频率;(3)在一个记录中的最少点数。11解:(1)TP而F10Hz TPsF10 最小纪录长度为 0.1s 1110310KHzT0.11 fs2fh fhfs5KHz2 允许处理的信号的最高频率为5KHz(2)fs TP0.11031000,又因N必须为2的整数幂T0.1 一个纪录中的最少点数为:N2101024(3)N
用直接I型及典范型结构实现以下系统函数
34.2z10.8z2H(z)20.6z10.4z2
121.52.1z10.4z21.52.1z0.4zH(z)12121(0.3z0.2z)10.3z0.2z解:H(z)
∵1anznn1m0NbznMmY(z)X(z)
∴a10.3,a20.24(z1)(z21.4z1)H(z)(z0.5)(z20.9z0.8)
2.用级联型结构实现以下系统函数b01.5,b12.1,b20.4
试问一共能构成几种级联型网络。11kz12kz2H(z)A121zzk1k2k解:
4(1z1)(11.4z1z2)112(10.5z)(10.9z0.8z)
∴ A4
111, 110.5 , 210 , 121.4 ,210 , 120.9 ,221 220.8
由此可得:采用二阶节实现,还考虑分子分母组合成二阶(一阶)基本节的方式,则有四种实现形式。
3.给出以下系统函数的并联型实现。
5.21.58z11.41z21.6z3H(z)112(10.5z)(10.9z0.8z)
解:对此系统函数进行因式分解并展成部分分式得:
5.21.58z11.41z21.6z3H(z)112(10.5z)(10.9z0.8z)
0.210.3z14110.5z10.9z10.8z2 G0 110.5 , 210,120.9 ,220.8
010.2 , 110,021 , 120.3
4.用横截型结构实现以下系统函数:
11H(z)1z116z112z11z11z126
解:
11H(z)(1z1)(16z1)(12z1)(1z1)(1z1)26
111122(1z12z1z)(1z6zz)(1z)26
1537(1z1z2)(1z1z26
2)(z11)8205220581z1zz3z4z531212 5.已知FIR滤波器的单位冲击响应为
0.3n(h(n)(n)N1n01)0.n72(2)n0.11(3n)0
试画出其级联型结构实现。
H(z)根据h(n)zn得:
220.z70.z3114
1H(z)10.z3z0.12)1z23
(10.z20.)(1z10.1z2 0.4而FIR级联型结构的模型公式为:
H(z)(0k1kz12kz2)k1N2
对照上式可得此题的参数为:
011 , 021, 110.2 , 120.1210.3 , 220.4
6.用频率抽样结构实现以下系统函数:
52z33z6H(z)1z1
抽样点数N = 6,修正半径r0.9。解;
因为N=6,所以根据公式可得:
H(z)2166(1rz)H0(z)H3(z)Hk(z)6k1(53z3)(1z3)H(z)1z1 (53z3)(1z1z2)故 H(k)H(Z)Z2k/N (53ejk)(1e因而 H(0)24,H(1)223j,H(2)0 H(3)2,H(4)0,H(5)223j
j3kej2k3)则 H0(z)H(0)241rz110.9z1H(3)2 H3(z)1rz110.9z1
0111z121求 : Hk(z)k1 时 :H1(z)2212zrcosrzN
012ReH(1)2Re[223j]411(2)(0.9)ReH(1)W613.643.6z1H1(z)10.9z10.81z2k2 时 :02120,H2(z)0 7.设某FIR数字滤波器的系统函数为:
1H(z)(13z15z23z3z4)5
试画出此滤波器的线性相位结构。解:由题中所给条件可知:
1331h(n)(n)(n1)(n2)(n3)(n4)5555
则 h(0)h(4)10.253 h(1)h(3)0.65 h(2)1N12 2即h(n)偶对称,对称中心在 n处,N 为奇数(N5)。8.设滤波器差分方程为:
y(n)x(n)x(n1)11y(n1)y(n2)34
⑴试用直接I型、典范型及一阶节的级联型、一阶节的并联型结构实现此差分方程。
⑵求系统的频率响应(幅度及相位)。
⑶设抽样频率为10kHz,输入正弦波幅度为5,频率为1kHz,试求稳态输出。解:
(1)直接Ⅰ型及直接Ⅱ:
根据 y(n)ak1Nky(nk)bx(nk)可得:kk0M
11a1 , a234;
b01 , b11
一阶节级联型:
1z1H(z)111z1z2341z1 11011101(1z)(1z)66
1z111
(10.7z)(10.36z)
一阶节并联型:
H(z)1z1(111011101z)(1z)66
17171010220220110111011z1z66
1.60.610.7z110.36z1
1z1(2)由题意可知 H(z)111z1z234 1ejH(e)1j12j1ee34 j(1cos)jsin11111cosco2sjsinsin23443
幅度为:
H(ej)
(1cos)2sin21111(1coscos2)2(sinsin2)23434
相位为:
sinargH(ej)arg)tg(1cos
11sinsin24tg(3arg)111cosco2s34
(3)输入正弦波为 : x(t)5sin(2t103)
3由 T210T12 可得:
又抽样频率为10kHz,即抽样周期为
13T0.1100.1ms31010
∴在x(t)的一个周期内,采样点数为10个,且在下一周期内的采样值与(0,2)间的采样值完全一样。所以我们可以将输入看为 周期为:T11103s1ms1000
5sin10x(n)5sin2103nT32104n1 5sinn(n0 ,1 ,5
由此看出,9)
00.2
根据公式可得此稳态输出为:
y(n)5H(ej0)cos0nargH(ej0)12.13cos0.2n51.6
4.试用N为组合数时的FFT算法求N12的结果(采并画出流图。1.如果一台通用计算机的速度为平均每次复乘需50 s 计算需要多少时间,用FFT运算需要多少时间。
每次复加5 s,用它来计算512点的DFT[x(n)],问直拉对于0nN,有解:依题意:N34r1r2,解: ⑴ 直接计算:
复乘所需时间: T61510N2 51065122 1.31072s
复加所需时间: T20.5106N(N1)0.5106512(5121)0.130816s TT1T21.441536s⑵用FFT计算:
复乘所需时间: T61510N2log2N 51065122log2512 0.01152s
复加所需时间: T20.5106Nlog2N 0.5106512log2512 0.002304s TT1T20.013824s
nn1r2n0,n10,1,2n00,1,2,3 同样: 令Nr2r1 对于频率变量k(0kN)有kkk10,1,2,31r1k0,k00,1,2x(n)x(n1r2n0)x(4n1n0)x(n1,n0)X(k)X(k1r1k0)X(3k1k0)X(k1,k0)11X(k)x(n)Wnk12n032 x(n(4n1n0)(3k1k01,n)0)W12n00n10