数字信号处理复习资料重点.docx
《数字信号处理复习资料重点.docx》由会员分享,可在线阅读,更多相关《数字信号处理复习资料重点.docx(13页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、1 .数字信号处理概述.离散时间信号2 .线性时不变系统.信号与系统的相互作用3 .离散傅里叶变换.快速傅里叶变换4 .数字滤波器概述.有限冲激响应滤波器5 .无限冲激响应滤波器收系统的频率响应是回波信号的傅里叶变换先求共粗,再乘上一个相移因子。时域:h(n) = Ks(no-n),即最佳接收系统的单位冲激响应式回波信号在时间轴上先折叠,再平移。频域分析最佳接收系统的幅频响应等于回波信号的幅频特性,在没有信号分量的频率上,噪声完全被抑制。 信号通过系统之后,相位都被抵消了,这样所有的不同频率的信号在0时刻的相位都是相同的,因而 在这个时刻点上的输出具有最大值。时域分析s0(n)=rss(n-n
2、o),系统输出信号实际上就是回波信号的自相关函数。最佳接收系统又称为相关接收系 统。一般情况下,噪声和信号是没有相似性的,因此,通过相关接收系统之后,很自然就尽可能抑制 了噪声。信号和自身是最相似的,就尽可能的放大了信号。在高斯白喔声背景下,错误概率最小准则下的最佳接收系统等效于输出信噪比最大准则下的最佳接收系 统。卷积表征的是一般意义上的信号与系统的相互作用,即任意的信号通过LTI系统(由单位冲激响应来表示), 都可以表示为卷积的关系。而相关表示的则是信号通过最佳接收系统这一特殊的系统。对确定性的信号,特别是非周期的确定性信号,常用能量谱来描述。面对随机信号,由于理论上持续 时间无限长,不满
3、足绝对可积与能量可积的条件,因此不存在傅里叶变换,所以常用功率谱来描述。周 期性的信号,也同样是不满足傅里叶变换的条件,也常用功率谱来描述。第五章离散傅里叶变换DFT离散傅里叶变换/2%X(k)=漆5离散傅里叶逆变换X(n)=初版 X(k)2rr/N)nkX(k)通常为复数,可用实部虚部的方式,也可以表示为幅度相位的方式。频率,幅度,初始相位这三个参数完全决定了一个复正弦信号。对离散时间信号来说,信号的频谱表示为信号的DTFT。比较DTFT和DFT的表达式,可以很容易的发 现:X(k) = X3 H 3=(27r/N)k可以发现,x(k)是X(*)的采样值,即DFT是DTFT的采样值,采样间隔
4、为2H/N.在3 = 2TTf, f表示信号的模拟频率(即线频率),fs为系统采样频率,3为对应的数字频率。现实中 的信号多是模拟信号,离散信号多是从模拟信号采样得到的。模拟频率的物理意义非常清晰(表示物体震 动快慢的物理量)。DFT结果幅度和相位的理解在X(k)=一片中,令Bn=x(n) ,7张。n,Bn为单位圆上的旋转向量, 旋转的速度为= 6)o- 2nk0/N ,则X(k)= Xn=o Bn 随着叩的增加,向量相加的结果的幅度会逐步减小,相位则会逐步增加。同时在DFT的结果中包含了 原始信号的初相信息。DFT具有:隐含的周期性和循环卷积的特性,与DTFT和Z变换不同。X(k)的幅频特性
5、和相频特性都隐含着是周期性了。X(k)和x(n)都是周期函数,周期分别为N和no 频域的采样同样等效于时域的周期延拓,并且延拓的周期为N。只要对x(n)进行DFT ,客观上就已经 将x(n)进行周期延拓了(N可以很大,主观上我们还没有意识至少。隐含着周期性条件下的卷积称为循环卷积,又称为圆周卷积。数学公式表示为:y(n) = h(n)0x(n) = 霹- 7m)Q , 0 n N-1 ,其中(5 - m)N表示的是周期性的x(n)移位m点之后,再取OW n 0对于任意给定的k=k0,都可以看作是信号x(m)通过系统%(m)在m=N时刻的输出。因此,x(n)的DFT结 果X(k),可以看作是信号
6、x(n)通过一组滤波器/1-(九)在n=N时刻的输出。令hk(n) = exp(j(2TT/N)kn) , n=l,2,3,4N ,由hk(n)的频率响应可知,当k取值不同时,hk(n)的幅频特性 形状都完全一样,都具有sine函数的形状,不同的只是频率中心发生变化,第ko个滤波器幻。5)的频率 中心位于3=2T!ko/N。(注意观察滤波器组的幅频特性)在信号处理中,频域分析基本上就等同于傅里叶分析,但就傅里叶分析来说,又有多种具体的方法,比 如对离散信号有离散时间傅里叶变换(DTFT),离散傅里叶变换(DFT 对连续信号有傅里叶级数(FS 1 傅里叶变换(FT I1 .傅里叶级数FS在满足一
7、定条件下,任何一个周期信号都可以分解为正弦信号的叠加。周期连续 非周期离 散.傅里叶变换FT即傅里叶积分。非周期连续 非周期连续.离散时间傅里叶变换DTFT对非周期的连续信号进行时域采样,得到非周期的离散信号,非周期离散信号的傅里叶分析称为 离散时间傅里叶变换。DTFT是周期的。非周期离散。周期连续.离散傅里叶变换DFT对周期离散信号的傅里叶分析称为离散傅里叶变换。周期离散 周期离散注:实际上,按照传统命名规则,离散傅里叶变换叫离散傅里叶级数或离散时间傅里叶级数更为妥当。 而对非周期离散信号的傅里叶分析才应该叫离散傅里叶变换。但本书中将离散傅里叶级数的分析隐含在 离散傅里叶变换中。在FS、FT
8、、DTFT、DFT中,时域的周期性对应着频域的离散性,时域的离散性对应着频域的周期性。 最重要的是DFT。因为前面三种都需要假定信号的时域或者频域是无限长的。计算机等数字设备只能处 理数字信号,也即要求无论是时域还是频域,都要是离散的。DFT 与 DTFTDFT可以看作是DTFT在频域上的等间隔采样的结果。DFT定义为刈叫间隔为g=2TI/N的等间隔采样, 因而不存在混叠。由于DFT不存在时域上的混叠,因此由X(k)可以不失真地恢复出Xw),数学表达式为x(蚂= X二i X,(爷5=中 n=0n=0n=0k=0.泄漏造成泄漏的原因主要是DFT所等效的N个滤波器组中,没有一个滤波器的中心频率和信
9、号的频 率相等,导致单频复正弦信号在滤波器主瓣内幅度有所抑制,信号能量有一部分散落到滤波器的 副瓣内。频谱泄漏是经常的,危害也很大,会导致信号的误判。分为主瓣泄漏和副瓣泄漏。实际的信号都被噪声污染了,副瓣的值很难参考。主瓣泄漏是因为复正弦信号的频率同时处于两个滤波器的主瓣范围内副瓣泄漏则是由于滤波 器存在副瓣。为了抑制频谱泄漏,不仅要降低滤波器的主瓣宽度,同时还希望降低副瓣的宽度。1 .加窗数据截短的过程称为加窗,数据截短会带来频谱的泄漏。直接对数据进行截短,等效于将信号 与脉冲信号相乘。时域上的突然变化,意味着频域上存在很高的频率分量。所以只要抑制了这些高 频分量,就能降低副瓣电平,从而有效
10、抑制副瓣的泄漏。窗函数起始处和结束处越平坦,对副瓣 的抑制应该会越好。频率分辨率通常与主瓣宽度的变化趋势一致,即主瓣越宽,频率分辨率越低。2 .分辨率主瓣半功率点之间的宽度。在DFT中,通常用2H/N来表示脉冲信号的数字频率分辨率,用1/T表 示模拟频率分辨率。脉冲信号的频率分辨率为脉冲宽度的倒数。为了提高频率分辨率,要尽量增加 信号的有效长度。当信号持续时间为无限长的时候,其频率分辨率趋于0。时间分辨率越好,则 对应的频率分辨率越差,这就是信号处理中的不确定性原理。3 .补零DFT可以看作是对离散信号连续频谱在频域的采样,只要频域采样点数M大于等于时域长度N ,则 可以减小栅栏效应,同时也满
11、足频域采样不混叠的要求。在频域上增加采样点数等效于在时域上补 零。从信息的角度,补零并不增加新的信息,自然也就谈不上改善分辨率。区分物理分辨朝口计算分辨率。矩形窗的物理分辨率为0.89F/N ,计算分辨率为2WN。两者非常 接近,但不完全相等。4 .增益信号处理系统的一个目标是尽量放大信号而抑制噪声,这个目标通常是用增益来表征的。增益G =安,分子分母分别表示输出输入信噪比。SNRinDFT能够很好地抑制噪声,其增益与DFT的点数N有关,G = lOloglON,对信号所在的窄带滤波器 来说,信号能完全通过,而噪声却只能通过1/N。第六章快速傅里叶变换FFT在工程实践中,DFT面临的主要问题是
12、运算效率低,很难用于实时处理。X(k)=一度*叱;%5)卬炉,W炉为自旋转因子DFT的运算量与M成正比,显然N越小越有利,因此就很自然地希望能将大点数的DFT分解为多个 小点数的DFT来计算,分解过程中充分利用旋转因子的对称性,周期性,可约性和特殊值。FFT算法主要分为两大类:按时间抽取方法(DIT)和按频率抽取方法(DIF)根据抽取长度,FFT又 分为基2基-4等算法,在一个DFT运算中还可以同时采用不同基的混合基算法,也称为分裂基算法。基-2 DIT FFT算法按时间抽取算法是指在时间上将信号长度逐步减小的算法,基-2算法则是指以2为基底降低信号长度的 过程。X(k) = X0(k) +
13、W1fxi(k) , k= 0,1,., N/2 -1X(k+N/2) = X0(k)-赚Xi (A) , k= 0,1,,N/2 -1一次分解后,运算量近似为直接计算的1/2.一点的DFT就是信号本身,两个一点的DFT或者一个两点的 DFT,正好构成一个基本的蝶形运算。因此,通常将大点数DFT按级分解,直到DFT的点数变为2.(即 将N点的DFT分解为N/2个2点的DFT ) 级的概念:M= log同址运算:原位运算,指的是在运算过程中,输入和输出占用相同的存储地址,从而可以节约存储空间, 进而节约存储空间。存放数据的存储单元始终为N个。输入信号不是按自然)11页序存储的看起来好像杂乱无章,
14、但实际上是按照倒位序排列的(偶数在上,奇 数在下I倒位序的实现:当n=(时,不必调换。否则,必须将x(n)与)调换。在DIF FFT (按频率抽取)中,输入x(n堤正常JII踮的输出X(k)是倒位序的。分裂基FFT算法,先是如基-2算法一般将x(n)分解为偶数点部分和奇数点部分。对偶数点部分按基-2算 法来分解,对奇数点部分是按基-4算法来分解。分裂基FFT算法同样也可以按时间抽取或者按频率抽取。 在满足基-4算法运算条件的前提下,基-4算法效率比基-2算法更高。在已知的N=2M的各种算法中, 分裂基算法的运算效率是最高的。基-2 ,基-4和分裂基算法都是常用的FFT算法。基-4算法虽然运算效
15、率高于基-2算法,但其所要求的 N=4M的条件比基-2算法要严苛。在不满足的情况下,补零的个数过多,反过来B斜氐了运算效率。分裂 基算法的运算效率也高于基-2算法,但编程的复杂度也比较高。因此,在实际中最常用的还是基-2算 法。旋转因子的生成在实际的FFT运算中,通过计算第=cos (等)一由皿2河0)来生成旋转因子。一种有效的方法是查表法。IDFT的快速计算IDFT : x =92屋(X”卬胪实信号的FFT(提高效率).两个实信号的FFT1 .一个2N点的实信号的FFT将实信号的FFT转换为复信号的FFT。在数字信号处理的工程实现中,往往用DSP芯片来实现具体的算法对每一种型号的DSP ,硬
16、件开发商 都会给出优化后的FFT算法源代码一般情况下直接调用就可以。但是,在某个具体的应用中,可能会遇 到大点数FFT的问题,因为芯片内存不足而导致FFT需要分解的问题。两类特殊的频率分析算法有时只需要计算某个或者某几个频点的X(k)值,而不需要计算所有的X(k)的值。也许有时候只需要分析 信号在某一个很小频段的频段值,而且要对这很小频段范围的频谱值进行比较精细的分析,此时也无 法直接用FFT算法。1. Goertzel 算法利用递推的方法来计算频谱值,在只需要计算很少几个频点的频谱值得情况下,Goertzel的效率高于 FFT算法。2. ChirpZ 算法利用卷积的方式来实现频谱的分析,期间
17、也用到了 FFT。当分析的信号频带很窄的时候,并且要求 的分辨率非常高的情况下,用ChirpZ算法比补零FFT效率更高。实际上,ChirpZ算法不仅是分析频域, 它还能分析Z平面上任意一段弧线的响应。第七章数字滤波器概述数字滤波器是一个按预定的有限精度算法实现的、将输入的数字信号转换为所需要的输出数字信号的线性 时不变系统LTL滤波器是是一种对信号有处理作用的器件或电路,其主要作用是让有用信号尽可能无衰减的通过,对无 用信号尽可能大地衰减。由电阻,电容和电感等元件组成的模拟滤波器广泛应用于各种信号处理系统 中,优点是处理速度快,处理带宽大,无需ADC及DAC器件;其缺点是稳定度及精度比较差,可
18、 重复性不强及抗干扰能力弱。数字滤波主要有两种用途:信号的分离和失真信号的恢复没有特别指明的情况下,数字滤波器与线性时不变系统是等价的,因此,描述LTI系统的工具均可用 于描述数字滤波器,这些描述工具包括:差分方程,单位冲激响应,传递函数(H(z)=Y(z)/X),频率响 应(幅频响应和相频响应)等。数字滤波器的分类.从单位冲激响应的角度分类有限冲激响应FIR和无限冲激响应HR.从频率响应的角度分类低通滤波器,高通滤波器、带通滤波器、带阻滤波器、全通滤波器(用于改善信号的相频响应)实际数字滤波器的性能参数截止频率:滤波器增益为0.707时的频率,即-3dB所对应的频率。通带边沿频率:滤波器增益
19、为I-6时的频率。截止频率通常要大于通带边沿频率。增益和衰减都表示滤波器对某个频率分量幅度上的影响,增益通常指对信号的放大,衰减通常对 信号的抑制。假定通带偏差距=0.1 ,阻带偏差演=0.05 ,如果理想增益为1 ,则通带边沿增益为1- 0.1=0.9 或 20log100.9=-0.915dB ,阻带边沿增益为 0.05 或 20logi00.05=-26.02dBo 用 dB 做单位时,衰减是 增益的负值。这样,也可以说是通带边沿的衰减为0.915dB ,阻带边沿的衰减为26.02dB。表征实际的数字滤波器的性能参数时,更重要的却是单位阶跃响应,即滤波器输入为单位阶跃信号时 滤波器的输出
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字信号 处理 复习资料 重点
限制150内