[发明专利]一种基于EEMD的HHT的太阳黑子面积周期特征分析的方法

权利要求书

1.一种基于EEMD的HHT的太阳黑子面积周期特征分析的方法，其特征在于：所述方法的具体步骤如下：

步骤1、对太阳黑子面积数据进行平滑处理，得到平滑后的数据；

步骤2：将平滑后的数据进行集成经验模式分解，分解后得到从高频到低频的一系列不同频率的本质模态函数集合；其中本质模态函数集合包括所有的本质模态函数分量及残差，每个本质模态函数分量代表着太阳黑子面积信号中的一个固有的振动模态；

步骤3：根据时间间隔确定采样频率，根据数据长度和采样频率得到频率分辨率，由奈奎斯特定理得出最大、最小频率值，对本质模态函数集合中的每一个IMF分量进行希尔伯特黄变换，得到其对应的瞬时频率谱；

步骤4：统计每一个瞬时频率谱中每个频率出现的次数，将统计的频率出现的次数的最大值变为1，并以该值为基准对所有频率的出现次数作归一化处理，得到对应的频率概率分布直方图；

步骤5：用高斯函数对频率概率分布直方图进行单高斯拟合，拟合出的高斯曲线的峰值的横坐标作为该本质模态函数的平均频率，取一个标准差的置信区间，得到该本质模态函数的频率区间；

步骤6：步骤5中得到的平均频率的倒数为平均周期，平均频率加一个标准差的和的倒数为该本质模态函数的最小周期，平均频率减去一个标准差的差的倒数为该本质模态函数的最大周期；

所述步骤3中，希尔伯特黄变换参数：采样频率f_s=1/150，频率分辨率f₀=f_s/N，由奈奎斯特定理，最大频率f_max=f_s /2，最小频率f_min=f₀；其中，N为平滑后数据的长度；

所述步骤2中，集成经验模式分解参数：添加噪声等级Nstd=0.2，添加噪声次数NR=300，最大本质模态函数个数numIMF=7，最大迭代筛选次数maxSift=10，得到IMF1-IMF6和残差。