[发明专利]基于小生境差分进化的多重信号分类谱峰搜索方法

权利要求书

1.一种基于小生境差分进化的多重信号分类谱峰搜索方法，其特征在于，所述搜索方法包括如下步骤：

(1)针对q个远场均匀横电磁波入射到M个阵元的阵列的情形，建立阵列接收数据模型X(t)＝A(θ)S(t)+N(t)；其中：

A(θ)＝[a(θ₁),a(θ₂),…,a(θ_q)]是天线阵列的空域导向矢量矩阵；

是第i个信号的导向矢量，i＝1,2,…,q；上角标T表示矩阵的转置；

φ_m(θ_i)是第i个信号s_i(t)入射到第m个阵元上时，相对于参考阵元的滞后相位；

S(t)＝[s₁(t) s₂(t) … s_q(t)]^T是波达方向为θ₁,θ₂,…,θ_q的入射信号矢量；

N(t)＝[n₁(t) n₂(t) … n_M(t)]^T是加性零均值白噪声矢量，功率为

(2)计算X(t)的协方差矩阵R＝E[X(t)X^H(t)]；

(3)对协方差矩阵R进行特征分解其中∑_S为信号特征值组成的对角阵，∑_N为噪声特征值组成的对角阵，得到M×q维信号子空间矩阵U_S和M×(M-q)维噪声子空间矩阵U_N，空间谱函数的q个极大值点为所需估计的q个远场入射信号的波达方向；

(4)针对于空间谱函数以基于小生境差分进化的搜索方法代替原有固定步长搜索，设置种群数目p_num及总迭代次数total，并对初始种群进行初始化，得到初始种群个体{⁰x₁,⁰x₂,…,⁰x_{p_num}}，其中左上角标识表示该种群样本是迭代第几次的样本；

(5)开始迭代过程，计算交叉概率CR(n)及缩放因子F(n)，其中n表示第n次迭代，n∈[0,total]；每一代的F(n)及CR(n)通过下式给出：

其中n是当前进化代数；F_min,F_max限定了缩放因子F(n)的取值范围满足F(n)∈[F_min,F_max]；CR_min,CR_max限定了交叉概率CR(n)的取值范围满足CR(n)∈[CR_min,CR_max]；

(6)通过自适应小生境识别方法，将种群{ⁿx₁,ⁿx₂,…,ⁿx_{p_num}}中的每一个体归类至各自的小生境中，并保留每个小生境的中心点；

(7)对每个小生境中的每一个种群个体进行差分进化操作；

(8)若进化后个体的适应度大于父代个体的适应度，则用进化后的种群个体代替父代个体；

(9)重复步骤(5)至(8)，直至迭代次数达到total；此时每个小生境的中心点x_c1,x_c2,…,x_cq即为估计的q个入射信号的波达方向。

2.如权利要求1所述的方法，其特征在于，所述步骤(1)中，第i个信号入射到天线阵列原点，设第m个阵元在坐标系中的位置坐标为第i个信号在第m个阵元处所产生的相对于参考阵元的时延为对应的相位差为λ为信号波长，c为光速。

3.如权利要求1所述的方法，其特征在于，所述步骤(2)中，通过有限长的观测数据来近似得到协方差矩阵R，根据最大似然原理，R的估计值满足其中L是采样快拍数。