[发明专利]基于差分进化算法的分数阶希尔伯特变换阶次优化方法

说明书

技术领域

本发明涉及信号分析处理技术领域，特别是涉及一种基于差分进化算法的分数阶希尔伯特变换阶次优化方法。

背景技术

希尔伯特变换在信号处理中是重要的分析方法，通常用于瞬时参数的计算和窄带信号的表示。它将一个实信号变换为对应的解析信号，使信号只包含正频率部分，从而降低信号的采样率。希尔伯特变换定义为幅度相等、相移±90°的空间滤波器，对输入信号进行希尔伯特变换可有效地提取其包络，然后对包络信号进行频域分析则可以获得特征频率。

分数阶希尔伯特变换(FHT)定义为不同分数阶次(-2～2)下幅度相等、相移-90°～90°的改良型空间滤波器，利用较优阶次对信号进行分数阶希尔伯特变换包络分析，可以获得更显著的特征频率。而分数阶希尔伯特变换的阶次优化性成为信号分析时评估特征频率显著性的主要依据。在信号分析处理技术的发展中，对信号特征频率显著性的要求越来越高，因此信号进行包络分析时对分数阶希尔伯特变换的阶次的优化性提出了比较苛刻的要求。

如图1所示，分数阶希尔伯特变换阶次优化的传统方法的基本思想是，设置步进精度使阶次从-2至2等间隔改变，然后信号分析对比后得到优化后的阶次。步进精度即间隔为a，取值范围为-2～2，a值越小精度越高，最终阶次的误差则越小。信号分析方法的步骤为先在当前阶次p下进行信号FHT，然后对信号包络进行峭度(K)计算，峭度最大的阶次为最终阶次。

假设设置步进间隔a＝0.01，首先读取输入的信号数据并初始化当前阶次p＝0和目前最终阶次pm＝0。然后进入for循环在当前阶次p下进行信号FHT和峭度分析，若当前阶次p比目前最终阶次pm更好，则令pm＝p；否则跳过此步骤使当前阶次p加一步进间隔a，进入下一阶次的循环中。若阶次溢出即p大于1则完成优化结束程序，pm为获得的FHT的优化后的阶次。最后利用最终阶次pm对振动信号进行分数阶希尔伯特变换和包络分析得到特征频率。

传统分数阶希尔伯特变换阶次优化方法，程序的运行时间取决于阶次精度的设置，寻求阶次的更优是以牺牲运行时间为代价的。造成这种现象的原因，是程序利用for循环将步进精度范围内的每一个阶次全部遍历一遍，这样步进精度设置的越高(a越小)，for循环的次数越多，阶次优化的时间则越长。并且步进精度每高一个数量级时，程序的运行时间也会相应的高出一个数量级。如想得到10^-6精度的优化阶次，则需要进行106次的计算。

发明内容

本发明针对现有分数阶希尔伯特变换阶次优化方法的程序运行时间长、效率低的技术问题，提供一种程序运行时间短、效率高的基于差分进化算法的分数阶希尔伯特变换阶次优化方法。

为此，本发明提供一种基于差分进化算法的分数阶希尔伯特变换阶次优化方法，包括以下步骤：

步骤1，读取需要处理的信号数据，存入待处理的寄存器中，选择需要优化的阶次参数，设置优化的收敛精度，即允许的最大误差；

步骤2，使用差分进化算法，设置合适的算法参数和迭代次数，初始化种群，循环进行变异操作、交叉操作、峭度分析和选择操作来优化分数阶阶次，直到满足设置的收敛精度或达到最大迭代次数；

步骤3，对信号进行频域分析处理，使用优化后的阶次对信号进行分数阶希尔伯特变换，得到变换后的包络信号，对包络信号进行快速傅里叶变换，得到优化后的频谱信号和特征频率。

优选地，步骤2中初始化种群是在变量允许范围内随机生成，其表达式如下：

其中，rand()是0至1范围内的随机数，和分别代表变量的上下边界值；G表示第G代种群，包含NP数量的个体，种群中的每一个个体代表问题的一个解，每个个体有Dim维，其表达式如下：

X_G＝x_i，G＝{x_i1，G，x_i2，G，…，x_iDim，G} (1)

其中，i＝1，2，…，NP，G＝1，2，…，Gm；Gm代表最大进化数。

优选地，步骤2中差分变异操作的表达式如下：

v_i，G+1＝x_r1，G+F×(x_r2，G-x_r3，G) (3)

其中r1，r2，r3∈[1，NP]，是随机抽取的整数，且满足i≠r1≠r2≠r3，v_i，G+1表示变异矢量，F是变异算子。

优选地，步骤2中交叉操作表达式如下：