一种基于子空间分解的多辐射源被动直接时差定位方法

1.一种基于子空间分解的多辐射源被动直接时差定位方法，其特征在于，包括以下步骤：步骤1：初始化设置，将辐射源位置坐标的搜索网格、观测网络中各观测站的位置坐标写入主机内存；步骤2：将各观测站的接收信号序列进行离散傅里叶变换后拼接得到总的观测样本向量，并利用多帧观测样本向量计算样本协方差矩阵并将其谱分解得到噪声子空间和辐射源个数；步骤3：利用噪声子空间与由位置坐标参数化的时差矩阵的正交性构造搜索目标函数；步骤4：在同频辐射源位置坐标的搜索网格上搜索目标函数峰值估计辐射源的位置坐标。 2.根据权利要求1所述的一种基于子空间分解的多辐射源被动直接时差定位方法，其特征在于，所述步骤1的具体执行过程如下：初始化设置，将同频辐射源位置坐标的搜索网格l k ∈Ω,k＝1,...,K和各观测站的位置坐标p n ,n＝1,...,N写入主机内存。 3.根据权利要求2所述的一种基于子空间分解的多辐射源被动直接时差定位方法，其特征在于，所述步骤2的具体执行过程如下：将各观测站的接收信号序列进行离散傅里叶变换后拼接得到总的观测样本向量 其中，是第n个观测站接收信号的第t帧离散频谱向量，f m ,m＝1,...,M为频点；利用L帧观测样本向量计算样本协方差矩阵 将 谱分解，并将特征值分为组内方差最小且组间均值差最大的两组特征值：其中， 和 分别为大特征值组为对角元素的对角矩阵和小特征值组为对角元素的对角矩阵；小特征值组对应的特征向量构成噪声子空间 大特征值组对应的特征向量构成信号子空间 辐射源个数P等于对角矩阵Λ S 阶数的 4.根据权利要求3所述的一种基于子空间分解的多辐射源被动直接时差定位方法，其特征在于，所述步骤3的具体执行过程如下：位置坐标l k 参数化的时差矩阵为：利用噪声子空间U N 与由位置坐标l k 参数化的时差矩阵Q n (l k )的正交性构造搜索目标函数：其中，λ min (·)表示求最小特征值；Ψ(l k )＝[Vec(G 1 Q 1 (l k )) Vec(G 2 Q 2 (l k )) … Vec(G N Q N (l k ))]，Vec(·)为向量化运算； 即将噪声子空间的转置 按每M列分块。 5.根据权利要求4所述的一种基于子空间分解的多辐射源被动直接时差定位方法，其特征在于，所述步骤4的具体执行过程如下：在同频辐射源位置坐标的搜索网格上搜索目标函数最大的P个峰值坐标：则确定的辐射源位置为：