课后习题及答案_第4章快速傅里叶变换--习题答案
- 格式:pdf
- 大小:73.63 KB
- 文档页数:5
数字信号处理教程课后习题答案目录第一章离散时间信号与系统第二章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 的分段处理。
第四章习题4.6求下列周期信号的基波角频率Ω和周期T解 ⑴角频率为Ω = IOO rad∕s,周期丁=盲=p÷ξ ⑵角频率为I fi=号■rad∕s,周期= 4 s(3) 角频率为Ω = 2 rad 倉,周期T = ~ = Tr S (4) 角频率为Q =兀rad∕ s,周期T=^ = 2 sΩ(5) 角频率为 Ω — rad∕s*周期 T=-^ = 8 s4 12⑹角频率为C =話rad∕s,周期T = -jy = 60 s4.7用直接计算傅里叶系数的方法, 求图4-15所示周期函数 的傅里叶系数(三角形式或指数形式)(1) e j100t(2) cos[,t - 3)](3) cos(2t) sin(4t) ⑷ cos(2 兀 t) +cos(3πt) +cos(5 兀 t)(5)π π cos( t) sin( t)2 4(6)JEJITEcos( t) cos( t) cos( t)2 35-2 -1 O 12 3 r(IJ)图4-15f>~ 十解 ⑴周期T = 4,1Ω = Y =亍r 则有H ,4⅛ - 1 ≤ r ≤ 4⅛+ 1/⑺=II∣07 4⅛ + 1 < r < 4⅛ + 3由此可得-Tu rt = ~∖ ' τ fit) cost nΩt)dt= -∣^∣ /(f)cos(^ψ^)df J- J —⅛ 乙-.:—2 I(2}周期丁=2・0 =年=兀,则有由此可得1 + e -jrhr2π( I - √ )所含有的频率分量)dr =2 J -[2『亍=Wl f(t)sm(ττΩt)dt =1 J -T2——SInnπ (才),= om 小山(竽)出ISin(Jrt) 9fm=! 0,2⅛ ≤ r ≤ 2⅛ + 12⅛ + 1 < r < 2⅛ + 2F ri ]ft1 Γl=TJV Cf)^dr =⅛J r ∣/(r)e-7iβ,dr — -7- Sin(^f)e -dr -I ZJV4.10利用奇偶性判断图4-18示各周期信号的傅里叶系数中扣 =O* ± 1 * + 2・・图 4-18解 (1)由旳⑺的波形可矩Λ<r) =√√-n =-∕l (f ⊂f)亠 IU Jr = f(t)cos( riΩt )df 则有丿 丁人 ,jj = 0.1,2,-[仇=0"[J =盘?=应丄=*" =QE=仇=仏=*八=0 则∕√r)的傅里叶级数中含有的频率分量为奇次余弦波亠 (2)由f 2(t)的波形可知则有— ■ ??f(t)s}n(tιΩt )d r ⅛ =A rz fl , J Tni JJO则f 2(t)的傅里叶级数中含有的频率分量为正弦波*(3)由 f 3(t)的波形可⅛l∕3<f) = f 3(~r)则有Γ⅛ = 0, n/(z)cos( fiΩt >d;(4)% 4召=亍即ΛG)的傅里叶级数中含有的频率分量为偶次余弦波* 由/<(0的波形可知,人⑺为奇谐函数■即fdι) =一 fZ 土 £)b 2 = h A = b 6 =・*・=0则有 U即人")的傅里叶级数中只含有奇次谐波•包括正弦波和余弦披"4-11 求u(t)的三角形式傅里叶系数。
第一章习题习题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 ===错误!未找到引用源。
第一章快速傅里叶变换(FFT )4.1 填空题(1)如果序列)(n x 是一长度为64点的有限长序列)630(≤≤n ,序列)(n h 是一长度为128点的有限长序列)1270(≤≤n ,记)()()(n h n x n y *=(线性卷积),则)(n y 为 点的序列,如果采用基FFT 2算法以快速卷积的方式实现线性卷积,则FFT 的点数至少为 点。
解:64+128-1=191点; 256(2)如果一台通用机算计的速度为:平均每次复乘需100s μ,每次复加需20s μ,今用来计算N=1024点的DFT )]([n x 。
问直接运算需( )时间,用FFT 运算需要( )时间。
解:①直接运算:需复数乘法2N 次,复数加法)(1-N N 次。
直接运算所用计算时间1T 为s s N N N T 80864.12512580864020110021==⨯-+⨯=μ)(② 基2FFT 运算:需复数乘法N N2log 2次,复数加法N N 2log 次。
用FFT 计算1024点DTF 所需计算时间2T 为s s N N N NT 7168.071680020log 100log 2222==⨯+⨯=μ。
(3)快速傅里叶变换是基于对离散傅里叶变换 和利用旋转因子k Nj e π2-的来减少计算量,其特点是 _______、_________和__________。
解:长度逐次变短;周期性;蝶形计算、原位计算、码位倒置 (4)N 点的FFT 的运算量为复乘 、复加 。
解:N NL N mF 2log 22==;N N NL aF 2log ==4.2 选择题1.在基2DIT —FFT 运算中通过不断地将长序列的DFT 分解成短序列的DFT ,最后达到2点DFT 来降低运算量。
若有一个64点的序列进行基2DIT —FFT 运算,需要分解 次,方能完成运算。
A.32 B.6 C.16 D. 8 解:B2.在基2 DIT —FFT 运算时,需要对输入序列进行倒序,若进行计算的序列点数N=16,倒序前信号点序号为8,则倒序后该信号点的序号为 。
第一章 傅里叶分析部份习题解答及参考答案[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.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 流程图,即迭代次数计算子程序。
信号与系统课后习题参考答案1试分别指出以下波形就是属于哪种信号?题图1-11-2试写出题1-1图中信号得函数表达式。
1-3已知信号与波形如题图1-3中所⽰,试作出下列各信号得波形图,并加以标注。
题图1-3⑴⑵⑶⑷⑸⑹⑺⑻⑼1-4已知信号与波形如题图1-4中所⽰,试作出下列各信号得波形图,并加以标注。
题图1-4⑴⑵⑶⑷⑸⑹⑺⑻⑼1-5已知信号得波形如题图1-5所⽰,试作出信号得波形图,并加以标注。
题图1-51-6试画出下列信号得波形图:⑴⑵⑶⑷1-7试画出下列信号得波形图:⑴⑵⑶⑷⑸⑹1-8试求出以下复变函数得模与幅⾓,并画出模与幅⾓得波形图。
⑴⑵⑶⑷1-9已知信号,求出下列信号,并画出它们得波形图。
1-10试作出下列波形得奇分量、偶分量与⾮零区间上得平均分量与交流分量。
题图1-101-11试求下列积分:⑴⑵⑶⑷⑸⑹1-12试求下列积分:⑴⑵⑴(均为常数)⑵⑶⑷⑸⑹⑺⑻1-14如题图1-14中已知⼀线性时不变系统当输⼊为时,响应为。
试做出当输⼊为时,响应得波形图。
题图1-14 1-15已知系统得信号流图如下,试写出各⾃系统得输⼊输出⽅程。
题图1-151-16已知系统⽅程如下,试分别画出她们得系统模拟框图。
⑴⑵⑶1-17已知⼀线性时不变系统⽆起始储能,当输⼊信号时,响应,试求出输⼊分别为与时得系统响应。
第⼆章习题2-1试计算下列各对信号得卷积积分:。
⑴(对与两种情况)⑵⑶⑷⑸⑹2-2试计算下列各对信号得卷积与:。
⑴(对与两种情况)⑵⑶⑷⑸⑹2-3试计算下图中各对信号得卷积积分:,并作出结果得图形。
题图2-32-4试计算下图中各对信号得卷积与:,并作出结果得图形。
题图2-42-5已知,试求:⑴⑵⑶2-7系统如题图2-7所⽰,试求系统得单位冲激响应。
已知其中各⼦系统得单位冲激响应分别为:题图2-72-8设已知LTI 系统得单位冲激响应,试求在激励作⽤下得零状态响应。
2-9⼀LTI 系统如题图2-9所⽰,由三个因果LTI ⼦系统级联⽽成,且已知系统得单位样值响应如图中。
1 信号分析与处理课后习题答案第五章快速傅里叶变换1.1.如果一台通用计算机的速度为平均每次复乘需要如果一台通用计算机的速度为平均每次复乘需要50us 50us,每次复加需要,每次复加需要10us 10us,,用来就散N=1024点的DFT DFT,问:,问:(1)直接计算需要多少时间?用FFT 计算呢?(2)照这样计算,用FFT 计算快速卷积对信号进行处理是,估计可实现实时处理的信号最高频率?解:分析:直接利用DFT 计算:复乘次数为N 2,复加次数为N(N-1);利用FFT 计算:复乘次数为20.5log N N ,复加次数为2log N N ;(1)直接DFT 计算:复乘所需时间2215010245052.4288T N us us s=´=´=复加所需时间2(1)101024(10241)1010.47552T N N us us s=-´=-´=所以总时间1262.90432DFT T T T s=+=FFT 计算:复乘所需时间3220.5log 500.51024log 1024500.256T N N us us s =´=´´´=复加所需时间422log 101024log 1024100.1024T N N us us s =´=´´=所以总时间为340.3584FFT T T T s =+=(2)假设计算两个N 长序列1()x n 和2()x n 的卷积计算过程为如下:第一步:求1()X k ,2()X k ;所需时间为2FFTT ´第二步:计算12()()()X k X k X k =·,共需要N 次复乘运算所需时间为501024500.0512To N us us s=´=´=第三步:计算(())IFFT X k ,所需时间为FFTT 所以总时间为230.35840.0512 1.1264FFT T T To s s s=´+=´+=容许计算信号频率为N/T=911.3Hz 2.2.设设x(n)x(n)是长度为是长度为2N 的有限长实序列,()X k 为x(n)x(n)的的2N 点得DFT DFT。
通信原理[张会⽣]课后习题答案解析思考题1-1 什么是通信?常见的通信⽅式有哪些?1-2 通信系统是如何分类的?1-3 何谓数字通信?数字通信的优缺点是什么?1-4 试画出模拟通信系统的模型,并简要说明各部分的作⽤。
1-5 试画出数字通信系统的⼀般模型,并简要说明各部分的作⽤。
1-6 衡量通信系统的主要性能指标是什么?对于数字通信具体⽤什么来表述?1-7 何谓码元速率?何谓信息速率?它们之间的关系如何?习题1-1 设英⽂字母E出现的概率=0.105,X出现的概率为=0.002,试求E和X的信息量各为多少?1-2 某信源的符号集由A、B、C、D、E、F组成,设每个符号独⽴出现,其概率分别为1/4、1/4、1/16、1/8、1/16、1/4,试求该信息源输出符号的平均信息量。
1-3 设⼀数字传输系统传送⼆进制信号,码元速率RB2=2400B,试求该系统的信息速率Rb2=?若该系统改为传送16进制信号,码元速率不变,则此时的系统信息速率为多少?1-4 已知某数字传输系统传送⼋进制信号,信息速率为3600b/s,试问码元速率应为多少?1-5 已知⼆进制信号的传输速率为4800b/s,试问变换成四进制和⼋进制数字信号时的传输速率各为多少(码元速率不变)?1-6 已知某系统的码元速率为3600kB,接收端在l⼩时内共收到1296个错误码元,试求系统的误码率=?1-7 已知某四进制数字信号传输系统的信息速率为2400b/s,接收端在0.5⼩时内共收到216个错误码元,试计算该系统=?l-8 在强⼲扰环境下,某电台在5分钟内共接收到正确信息量为355Mb,假定系统信息速率为1200kb/s。
(l)试问系统误信率=?(2)若具体指出系统所传数字信号为四进制信号,值是否改变?为什么?(3)若假定信号为四进制信号,系统传输速率为1200kB,则=?习题答案第⼀章习题答案1-1 解:1-2 解:1-3 解:1-4 解:1-5 解:1-6 解:1-7 解:1-8 解:思考题2-1 什么是狭义信道?什么是⼴义信道?(答案)2-2 在⼴义信道中,什么是调制信道?什么是编码信道?2-3 试画出调制信道模型和⼆进制⽆记忆编码信道模型。
傅里叶光学教程(黄婉云)课后习题解答第一章数学基础知识在信息光学中,有一些广泛使用的函数,包括脉冲函数、梳状函数等,用于描述各种物理量。
另外还有一些重要的数学运算,如卷积、相关、傅里叶变换等,用于讨论和分析各种物理过程。
本章主要介绍这些函数及计算方法。
1.1 常用函数1. 阶跃函数(Step function)x,00,,a,xx1,一维: step()0,,,aa2,x,10,,a,a的正负决定阶跃函数的取向,阶跃函数作用如同一个开关,可在某点开启或关闭一个函数。
2. 符号函数(Sign function)x,10,,a,xx, sgn()00,,,aa,x,,,10,a,a的正负决定符号函数的取向,符号函数用来改变一个变量或函数的正负。
xxx,,sgn11阶跃函数与符号函数的关系: stepxx()sgn(),,223. 矩形函数(Rectangle function)a,1x,x,rect(),一维: 2,a,0others,表示函数以0为中心,宽度为a(a>0),高度为1的矩形。
在时间域,矩形函数可以描写照相机快门;在空间域,矩形函数可以描写无限大不透明屏上单缝的透过率,故被称为门函数。
并且矩形函数可以作为截取函数。
2xyxy二维: rectrectrect(,)()(),,ababab,1,,xy,,,, 22,,0,others,二维矩形函数可以描写无限大不透明屏上矩形孔的透过率。
矩形函数可以移位和改变比例以及高度:a,hxx,,,xx,,00 hrect(),2,a,0,others,4. 三角函数(Triangle function),xx1,,xa,一维: tri(),a,a,0others,xyxx二维: ,,,,(,)()()abab,xxxy(1)(1),,1,,,, ,abab,,0,others,可用来表示一个光瞳为矩形的非相干成像系统的光学传递函数。
5. Sinc 函数(Sinc function)xxasin(),,,sin(),0ca一维: axa,xxc,,0sin()1时, axxnac,,,时,sin()0 axyxy二维: sin(,)sin()sin()ccc, abab3可用来描述单缝和矩孔的夫琅和费衍射振幅分布,其平方表示衍射图样6. 高斯函数(Gaussian function)x2一维: Gausxa()exp[/],,,,,a高斯函数也称为正态分布函数。
第一章数字信号处理概述简答题: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 。