课后习题及答案第4章快速傅里叶变换习题答案.pdf
傅里叶变换习题及答案
傅里叶变换习题及答案傅里叶变换习题及答案傅里叶变换是一种重要的数学工具，广泛应用于信号处理、图像处理、物理学等领域。
它能够将一个函数表示为一系列正弦和余弦函数的和，从而将时域上的函数转换为频域上的函数。
为了帮助读者更好地理解和掌握傅里叶变换，本文将介绍一些傅里叶变换的习题，并提供相应的答案。
1. 问题：计算函数 f(t) = 2cos(3t) + 4sin(5t) 的傅里叶变换。
解答：根据傅里叶变换的定义，我们可以将 f(t) 表示为一系列正弦和余弦函数的和。
首先，我们需要计算 f(t) 的频谱。
根据欧拉公式，我们可以将 cos(3t) 和sin(5t) 表示为指数形式。
cos(3t) = (e^(3it) + e^(-3it)) / 2sin(5t) = (e^(5it) - e^(-5it)) / 2i将上述结果代入 f(t) 的表达式中，得到：f(t) = 2((e^(3it) + e^(-3it)) / 2) + 4((e^(5it) - e^(-5it)) / 2i)= e^(3it) + e^(-3it) + 2i(e^(5it) - e^(-5it))接下来，我们需要计算 f(t) 的傅里叶变换F(ω)。
根据傅里叶变换的定义，可以得到：F(ω) = ∫[(-∞,∞)] f(t)e^(-iωt) dt将 f(t) 的表达式代入上述公式中，并进行积分计算，得到：F(ω) = ∫[(-∞,∞)] (e^(3it) + e^(-3it) + 2i(e^(5it) - e^(-5it)))e^(-iωt) dt通过对每一项进行积分计算，最终得到F(ω) 的表达式为：F(ω) = π(δ(ω - 3) + δ(ω + 3)) + 2πi(δ(ω - 5) - δ(ω + 5))其中，δ(x) 表示Dirac δ 函数。
2. 问题：计算函数 f(t) = e^(-2t)u(t) 的傅里叶变换。
解答：首先，我们需要将 f(t) 表示为指数形式。
数字信号处理答案(第三版)清华大学
数字信号处理教程课后习题答案目录第一章离散时间信号与系统第二章Z变换第三章离散傅立叶变换第四章快速傅立叶变换第五章数字滤波器的基本结构第六章无限长单位冲激响应（IIR）数字滤波器的设计方法第七章有限长单位冲激响应（FIR）数字滤波器的设计方法第八章数字信号处理中有限字长效应第一章 离散时间信号与系统1 .直接计算下面两个序列的卷积和)n (h *)n (x )n (y =请用公式表示。
分析：①注意卷积和公式中求和式中是哑变量m （ n 看作参量）， 结果)(n y 中变量是 n ,; )()()()()(∑∑∞-∞=∞-∞=-=-=m m m n x m h m n h m x n y ②分为四步 （1）翻褶（ -m ），（2）移位（ n ）,（3）相乘，; )( )( 4n y n n y n 值的，如此可求得所有值的）相加，求得一个（③ 围的不同的不同时间段上求和范一定要注意某些题中在 n00 , 01()0 , ,()0,n n n a n N h n n n n x n n n β-⎧≤≤-=⎨⎩⎧≤⎪=⎨<⎪⎩其他如此题所示，因而要分段求解。
)(5.0)(,)1(2 )()4()(5.0)(,)2( )()3()()(,)( )()2()()(,)( )()1(3435n u n h n u n x n R n h n n x n R n h n R n x n R n h n n x n n n =--==-=====δδ2 .已知线性移不变系统的输入为)n (x ,系统的单位抽样响应 为)n (h ,试求系统的输出)n (y ,并画图。
分析：①如果是因果序列)(n y 可表示成)(n y ={)0(y ，)1(y ，)2(y ……}，例如小题（2）为)(n y ={1，2，3，3，2，1} ;②)()(*)( , )()(*)(m n x n x m n n x n x n -=-=δδ ;③卷积和求解时，n 的分段处理。
数字信号处理课后答案第3和4章fb
[e 2
n0
N 1
j 0 n
e
]e
j N j N 1 e 0 1 e 0 2π 2π j( 0 k) j( 0 k) 2 N N 1 e 1 e
第３章
离散傅里叶变换(DFT)及其快速算法 (FFT)
解法二
因为
由DFT共轭对称性可得同样结果。
离散傅里叶变换(DFT)及其快速算法 (FFT)
(6)
2π kn X (k ) cos mn W N N n0
N 1
N 1
1 2
j
2π N
mn
(e
e
-j
2π N
mn
-j
2π N
kn
)e
n0
1 2
N 1
j
2π N
(mk )n
e
1 2
n0
N 1
第３章
离散傅里叶变换(DFT)及其快速算法 (FFT)
教材第3章习题与上机题解答
1． 计算以下序列的N点DFT， 在变换区间0≤n≤N－1内， 序列定义为 (1) x(n)=1 (2) x(n)=δ(n) (3) x(n)=δ(n－n0) 0<n0<N (4) x(n)=Rm(n) 0<m<N
1 2j
j(
2π N
mn )
[e
e
j(
2π N
mn )
]
2π sin mn N
n=0, 1, …, N－1
第３章
离散傅里叶变换(DFT)及其快速算法 (FFT)
3． 已知长度为N=10的两个有限长序列：
课后习题及答案第4章快速傅里叶变换习题答案.pdf
和共轭反对称分量， 即
F(k)=X(k)+jY(k)=Fep(k)+Fop(k) 计算一次 N 点 IFFT 得到
f(n)=IFFT［F(k)］=Re［f(n)］+j Im［f(n)］ 由 DFT 的共轭对称性可知
Re［f(n)］=IDFT［Fep(k)］=IDFT［X(k)］=x(n) j Im［f(n)］=IDFT［Fop(k)］=IDFT［jY(k)］=jy(n)
X (k + N ) = X1(k) −W2kN X 2 (k)
k = 0,1,L, N −1
由上式可解出
X1(k)
=
1 2
[
X
(k)
+
X
(k
+
N )]
X
2
(k)
=
1 2
[X
(k)
+
X
(k
+
N
)]W2−Nk
k = 0,1, 2,L, N −1
由以上分析可得出运算过程如下：
（1）由 X(k)计算出 X1(k)和 X2(k):
Xk=conj(Xk)；
%对 Xk 取复共轭
xn=conj(fft(Xk， N))/N； %按照所给算法公式计算 IFFT
分别对单位脉冲序列、 长度为 8 的矩形序列和三角序列进行 FFT， 并调
用函数 ifft46 计算 IFFT 变换， 验证函数 ifft46 的程序 ex406.m 如下：
%程序 ex406.m
Tc = 2TF +1024 次复数乘计算时间 = 2 × 0.1536×10−3 +10×10−9 ×1024
= 0.317 44 ms 可实时处理的信号最高频率 fmax 为
《数字信号处理》课后答案
数字信号处理课后答案 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 分别表示系统输入和输出，判断系统是否是线性非时变的。
第四章 快速傅里叶变换
运 算
X[3] 2WN0 3WN3 3WN6 2WN9 1 j
复数加法 N(N-1)
复数乘法 N 2
效 率
?
1.将长序列DFT分解为短序列的DFT 由于DFT的运算量与N2成正比，N越小，
DFT的运算量越少。所以如果把一个长序列 的DFT能分解为短序列的DFT的组合，则显 然可取得减少运算工作量的效果。
N 1
X (k) DFT[x(n)]
x(n)W
kn N
n0
N 1
N 1
x(n)WNkn＋ x(n)WNkn
n0 n为偶数
n0 n 为奇数
N / 21
N / 21
＝ x(2r)WN2rk x(2r 1)WN(2r1)k
r 0
r 0
N / 21
N / 21
＝ x1 (r)WN2rk WNk x2 (r)WN2rk
X 3 (k) WNk / 2 X 4 (k)
k 0,1, , N 1 4
将系数统一为以Ｎ为周期，即
，可得
Wk N/2
WN2k
X X
1 1
(k (k
) X 3 (k) WN2k X 4 (k) N / 4) X 3 (k) WN2k
X
4
(k
)
k 0,1, , N 1 4
同样，对也可进行类似的分解。一直分解
第4章 快速傅里叶变换
4.1 快速傅里叶变换
4点序列{2，3，3，2} DFT的计算复杂度 如
N 1
X [k] x[k]WNkn
何 提
n0
高
X[0] 2WN0 3WN0 3WN0 2WN0 10
DFT
X[1] 2WN0 3WN1 3WN2 2WN3 1 j
(完整版)通信原理第七版课后答案解析樊昌信
第一章习题习题1.1 在英文字母中E 出现的概率最大，等于0.105，试求其信息量。
解：E 的信息量：()()b 25.3105.0log E log E 1log 222E =-=-==P P I习题1.2 某信息源由A ，B ，C ，D 四个符号组成，设每个符号独立出现，其出现的概率分别为1/4，1/4，3/16，5/16。
试求该信息源中每个符号的信息量。
解：b A P A P I A 241log )(log )(1log 222=-=-==b I B 415.2163log 2=-= b I C 415.2163log 2=-= b I D 678.1165log 2=-= 习题1.3 某信息源由A ，B ，C ，D 四个符号组成，这些符号分别用二进制码组00，01，10，11表示。
若每个二进制码元用宽度为5ms 的脉冲传输，试分别求出在下列条件下的平均信息速率。
（1） 这四个符号等概率出现； （2）这四个符号出现概率如习题1.2所示。
解：（1）一个字母对应两个二进制脉冲，属于四进制符号，故一个字母的持续时间为2×5ms。
传送字母的符号速率为Bd 100105213B =⨯⨯=-R 等概时的平均信息速率为b 2004log log 2B 2B b ===R M R R（2）平均信息量为符号比特977.1516log 165316log 1634log 414log 412222=+++=H 则平均信息速率为 b 7.197977.1100B b =⨯==H R R习题1.4 试问上题中的码元速率是多少？解：311200 Bd 5*10B B R T -===错误!未找到引用源。
习题1.5 设一个信息源由64个不同的符号组成，其中16个符号的出现概率均为1/32，其余48个符号出现的概率为1/96，若此信息源每秒发出1000个独立的符号，试求该信息源的平均信息速率。
解：该信息源的熵为96log 961*4832log 321*16)(log )()(log )()(22264121+=-=-=∑∑==i i i i M i i x P x P x P x P X H=5.79比特/符号因此，该信息源的平均信息速率 1000*5.795790 b/s b R mH ===错误!未找到引用源。
数字信号处理第4章习题解答教材
DFT [x2 (n)]
DFT {Im[ w(n )]}
1 j Wop (k )
1 2j
[W
((k )) N
W
* (( N
k )) N
]RN
(k)
解：由题意 X k DFT xn，Y k DFT y n 构造序列 Z k X k jY k 对Z k 作一次N点IFFT可得序列z n z(n) IDFT Z k
Re[w(n)] j Im[w(n)]
Wep (k) Wop (k)
由x1(n) Re[w(n)]得
X1(k) DFT[x1(n)] DFT{Re[w(n)]} Wep (k)
1 2
[W
((k
))
N
W *((N
k ))N
]RN
(k)
由x2 (n) Im[w(n)]得
X 2 (k )
(2) 按频率抽取的基-2FFT流图
同样共有L = 4级蝶形运算，每级N / 2 = 8个蝶形运算
基本蝶形是DIT 蝶形的转置
X m1(k )
X m1( j)
WNr
-1
X m (k ) Xm( j)
每个蝶形的两节点距离为2Lm ，即从第一级到 第四级两节点距离分别为8，4，2，1。
系数WNr的确定：r (k )2 2m1 即k的二进制左移m 1位补零
3. N=16 时，画出基 -2 按时间抽取法及按频率抽取法 的 FFT 流图（时间抽取采用输入倒位序，输出自然数 顺序，频率抽取采用输入自然顺序，输出倒位序）。
解： 自然序
倒位序
0 0000 0000 0 1 0001 1000 8 2 0010 0100 4 3 0011 1100 12 4 0100 0010 2 5 0101 1010 10 6 0110 0110 6 7 0111 1110 14
《数字通信基础》第4章答案
4-1 已知线性调制信号表示式如下： (1) cos t cos c t (2) (1 0.5sin t ) cos c t 式中， c 6 。试分别划出它们的波形图和频谱图。 [解] (1) 波形图：
频谱图：
S()
/2
c
c
c +
c
4-7 设某信道具有均匀的双边噪声功率谱密度 Pn ( f ) =0.5103W/Hz， 在该信道中传输抑制
载波的双边带信号，并设调制信号 m(t ) 的频带限制在 5kHz，而载波为 100kHz，已调信号 的功率为 10kW。若接收机的输入信号在加至解调器之前，先经过带宽为 10kHz 的一理想带 通滤波器滤波，试问： (1) 该理想带通滤波器的中心频率为多大？ (2) 解调器输入端的信噪功率比为多少？ (3) 解调器输出端的信噪功率比为多少？ (4) 求出解调器输出端的噪声功率谱密度，并用图形表示出来。 [解] (1) 该理想带通滤波器的中心频率为 100kHz。 (2)
Si 10 103 (W)， N i n0 B 2 0.5 103 10 103 10 (W)。所以，
Si 10000 1000 。 Ni 10
(3) 因为抑制载波的双边带调制的信噪比增益 G 2 ，所以
So S G i 2 1000 2000 。 No Ni
接收机输入信号功率
1 2 1 m (t ) Pm (t )df 2 2 fm n f 1 2 m df 0 2 2 fm 1 nm f m 4 1 （2）输出信号为 mout (t ) m(t ) ，所以输出信号功率 2 1 2 1 PSout m (t ) nm f m 4 8 Psin
傅立叶光学习题解答及参考答案
第一章 傅里叶分析部份习题解答及参考答案[1-1] 试分别写出图X1－1中所示图形的函数表达式。
图X1-1 习题[1-1]各函数图形解：(a)−∧L x x a 0 (b) () ∧−−L x b a L x a 2rect(c) ()x L x a sgn 2rect (d) x L x cos 2rect[1-2] 试证明下列各式。
(1) += 21comb 21comb x x- (2) ()()x i e x x x πcomb comb 2comb +=(3)()()()x x N x N ππsin sin lim comb ∞→= (4) ()()xx x πωδωsin lim ∞→=(5)()()∫∞∞−=ωωπδd cos 21x x (6)()ωπδωd 21∫∞∞−±=x i e x解：(1)原式左端∑∑∞−∞=∞−∞=+−−=−−=m n m x n x 12121δδ 令()1−=m n=−+=∑∞−∞=m m x 21δ右端 (2)()∑∑∞−∞=∞−∞=−=−= n n n x n x x 2222comb δδ n 2只取偶数()()∑∞−∞=−=m m x x δcomb()()πδδππm m x e m x e x m im m x i cos 2comb ∑∑∞−∞=∞−∞=−=−=当=m 奇数时，()()0comb comb =+xi ex x π；当=m 偶数时，令n m 2=，则12 cos =x π，并且有： ()()()∑∞−∞=−=+n n x x x 22e comb comb xi δπ 得证。
(3)由公式(1-8-7)知：()∑∞−∞=−=n nxex π2i comb上式可视为等比级数求和，其前N 项之和为：()()()()()x Nx e e e e e e e e q q a S x i x i x i Nx i Nx i Nx i x i Nx i N N ππππππππππsin sin 1111221=−−=−−=−−=−−−−−− 所以 ()()()x Nx S x N N N ππsin sin limlim comb ∞→∞→==得证。
程佩青《数字信号处理教程》(第4版)(课后习题详解 快速傅里叶变换(FFT))
4.2 课后习题详解4-1 如果一台通用计算机的速度为平均每次复乘40ns ，每次复加5ns ，用它来计算512点的DFT[x （n ）]，问直接计算需要多少时问?用FFT 运算需要多少时间?若做128点快速卷积运算，问最低抽样频率应是多少?解：①直接利用DFT 计算：复乘次数为N 2，复加次数为N （N-1）。
②利用FFT计算：复乘次数为，复加次数为N㏒2N 。
（1）直接计算复乘所需时间复加所需时间所以（2）用FFT 计算复乘所需时间复加所需时间所以4-2 N ＝16时，画出基-2按频率抽选法的FFT 流图采用输入自然顺序，输出倒位序），统计所需乘法次数（乘±1，乘±j 都不计在内）。
根据任一种流图确定序列x （n ）＝4cos （n π／2）（0≤n ≤15）的DFT 。
解：按频率抽取法的FFT 流图中的复数乘法出现在减法之后，其运算量为复数乘法：；复数加法：；由于N ＝16，有，，，不需要乘法。
按频率抽取，见图4-1（a ）。
图4-1（a ）运算量：复数乘法：由于，，，不需要乘法。
由图P4.2（a ）可知，共有的个数为1＋2＋4＋8＝15有的个数为1＋2＋4＝7所以总的乘法次数为32-15-7＝10（个）复数加法：举例：对序列x （n ）＝4cos （n π／2）（0≤n ≤15）可表示为由于N ＝16，可采用P4.2（b ）的流图。
设Xi （k ）＝（i ＝1，2，3，4）分别为第i 级蝶形结构的输出序列，则由P4.2（b ）的流图可知由于采用的是顺序输入、逆序输出的结构，因此输出X （k ）与X 4（k ）为逆序关系，即，为k 二进制逆序值由此可知，x （n ）的DFT 为X （4）＝X 4（2）＝32，X （12）＝X 4（3）＝12图4-1（b ）4-3 用MATLAB 或C 语言编制以下几个子程序。
（1）蝶形结运算子程序；（2）求二进制倒位序子程序；（3）基-2 DIT FFT 流程图，即迭代次数计算子程序。
2.快速傅里叶变换
快速Fourier变换()()()()10210,1,2,,Four 1,ie 1r 2jk i N N j ex k X j k N i N -===-=-∑π离散逆变换：()()()()()()()()2100,1,2,,Fourier 0,1,2,,1,111N n nj i N ex x x x N X j x n j N i -⎛⎫⎪⎝⎭=- -==-=-∑π离散变换给号序列::定信FFT-Fast Fourier Transform Fourie FFT FFT ,r FFT.逆变换的形式与完全相同结论平行地用到逆（）：快速变换优化()()()()()()()()()()()1111000001010000112222()()1nji N mn j mn j nj n n mj j N N n j n j n j n j n j m N N N N N e W W a W W W b W W W π+-+===⇓⇓=-010101012,,0,1,,10,1,,1,0,10,12,0,1,,1N m j n j n N j m j mj j j n n n n n m ==-=-⎧=+⎨=⎩=⎧=+⎨=-⎩设将和分别改写成22221(3)1,i m N m i m i N m N N N N N W e W W e W W eW ---====-===πππ记则计算复杂度降低！()()()()()()()()()()()100010101101001001100000111000,2,1,0,1,0,1;0,1,6271m n j n j N m n n j n m n j m n z n j W x n n W X j j z n j n j j m x n n W 对原信号偶数和奇数序列信号傅里叶变换-=-==⎧=+⎪⎪⎪⎪⎨⎪⎪=-⎪⎪===-⎩+∑∑∑()()()()()()()()()()()()()()0010000110010101121101011001110000100,2141Fourie 5r ,22nj N i N n m n j n j n j m N n n m n j n j N n n m j n j m X j X j j x n ex n n W W j j n n n W x n n W --=-=-======+=+=+⎡⎤+⎢⎦-⎥=⎣-∑∑∑∑∑π根据离散变换公式：进行一次分解快速Fourier 变换()()()()()()()()()()()()()()0010110101101010010011110011000010012,2221,1,,0,1,0,1;0,16,1m n j n j N m n m m n j n j m m n n n j n z n j W x n n W x n W x n W X j j z n j n j j m X j X j -=--====++===-===-∑∑∑∑快速Fourier 变换()()()()()()()()()001101001001001100006,2*,17,m n n n j n j N m j n z n j x n n X j j z n W W j -==⎧=+⎪⎪⎨⎪=-⎪⎩∑∑()(){}(){}()()()()()()()()()()()001000002200.601:,:111111:1,:17.66222n j N n j n j N N n j m W Wa n m m m m m m n m m m m m mb W N W N m N 因子和都是事先算好存储在计算机内的复杂度分奇偶项分析：时，次复数乘法，次复数加法乘法加法时，次复数乘法，次复数加法乘法加法复，，实际处理数据与计算复杂度杂度分析：2复数加法分析一次优化：总共需要做次复数乘法和次复数加法==+-⨯⨯-=+-⨯+⨯-≠0100,1,0,10,1,1n j j m ===-快速Fourier变换进行一次分解3N 为例分析FFT计算过程：232N =为例分析FFT计算过程：41+22+141242+24+1824⨯⨯⨯⨯⨯⨯FFT 计算量:=次复数乘法,=次复数加法2iN N W e 离散傅里叶逆变换:8个信号序列的离散傅里叶逆变换其中π=2iNN W e 离散傅里叶逆变换:8个信号序列的离散傅里叶逆变换其中π=()22=2,2FFT 1log 22log k N k kN N N kN N N==复执行一个以为底的完整数乘法：复数加法：的计算复杂度次优化：2log 0,N N N由于→→∞210FFT ,.,=2=1024()FFT 200!N N 与离散傅立叶变换需要次运算的在数量级上有了重大改进节省的工作量相当惊人比如对对实际问题来讲,这仅是一个很小的数字原算法的复数乘法次数就超过的倍()()()FFT 322FFT 21FFT ,)kk m N N ===计算直到为止的过程称为以为底的变换信号个数信号个数任何一个大于的自然数都可以作为变换的底在同一计算过程中还可以混合使用以3为底的换多个底数变()()()()()()()()()00101010100100110000010,2*,1,0,1,0,1;0,9,811m n j n j N m n n j n z n j W x n n W X j j z n j n j j m -==⎧=+⎪⎪⎪⎪=-⎨⎪⎪===-⎪⎪⎩∑∑James W. Cooley John W. Tukey 美国普林斯顿大学教授IBM 研究员James W. Cooley,John W. Tukey.An Algorithm for the Machine Calculation of Complex Fourier Series. Math. Comp.19(1965), 297-301.。
4快速傅里叶变换
X(N 2
k)
X1(k) WNk
X 2 (k )
(后一半)
例如 N=8 时的DFT,可以分解为两个N/2=4点的DFT.具 体方法如下:
(1)n为偶数时,即: x(0), x(2), x(4), x(6);
分别记作:
x1(0) x(0),
x1(1) x(2),
x1(2) x(4),
x1(3) x(6);
X (k) X1(k) WNk X 2 (k)
前一半
(k
0,1,,
N 2
1)
X (k) X1(k) WNk X 2 (k) 后一半
(k
N 2
,,
N
1)
实现上式运算的流图称作蝶形运算（N/2个蝶形)
X1(k) 1
X 2(k)
WNk
(前一半)
1
X (k) X1(k) WNk X 2 (k)
1
1
-1
由于
W W r
(
k
N 2
)
rk
N
N
（周期性）,所以:
2
2
N
X1( 2
k)
x (r)W N
2
1
1
r
(
N 2
k
)
N
2
r 0
N 2
1
x1(r)WNrk
2
r 0
X1(k)
同理,
X2(
N 2
k)
X 2 (k )
这就是说,X1(k),X2(k)的后一半,分别等于其前一
半的值。
又由于
W(
N 2
k
)
N
N
W2 N
N/4点 DFT
第4章课后习题答案及讲解
4-2 根据图P4—1所示的调制信号波形，试画出DSB及AM信号的波形图，并比较它们分别通过包络检波器后的波形差别。
解：DSB信号及包络检波后输出AM信号及包络栓波后输出由此可见，对DSB信号采用包络检波法不能正确还原基带信号。
4-3已知调制信号m(t)=cos(2000πt)+cos(4000πt)载波为cos104πt，进行单边带调制，试确定该单边带信号的表示式，并画出频谱图。
解：DSB信号为:S DSB(t)= [cos(2000πt)+ cos(4000πt)] cos104πt= 1/2[cos(12000πt)+cos(8000πt)] +1/2[cos(14000πt)+cos(6000πt)]SSB 信号为：上边带S SSB (t)= 1/2·cos(12000πt)+ 1/2·cos(14000πt)-8000π 0 6000π ω 下边带S SSB (t)= 1/2·cos(6000πt)+ 1/2·cos(8000πt)-14000π 012000π ω4-6 某调制系统如图P4-4所示。
为了在输出端同时分别得到f 1(t)及f 2(t)，试确定接收端的c 1(t)和c 2(t)。
解：该调制系统采用相干解调，设c1(t)=cos(ω1t+φ1)则接收端相乘器输出r1(t)=[f1(t) cosω0t + f2(t) sinω0t] cos(ω1t+φ1)= f1(t) cosω0t cos(ω1t+φ1) + f2(t) sinω0t cos(ω1t+φ1)=1/2 f1(t) [ cos(ω0t+ω1t+φ1)+ cos(ω0t- ω1t- φ1)]+1/2 f2(t) [ sin(ω0t+ω1t+φ1)+ sin(ω0t- ω1t- φ1)]若要经过低通滤波器后得到f1(t)，应有ω1=ω0，φ1=0，即c1(t)= cosω0t同理可得c2(t)= sinω0t思考题：4-11 什么是频分复用？答：频分复用(Frequency Division Multiplexing) 是按频率分割多路信号的方法，即将信道的可用频带分成若干互不交叠的频段，每路信号占据其中的一个频段。
2021年信号与线性系统分析(第四版)习题答案
1-1画出下列各信号的波形【式中)()(t t t r ε=】为斜升函数。
（2）∞<<−∞=−t et f t,)( （3）)()sin()(t t t f επ=（4）)(sin )(t t f ε= （5）)(sin )(t r t f = （7）)(2)(k t f kε= （10）)(])1(1[)(k k f kε−+=解：各信号波形为 （2）∞<<−∞=−t et f t,)(（3）)()sin()(t t t f επ=（4）)fεt=(sin)(t （5）)tf=r(sin)(t（7）)tf kε(k=(2)（10）)f kεk−=(k+]()1(1[)1-2 画出下列各信号的波形[式中)()(t t t r ε=为斜升函数]。
（1）)2()1(3)1(2)(−+−−+=t t t t f εεε （2）)2()1(2)()(−+−−=t r t r t r t f （5）)2()2()(t t r t f −=ε （8）)]5()([)(−−=k k k k f εε （11）)]7()()[6sin()(−−=k k k k f εεπ（12）)]()3([2)(k k k f k−−−=εε 解：各信号波形为（1）)2()1(3)1(2)(−+−−+=t t t t f εεε（2）)2()1(2)()(−+−−=t rt rt rtf（5）)2()2()(ttrtf−=ε（8）)]5()([)(−−=k k k k f εε（11）)]7()()[6sin()(−−=k k k k f εεπ（12）)]()3([2)(kkkf k−−−=εε1-3 写出图1-3所示各波形的表达式。
1-4 写出图1-4所示各序列的闭合形式表达式。
1-5 判别下列各序列是否为周期性的。
如果是，确定其周期。
（2）)63cos()443cos()(2ππππ+++=k k k f （5）)sin(2cos 3)(5t t t f π+=解：1-6 已知信号)(t f 的波形如图1-5所示，画出下列各函数的波形。
吴大正《信号与线性系统分析》(第4版)章节题库(傅里叶变换和系统的频域分析)【圣才出品】
第 4 章 傅里叶变换和系统的频域分析
一、选择题 1．图 4-1 所示系统由两个 LTI 子系统组成，已知子系统 H1 和 H2 的群时延分别为 τ1 和 τ2，则整个系统的群时延 τ 为（ ）。
图 4-1 A．τ1＋τ2 B．τ1－τ2 C．τ1·τ2 D．max（τ1,τ2） 【答案】A
9．如图 4-2 所示信号 f1（t）的傅里叶发换 F1（jω）已知，求信号 f2（t）的傅里叶发 换为（ ）。
图 4-2
【答案】A
【解析】由题意知， f2 (t) f1(t t0 ) 。由于 f2（t）＝f1（－（t＋t0）），根据傅里叶 发换的反转性质和时秱性质可知， F2 ( j) F1( j)e jt0 。
4．设 f（t）的频谱函数为 F（jω），则
的频谱函数等于（ ）。
【答案】D
2 / 150
圣才电子书 十万种考研考证电子书、题库视频学习平台
【解析】
可写为 f[－1/2（t－6）],根据傅里叶发换的尺度发换性质，
x（at）
|
1 a
|
[x(w
/
a）]，得
f[－1/2（t）]
A．x（t）＝－4Sa[2π（t－3）]
B．x（t）＝4Sa[2π（t＋3）]
C．x（t）＝－2Sa[2π（t－3）]
D．x（t）＝2Sa[2π（t＋3）]
【答案】A
【解析】常用的傅里叶发换对
Sa(ct)
c
G2c
()
令c 2 ，则有 4Sa(2t) 2G4 ()
ቤተ መጻሕፍቲ ባይዱ
再由傅里叶发换的时秱性质，有
4Sa[2 (t 3)] 2G4 ()e j3
数字信号处理教程 程佩青 课后题答案
第一章 离散时间信号与系统2.任意序列x(n)与δ(n)线性卷积都等于序列本身x(n)，与δ(n-n 0)卷积x(n- n 0),所以（1）结果为h(n) (3)结果h(n-2) （2（4）3 .已知 10,)1()(<<--=-a n u a n h n，通过直接计算卷积和的办法，试确定单位抽样响应为 )(n h 的线性移不变系统的阶跃响应。
4. 判断下列每个序列是否是周期性的,若是周期性的,试确定其周期：)6()( )( )n 313si n()( )()873cos()( )(ππππ-==-=n j e n x c A n x b n A n x a分析：序列为)cos()(0ψω+=n A n x 或)sin()(0ψω+=n A n x 时，不一定是周期序列，nmm m n n y n - - -∞ = - ⋅ = = ≥ ∑ 2 31 2 5 . 0 ) ( 01当 3 4n m nm m n n y n 2 2 5 . 0 ) ( 1⋅ = = - ≤ ∑ -∞ = - 当 aa a n y n a a an y n n h n x n y a n u a n h n u n x m m nnm mn -==->-==-≤=<<--==∑∑--∞=---∞=--1)(11)(1)(*)()(10,)1()()()(:1时当时当解①当=0/2ωπ整数，则周期为0/2ωπ；②；为为互素的整数）则周期、（有理数当 , 2 0Q Q P QP =ωπ ③当=0/2ωπ无理数 ，则)(n x 不是周期序列。
解：（1）0142/3πω=，周期为14 （2）062/13πω=，周期为6 （2）02/12πωπ=，不是周期的 7.（1）[][]12121212()()()()()()[()()]()()()()[()][()]T x n g n x n T ax n bx n g n ax n bx n g n ax n g n bx n aT x n bT x n =+=+=⨯+⨯=+所以是线性的T[x(n-m)]=g(n)x(n-m) y(n-m)=g(n-m)x(n-m) 两者不相等，所以是移变的y(n)=g(n)x(n) y 和x 括号内相等，所以是因果的。
数字信号处理习题集(附答案)
第一章数字信号处理概述简答题：1．在A/D变换之前和D/A变换之后都要让信号通过一个低通滤波器，它们分别起什么作用？答：在A/D变化之前为了限制信号的最高频率，使其满足当采样频率一定时，采样频率应大于等于信号最高频率2倍的条件。
此滤波器亦称为“抗混叠”滤波器。
在D/A变换之后为了滤除高频延拓谱，以便把抽样保持的阶梯形输出波平滑化，故又称之为“平滑”滤波器。
判断说明题：2．模拟信号也可以与数字信号一样在计算机上进行数字信号处理，自己要增加一道采样的工序就可以了。
（）答：错。
需要增加采样和量化两道工序。
3．一个模拟信号处理系统总可以转换成功能相同的数字系统，然后基于数字信号处理理论，对信号进行等效的数字处理。
（）答：受采样频率、有限字长效应的约束，与模拟信号处理系统完全等效的数字系统未必一定能找到。
因此数字信号处理系统的分析方法是先对抽样信号及系统进行分析，再考虑幅度量化及实现过程中有限字长所造成的影响。
故离散时间信号和系统理论是数字信号处理的理论基础。
第二章 离散时间信号与系统分析基础一、连续时间信号取样与取样定理计算题：1．过滤限带的模拟数据时，常采用数字滤波器，如图所示，图中T 表示采样周期（假设T 足够小，足以防止混叠效应），把从)()(t y t x 到的整个系统等效为一个模拟滤波器。
（a ） 如果kHz rad n h 101,8)(=π截止于，求整个系统的截止频率。
（b ） 对于kHz T 201=，重复（a ）的计算。
解 （a ）因为当0)(8=≥ωπωj e H rad 时，在数 — 模变换中)(1)(1)(Tj X Tj X Te Y a a j ωω=Ω=所以)(n h 得截止频率8πω=c 对应于模拟信号的角频率c Ω为8π=ΩT c因此 Hz Tf c c 6251612==Ω=π 由于最后一级的低通滤波器的截止频率为Tπ，因此对T8π没有影响，故整个系统的截止频率由)(ωj e H 决定，是625Hz 。
课后习题答案及讲解
第4章课后习题答案及讲解(总6页)--本页仅作为文档封面，使用时请直接删除即可----内页可以根据需求调整合适字体及大小--4-2 根据图P4—1所示的调制信号波形，试画出DSB及AM信号的波形图，并比较它们分别通过包络检波器后的波形差别。
解：DSB信号及包络检波后输出AM信号及包络栓波后输出由此可见，对DSB信号采用包络检波法不能正确还原基带信号。
4-3已知调制信号m(t)=cos(2000πt)+cos(4000πt)载波为cos104πt，进行单边带调制，试确定该单边带信号的表示式，并画出频谱图。
解：DSB信号为:S DSB(t)= [cos(2000πt)+ cos(4000πt)] cos104πt= 1/2[cos(12000πt)+cos(8000πt)]+1/2[cos(14000πt)+cos(6000πt)]SSB信号为：上边带S SSB(t)= 1/2·cos(12000πt)+ 1/2·cos(14000πt)-8000π 0 6000πω下边带S SSB(t)= 1/2·cos(6000πt)+ 1/2·cos(8000πt)-14000π 0 12000πω4-6 某调制系统如图P4-4所示。
为了在输出端同时分别得到f1(t)及f2(t)，试确定接收端的c1(t)和c2(t)。
（发送（接收解：该调制系统采用相干解调，设c1(t)=cos(ω1t+φ1)则接收端相乘器输出r1(t)=[f1(t) cosω0t + f2(t) sinω0t] cos(ω1t+φ1)= f1(t) cosω0t cos(ω1t+φ1) + f2(t) sinω0t cos(ω1t+φ1)=1/2 f1(t) [ cos(ω0t+ω1t+φ1)+ cos(ω0t- ω1t- φ1)]+1/2 f2(t) [ sin(ω0t+ω1t+φ1)+ sin(ω0t- ω1t- φ1)]若要经过低通滤波器后得到f1(t)，应有ω1=ω0，φ1=0，即c1(t)= cosω0t同理可得c2(t)= sinω0t思考题：4-11 什么是频分复用答：频分复用(Frequency Division Multiplexing) 是按频率分割多路信号的方法，即将信道的可用频带分成若干互不交叠的频段，每路信号占据其中的一个频段。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Tc = 2TF +1024 次复数乘计算时间 = 2 × 0.1536×10−3 +10×10−9 ×1024
= 0.317 44 ms 可实时处理的信号最高频率 fmax 为
f
max
≤
1 2
Fs
=
1 2
·
由此可见，
1024 = 用TDc SP
1 ·3.1158 MHz=1.6129 MHz 2 专用单片机可大大提高信号处理速度。
令
y(n)=x1(n)+jx2(n)
Y(k)=DFT［y(n)］
k=0, 1, …, N－1
则
X1(k)
=
DFT[ x1 (n)]
=
Yep (k )
=
1 2
[Y
(k)
+Y
*(N
−
k )]
jX 2 (k)
=
DFT[ jx2 (n)]
=
Yep (k )
=
1 2
[Y
(k)
−Y*(N
−
k )]
2N 点 DFT［x（n）］=X（k）可由 X1(k)和 X2(k)得到
x2(n)=x(2n+1)
n=0, 1, …, N－1
2
根据 DIT-FFT 的思想， 只要求得 x1(n)和 x2(n)的 N 点 DFT， 再经过简单 的一级蝶形运算就可得到 x(n)的 2N 点 DFT。因为 x1(n)和 x2(n)均为实序列， 所 以根据 DFT 的共轭对称性， 可用一次 N 点 FFT 求得 X1(k)和 X2(k)。 具体方法 如下：
（2） 与（1）相同， 设
x1(n)=x(2n)
n=0, 1, …, N－1
x2(n)=x(2n+1)
n=0, 1, …, N－1
X1(k)=DFT［x1(n)］
k=0, 1, …, N－1
X2(k)=DFT［x2(n)］
k=0, 1, …, N－1
则应满足关系式
X (k) = X1(k) + W2kN X 2 (k)
积的计算时间 Tc 约为 Tc = 2TF +1024次复数乘计算时间 = 71680 µs + 4×1024 µs = 65536 µs
所以， 每秒钟处理的采样点数（即采样速率）
Fs
<
1024 65536 ×10−6
= 15
625
次/秒
由采样定理知， 可实时处理的信号最5 2
X (k) = X1(k) + W2kN X 2 (k)
X (k + N ) = X1(k) −W2kN X 2 (k)
k = 0,1,L, N −1
这样， 通过一次 N 点 IFFT 计算就完成了计算 2N 点 DFT。 当然还要进行由
Y(k)求 X1(k)、 X2(k)和 X(k)的运算(运算量相对很少)。
用 FFT 计算 1024 点 DFT 所需计算时间 TF 为
TF
=
5 ×10−6
×
N 2
lbN
+
NlbN
×10−6
= 5×10−6 × 1024 ×10 +1024×10×10−6 2
= 30.72 ms
快速卷积时， 需要计算一次 N 点 FFT（考虑到 H(k)= DFT［h(n)］已计算好存
入内存）、 N 次频域复数乘法和一次 N 点 IFFT。 所以， 计算 1024 点快速卷
所以，
DSP
在数字信号处理领域得到广泛应用。机器周期小于 1 ns 的 DSP 产品已上市， 其
处理速度更高。
3. 解: 因为 x(n)和 y(n)均为实序列， 所以， X(k)和 Y(n)为共轭对称序列，
jY(k)为共轭反对称序列。 可令 X(k)和 jY(k)分别作为复序列 F(k)的共轭对称分量
故 x(n) = 1 [ f (n) + f ∗(n)] 2
y(n) = 1 [ f (n) − f ∗ (n)] 2j
4. 解: 本题的解题思路就是 DIT-FFT 思想。
（1）在时域分别抽取偶数和奇数点 x(n)，得到两个 N 点实序列 x1(n)和 x2(n)：
x1(n)=x(2n)
n=0, 1, …, N－1
3
X1(k)
=
1 2
[
X
(k)
+
X
(k
+
N
)]
X2
(k)
=
1 2
[X
(k)
+
X
(k
+
N
)]W2−Nk
② 由 X1(k)和 X2(k)构成 N 点频域序列 Y(k):
Y(k)=X1(k)+jX2(k)=Yep(k)+Yop(k)
其中， Yep(k)=X1(k), Yop(k)=jX2(k)， 进行 N 点 IFFT， 得到
X (k + N ) = X1(k) −W2kN X 2 (k)
k = 0,1,L, N −1
由上式可解出
X1(k)
=
1 2
[
X
(k)
+
X
(k
+
N )]
X
2
(k)
=
1 2
[X
(k)
+
X
(k
+
N
)]W2−Nk
k = 0,1, 2,L, N −1
由以上分析可得出运算过程如下：
（1）由 X(k)计算出 X1(k)和 X2(k):
1
TD=10×10－9×10242+10×10－9×1 047 552=20.961 28 ms
用 FFT 计算 1024 点 DFT 所需计算时间 TF 为
TF
= 10 ×10−9
×
N 2
l
bN
+10 ×10−9
×
N
l bN
= 10−8 × 1024 ×10 +10−8 ×1024×10 2
= 0.1536 ms 快速卷积计算时间 Tc 约为
和共轭反对称分量， 即
F(k)=X(k)+jY(k)=Fep(k)+Fop(k) 计算一次 N 点 IFFT 得到
f(n)=IFFT［F(k)］=Re［f(n)］+j Im［f(n)］ 由 DFT 的共轭对称性可知
Re［f(n)］=IDFT［Fep(k)］=IDFT［X(k)］=x(n) j Im［f(n)］=IDFT［Fop(k)］=IDFT［jY(k)］=jy(n)
第 4 章 快速傅里叶变换（FFT）
习题答案
1. 解: 当 N=1024=210 时， 直接计算 DFT 的复数乘法运算次数为
N2=1024×1024=1 048 576 次
复数加法运算次数为
N（N－1）=1024×1023=1 047 552 次
直接计算所用计算时间 TD 为
TD=4×10－6×10242+1 047 552×10－6=5.241 856 s
= 7.8125
kHz
应当说明， 实际实现时， fmax 还要小一些。 这是由于实际中要求采样频 率高于奈奎斯特速率， 而且在采用重叠相加法时， 重叠部分要计算两次。 重 叠部分长度与 h(n)长度有关， 而且还有存取数据和指令周期等消耗的时间。 2. 解: 与第 1 题同理。
直接计算 1024 点 DFT 所需计算时间 TD 为