第4章模拟信号的数字处理DigitalProcessingofContinuous-TimeSignalsSCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.22020/1/16本章主要内容模拟信号与数字信号的相互转换模拟信号的数字处理模拟系统的数字模仿模拟信号的频谱分析SCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.32020/1/164.1模拟信号数字处理原理方框图作为数字信号处理的第一步,要将现实中连续时间信号变成数字信号。如何做到无损变换预滤波平滑滤波数模转换DAC数字信号处理量化编码采样()xt()xn[]xn[]yn()yt()yn模数转换ADC4.2模拟信号与数字信号的相互转换SCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.52020/1/16模拟信号与数字信号的相互转换目标:无损变换抽样:就是利用周期性抽样脉冲序列p(t),从连续信号xa(t)中抽取一系列的离散值,得到抽样信号(或称抽样数据信号)即离散时间信号,以表示。抽样是模拟信号数字化的第一环节,再经幅度量化编码后即得到数字信号x(n)。)(ˆtxaSCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.62020/1/164.2.1时域采样定理对模拟信号的采样()axtˆ()axt()pt()axtt()pttTˆ()axttS()axt()axt()axt()ptˆ()axttttTSCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.72020/1/16理想采样信号的频谱是以为周期进行延拓并累加而得。如何恢复原信号的频谱?加低通滤波器,传输函数为()aXj0cc()Pj0ss()aXjcssc2sˆ()aXj()aXjs2()02ssTGj11ˆ()()()()2aaasmXjXjPjXjjmTˆ()aXj时域采样定理SCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.82020/1/16理想采样的恢复ˆ()aXjcssc2s()aXjcc()Gj采样信号的频谱混叠()aXj0cc()Pj0ssˆ()aXjssc2s频谱混叠的原因:2scSCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.92020/1/16奈奎斯特采样定理模拟信号为带限信号;最高角频率为如采样频率满足则理想采样信号通过增益为T,截止频率为的低通滤波器,可唯一地恢复原模拟信号采样频率的确定防止混叠并降低运算量为原则加预滤波器2sccˆ()axt()axt2s()Gj()axt(34)scSCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.102020/1/16时域离散信号与理想采样信号之间的关系理想采样信号采样序列频谱之间的关系?ˆ()()()aanxtxnTtnT()()()aatnTxnxtxnT()axttˆ()axtt()xnnSCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.112020/1/16频谱时域离散信号ˆ()()()aanxtxnTtnT()()()aatnTxnxtxnT()xnnˆ()axttˆ()aXjcssc2s()jXecT22cT1ˆ()()aasmXjXjjmT12()()jakkXeXjTT22ssTTTT理想采样信号理想采样信号与时域离散信号的频谱关系SCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.122020/1/164.2.3带通信号的采样带通信号中心频率最低频率最高频率带宽特点:中心频率较高,带宽较窄,采样频率?lff()aXjf0hf0flfhf0flfhfBhlfff0BffSCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.132020/1/16如最高频率是带宽的整数倍,采样频率如下,不会造成频谱混叠:理想采样信号的频谱为原模拟信号的频谱以采样频率为周期进行周期延拓。例K=4hf2sBFfhBfKfsFfˆ()aXjf0hflflfhff()aXjf0hflflfhfBfSCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.142020/1/16带通信号的恢复理想采样信号通过一理想带通滤波器()0lhpTfffHjf其它sFfˆ()aXjf0hflflfhfSCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.152020/1/16降低最低频率到,带宽展宽为小于的最大整数例解:hBfKf/hBrfflfBhlfffhBfrfrr4kHz32kHz83hBhBffrff不为整数2=2kHz24kHzBhsBrffrFf问题:如果SCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.162020/1/16sFfˆ()aXjf0hflfhff()aXjf0hflflfhflfSCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.172020/1/164.2.5将数字信号转换为模拟信号预滤波平滑滤波数模转换DAC数字信号处理量化编码采样()xt()xn[]xn[]yn()yt()yn理想恢复(略)SCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.182020/1/164.3对数字信号处理部分的设计考虑见书上例题P103SCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.192020/1/164.4线性模拟系统的数字模拟用数字系统模拟线性模拟系统的外部特性A/D()xn()yn()axt()jHeD/A()ayt()axt()Hj()aytSCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.202020/1/16数字系统单位脉冲响应与模拟系统单位冲激响应之间的关系为结论:采样频率满足采样定理D/A为理想恢复信号与系统均为带限信号与带限系统可以用数字系统模拟线性模拟系统,进行分析与设计()()()aatnThnhthnT()axt()Hj()aytA/D()xn()yn()axt()jHeD/A()aytSCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.212020/1/164.5模拟信号的频谱分析问题:(1)模拟信号的频谱分析不适于用计算机实现,一般是将模拟信号转换为时域离散信号,用DFT(FFT)进行频谱分析。需解决哪些问题?时域抽样时域截断(2)用DFT(FFT)进行模拟信号的频谱分析,误差如何?SCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.222020/1/16用DFT进行谱分析的基本思路()()aaxtXj连续信号连续频谱DTFT()()()jaxnxnTXe采样连续周期频谱DFT()()xnXk在频率区间[0,2]上N电等间隔采样SCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.232020/1/16模拟信号的傅里叶变换设连续信号的持续时间为Tp,最高频率为fc,令t=nT,N=Tp/T,dt=T,Fs=1/T2fc,得Xa(jf)的零阶近似:频域等间隔采样,在区间采样N点,采用间隔为F(频谱分辨率)()axt2()FT[()]()jftaaaXjfxtxtedt[0,]sF11spFFNNTT120()()NjkFnTanXjkFTxnTe210()NjknNanTxnTe120()()NjfnTaanXjfTxnTeSCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.242020/1/16令,则类似地210()()NjknNanXjkFTxnTe()(),()()aXkXjkFxnxnT210()()NjknNnXkTxneDFT01[()]TxnkN1IDFT[(01)]()XknTxnNSCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.252020/1/16fc-fcffktttxa(t)x(n)xa(nT)fs/2-fs/20N-1TpNTTFXa(jf)X(jf)Xa(k)2()FT[()]()jftaaaXjfxtxtedt120()()NjfnTanXjfTxnTe210()()DFT[()]NjknNanXkTxneTxnSCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.262020/1/16重要参数的选择minminminminminminmax//1/spsNFFTNFNTFminmax221/scscFfFfT时域抽样:1/pFT频域抽样(频率分辨率):cf信号最高频谱:SCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.272020/1/16解:最小记录时间最大采样间隔最少采样点数谱分析范围如何提高谱分辨率?例4.5.1对模拟信号进行频谱分析,要求谱分辨率信号最高频率计算:最小记录时间、最大采样间隔、最少采样点数及谱分析范围;如何提高谱分辨率?minpT10FHz2.5cfkHzminmin/22500/10500sNFFmaxTminNmin1/1/100.1pTFsmaxmin1/sTF31/(22500)0.210smin1/pFTFminpTminNmin0~sFSCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.282020/1/16DFT与连续时间信号傅氏变换间的相对数值关系DFT[()]()(),01kTxnXkXfkN21,,22skkpskkkksfkfkFkkNNTffkfTN1IDFT[()](),01XxnNTnkSCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.292020/1/164.5.2用DFT(FFT)对模拟信号进行谱分析的误差误差来源1.频谱混叠时域采样时,应满足采样定理的条件。为避免频谱混叠,可以作预滤波。()aXj0cc()Pj0ssˆ()aXjssc2s预滤波平滑滤波数模转换DAC数字信号处理量化编码采样()xt()xn[]xn[]yn()yt()ynSCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.302020/1/162.截断效应如x(n)为无限长序列,必须截短成有限长序列,因截断引起的误差现象称为截断效应矩形窗傅里叶变换()()()NNxnxnRn()NRn()1()DTFT[()]()()21()()2jjjNNfNjjfNXexnXeReXeRed1()2()DTFT[()]sin(/2)()DTFT[()]()sin(/2)jNjjjfNNfNXexnNReRneReSCHOOLOFPHYSICSANDTECHNOLOGYN.N.U.312020/1/16主瓣宽度:有许多旁瓣,截断后信号的频谱与原信号的频谱不同1()2sin(/2)()DTFT[()]()sin(/2)NjjjfNNfNNReRneResin(/2)