舒曼谐振信号是在地球和电离层谐振腔内传播的电磁波.doc

舒曼谐振信号处理 舒曼谐振信号是在地球和电离层谐振腔内传播的电磁波，主要由全球雷电激励产生。舒曼谐振对于全球雷电分布、空间天气以及电离层底部探测等研究都有重要的意义。论文针对舒曼谐振观测信号微弱，且易受观测环境噪声和干扰的影响等问题，进行了如下研究: 论文在传统方法的基础上，首先采用主分量分析（PCA）方法，对时域信号进行预处理，以提高信号的信噪比。在此基础上采用了周期图谱法去除部分噪声和干扰，进而采用包络滤波法代替文献中洛仑兹拟合来平滑信号的幅度谱，获得了较为理想的平滑曲线。 论文采用永胜观测站的观测数据，分析了舒曼谐振各阶信号的幅度和中心频率的日变化规律。结果表明，经过本文方法处理后提高了舒曼谐振信号的信噪比，比单纯的周期谱图法更能反映实际信号的变化规律。 一．舒曼谐振数据的获取 数据于2010年--2011年收集于云南永胜观测站（经度100.76，纬度26.71），该观测站检测三个方向的数据组成成分（南北向磁场，东西向磁场和垂直向电场。在本次研究中，仅研究南北向和东西向的磁场数据，该数据主要受来自东南亚雷暴源源和非洲雷暴源的影响，数据采样频率是100Hz，处理时长为10s，即1000点。 舒曼谐振数据的时域波形 舒曼谐振原始信号的时域波形，如图2-1所示。舒曼谐振的磁感应强度在1pT以下，非常微弱。信号的时域分辨率为0.01s，从图2-1中可以看出，舒曼谐振信号的周期性很差，幅度在0.015pT左右。 图 2-1原始信号时域波形 图 2-2原始信号时域波形 舒曼谐振原始信号的频域波形，如图2-2所示。对舒曼谐振10s的数据（1000点）作FFT，取幅度谱，可以看出，舒曼谐振数据包含丰富的频率分量，噪声和干扰比较强，信号的谱峰不明显，无法直接进行参数的提取。因此，需要进行后续的处理。 二．舒曼谐振信号的处理 本章分别介绍了对原始舒曼谐振信号进行的处理，包括：一、削弱噪声和干扰的处理，即周期图谱法、主分量分析、中值滤波；二、用来拟合幅度谱的方法：洛仑兹拟合、包络滤波。 2.1 周期图谱法 本处理中，将30组数据的幅度谱叠加，每组数据时长10s(1000点)，即得到叠加后的时长为5min的数据的幅度谱。 经过周期谱图法处理后的舒曼谐振数据的幅度谱如图2-3所示，其中蓝色曲线代表原始SR信号的幅度谱，红色曲线代表经周期谱图法SR信号的幅度谱，从图中可以看出，经周期谱图法处理后，很多强干扰和噪声明显降低，前四阶谐振峰可以比较明显地看出。 图2-3 周期图谱法处理后频谱 2.2 主分量分析 在PCA处理过程中，构造的Hankel矩阵大小为901×100，对其求协方差矩阵，然后进行特征值分解，求得的协方差矩阵的100个特征值的分布。根据曲线发生突变的点的位置大约在前30个点处，取,即取前30个较大的特征值来重构数据阵。 周期谱图法处理后的幅度谱与先经过周期谱图法再进行PCA处理后的舒曼谐振信号幅度谱的对比，如图2-4，可以很明显地看出，相较于仅经过周期谱图法处理的信号幅度谱，经过PCA处理以后，各阶谱峰更加显著，信噪比也得到了较大的改善，更有利于后续的参数提取。这是因为PCA处理能去除能量比较小的噪声和干扰，使得谷值明显降低，同时能很好地保留各阶谐振峰，因此，得到了信噪比显著改善，峰谷比十分显著的频谱。 图2-4 周期谱图后PCA处理的SR信号幅度谱 2.3 洛仑兹拟合 本文在进行洛仑兹拟合时，采用的是分段拟合的方法，将经过前期处理后的幅度谱分成四段，第一段在8Hz谱峰左右，第二段在14Hz谱峰左右，第三段在20Hz谱峰左右，第四段在26Hz谱峰左右。洛仑兹拟合前后的波形图如图2-5所示，可见，洛仑兹拟合曲线基本可以反映舒曼谐振信号的幅度谱，得到了光滑的四个谱峰。但是，在第二阶谐振峰附近，由于频率为15Hz左右的强干扰的存在，导致谱峰向右偏移，这样会导致中心频率的提取不准确。因此，洛仑兹拟合的方法虽然可以得到光滑清晰的各阶谱峰，但是，对于舒曼谐振信号各分量各阶模式的幅度和中心频率的提取存在着一定的偏差。 图2-5 洛仑兹拟合前后SR信号幅度谱 2.4 包络滤波 对经过PCA及周期谱图法处理的信号的幅度谱做FFT，求得的幅度谱如图3-8所示，其中低频分量（5Hz）代表着幅度谱的包络，而图中高频部分（幅度很小）是幅度谱中的毛刺干扰。通过设计一个巴特沃斯IIR型低通滤波器，截止频率在4Hz左右，来去除高频干扰。 经过包络滤波处理的幅度谱如图2-6所示，对比处理前后的曲线，可见，经过包络滤波曲线变得平滑，且与图2-5中洛仑兹拟合的曲线对比，其峰值的位置不受干扰的影响，因此对于后续的中心频率的提取更为准确。 对比洛仑兹拟合和包络滤波得到的幅度谱曲线，可以看出，