[发明专利]一种基于稀疏重构最优化算法的信号包络线提取方法

摘要

本发明公开了一种基于稀疏重构最优化算法的信号包络线提取方法，其采用了基于粒子群稀疏重构最优化算法，将改变DCT基的频带宽度的变化因子作为基于粒子群稀疏重构最优化算法中的每个粒子的属性，通过多次迭代获得上包络线全局最优位置和下包络线全局最优位置，对应作为用于改变粒子确定的上包络线获取过程所需的DCT基的频带宽度的变化因子和用于改变粒子确定的下包络线获取过程所需的DCT基的频带宽度的变化因子，进而得到最适合上包络线变化趋势的最佳DCT基和最适合下包络线变化趋势的最佳DCT基，最终自适应地提取出最佳上包络线和最佳下包络线，提高了上包络线和下包络线的提取精度，同时提高了上包络线和下包络线提取的抗噪性能。

主权项

1.一种基于稀疏重构最优化算法的信号包络线提取方法，其特征在于包括以下步骤：步骤1)：将所要提取的包络线信号记为x(t)，将x(t)表示为列向量形式，x(t)＝[x1(t),x2(t),…,xN‑1(t),xN(t)]T；其中，t表示采样时间，t的单位为秒，0≤t≤3，N表示采样点的总个数，N为正整数，N≥250，x1(t)、x2(t)、xN‑1(t)、xN(t)对应表示x(t)中的第1个采样值、第2个采样值、第N‑1个采样值、第N个采样值，符号“[]”为向量表示符号，[x1(t),x2(t),…,xN‑1(t),xN(t)]T为[x1(t),x2(t),…,xN‑1(t),xN(t)]的转置；步骤2)：从x(t)中提取出所有极大值点和所有极小值点；然后将从x(t)中提取出的所有极大值点构成的列向量记为y_d(t)，并将从x(t)中提取出的所有极小值点构成的列向量记为y_s(t)，其中，x_d,1、x_d,2和对应表示从x(t)中提取出的第1个极大值点、第2个极大值点和第K₁个极大值点，K₁表示从x(t)中提取出的极大值点的总个数，K₁为正整数，1≤K₁<N，为的转置，x_s,1、x_s,2和对应表示从x(t)中提取出的第1个极小值点、第2个极小值点和第K₂个极小值点，K₂表示从x(t)中提取出的极小值点的总个数，K₂为正整数，1≤K₂<N，为的转置；步骤3)：基于粒子群稀疏重构最优化算法，提取x(t)的上包络线和下包络线，具体过程为：3)_1、令α表示粒子群的迭代次数，令αmax表示粒子群的最大迭代次数；设粒子群中有L个粒子；其中，α为正整数，α的初始值为1，αmax为正整数，30≤αmax≤200，L为正整数，1并计算当前粒子所对应的x(t)的下包络线，记为3)_8、构建一个维数为(N‑1)×N的一阶微分矩阵，记为D，接着根据D和确定当前粒子所对应的用于体现的变化程度的上包络线平稳度函数，记为E_d,i(t)，并根据D和确定当前粒子所对应的用于体现的变化程度的下包络线平稳度函数，记为E_s,i(t)，其中，符号“|| ||₂”为求矩阵的2‑范数符号；3)_9、计算当前粒子的上包络线适应度函数，记为f_d,i，并计算当前粒子的下包络线适应度函数，记为其中，λ为正则化参数，λ>0，符号“|| ||₁”为求矩阵的1‑范数符号；3)_10、在第α次迭代过程中，将粒子群中下一个待处理的粒子作为当前粒子，然后返回步骤3)_3继续执行，直至粒子群中的所有粒子处理完毕；3)_11、在第α次迭代过程中，根据L个上包络线适应度函数，从L个粒子中找出使上包络线平稳度函数最小的粒子，将找出的粒子的位置作为粒子群中的所有粒子的上包络线全局最优位置，记为gd,best；并针对粒子群中的每个粒子，从第1次迭代过程至第α次迭代过程中找出使上包络线平稳度函数最小的迭代序号，将找出的迭代序号对应的迭代过程中该粒子的位置作为该粒子的上包络线局部最优位置，将粒子群中的第i个粒子的上包络线局部最优位置记为pd,i,best；同样，根据L个下包络线适应度函数，从L个粒子中找出使下包络线平稳度函数最大的粒子，将找出的粒子的位置作为粒子群中的所有粒子的下包络线全局最优位置，记为gs,best；并针对粒子群中的每个粒子，从第1次迭代过程至第α次迭代过程中找出使下包络线平稳度函数最大的迭代序号，将找出的迭代序号对应的迭代过程中该粒子的位置作为该粒子的下包络线局部最优位置，将粒子群中的第i个粒子的下包络线局部最优位置记为ps,i,best；3)_12、判断α<αmax是否成立，如果α<αmax成立，则更新上包络线获取过程中粒子群中的每个粒子的速度和位置、下包络线获取过程中粒子群中的每个粒子的速度和位置，将上包络线获取过程中粒子群中的第i个粒子更新后的速度记为vd,i,back，vd,i,back＝w×vd,i,pre+c1×r1×(pd,i,best,pre‑pd,i,pre)+c2×r2×(gd,best,pre‑pd,i,best,pre)，将上包络线获取过程中粒子群中的第i个粒子更新后的位置记为pd,i,back，pd,i,back＝pd,i,pre+vd,i,pre，将下包络线获取过程中粒子群中的第i个粒子更新后的速度记为vs,i,back，vs,i,back＝w×vs,i,pre+c1×r1×(ps,i,best,pre‑ps,i,pre)+c2×r2×(gs,best,pre‑ps,i,best,pre)，将下包络线获取过程中粒子群中的第i个粒子更新后的位置记为ps,i,back，ps,i,back＝ps,i,pre+vs,i,pre；然后将上包络线获取过程中粒子群中的每个粒子更新后的速度和位置对应作为下一次迭代过程中上包络线获取过程中该粒子的速度和位置，将下包络线获取过程中粒子群中的每个粒子更新后的速度和位置对应作为下一次迭代过程中下包络线获取过程中该粒子的速度和位置，接着令α＝α+1；再返回步骤3)_2继续执行；其中，w为惯性因子，w的取值范围为[0,1]，vd,i,pre表示上包络线获取过程中粒子群中的第i个粒子更新前的速度，c1和c2均为学习因子，r1和r2均为区间[0,1]内的随机数，pd,i,best,pre表示上包络线获取过程中粒子群中的第i个粒子更新前的上包络线局部最优位置，pd,i,pre表示上包络线获取过程中粒子群中的第i个粒子更新前的位置，gd,best,pre表示更新前粒子群中的所有粒子的上包络线全局最优位置，vs,i,pre表示下包络线获取过程中粒子群中的第i个粒子更新前的速度，ps,i,best,pre表示下包络线获取过程中粒子群中的第i个粒子更新前的下包络线局部最优位置，ps,i,pre表示下包络线获取过程中粒子群中的第i个粒子更新前的位置，gs,best,pre表示更新前粒子群中的所有粒子的下包络线全局最优位置，α＝α+1中的“＝”为赋值符号；如果α<αmax不成立，则结束迭代过程，然后将第α次迭代过程中获得的gd,best赋值给用于改变一个粒子确定的上包络线获取过程所需的DCT基的频带宽度的变化因子，将第α次迭代过程中获得的gs,best赋值给用于改变一个粒子确定的下包络线获取过程所需的DCT基的频带宽度的变化因子，再按照步骤3)_4至步骤3)_7的过程，以相同的方式执行一次，获得x(t)的最终上包络线即为最佳上包络线和x(t)的最终下包络线即为最佳下包络线。