《第二章_离散信号频谱的窗谱校正方法99841》由会员分享，可在线阅读，更多相关《第二章_离散信号频谱的窗谱校正方法99841（48页珍藏版）》请在金锄头文库上搜索。
1、第二章_离散信号频谱的窗谱校正方法99841第二章_离散信号频谱的窗谱校正方法.txt不要放弃自己！ -（妈妈曾经这样对我说，转身出门的一刹那，我泪流满面，却不想让任何人看见！） 看到这一句小编也心有感触，想起当初离家前往几千里外的地方的时候，妈妈也说过类似的话，但是身为男儿，必须创出一片天，才能报答父母的养育之恩！第二章_离散信号频谱的窗谱校正方法.txt52每个人都一条抛物线，天赋决定其开口，而最高点则需后天的努力。没有秋日落叶的飘零，何来新春绿芽的饿明丽？只有懂得失去，才会重新拥有。华中理工大学博士学位论文 第二章离散信号频谱的窗谱校正方法 -基本理论 2.1 引言利用DFT可以对离散信
2、号进行频谱分析，但是计算工作量相当大，因此，在快速算法没用发明之前，DFT并没用多大的实际意义。直到1965年，Cooley-Tukey在计算数学杂志上首先提出FFT算法之后， DFT才得到广泛的应用。这一快速算法的出现对数字信号分析领域的发展起到了极大的推动作用。从此以后，它作为频谱分析的基础得到了广泛的应用75,76,77,78。由于计算机只能对信号的有限多个样本进行计算，信号的 FFT谱分析也只能在时域信号的有限区间内进行，这就不可避免地存在由于时域截断(加矩形窗)而产生泄漏61，使谱峰值减小，精度降低，求得的信号相位更是面目全非。在数字信号处理中，由 DFT或 FFT得到的幅值谱是离散
3、谱，是信号与窗函数频谱卷积后，按频率分辨率 fs = fs / N ( fs为信号采样频率，N为分析信号样本长度)等间隔频域抽样的结果（如图 21所示）78。 A幅值f图2-1 频谱抽样的离散谱线如果周期信号的频率正好表2-1 离散频谱幅值、相位和频率误差表落在某一谱线上，经FFT后得矩形窗 Hanning窗 Hamming窗幅值误差（） 0-36.4 0-15.3 0-18.3 相位误差（0） 90 90 90频率误差（Hz） 05fs. 05fs. 05fs. 到的频率、幅值和相位是准确的。在一般情况下，信号频率落于两条相邻谱线之间，由于谱线不在主瓣中心，由峰值谱线反映的频率和幅值都不准1
4、0华中理工大学博士学位论文 确，相位误差更大。从理论上分析，加矩形窗时，最大误差可达 36 4.%，即使加其他窗时，也不能完全消除这一影响，在加Hanning窗时，只进行幅值恢复时的最大幅值误差仍高达15 3.%，相位误差将更大，表2-1是离散频谱只进行幅值恢复，不进行其他处理时幅值、相位和频率误差141。 2.2 单频谐波的频谱分析误差产生原因无限长信号 x()t（如振动信号、噪声信号等）的频谱分析所采用的方法为对信号进行截取，然后再对截取得到的有限长度信号进行频谱分析。窗函数 wt()的作用就相当于对无限长的信号开一窗口，从窗口中取出一段数据，从而完成信号的截取。窗函数都是选择实偶函数，并
5、在时域上将窗函数的中心放于被分析的那段信号的中心。加窗信号的傅氏变换为： Fx. () ()() 2 dt.(2.1) () twtxtwteTTjft=.() .其中，wtT()由对称窗wt()在时间上平移T/2得到，即 T() t.T/) wt=w( 2 .(2.2) 设wt()的傅氏变换（如图2-2a） Fwt =W( f( ) ) .(2.3) 根据傅氏变换的奇偶性质，当wt()是实偶函数时，Wf()此时也为实偶函数。又由傅氏变换的时移特性可知（如图2-2b）， Fw t =Wfe () jfT T( ) .(2.4) 设有一周期信号 x()t=ACos(2.f0 .t+.),则其傅氏
6、变换结果为（如图2-2c）： A .j. Aj.() 0 .f0 ) .(2.5) Xf =.e .(f +f) +.e .(f 22 根据卷积定理，加窗后的谐波信号 x()twt.T()的傅氏变换可表示为(图2-2d)： Xf () =Fxt w t T( ) =Fxt () Fw t ( ) T () . T （这里“*”表示卷积） AA.jj . . . j fT = .e .( + ) +.e .( .Wfe ff ff ) () 0022 A. + )+. AjTf f 0).jTf f . . (0 (=.Wf fe + ) +. ( . )( Wf fe . .(2.6) 002
7、2 由此可知，在加窗信号的傅氏分析中，当 ff0时，将存在泄漏情况。此时的幅值及相位分别为： AY=.Wf f ( . 0 ) .(2.7) 2 =. . Tf( .f +.(2.8) 0) 11华中理工大学博士学位论文 对窗长度T = N /fs 作归一化处理，则T = 1，且令f = f . f0代入上面两式可得： A (. ) =.f +. .(2.10) Y = 2 .Wf .(2.9) W ( f )w () Ttt f-T/2 T/2 1.1/0.2/2/1/时域波形窗谱模函数 (a)对称矩形窗函数的时域波形和频谱模函数 wtT ()1 t 0T T .2/.1/0 2/1/W()
8、fTf时域波形窗谱模函数 (b)实际窗函数的时域波形和频谱模函数 T0 Ax(t) t X ( f )A/2A/2 f -f 0 f00 时域波形傅里叶变换模函数 (c)单频率谐波的时域波形和频谱模函数 A xtT( ) 0 T t Y n yKyK.1 k k+2k+1k-1k-2 f00 时域波形离散频谱模函数 (d)单频率谐波离散频谱模函数图2-2 单频率谐波离散频谱的误差产生原因12华中理工大学博士学位论文 显然，当 f=f0时，Y= A ，=.不存在泄漏情况，得到的幅值、相位和频率都是2准确无误差的。在大多数情况下，当 ff0时。由加窗信号的傅氏分析得到的频率 f、幅值Y和相位并不是
9、真实值，且有旁瓣产生，这就是所谓的离散频谱的栅栏效应、梳状效应、能量泄漏和假频等（如图 2-2d所示）。当信号真实频率位于两个相邻离散谱线中间时，即 fK.1 =f0 .fs/2 ， fK =f0 +fs/2 （这里 fs为频率分辩率）时，求得的信号幅值、相位和频率的误差最大。 2.3 离散频谱信号的窗谱校正方法假设加窗信号的频谱主瓣中心为 f（即为信号的真实频率），信号幅值为 A0；加窗信号FFT结果的频谱中，最高的频谱频率(0) 为 f，高度为Y；次高谱线频率为 f2，高度为 Y2。显然， f和 f相差仅为一个频率分辨率(1) ，对此归(1) 一化后，即有 f1 =f2 1。当 f1 f2
10、时取“+(1) ”号(2) ；当 f1 f2时取“”号。加窗信号的频谱由窗频谱对信号真实频谱调制而得到，因而可利用窗函数的频谱图形对傅氏分析结果进行校正，以求出真实的频率、幅值和相位6567。 (1) 频率校正 0 WWxyf+1 f f+1f0 yyxyf+1 f f+1ff0 A图2-3 窗函数的频谱函数图2-4 频率和幅值的校正频率校正即求出信号幅值谱图上窗谱主瓣中心所处的横坐标 f0。设窗函数的频谱函数为 Wfff，其窗谱函数为()， W()对称于 Y轴（如图 2-3所示）。对于任一 Wf()，相应的信号离散频谱幅值为 y；对于处于Y轴右侧的 f+1谱线，其窗谱函数为Wf1，相应的信号离散频谱为 yf+1(1、 yf和 yf+1(+)(+)幅值(f) ；由Wf)、Wf可求出 f的值，即求出了频率修正量 ffff=.0。由于W()的函数表达式为已知，可构成一函数： Wf) yf(mFf = () = .
