低频电磁信号的频率细化技术研究.docVIP

汽轮机震动信号的频率细化技术研究 摘要：在低频电磁信号中，为了得到某一窄带内频谱的精细结构，提高频谱的频率分辨率，满足良好的实时性要求，本文提出了基于快速傅立叶变换算法的频率细化技术。通过复调制频率细化，将某一局部频谱放大，便于对频谱进行快速、）方法 1．1 复调制（）细化方法的基本原理 复调制方法(以下简称)能以指定的、足够高的采样频率分析频率轴上任一窄带内信号的频谱结构。在序列变换点数相同的情况下, 能获得更高的频率分辨率;或者,在相同的频率分辨率下, 比基带需要更少的傅立叶变换点数。 图1 原理框图 Fig.1 The principle block diagram of 设模拟信号为，经抗混滤波、转换后得到采样时间序列 ，为采样频率，为细化分析频带，为需要细化的频带中心频率, 为细化倍数, 为分析的点数, 具体算法过程可归纳为以下5个步骤： （1）复调制移频 所谓复调制移频就是将频域坐标向左移或向右移，使得被观察频段的起点为频域坐标的零频位置。这里对离散信号用进行复调制，把需要细化的频带中心频率移至频率轴原点，得到?: （1） （2）数字低通滤波 为保证重新采样后不发生频谱混叠,必须进行抗混叠滤波, 滤出所需分析频段信号。设频率细化倍数为，则低通滤波器的截止频率。 （3）重新采样 信号经过移频和低通滤波后,分析信号频带变窄,因而可以得到较低的采样频率进行重采样, 比原采样频率降低了,即对原采样点每隔点在抽样一次。 （4）复处理 对重新采样后的点复序列进行复处理，得到条谱线，其频率分辨率，分辨率提高了倍。 (5) 频率调整 将上述谱线移至实际频率处即得到细化后的频谱。设第(4)步得到的频谱为,细化后的频谱为。 （2） 1．2 流程图 根据复调制细化的原理，我们得到复调制细化的流程图： 图2 流程图 Fig.2 Flow chart 在流程图的基础上对复调制细化进行编程，得到其实现的程序。然后通过matlab对信号进行仿真，同时，对仿真结果进行分析。 2 复调制频率细化的matlab实现 2．1 模拟低频电磁信号 采用一个模拟信号，针对频率细化方法进行仿真模拟。设信号由频率分别为和,振幅分别为和的两个正弦信号组成。 （3） 采样频率为，采样点数为1024，放大倍数为10。对于信号进行采样，则频谱的分辨率为，而.所以要把这两个频率分辨开来就需进行频率细化，将频谱在的范围内细化倍。 2．2 对模拟信号进行复调制细化仿真 （4） 在matlab中，对信号进行仿真分析，结果如下所示： 2．2．1 信号的时域图和频谱图 图3 信号的时域波形图 图4 信号的频谱图 Fig.3 Time domain wavefom diagram of the signal Fig.4 spectrum diagram of the signal 从图4中可以看出两个频率的频谱是重叠在一起的，无法区分开。 2．2．2 对信号的频谱进行移频，移到零频位置对其重新抽样 图5 移到零频位置的频谱图 图6 重新采样后的频谱图 Fig.5 Move to the location of fig.6 Re-sampled spectrum diagram zero frequency of the spectrum diagram 从图6中我们可以看到，经重新采样后，对其进行快速傅里叶变换，在零频位置已经将两个混叠在一起的两个频率的频谱区分开来，达到细化的目的。 2．2．3 将频谱恢复到原始位置： 图7 恢复原位置的细化谱 Fig.7 Restore the original location zoom spectrum 图3、4、5、6、7反映了复调制频率细化的细化过程，通过仿真过程，我们看到了频率细化的有效、可行。将两个混叠的频谱区分开来，而且过程非常的清晰明了。 2．3 线性调频的仿真结果 对同一信号进行线性调频变换，其细化结果如图所示： 图8 信号的频谱图 图9 细化频谱图 Fig.8 spectrum diagram of the signal fig.9 The zoom spectrum diagram of 如图8、9所示，经过细化后，这两个频率的频谱依然重叠在一起。