短记录加汉宁窗的频谱校正.pdfVIP

维普资讯 振 动 与 冲 击 第27卷第4期 JOURNALOFVIBRA~ON AND SHOCK 短 记 录 加 汉 宁 窗 的 频 谱 校 正 陈奎孚 ， 王建立 ， 张森文 (1．中国农业大学 理学院74#，北京 100083；2．邢台职业技术学院，河北邢台 054035； 3．暨南大学 应用力学研究所，广州 510632) 摘 要 ：基于解析单频模型建立的频谱校正方法难以应用于短记录的频谱校正，因为这时正频率附近的频谱显著 偏离解析单频模型，背后的原因是实数信号中负频率分量的泄漏干扰。针对短记录加汉宁窗，给出了一种显式频谱校正 方法，它利用局部谱峰附近的三条谱线。采用仿真手段对提出的方法进行了考核。不同仿真样本所含的信号周期数 (CiR)从 0．05变化到 5，步长为 0．O1；相位从 0。变化到 179。，步长为 1。。结果表明：1)当CiR1时，最大频率误差小于 10 。△甜(△甜是快速傅里叶变换的频率分辨率)，而幅值相对误差和相位误差的上限分别不超过 10 。和 10 度；2)当 CiR1时，误差总体趋势随频率下降而增加，但即使对CiR：0．05，频率误差也不超过4×10 △甜；3)精度最高的条件仍 然是整周期采样。 关键词：短记录；频谱 ；快速傅里叶变换(FFT)；汉宁窗 中图分类号：TN911．6：0329 文献标识码：A 快速傅里叶变换(FastFourierTrans~rm：FVr)是处 这里 。，A。，。和B。分别为信号的圆频率、幅度、相位 理周期信号的经典方法，但是 由于栅栏效应，直接从 和直流偏移等 四个参数。假定模型(1)在 [0，T]上截 FFr谱线读 出的信号参数的误差有时可能很大。例 出，而[0，T]之外为零 。 如：对于矩形窗，在极端情形下，幅值误差可达 36％，相 下文需要用到记录中的周期数 (CyclesinRecord： 位误差可达90。 。为了降低误差，目前已经发展了 CiR)这个概念。CiR在数学上等于 。／(21T)= 。／ 形形色色的校正方法，如：幅值 比值法 一、相位差 (△ )，这里△ =21T／是FFT的频率分辨率。在 FFT 法 、双窗法 等。对解析单频 (AnMytiealSingle 之前要加如下所示的汉宁窗 Tone：AST)信号，或频率成分非密集情形 ，这些校正方 法都具有优异的性能。但是对于低频成分，其校正精 W((￡￡))=：f{ 2 0…≤￡一≤ ((2)) 度随信号频率的变化而变化 J，特别是当信号频率非 0 t或者 tT 常低或者接近奈奎斯特频率时，这些校正方法的精度 (￡)加上述汉宁窗之后的傅里叶变换为 很不理想 。对4项余弦窗的频谱插值，钱吴等则明 ()=Aoexp(～jo) ( + o)／2+ 确指出采集的样本至少要包含二十几个周期 ¨。 Aoexp( o)W( 一 o)／2+BoW() (3) 频率的高低与记录的长短有关。当采用 FFYr对短 这里W()是 (t)的频谱，它可表示为 记录作谱分析时，除了FFYr栅栏效应会引起误差外，信 ()=D()／P() (4) 号中的负频率分量泄漏也使真实频率附近的谱值不满 其 中分子和分母分别为 足AST模型，从而导致相应的校正方法失效。本文将 )：吾sin()exp(一譬)， 建立短记录的加汉宁窗的频谱模型，利用FFT频域采 样的特殊性，构造一个无须迭代 的显式校正方法。通