[发明专利]一种基于改进鸡群算法的盲源分离方法

权利要求书

1.一种基于改进鸡群算法的盲源分离方法，其特征在于，包括以下步骤：

(1)采集源信号S(t)＝[s₁(t),s₂(t),...,s_M(t)]^T，随机产生非奇异混合矩阵A，混合矩阵A对源信号S(t)进行线性混合得到观测信号X(t)＝[x₁(t),x₂(t),...,x_M(t)]^T，X(t)＝AS(t)；其中s_M(t)是源信号S(t)的第M个分量，x_M(t)是观测信号X(t)的第M个分量，t为时间序列，上标T表示共轭转置，M为正整数，A是M×M维矩阵；

(2)对步骤(1)中得到的观测信号X(t)进行中心化和白化后得到预处理信号Z(t)，中心化和白化为现有成熟技术，具体原理此处不再赘述；

(3)随机产生分离矩阵作为改进鸡群算法初始粒子,根据得到的预处理信号Z(t)，利用改进鸡群算法迭代更新得到最优解即最佳分离矩阵W，将X(t)送至W得到最佳分离信号Y(t)＝[y₁(t),y₂(t),...,y_M(t)]^T，Y(t)＝WX(t)，完成混合信号的分离。

2.如权利要求1所述的一种基于改进鸡群算法的盲源分离方法，其特征在于，所述步骤(3)中利用改进鸡群算法得到最优解即最佳分离矩阵W的具体步骤如下：

(3.1)初始化鸡群，设置最大迭代次数t₁＝M₁，鸡群粒子数N＝100，随机产生分离矩阵作为鸡群粒子，定义公鸡粒子个数N_R＝0.15N,母鸡粒子个数N_H＝0.7N，小鸡粒子个数N_C＝0.25N，妈妈母鸡粒子个数N_M＝0.5N_H；

(3.2)设置适应度函数fitness，将预处理后的信号Z(t)送至随机产生的分离矩阵(鸡群粒子)得到初始分离信号，对初始分离信号进行中心化、白化操作，代入适应度函数fitness计算鸡群粒子的适应度值，设定粒子当前最好位置和鸡群全局最好位置，鸡群算法迭代次数t₁＝1；

(3.3)如果t₁/G＝1(即当前为第一代)，从小到大排序适应度值并根据公鸡、母鸡和小鸡粒子个数确定公鸡、母鸡和小鸡的划分，建立鸡群等级制度，将鸡群分成数个子群并确定母鸡粒子和小鸡粒子的对应母子关系(每个子群中有一个公鸡粒子和若干母鸡粒子和小鸡粒子构成)，其中，G表示开始更新等级制度、支配关系和母子关系的代数，G＝10；

(3.4)根据公式(1)：

x_ir,j(t₁+1)＝x_ir,j(t₁)·(1+Φ(0,σ²))……………………………(1)

来更新公鸡粒子位置，其中，x_ir,j(t₁),x_ir,j(t₁+1)分别代表公鸡粒子ir在第t₁次和t₁+1次迭代中于第j维空间所处的位置；Φ(0,σ²)是一个方差为σ²的高斯分布，σ²表达式为：

${\mathrm{\σ}}^2=\left\{\begin{array}{cc}1& if{f}_{ir}\≤f_{kr}\\ \exp \left(\frac{f_{kr}-f_{ir}}{\left|f_{ir}\right|+\mathrm{\ϵ}}\right)& else\end{array}\right.,kr\∈\[1,N_R\],kr\≠ir\mathrm{...}\left(2\right)$

其中，f_ir和f_kr代表公鸡粒子ir和公鸡粒子kr的适应度值，ε是一个极小的常数，用来保证分母有意义，N_R为整个鸡群公鸡粒子数目，kr为所有公鸡粒子中除去ir后的任一个体，当公鸡粒子ir的适应度值好于公鸡粒子kr的适应度值，方差σ²为1，公鸡粒子ir的搜索空间变大，反之σ²降低，公鸡粒子ir的搜索空间缩小；

(3.5)母鸡粒子将跟随其子群的公鸡粒子进行搜索，同时也跟随其他子群的公鸡粒子进行搜索，根据公式(3)

x_ih,j(t₁+1)＝x_ih,j(t₁)+C₁·θ·(x_r1,j(t₁)-x_ih,j(t₁))+C₂·θ·(x_r2,j(t₁)-x_ih,j(t₁))……(3)

更新母鸡粒子位置，其中，x_r1,j(t₁),x_r2,j(t₁)分别代表母鸡粒子所属子群的公鸡粒子和其他子群公鸡粒子的位置信息，θ为0到1之间均匀分布的随机数，C₁和C₂分别代表母鸡粒子参考自身子群和其他子群权重，根据公式(4)、(5)

C₁＝exp((f_ih-f_r1)/(abs(f_ir+ε))…………………………(4)

C₂＝exp((f_r2-f_ir))…………………………(5)

得到，其中，f_ih和f_r1分别代表母鸡粒子ih和所属子群公鸡粒子r₁的适应度值，f_r2代表随机选取的其他子群公鸡粒子的适应度值；

(3.6)小鸡粒子不仅跟随其子群的妈妈母鸡粒子进行搜索，同时向所在子群公鸡学习，根据式(6)

x_ic,j(t₁+1)＝w·x_ic,j(t₁)+F·(x_m,j(t₁)-x_ic,j(t₁))+C·(x_r,j(t₁)-x_ic,j(t₁))…………(6)

更新小鸡位置，其中，x_m,j(t₁)代表小鸡粒子所跟随妈妈母鸡粒子的位置信息，x_r,j(t₁)代表妈妈母鸡粒子自身所在子群中的公鸡粒子位置信息，C为学习因子，取值0.5，表示小鸡粒子向自身所在子群中公鸡粒子学习的程度，w为小鸡粒子自身惯性权重，这里采用自适应惯性权重，通过公式(7)

得到粒子自身自适应权重，其中，wmax是惯性权重的最大值，wmin是惯性权重的最小值，τ取50，F为跟随系数，表示小鸡粒子跟随妈妈母鸡粒子寻找食物；

(3.7)利用公式(1)～(7)更新公鸡、母鸡和小鸡粒子位置后计算每个粒子的适应度值，更新鸡群的粒子当前最好位置和鸡群全局最好位置；

(3.8)t₁＝t₁+1，若达到迭代次数，则停止迭代，得到最优位置(最优解)，即最佳分离矩阵W，否则转到步骤(3.3)。

3.如权利要求2所述的一种基于改进鸡群算法的盲源分离方法，其特征在于，所述改进鸡群算法适应度函数fitness为分离信号负熵与峭度绝对值的加权平均，具体实现为：

fitness＝-(0.4×fitness1+0.6×fitness2|/2)…………………(8)

式(8)中，表示信号负熵，其中y_i表示第i个分离信号；表示信号峭度的绝对值，其中kurt(y_i)为第i个分离信号的峭度，适应度函数值越小表明分离效果越好。