数字信号处理第三章作业评讲
数字信号处理 Chapter03答案
11
3.2 Properties of the z-Transform
Ex. ( linearity) x(n) = [3(2n) – 4(3n)] u(n) 3 4 – 1 – 2z –1 1 – 3z –1
X(z) =
ROC: |z| > 3
12
3.2 Properties of the z-Transform
z = re
jθ
=
n =−∞
∑ x ( n )r
−n
∞
− n − jθ n
e
X ( z) ≤
n =−∞
∞
∑
−1
x (n) r
+∑
n=0
∞
∞
x ( n) rn
x (n) rn
≤ ∑ x ( −n ) r + ∑
n n =1 n =0
7
3.1 The z-Transform
3.1.1 The Direct z-Transform
3.2 Properties of the z-Transform
X(z) = ∑ x(n) z – n
14
3.2 Properties of the z-Transform
X(z) = ∑ x(n) z – n
15
X(z) = ∑ x(n) z – n
16
3.3 Rational z-Transforms
1 2 −1 1 2 2
X ( z ) = 1+ z + (
X ( z) = 1 1− z
1 2 −1
)
z + .... + (
−2
1 n 2
)
z −n
《数字信号处理教程》(第三版)第三章
N 1 N 1
~ km ~ kn x1 (m ) WN X 2 (k ) WN n 0 m 0
1 N
N 1
N 1 m 0
~ ~ (m ) X (k ) W ( m n)k x1 2 N
n 0
N 1
~ ~ x1 ( m ) x2 ( n m )
域是连续周期的;若此时我们对频域的连续信号抽样,
人为的使其离散化,这样,频域的离散又导致时域的周 期化。于是有:
时域离散、周期
频域周期、离散
3.3 周期序列的离散傅立叶级数(DFS)
~ 设x (n)是周期为N的一个周期序列 ~ ~ x ( n) x (n rN ) ,r为任意整数
注:不论是离散的,还是连续的周期序列,均可用傅立叶级数 表示。离散的周期序列用离散傅立叶级数表示。(任一个周
第三章
离散傅立叶变换
理解傅里叶变换的几种形式
了解周期序列的傅里叶级数及性质,掌握周期卷 积过程
理解离散傅里叶变换及性质,掌握圆周移位、共 轭对称性,掌握圆周卷积、线性卷积及两者之间 的关系
了解频域抽样理论
理解频谱分析过程
了解序列的抽取与插值过程
3.2Leabharlann 傅立叶变换的几种可能形式傅立叶变换
~ x (n) x(n模N ) x(( n)) N
其中,(n模N)或((n))N数学上表示“n对N取余数或取模值”。
~ ~ 和 ~ 所对应的x(n)。 例: ( n)的周期为N=9,求 x ( 25) x ( 5) x
~ x (25) x(25模9) x(( 25))9 x(7) ~ x (5) x(5模9) x(( 5)) x(4)
数字信号处理[方勇]第三章习题答案解析
3-1 画出)5.01)(25.01()264.524.14)(379.02()(211211------+--+--=z zz z z z z H 级联型网络结构。
解:243-2 画出112112(23)(465)()(17)(18)z z z H z z z z --------+=--+级联型网络结构。
解:()x n ()y n 243-3 已知某三阶数字滤波器的系统函数为1211252333()111(1)(1)322z z H z z z z -----++=-++,试画出其并联型网络结构。
解:将系统函数()H z 表达为实系数一阶,二阶子系统之和,即:()H z 11122111111322z z z z ----+=+-++ 由上式可以画出并联型结构如题3-3图所示:)题3-3图3-4 已知一FIR 滤波器的系统函数为121()(10.70.5)(12)H z z z z ---=-++,画出该FIR滤波器的线性相位结构。
解: 因为121123()(10.70.5)(12)1 1.30.9H z z z z z z z ------=-++=+-+,所以由第二类线性相位结构画出该滤波器的线性相位结构,如题3-4图所示:()x n 1-1-1z -题3-4图3-5 已知一个FIR 系统的转移函数为:12345()1 1.25 2.75 2.75 1.23H z z z z z z -----=+--++求用级联形式实现的结构流图并用MATLAB 画出其零点分布及其频率响应曲线。
解: 由转移函数可知,6=N ,且)(n h 偶对称,故为线性相位系统,共有5个零点,为5阶系统,因而必存在一个一阶系统,即1±=z 为系统的零点。
而最高阶5-z 的系数为+1,所以1-=z 为其零点。
)(z H 中包含11-+z 项。
所以:11()()(1)H z H z z -=+。
1()H z 为一四阶子系统,设12341()1H z bz czbz z ----=++++,代入等式,两边相等求得12341()10.2530.25H z z zz z ----=+-++,得出系统全部零点,如图3-5(b )所示。
数字信号处理第三版(姚天任、江太辉) 答案 第三章
3.1 图 P3.1 所示的序列 x(n) 是周期为 4 的周期性序列。请确定其傅里叶级数的系数 X (k) 。
∑ ∑ ∑ 解: X (k)
=
N −1
x(n)WNnk
=
N −1
x(−n)WNnk
=
−( N −1)
x(n)WN−nk
=
X (−k)
解:图 P3.5_1 所示的是计算这两个序列的周期卷积 x3 (n) 的过程,可以看出,x3 (n) 是 x1 (n) 延时 1 的结果, 即 x3(n) = x1(n −1) 。
3.6 计算下列序列的N点DFT:
(1) x(n) = δ (n)
(2) x(n) = δ [(n − n0 )]N * RN (n), 0 < n0 < N
总计需要时间: (105 + 21)s = 126s
用 FFT 计算 DFT:
复数乘法:
N 2
log2
N
=
5120次, 5120 ×100μ s
≈
0.512s
复数加法: N log2 N = 10240次,10240× 20μs ≈ 0.2048s
总计需要时间: (0.512 + 0.2048)s = 0.7168s
(2) x2 (n) = x ⎡⎣(2 − n)⎤⎦4 R4 (n)
解: x1(n) 和 x2 (n) 的图形如图 P3.7_1 所示:
3.8 图 P3.8 表示一个 4 点序列 x(n) 。 (1)绘出 x(n) 与 x(n) 的线性卷积结果的图形。 (2)绘出 x(n) 与 x(n) 的 4 点循环卷积结果的图形。 (3)绘出 x(n) 与 x(n) 的 8 点循环卷积结果的图形,并将结果与(1)比较,说明线性卷积与循环卷
数字信号处理(方勇)第三章习题答案
数字信号处理(方勇)第三章习题答案3-1 画出)5.01)(25.01()264.524.14)(379.02()(211211------+--+--=z zz z z z z H 级联型网络结构。
解:23-2 画出112112(23)(465)()(17)(18)z z z H z z z z --------+=--+级联型网络结构。
解:()x n ()y n 243-3 已知某三阶数字滤波器的系统函数为1211252333()111(1)(1)322z z H z z z z -----++=-++,试画出其并联型网络结构。
解:将系统函数()H z 表达为实系数一阶,二阶子系统之和,即:()H z 11122111111322z z z z ----+=+-++ 由上式可以画出并联型结构如题3-3图所示:)题3-3图3-4 已知一FIR 滤波器的系统函数为121()(10.70.5)(12)H z z z z ---=-++,画出该FIR 滤波器的线性相位结构。
解: 因为121123()(10.70.5)(12)1 1.30.9H z zz z z z z ------=-++=+-+,所以由第二类线性相位结构画出该滤波器的线性相位结构,如题3-4图所示:()x n 1-1-1z -题3-4图3-5 已知一个FIR 系统的转移函数为:12345()1 1.25 2.75 2.75 1.23H z z z z z z -----=+--++求用级联形式实现的结构流图并用MATLAB 画出其零点分布及其频率响应曲线。
解: 由转移函数可知,6=N ,且)(n h 偶对称,故为线性相位系统,共有5个零点,为5阶系统,因而必存在一个一阶系统,即1±=z 为系统的零点。
而最高阶5-z 的系数为+1,所以1-=z 为其零点。
)(z H 中包含11-+z 项。
所以:11()()(1)H z H z z -=+。
数字信号处理课后答案+第3章(高西全丁美玉第三版)
所以
kn DFT[ X (n)] = ∑ X (n)W N n =0 N −1
N −1 mn kn = ∑ ∑ x(m)W N W N n =0 m =0
N −1
n = ∑ x ( m)∑ W N ( m + k ) m =0 n =0
X (k ) − W X (k ) = ∑ WNkm − ( N − 1)
k N m =1
kn = ∑ W N − 1 − ( N − 1) = − N n =0 N −1
N −1
所以,
X (k ) =
−N , k ≠ 0 ,即 k 1 − WN N ( N − 1) k =0 2 X (k ) = −N k = 1, 2, ⋯, N − 1 k 1 − WN
=
1− e
−j
2π (m−k ) N N 2π (m−k ) N
1− e
−j
N = 0
k =m k≠m
0≤k≤N-1
(6) X (k ) = ∑ cos
n =0
N −1
1 2π kn mn ⋅ WN = (e 2 N n =0
∑
N −1
j
2π mn N
+e
-j
2π 2π mn - j kn N )e N
j
2π mn N ,
0<m< N
2π x(n) = cos mn , 0 < m < N N
(7) (8) (9)
x(n)=ejω0nRN(n) x(n)=sin(ω0n)RN(n) x(n)=cos(ω0n)RN(N)
第三章数字信号处理课后答案刘顺兰版
第三章数字信号处理课后答案刘顺兰版第三章部分习题解答(数字信号处理(第二版),刘顺兰,版权归作者所有,未经许可,不得在互联网传播)3.1如果一台通用计算机的速度为平均每次复乘需100μs ,每次复加需20μs ,今用来计算N=1024点的)]([n x DFT ,问用直接运算需要多少时间,用FFT 运算需要多少时间?解:∑?=====101010,21024,)()(N n nk N M N Wn x k X直接运算所需的总时间为s N N s N T d μμ20)1(1002×?+×=秒分62126201023102410010242=≈××+×=s s s μμFFT 运算所需总时间为s NM s M N T F μμ201002×+×=s s s 717.02010102410010102421=××+×××=μμ3.2在基-2FFT 算法中,最后一级或开始一级运算的系数10==N p N W W ,即可以不做乘法运算。
问(1)乘法可节省多少次,所占百分比为多少?解:可节省2N 次,所占百分比为 %100log 1%100log 2222×=×N N N N 如 8=N 则为%3.33%10031≈×3.11以20kHz 的采样率对最高频率10kHz 的带限信号()a x t 采样,然后计算)(n x 的1000N =个采样点的DFT ,即210()()N j nk N n X k x n eπ??==∑,1000N =.(1)试求频谱采样点之间的频率间隔是多少?(2)在()X k 中,200k =对应的模拟频率是多少?(3)在()X k 中,700k =对应的模拟频率是多少?解:(1)频谱采样点之间的频率间隔为:20000201000s f f Hz N Δ=== (2)200k =对应的模拟频率为20000200400041000s k f f k Hz kHz N ==×== (3)因700k =大于N/2,故其对应的模拟频率为 20000()300600061000s k f f N k Hz kHz N =?=×== 3.12 对一个连续时间信号)(t x α采样1s 得到一个4096个采样点的序列:(1) 若采样后没有发生频谱混叠,)(t x α的最高频率是多少?(2) 若计算采样信号的4096点DFT,DFT 系数之间的频率间隔是多少Hz?(3) 假定我们仅仅对Hz f Hz 300200≤≤频率范围所对应的DFT 采样点感兴趣,若直接用DFT,要计算这些值需要多少次复乘?若用按时间抽取FFT 则需要多少次?解:(1)由题意可知:4096s f Hz =,故)(t x α的最高频率/22048h s f f Hz == (2)409614096s f f Hz N Δ=== (3)直接用DFT 计算,所需要的复乘次数为(3002001)1014096413696d M N =?+=×=若用按时间抽取FFT 则需要的复乘次数为10log 204812245762F N M N ==×= 3.17若给定两个实序列)(1n x 、)(2n x ,令:)()()(21n jx n x n g +=,)(kG 为其傅里叶变换,可以利用快速傅里叶变换来实现快速运算,试利用傅里叶变换的性质求出用)(k G 表示的)(1n x 、)(2n x 的离散傅里叶变换)(1k X 、)(2k X 。
数字信号处理课后第三章习题答案
1 e j 0 N
2 j(0 k ) N 1 e
k 0, 1, , N 1
(8) 解法一
直接计算:
1 j 0 n x8 (n) sin(0 n) RN (n) [e e j 0 n ] R N ( n ) 2j
X 8 (n)
n 0
N 1
kn x8 (n)WN
k 0, 1, , N 1
(4)
X (k ) WNkn
n 0
m1
π j ( m1) k 1 WNkm N e 1 WNk
π sin mk N R (k ) N π sin k N
第3章
离散傅里叶变换(DFT)及其快速算法 (FFT)
所以
DFT[ X (n)] X (n)W
n 0
N 1
N 1
kn N
N 1 mn kn x(m)WN WN n 0 m 0
N 1
n ( m k ) x(m)WN m 0 n 0
N 1
第3章
由于
离散傅里叶变换(DFT)及其快速算法 (FFT)
第3章
离散傅里叶变换(DFT)及其快速算法 (FFT)
(10) 解法一
X (k )
n 0
N 1
kn nWN
k 0, 1, , N 1
上式直接计算较难, 可根据循环移位性质来求解X(k)。 因为x(n)=nRN(n), 所
以
x(n)-x((n-1))NRN(n)+Nδ(n)=RN(n) 等式两边进行DFT, 得到
1 [e j0 n e j0 n ] e 2 j n 0
数字信号处理第三章习题作业答案
1 e 当 k 2, 4, 6,... 时,X 1 (k ) 0
序列3:
x3 (n) x1 (n) x1 (n 4)
根据序列移位性质可知
X 3 (k ) X1 ( k ) e j k X1 ( k ) (1 e j k )
即 x(n) 是以 n 0 对称轴的奇对称
故这三个序列都不满足这个条件
(3)由于是8点周期序列,其DFS:
nk X (k ) x(n )WN x (n )e n 0 n 0 N 1 7 j 2 nk 8
序列1:
X 1 (k ) e
n 0
3
y 解: 序列 x(n) 的点数为 N1 6 , (n) 的点数为 N 2 15, 故 x(n) y (n) 的点数应为
N N1 N 2 1 20
是线性卷积以15为周期周期延拓后取主值序列 19( N 1) 0
15 ( L)
又 f (n) 为 x(n) 与 y (n) 的15点的圆周卷积,即L=15。
第三章习题讲解
n 1, 0 n 4 h(n) R4 (n 2) 3.设 x(n) 其他n 0, h 令 x(n) x((n))6 , ( n) h((n)) 6 ,
试求 x(n) 与 h (n) 的周期卷积并作图。
解:
y ( n ) x ( m )h ( n m )
4 ( L N 1)
15 ( L)
34 ( L N 1)
混叠点数为N-L=20-15=5 n 0 ~ n 4( N L 1) 故 f (n)中只有 n 5到 n 14的点对应于 x(n) y (n)
数字信号处理 刘顺兰第三章完整版习题解答
数字信号处理刘顺兰第三章完整版习题解答一、题目解答1. 题目利用时域抽样、频域抽样、零填充、插值法等,实现信号的变换。
1.1 时域抽样时域抽样是指将一个连续时间信号在时间轴上的等间隔位置上进行采样,可以得到一个离散时间信号。
时域抽样的原理是,将时间轴上的信号按照一定的时间间隔进行采样,每个采样点的振幅值就是该点对应的连续时间信号的振幅值。
时域抽样可以通过以下步骤进行实现:1.假设连续时间信号为x(t),采样频率为Fs(采样频率是指每秒采样的次数),采样间隔为Ts(采样间隔是指相邻两个采样点之间的时间间隔)。
2.根据采样频率和采样间隔,计算出采样点数N:N =Fs * T,其中T为采样时长。
为Ts。
4.在每段的中点位置进行采样,得到N个采样点。
5.将N个采样点按照时域顺序排列,即可得到离散时间信号。
1.2 频域抽样频域抽样是指将一个连续频谱信号在频率轴上的等间隔位置上进行采样,可以得到一个离散频谱信号。
频域抽样的原理是,将频率轴上的信号按照一定的频率间隔进行采样,每个采样频率点上的能量值就是该频率点对应的连续频谱信号的能量值。
频域抽样可以通过以下步骤进行实现:1.假设连续频谱信号为X(f),采样频率为Fs(采样频率是指每秒采样的次数),采样间隔为Δf(采样间隔是指相邻两个采样频率点之间的频率间隔)。
2.根据采样频率和采样间隔,计算出采样点数N:N =Fs / Δf,其中Δf为采样频率点之间的频率间隔。
为Δf。
4.在每段的中点位置进行采样,得到N个采样频率点。
5.将N个采样频率点按照频域顺序排列,即可得到离散频谱信号。
1.3 零填充零填充是指在信号的末尾添加一些零值样本,使得信号的长度变长。
零填充的原理是,通过增加信号的长度,可以在时域和频域上提高信号的分辨率,从而更精确地观察信号的特征。
零填充可以通过以下步骤进行实现:1.假设原始信号为x(n),长度为N。
2.计算需要填充的长度L,L > 0。
数字信号处理答案第三章
= = =
0 0 1 j 2πn e 10 , n = 1, 2, . . . , k. 2
3.3
(a) X1 (z ) = = = = The ROC is (b)
1 3 ∞ 0
1 1 ( )n z −n − 1 ( )n z −n + 3 2 n=−∞ n=0 1
1 −1 1− 3 z
+ +
1 ( )n z n − 1 2 n=0 1 − 1, 1− 1 2z −1 2 z)
∞
1
1−1 −1 3zFra bibliotek(1 −
5 6 1 −1 )(1 3z
< |z | < 2. X2 (z ) = = = 1 ( )n z −n − 2n z −n 3 n=0 n=0 1 1−
1 −1 3z ∞ ∞
nan cosw0 nz −n nan ejw0 n + e−jw0 n −n z 2 60
© 2007 Pearson Education, Inc., Upper Saddle River, NJ. All rights reserved. This material is protected under all copyright laws as they currently exist. No portion of this material may be reproduced, in any form or by any means, without permission in writing from the publisher. For the exclusive use of adopters of the book Digital Signal Processing, Fourth Edition, by John G. Proakis and Dimitris G. Manolakis. ISBN 0-13-187374-1.
数字信号处理 答案 第三章
解: x1 ( n) 和 x2 (n) 的图形如图 P3.7_1 所示:
3.8 图 P3.8 表示一个 4 点序列 x( n) 。 (1)绘出 x( n) 与 x( n) 的线性卷积结果的图形。 (2)绘出 x( n) 与 x( n) 的 4 点循环卷积结果的图形。 (3)绘出 x( n) 与 x( n) 的 8 点循环卷积结果的图形,并将结果与(1)比较,说明线性卷积与循环卷 积之间的关系。
j [(2π k /10) + (π /10)]
={
3.7
N ,k=m或 k=−m 2 0,其 他
图 P3.7 表示的是一个有限长序列 x( n) ,画出 x1 ( n) 和 x2 (n) 的图形。 (1) x1 ( n) = x ⎡ ⎣( n − 2 ) ⎤ ⎦ 4 R4 (n)
(2) x2 ( n) = x ⎡ ⎣( 2 − n ) ⎤ ⎦ 4 R4 (n)
解: (1) X ( k )
= ∑ δ (n)WNnk = δ (0) = 1, 0 ≤ k ≤ N − 1
n=0
N −1
(2) X ( k ) =
∑ δ [(n − n )]
n =0 0
N −1
N
RN (n)WNnk = WNn0 k , 0 ≤ k ≤ N − 1
(3) (4)
X (k ) = ∑ a W
− jω N
−j
N ω 2
j
N ω 2
−j
N ω 2
⎛N ⎞ sin ⎜ ω ⎟ N −1 ) ⎝ 2 ⎠ e− j 2 ω = sin
ω
2
⎛N ⎞ sin ⎜ ω ⎟ ⎝ 2 ⎠ , ϕ (ω ) = − N − 1 ω | X (e jω ) |= ω 2 sin 2
数字信号处理(方勇)第三章习题答案
3-1 画出)5.01)(25.01()264.524.14)(379.02()(211211------+--+--=z zz z z z z H 级联型网络结构。
解:243-2 画出112112(23)(465)()(17)(18)z z z H z z z z --------+=--+级联型网络结构。
解:()x n ()y n 243-3 已知某三阶数字滤波器的系统函数为1211252333()111(1)(1)322z z H z z z z -----++=-++,试画出其并联型网络结构。
解:将系统函数()H z 表达为实系数一阶,二阶子系统之和,即:()H z 11122111111322z z z z ----+=+-++ 由上式可以画出并联型结构如题3-3图所示:)题3-3图3-4 已知一FIR 滤波器的系统函数为121()(10.70.5)(12)H z zz z ---=-++,画出该FIR滤波器的线性相位结构。
解: 因为121123()(10.70.5)(12)1 1.30.9H z zz z z z z ------=-++=+-+,所以由第二类线性相位结构画出该滤波器的线性相位结构,如题3-4图所示:()x n 1-1-1z -题3-4图3-5 已知一个FIR 系统的转移函数为:12345()1 1.25 2.75 2.75 1.23H z z z z z z -----=+--++求用级联形式实现的结构流图并用MATLAB 画出其零点分布及其频率响应曲线。
解: 由转移函数可知,6=N ,且)(n h 偶对称,故为线性相位系统,共有5个零点,为5阶系统,因而必存在一个一阶系统,即1±=z 为系统的零点。
而最高阶5-z 的系数为+1,所以1-=z 为其零点。
)(z H 中包含11-+z 项。
所以:11()()(1)H z H z z -=+。
1()H z 为一四阶子系统,设12341()1H z bzcz bz z ----=++++,代入等式,两边相等求得12341()10.2530.25H z zz z z ----=+-++,得出系统全部零点,如图3-5(b )所示。
数字信号处理课后答案第3和4章
用DFT/FFT对信号进行谱分析的误差表现在三个方面, 即混叠现象、 栅栏效应和截断效应。 截断效应包括泄漏和 谱间干扰。
第3章 离散傅里叶变换(DFT)及其快速算法
第3章 离散傅里叶变换(DFT)及其快速算法
(FFT)
xN(n)=IDFT[X(k)]为x(n)的周期延拓序列(以N为延拓周期) 的主值序列。 以后这一结论可以直接引用。
[例3.4.2] 已知 x(n)=R8(n), X(ejω)=FT[x(n)]
对X(ejω)采样得到X(k),
X(k)X(ej)|2πk, k0,1, ,5 6
第3章 离散傅里叶变换(DFT)及其快速算法
(FFT)
当然, 截取信号的长度要足够长。 但如果截取的长度 不够长, 而依靠在所截取的序列尾部加零点, 增加变换区 间长度, 也不会提高分辨率。 例如, 分析周期序列的频谱, 只观察了一个周期的1/4长度, 用这些数据进行DFT, 再通 过尾部增加零点, 加大DFT的变换区间N, 也不能分辨出是 周期序列, 更不能得到周期序列的精确频率。
令m=N-1-n, 则上式可写成
0
N1
X(k) x(m )W N k(n1) x(m )W N km
m N1
m 0
W N k(N 1 )X ( (k)N )R N (k)
第3章 离散傅里叶变换(DFT)及其快速算法
(FFT)
当 k N 时(N为偶数), 2
因为
X N 2 W N N 2(N 1 )X N 2 NW N N 2(N 1 )X N 2
数字信号处理教程习题分析与解答(程佩青)第三章
第三章 离散傅立叶变换1.如下图,序列x(n)是周期为6的周期性序列,试求其傅立叶级数的系数。
∑∑=-===56265)(~)(~)(X ~:n nkj nkn e n x W n x k π解kj k j k j kj kj e e e e e 562462362262621068101214πππππ-----+++++=计算求得:。
339)5(~; 33)4(~ ; 0)3(~; 33)2(~;339)1(~;60)0(~j X j X X j X j X X +=-==+=-==。
并作图表示试求设)(~),(~)(~ .))(()(~),()(.264k X n x k X n x n x n R n x == ∑∑=-===56265)(~)(~)(~:n nk jnk n en x W n x k X π解k j kj k j e e e πππ---+++=3231。
计算求得: 3)5(~; 1)4(~ ; 0)3(~ ;1)2(~; 3)1(~ ; 4)0(~j X X X X j X X ====-==。
的周期卷积并作图与试求令其它,设 )(~)(~,))(()(~,))(()(~,)2()(,040,1)(.3464n h n x n h n h n x n x n R n h nn n n x ==-=⎩⎨⎧≤≤+=解:在一个周期内的计算值)(~)(~*)(~)(~m n h n h n x n y -==)(~)(~*)(~)(~m n h n h n x n y -==等各序列。
试画出所示如图已知)())((),())3((,))(()())((),())((,))((,13)(.47755633665n R n x n R n x n x n R n x n R n x n x P n x ----)()()5()(x(n)(4)N n 0 ),n -(n )()3()()()2()()(cos )()1()(52000n R n n x n nR n x n R a n x n R n a n x DFT N N N N n N ==<<===δω闭合形式表达式点试求以下有限长序列的])21sin()2sin()21sin()2sin([21])()()()([21)(]1111[)(][)(])([)()(cos )()()(cos )(:0)2(21020)2(2102)2(21)2(21)2(21222)2(21)2(21)2(21222)()(211)(10)(2110211000000000000000000002002002022002ϖπϖϖπϖωωϖπϖϖπϖϖπϖπϖπϖϖϖϖπϖπϖπϖϖϖωωωωωωωωππππππ-++⋅=--+--=--+--=+=+===---+---------+-++-----+---=---=+--=---=-∑∑∑∑k N e N e k N e N e a e ee e e e eeeee e a k R ee ee a k R eea k R e e e a k R en a k X n R n a n x k N j N j k Nj Njk Nj k Nj k Nj NjNjN jk Nj k N j k Nj NjNjNjN k j N j k j N j N N n nj N n nk j N N n nkj n j n j N n N nkj N N N N N N N 解)(111121)(21)()(21)()(cos )( )()(cos )( ) 1 (:)2()2(10)2(10)2(1020010200000k R e ee e a k R e e a k R e e e a k Re n a k X n R n a n x N k N j N j k Nj Nj N N n nN j N n nk N j N N n nk N j n j n j N n N nk N j N ⎥⎥⎦⎤⎢⎢⎣⎡--+--=⎥⎦⎤⎢⎣⎡+=⎥⎦⎤⎢⎣⎡+===--+---=---=+--=---=-∑∑∑∑ωπωωπωωπωππωωπωω解⎥⎥⎥⎦⎤---⎢⎢⎢⎣⎡--=------+-++---)()()()(21)2(21)2(21)2(21222)2(21)2(21)2(212220000000ωπωπωπωωωωπωπωπωωωk Nj k N j k N j N jN jNjk N j k N j k N j Nj N j N j e e e eeee e e e e e a⎥⎥⎥⎥⎦⎤--⎢⎢⎢⎢⎣⎡+⋅=--+--)21sin()2sin()21sin()2sin(210)2(21020)2(21020000ωπωωπωωπωωπωk N e Nek N e Ne a k Nj Njk N j N jk Nj N N n nk NjnN n aea eak X n R a n x ππ210211)()()((2)--=---===∑)( )()( )()()( 0,)()( (3)02102010200k Re k Re n n k Re n x k X N n n n n x Nk n N j NN n nk N j N n NnkN j πππδδ--=--=-=-==<<-=∑∑)(1)( 11)1()())1(()(])1)2( 2[)1( 32()1)(()()()()( )()(411)1(32)1(321)1(110)1(1k R W Nk X N W W N k R W N k R N W N W W W N W W W nW nWW k X k R nW k X W k R nW k X n nR n x N kNkNkN N N n nk N N k N N k N k N k N N kN k N k N N n kn N N n nk Nk NN n N k n N k NN n N nkN N --=∴-=--+--=+--=-+-+++--++++=-=-==∴=∑∑∑∑∑-=---=+-=-=+-=∙∙∙∙∙∙)(kNN N n nkNN W Nk X n nR n x W n k X n R n n x --===∴=∑-=1)()()()4()( )()(5111022,则小题的结论根据第)(221111122)1(232)1(23210)1(2121)1(2)1()2()(12)2()(2)2(2)2()12()1(]1()2(4[)1(94)1)(()(k N kN kNN n nk NN n nkNk N N kN k N k N N k N k N k N N n kn NN n nk NkNN n kn N k NW N W N N k X W NN N k X N N nWN N W n N N W N W W W N W W W W n W nW k X W n k X W ---=∴----=+--=+--=-+--=-+-+++--++++=-=-=∑∑∑∑∑-=-=---=+-=-=+∙∙∙∙∙∙)∙∙∙±±±===∑-=,6,4,2,0)(~)3(?])0([)()2(?)()1(:;)(~1)(~).(~.61)/2(k k x X k X k X ek X Nn x n x N k nkN j 哪些序列列能做到成虚数外除时间原点使所有的哪些序列能够通过选择成为实数时间原点使所有的哪些序列能够通过选择问傅里叶级数这些序列可以表示成列如图画出了几个周期序π条件。
《数字信号处理》朱金秀第三章习题及参考答案
第三章习题答案 3.1 (1)非周期(2)N=1 (3)N=10 (4)N=4 (5)N=20 3.2 02s f f ωπ=,1s sf T = (1)0153,2f ωπ==;0.3s T =,05f π= (2)010,25f ωπ==;0.3s T =,0503f =(3)0,0.55f πω==;0.3s T =,013f =(4)03.5,8.75f ωπ==;0.3s T =,0356f =(5) ()()()(){}0.20.210.20.20.20.2(0.2)(0.2)1cos(0.2)()2130.6cos(0.2)() 1.8()0.6()211.80.6()0.6()2110.910.610.6j n j n n n j n j n n nj n j n j j n e e F n u n F e e u n F e u n F e u n ee ππππππωπωπππ-+-----+=+⎡⎤⎡⎤-=-•+-⎢⎥⎣⎦⎣⎦⎡⎤⎡⎤=-•-+-⎣⎦⎣⎦⎛⎫=-+ ⎪++⎝⎭3.3 function [X]=myDTFT(x, n, w)% 计算DTFT% [X]=myDTFT(x, n, w) %X=输出的DTFT 数组 %x=输入的有限长序列 %n=样本位置行向量 %w=频率点位置行向量 X=x*exp(-j*n ’*w)3.4 (1) 7()10.3j j X e eωω-=- (2)20.51()(10.5)10.5j j j j e X e e e ωωωω---=---(3)2()0.80.1610.4j j j e X e e ωωω--=⨯⨯-(4)112210.920.9()(10.9)10.9(10.9)j j j j j j e e X e e e e ωωωωωω-----⨯-⨯=-=---3.5(1) 23456()642246j j j j j j j X e e e e e e e ωωωωωωω------=++++++(2)234567()642246j j j j j j j j X e e e e e e e e ωωωωωωωω-------=+++++++ (3)234567()642246j j j j j j j j X e e e e e e e e ωωωωωωωω-------=+++---- (4)235678()642246j j j j j j j j X e e e e e e e e ωωωωωωωω-------=+++----3.6 00()()11()211j j j A X e ae ae ωωωωω---+⎡⎤=+⎢⎥--⎣⎦3.7 N=5,()5611()11j j j j j j e ee X e e e ωωωωωω----=+--N=25,()252611()11j j j j j j e e eX e e e ωωωωωω----=+-- N=100,()10010111()11j j j j j j e ee X e e e ωωωωωω----=+-- N=5,》n = -5:5; x =ones(1,11); % x(n)k = -500:499; w = (pi/500)*k; % [-pi, pi] X =1/11* x*exp(-j*pi/500*n'*k); % DTFT magX = abs(X); angX = angle(X); realX = real(X); imagX = imag(X); subplot(2,2,1); plot(w/pi,magX); gridxlabel('以pi 为单位的频率'); title('幅度部分'); ylabel('幅值') subplot(2,2,2); plot(w/pi,angX); gridxlabel('以pi 为单位的频率'); title('相位部分'); ylabel('弧度')-1-0.500.5100.51以pi 为单位的频率幅度部分幅值-1-0.500.51-4-2024以pi 为单位的频率相位部分弧度N=25,>> n = -25:25; x =ones(1,51); % x(n)k = -500:499; w = (pi/500)*k; % [-pi, pi] X =1/51* x*exp(-j*pi/500*n'*k); % DTFT magX = abs(X); angX = angle(X); realX = real(X); imagX = imag(X); subplot(2,2,1); plot(w/pi,magX); gridxlabel('以pi 为单位的频率'); title('幅度部分'); ylabel('幅值') subplot(2,2,2); plot(w/pi,angX); gridxlabel('以pi 为单位的频率'); title('相位部分'); ylabel('弧度')-1-0.8-0.6-0.4-0.200.20.40.60.81以pi 为单位的频率相位部分弧度-1-0.8-0.6-0.4-0.200.20.40.60.81以pi 为单位的频率幅度部分幅值N=100,>> n = -100:100; x =ones(1,201); % x(n)k = -500:499; w = (pi/500)*k; % [-pi, pi] X =1/201* x*exp(-j*pi/500*n'*k); % DTFT magX = abs(X); angX = angle(X); realX = real(X); imagX = imag(X); subplot(2,2,1); plot(w/pi,magX); gridxlabel('以pi 为单位的频率'); title('幅度部分'); ylabel('幅值') subplot(2,2,2); plot(w/pi,angX); gridxlabel('以pi 为单位的频率'); title('相位部分'); ylabel('弧度')-1-0.500.5100.51以pi 为单位的频率幅度部分幅值-1-0.500.51-4-2024以pi 为单位的频率相位部分弧度随着N 的增大,DTFT 的幅度特性主瓣越尖锐,旁瓣越小,越接近于1)(=n x 的DTFT 特性。
数字信号处理课后答案+第3章DFT+FFT
ej(02N πk)(N21)sin( 02N πk)N 2 sin(02N πk)/2
k0,1, ,N1
或
1ej0N X7(k)1ej(02 N k)
(8) 解法一 直接计算:
k0,1, ,N1
x 8 (n ) si0 n n )R (N (n ) 2 1 j[ej 0 n e j 0 n ]R N (n )
即 X 8 (k ) jX 7 o (k ) j1 2 [X 7 (k ) X 7 * (N k )]
结果与解法一所得结果相同。 此题验证了共轭对称性。
(9) 解法一 直接计算:
x9(n )co0 n s )R N ((n )1 2[ej 0 n e j 0 n]
N1
X9(k) x9(n)WNkn
x 7 ( n ) e j 0 n R N ( n ) [c 0 n ) o j ss i 0 n ) ( n R N ] ( n ( )
所以
x 8 (n ) si0 n )R N ( (n ) Im x 7 (n )[ ]
所以 D [ jx 8 ( F n ) D ] T [ j IF x m 7 ( n ) T ] X [ 7 o ] ( k )
(7) x(n)=ejω0nRN(n) (8) x(n)=sin(ω0n)RN(n) (9) x(n)=cos(ω0n)RN(N) (10) x(n)=nRN(n) 解: (1)
X(k)N n011WN knN n01ej2N πkn1 1 ee jj2 2N N π πkkN N
N k0 0 k1,2,,N1
(5)
X (k)
N 1 j 2π mn
eN
W Nkn
N 1 j 2π (mk )n
数字信号处理第三版西安科大出版高西全丁玉美课后答案第3与4章
x 6 ( n ) ID [X ( k F ),]n T 0 ,1 ,2 , ,5
第3章 离散傅里叶变换(DFT)及其快速算法
(FFT)
解:直接根据频域采样概念得到
x6(n ) x(n 6 l)R 6(n )R 6(n )R 2(n ) l
[例3.4.3] 令X(k)表示x(n)的N点DFT, 分别证明: (1) 如果x(n)满足关系式
yc(1)
x(1)
x(0)
x(L1)
x(2)
h(1)
yc(2)
x(2)
x(1)
x(0) x(3) h(2)
yc(L1) x(L1) x(L2) x(L3) x(0)h(L1)
第3章 离散傅里叶变换(DFT)及其快速算法
(FFT)
循环卷积定理: 若
yc(n)=h(n) L x(n) 则
~xN(n) x(niN) n
会发生时域混叠, xN(n)≠x(n)。
通过频率域采样得到频域离散序列xN(k), 再对xN(k)进行 IDFT得到的序列xN(n)应是原序列x(n)以采样点数N为周期进行 周期化后的主值区序列, 这一概念非常重要。
第3章 离散傅里叶变换(DFT)及其快速算法
(FFT)
(FFT)
3.1.2 重要公式
1) 定义
N1
X(k)DF [x(T n)N ] x(n)W N k n k=0, 1, …, N-1 n0
x(n)ID[X F(kT )N ]N 1N k 0 1X(k)W N kn
2) 隐含周期性
k=0, 1, …, N-1
N 1
N 1
X (k m ) N x (n ) W N (k m )n N x (n ) W N k nX (k )
《数字信号处理》第4版第三章作业评讲
X (3) X * (5) 1 x(0) (10 8 6 4) 2 8
(2) X (0)
x(n)W80n x(n) 6
n 0 n 0
7
7
3.19 解:
解 : (1)
1 而F 10Hz F 最小纪录长度为 0.1s
TP
TP
3.21 解: (1) Y (k ) W83k X (k ) 根据 DFT 的圆周移位性质
y(n) x3 (n) ,即将 x(n) 循环右移 3 位
x(n) 3,0,5,0, 4,0,0,0 y (n) 0,0,0,3,0,5,0, 4
(2)根据 DFT 的对称性质
在一个周期内的计算值35nwnw小题的结论根据第1000101010khzkhz一个纪录中的最少点数的整数幂必须为又因频率为允许处理的信号的最高最小纪录长度为根据dft的圆周移位性质2根据dft的对称性质
3.2 解:
~ 解 : X (k )
n 0
x ( n )W6nk ~
j k
5
1 e 3 e jk ~ ~ ~ 计算求得:X (0) 4 ; X (1) j 3 ; X (2) 1 ; ~ ~ ~ X (3) 0 ; X (4) 1 ; X (5) j 3 。
N
(k )
N
N 1 n 0
( n n0 ) e
j 2 n0k N R N
j 2 nk N R
(k )
e
(k )
解: (4)x( n) nR N ( n) X (k )
k WN X (k )
nW
n0 N 1 n0
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
2.解:
∑∑=-===
5062650)(~)(~)(~ :n nk j nk n e n x W n x k X π解
k j k j k j e e e πππ---+++=3231。
计算求得: 3)5(~ ; 1)4(~ ; 0)3(~ ; 1)2(~ ; 3)1(~ ; 4)0(~j X X X X j X X ====-== 3.解:
解:在一个周期内的计算值
5.解:
(1)
)
(~)(~*)(~)(~m n h n h n x n y -==)
(~)(~*)(~)(~m n h n h n x n y -=
=
k N j N N n nk N j n N n ae a e a k X n R a n x π1
0211)()
()((2)--=---===∑
)
( )()( )()()( 0,)()( (3)021********k R e k R e n n k R e n x k X N
n n n n x N k n N j N N n nk j N n N nk N j πππδδ--=--=-=-=
=
<<-=∑∑ )(1)( 11)1()())1(()
(])1)2( 2[)1( 32()1)(()
()()()( )
()(411)1(32)1(3210)1(101
0)1(10k R W N k X N W W N k R W N k R N W N W W W N W W W nW nW W k X k R nW k X W k R nW k X n nR n x N k N k N k N N N n nk N N k N N k N k N k N N k N k N k N N n k n N
N n nk N k N N n N k n N k N N n N nk N N --=∴-=--+
--=+
--=-+-+
++--+
+++=-=-==
∴=∑∑∑
∑∑-=---=+-=-=+-=∙∙∙∙∙∙)(解:
N N N n nk N N W N k X n nR n x W n k X n R n n x --=
==
∴=∑-=1)()()()4()( )()(5111022,则小题的结论
根据第)(解: 22
1111122)1(232)1(2321
0)1(210210)1(2
)1()2()(12)2()
(2)2(2)2()12()1(]
1()2(4[)1(94)1)(()(k N k N k N N n nk N
N n nk N
k N N k N
k N k N N k N k N k N N n k n N
N n nk N k N N n k n N
k N W N W N N k X W N N N k X N N nW
N N W
n N N W N W W W N W W W W n W n W k X W n k X W ---=∴--
--=+--=+--=-+--=-+-+++--+
+++=-=-=∑∑∑∑∑-=-=---=+-=-=+∙∙∙∙∙∙)
7.解:
(){18,21,24,9,12,15}
05y n n =≤≤
10.解:
直接利用圆周卷积公式求解
(){0,4,2,10,10,8,4}
06
f n n =-----≤≤
11.解:
)())(()()(650621n R m n x m x n y m ⎥⎦
⎤⎢⎣⎡-=∑=
()[][]相等。
与倍时,的整数为不一定为零),而当个其他的数值
的每两个值之间插入相当于在的周期为倍的的抽样点数是在一个周期内解)()(()1()(),)(()()(,)1,1,0()()( )()()()( 10)()( :1021010102r
k X k Y l r k r k X r N k Y r k X k Y N l lr k r k X e n x W n x W n y n y DFT k Y N k e n x n x DFT k X N n k n πj nk rN rN n N n nk rN N n nk j -∴-======
=-≤≤=
=∙∙∙∑∑∑∑-=--=-=-=-π
14.解:
10242: 2,1000101
.01.0 )3( 5 52
1 2 10101
.011 )2( 1.0 10
1 101 )1( :1033
==∴=⨯=≥∴=<∴>=⨯==∴≥∴≤=∙∙∙∙
∙∙N N T T N KHz
KHz f f f f KHz T f s
s T Hz F F T P s h h s s P P 为一个纪录中的最少点数的整数幂必须为又因频率为允许处理的信号的最高最小纪录长度为而解
22.()(2)(21)()() 0,1...1y n x n jx n g n jh n n N =++=+=-
*1()[()()]21()[()*()]2G k Y k Y N k H k Y k Y N k j
=+-=-- 22()()()
()()()0,1,2 (1)
k N k N X k G k W H k X k N G k W H k k N =++=-=-。