信号处理习题答案VIP

数字信号处理习题解答第二章数据采集技术基础有一个理想采样系统,其采样角频率Ωs=6π,采样后经理想低通滤波器HajΩ还原,其中现有两个输入,x1t=cos2πt,x2t=cos5πt;试问输出信号y1t,y2t有无失真为什么

分析：要想时域采样后能不失真地还原出原信号,则采样角频率Ωs必须大于等于信号谱最高角频率Ωh的2倍,即满足Ωs≥2Ωh;解：已知采样角频率Ωs=6π,则由香农采样定理,可得因为x1t=cos2πt,而频谱中最高角频率,所以y1t无失真；因为x2t=cos5πt,而频谱中最高角频率,所以y2t失真;设模拟信号xt=3cos2000πt+5sin6000πt+10cos12000πt,求：该信号的最小采样频率；若采样频率fs=5000Hz,其采样后的输出信号；分析：利用信号的采样定理及采样公式来求解;eq\o\ac○,1采样定理采样后信号不失真的条件为：信号的采样频率fs不小于其最高频率fm的两倍,即fs≥2fmeq\o\ac○,2采样公式解：1在模拟信号中含有的频率成分是f1=1000Hz,f2=3000Hz,f3=6000Hz∴信号的最高频率fm=6000Hz由采样定理fs≥2fm,得信号的最小采样频率fs=2fm=12kHz2由于采样频率fs=5kHz,则采样后的输出信号说明：由上式可见,采样后的信号中只出现1kHz和2kHz的频率成分,即若由理想内插函数将此采样信号恢复成模拟信号,则恢复后的模拟信号可见,恢复后的模拟信号yt不同于原模拟信号xt,存在失真,这是由于采样频率不满足采样定理的要求,而产生混叠的结果;第三章傅里叶分析傅里叶变换概述习题设序列xn=δn-m,求其频谱Xejω,并讨论其幅频和相频响应分析：求解序列的频谱有两种方法：eq\o\ac○,1先求序列的z变换Xz,再求频谱,即Xejω为单位圆上的z变换；eq\o\ac○,2直接求序列的傅里叶变换解：对序列xn先进行z变换,再求频谱,得则若系统的单位采样响应hn=xn,则系统的频率响应故其幅频和相频响应如图分别为幅频响应相频响应ωωφωωHejω1由图可见,该系统的频率响应具有单位幅值以及线性相位的特点;设xn的傅里叶变换为Xejω,试利用Xejω表示下列序列的傅里叶变换：分析：利用序列翻褶后的时移性质和线性性质来求解,即,解：1由于,,则故2由于故设Xejω是如图所示的信号xn的傅里叶变换,不必求出Xejω,试完成下列计算：分析：利用序列傅里叶变换的定义以及帕塞瓦定理来求解;序列的傅里叶变换公式为：正变换反变换帕塞瓦定理解：1由傅里叶正变换公式可知ω=0,则2由于ej0=1,则由傅里叶反变换公式可知n=0,故由帕塞瓦定理,得周期序列的离散傅里叶级数DFS如图所示,序列xn是周期为6的周期性序列,试求其傅里叶级数的系数;分析：利用DFS的定义求解,即,其中k=0~N-1解：已知N=6,则由DFS的定义得对上式依次取k=0~5,计算求得设,令,,试求与的周期卷积;分析：可以利用列表法求解,直观方便;由于eq\o\ac○,只要将列表中对应于某个n的一行中的值和第一行中与之对应的值相乘,然后再将所有乘积结果相加,就得到此n的值解：注意：本题需要利用下一节中有限序列与周期序列的关系以及序列循环移位的概念;在一个周期N=6内的计算卷积值eq\o\ac○,则与的周期卷积值n=0~5如下表所示：离散傅里叶变换DFT已知xn如图所示,为{1,1,3,2},试画出序列x-n5,x-n6R6n,xn3R3n,xn6,xn-35R5n和xn7R7n的略图;分析：此题需注意周期延拓的数值,也就是xnN中N的数值;如果N比序列的点数多,则需补零；如果N比序列的点数少,则需将序列按N为周期进行周期延拓,造成混叠相加形成新的序列;解：各序列的略图如图所示;试求下列有限长序列的N点离散傅里叶变换闭合形式表达式：分析：利用有限长序列的DFT的定义,即解：1因为,所以2因为,所以3由,得注意：为了便于求解,必须利用代数简化法消除掉上式中的变量n;所以4注意：本题可利用上题的结论来进行化简;由,则根据第3小题的结论：若则与上题同理,得所以试画出图示的两个有限长序列的六点循环卷积;分析：本题可以直接利用循环卷积的公式求解,也可以利用循环移位的概念来求解,即：有限长序列xn左移mm为正整数位的循环移位定义为且移位时,在主值区间n=0~N-1内,当某序列值从区间的一端移出时,与它同值的序列值又从区间的另一端移入;解：由循环卷积的定义,可知eq\o\ac○,6eq\o\ac○,eq\o\ac○,则根据循环移位的概念,将序列x1n循环右移3个单位后乘以3并取其主值序列n=0~5即可,其结果如图所示;如图所示的5点序列xn,试画出：xnxnxneq\o\ac○,5xnxneq\o\ac○,10xn分析：本题可由图解法来计算循环卷积,并利用循环卷积来求解线性卷积;同时应注意循环卷积代替线性卷积的条件：设两个有限长序列xn、hn的点数分别为N和M,其循环卷积的长度为L,则要用循环卷积代替线性卷积的条件是：循环卷积的长度L必须不小于线性卷积的长度N+M-1,即L≥N+M-1否则,在循环卷积周期延拓时会产生混叠;解：由于xn是5点序列,所以xnxn是5+5-1=9点序列,因此,xneq\o\ac○,10xn的前9个点n=0,1,…,8就是xnxn值,后一个点n=9为零,因为L点循环卷积等于线性卷积结果的L点周期延拓、混叠相加后的主值区间内的序列L可以是任意整数值;其运算结果分别如图a、b、c所示;已知两个有限长序列为试作图表示xn,yn以及fn=xneq\o\ac○,7yn;分析：直接利用循环卷积公式或图解法求解;解：其结果如图所示;习题已知xn是N点有限长序列,且Xk=DFTxn;现将它补零扩展成长度为rN点的有限长序列yn,即试求rN点DFTyn与Xk的关系;分析：利用DFT定义求解;yn是rN点序列,因而结果相当于在频域序列进行插值;解：由可得所以在一个周期内,Yk的采样点数是Xk的r倍Yk的周期为rN,相当于在Xk的每两个值之间插入r-1个其它的数值不一定为零,而当k为r的整数l倍时,Yk与相等;习题频谱分析的模拟信号以8kHz被采样,计算了512个采样点的DFT,试确定频谱采样之间的间隔,并证明你的回答;分析：利用频域采样间隔F0和时域采样频率fs以及采样点数N的关系fs=NF0;证：由得其中Ωs是以角频率为变量的频谱周期,Ω0是频谱采样之间的频谱间隔;又则对于本题有fs=8kHz,N=512所以习题设有一个频谱分析用的信号处理器,采样点数必须为2的整数幂,假定没有采用任何特殊数据处理措施,要求频率分辨力≤10Hz,如果采用的采样时间间隔为,试确定：最小记录长度；所允许处理信号的最高频率；在一个记录中的最小点数;分析：采样间隔T和采样频率fs满足fs=1/T,记录长度T0和频域分辨力F0的关系为T0=1/F0,采样定理为fs≥2fhfh为信号最高频率分量,一个记录中最少的采样总数N满足解：1因为T0=1/F0,而F0≤10Hz,所以即最小记录长度为;2因为,而fs≥2fh所以即允许处理信号的最高频率为5kHz;3又因N必须为2的整数幂所以一个记录中的最少点数为N=210=1024;快速傅里叶变换FFT如果一台通用计算机的速度为平均每次复乘5μs,每次复加μs,用它来计算512点的DFTxn,问直接计算需要多少时间,用FFT运算需要多少时间

分析：eq\o\ac○,1直接利用DFT计算：复乘次数为N2,复加次数为NN-1；eq\o\ac○,2利用FFT计算：复乘次数为,复加次数为；解：1直接计算复乘所需时间复加所需时间所以2用FFT计算复乘所需时间复加所需时间所以已知Xk,Yk是两个N点实序列xn,yn的DFT值,今需要从Xk,Yk求xn,yn的值,为了提高运算效率,试用一个N点IFFT运算一次完成;分析：我们来组成一个新的序列Xk+jYk序列,则有它的实部即为实序列xn,虚部即为实序列yn;解：依据题意,可知取序列对Zk作N点IFFT可得序列zn;又根据DFT线性性质由原题意可知,xn,yn都是实序列;再根据zn=xn+jyn,可得习题,N=16时,画出基-2按时间抽取法DIT及按频率抽取法DIF的FFT流图时间抽取采用输入倒位序,输出自然数顺序,频率抽取采用输入自然顺序,输出倒位序;分析：eq\o\ac○,1DIF法与DIT法的异同：不同点：DIT与DIF的基本蝶形图不同,DIF的复数乘法出现在减法之后,DIT的复数乘法出现在减法之前；相同点：DIT与DIF的运算量是相同的；eq\o\ac○,2DIF法与DIT法的关系：它们的基本蝶形是互为转置的;解：1按时间抽取DIT如图所示2按频率抽取DIF如图所示课堂思考题若是因果稳定序列,求证：证：设则由时域卷积定理,得即令上式的左右两边n=0,得又傅里叶反变换公式,得,则,所以课堂思考题在N=16时按时间抽取的基-2FFT算法中,若输入序列xn采用倒位序,输出序列Xk采用自然数顺序,试写出输入序列xn的排列顺序,并简述理由;答：N=16的基-2FFT算法中,输入序列xn倒位序排列顺序为x0、x8、x4、x12、x2、x10、x6、x14、x1、x9、x5、x13、x3、x11、x7、x15;其倒位序排序规则如表所示：自然顺序n自然顺序二进制数倒位序二进制数倒位序顺序0000000000100011000820010010043001111001240100001025010110101060110011067011111101481000000119100110019101010010151110111101131211000011313110110111114111001117151111111115第五章时域分析随机相位正弦波式中,x0,ω均为常数,φ在0~2π内随机取值,试求其自相关函数并作图;分析：利用自相关函数的定义求解,即解：由自相关函数的定义式,得可见,该随机相位正弦波的自相关函数只与角频率ω有关,而不含相位信息,这表明：正弦函数的自相关函数为失去了相位信息的同频率余弦函数;其自相关函数图形如图所示;x02/2x02/2Rxxτττ两个随机相位正弦波式中,A0,B0,ω,φ均为常数,θ在0~2π内的取值概率相同,即满足试求其互相关函数并作图;分析：利用互相关函数的定义求解,即解：由互相关函数的定义式,得可见,两个正弦函数的互相关函数仍为同频率的余弦函数,其最大峰值出现在τ=φ/ω处;RxyRxyττφ/ω第六章数字滤波器设计已知模拟滤波器的模方函数求模拟滤波器的传递函数;分析：利用模拟滤波器的模方函数|HjΩ|2与其传递函数Hs之间的关系式求解,即解：将s=jΩ,即Ω2=s2代入|HjΩ|2,得可见,系统有四个极点s1,2=±3,s3,4=±4和两对零点z1,2=±j2;为了得到一个稳定的滤波系统,则将左半平面的极点分配给Hs；并取虚轴上的一对共轭零点作为Hs的零点,以保证Hs收敛,故模拟滤波器的传递函数为试设计一个巴特沃思BW低通模拟滤波器,使滤波器的幅度响应在通带截止频率105rad/s处的衰减不大于3dB,在阻带截止频率4×105rad/s处的衰减不小于35dB;分析：按照§中所述的巴特沃思低通滤波器的设计过程来实现;eq\o\ac○,1先确定滤波器的阶数N由于公式1则滤波器的阶数公式2注意：N为正整数且截止频率公式3eq\o\ac○,2求解位于左半S平面上的极点公式4eq\o\ac○,3确定N阶巴特沃斯低通滤波器的传递函数公式5解：eq\o\ac○,1先确定滤波器的阶数N由题意可知,Ωp=105rad/s时,通带最大衰减αp=3dBΩs=4×105rad/s时,阻带最小衰减αs=35dB则代入公式1,求得参数γ和λ将参数γ、λ、Ωp和Ωs代入公式2,则滤波器的阶数将参数N、γ和Ωp代入公式3,可得截止频率eq\o\ac○,2求解位于左半S平面上的极点将参数Ωc和N代入公式4,得极点即eq\o\ac○,3确定巴特沃斯低通滤波器的传递函数Hs将参数N、Ωc和sk代入公式5,得巴特沃斯低通滤波器的传递函数式中Ωc=105rad/s试导出二阶巴特沃斯低通滤波器的系统函数,设Ωc=3rad/s;分析：本题利用模方函数求出其左半S平面极点,而求得系统函数;N阶巴特沃斯低通滤波器的模方函数定义为在上式中代入jΩ=s,可得而HsH-s在左半S平面的极点即为Hs的极点,因此其中,k0由来确定;注意：可以证明,系数k0=ΩcN;解：对于二阶N=2巴特沃斯低通滤波器,其模方函数为令jΩ=s,则有各极点满足则k=1,2时,所得的sk位于左半S平面,即为Hs的极点由以上两个极点构成的系统函数为代入条件,可得k0=9注：k0=Ωc2,故二阶巴特沃斯低通滤波器的系统函数试导出三阶巴特沃斯低通滤波器的系统函数,设Ωc=2rad/s;分析：与习题6.3同理,利用模方函数求出其左半S平面极点,而求得系统函数;解：对于三阶N=3巴特沃斯低通滤波器,其模方函数为令jΩ=s,则有各极点满足不难得知,当k=1,2,3时,相应的极点sk均位于左半S平面;则滤波器的系统函数Hs的极点因此,三阶巴特沃斯低通滤波器的系统函数为设模拟滤波器的系统函数为试利用冲激响应不变法,设计IIR数字低通滤波器;确定AF的传函Hs给定指标确定AF的传函Hs给定指标求解AF的单位冲激响应ht取拉氏反变换获得DF的单位冲激响应序列hn采样获得DF的传函HzZ变换令t=nT解：将Hs展开成部分分式,得对Hs取拉氏反变换,得对ht作周期为T的等间隔采样,得对hn取Z变换,得IIR数字低通滤波器的系统函数为设有一模拟滤波器采样周期T=2,试用双线性变换法将它转变为数字系统函数Hz;分析：双线性变换法是模拟系统函数的S平面和数字系统函数的Z平面之间是一一对应的关系,消除了频谱的混叠现象,其变换关系为解：将T=2代入变换公式,可得则数字系统函数用双线性变换法设计一个数字低通滤波器,采样频率fs=,截止频率fc=400Hz;分析：按照§中所述的采用双线性变换法的设计过程来实现;eq\o\ac○,1利用关系式ω=TΩ将给定的模拟域频率指标转化为数字域频率指标eq\o\ac○,2利用如下的预畸变补偿公式将数字域频率指标变换为补偿后的模拟域频率指标eq\o\ac○,3按补偿后的模拟域频率指标设计三阶巴特沃斯模拟滤波器Hs参见例6.2.4eq\o\ac○,4利用双线性变换公式,将模拟滤波器Hs变换为数字滤波器Hz,即T为采样周期解：此数字滤波器的截止频率由预畸变补偿,得相应的模拟滤波器的截止频率由习题6.4可知,三阶巴特沃斯模拟滤波器的系统函数其中,滤波器的系统函数Hs的极点故有将双线性变换公式和代入,可得三阶巴特沃斯数字滤波器的系统函数请选择合适的窗函数及窗宽N来设计一个线性相位低通滤波器要求其阻带最小衰减为45dB,过渡带宽为8π/51,试求出hn设截止频率ωc=π;分析：本题是真正实用的设计题,从中可以看到阻带衰减影响窗形状的选择当然用凯塞窗则可改变β来满足阻带衰减的要求,而窗宽N的选择则影响过渡带宽;按照§中所述的线性相位FIR数字低通滤波器的设计步骤来实现;eq\o\ac○,1给定