[发明专利]基于多基地雷达的杂波子空间下双门限恒虚警检测方法

权利要求书

1.一种基于多基地雷达的杂波子空间下双门限恒虚警检测方法，包括如下步骤：

(1)对雷达回波信号进行采集，得到N×L×(K+1)维原始数据矩阵：X＝[X₁,X₂,...,X_n,...,X_N]，X_n表示第n个局部雷达站中的原始数据，n＝1,2,...,N，N为局部雷达站个数，L为一个非相干积累周期中所包含的脉冲个数，K代表检测单元附近参考单元个数；

(2)对数据矩阵X中的元素进行广义似然比检测，得到N×L×(K+1)维经广义似然比检测后的局部检验统计量Z＝[Z₁,Z₂,...,Z_n,...,Z_N]，Z_n表示Z中的第n个数据元素；

(3)计算第一门限τ₁：

根据脉冲个数L、参考单元个数K、局部雷达站干扰个数D及局部虚警概率p_fa计算第一门限τ₁：

τ₁＝-1/(K-L+D+1)×log(p_fa)，

其中，log(·)表示对数运算；

(4)将局部检验统计量Z中的每个元素Z_n与第一门限τ₁相比较：

若Z_n≥τ₁，则将数据Z_n传输到融合中心，记传输到融合中心的数据为B＝[B₁,…,B_m,…,B_M]，B_m为B中第m个元素，m＝1,2,...,M，M为传输到融合中心的数据个数；

若Z_n＜τ₁，则数据Z_n不传输；

(5)根据传输到融合中心的数据个数m、脉冲个数L，参考单元个数K，干扰子空间个数D，局部虚警概率p_fa以及全局虚警概率P_fa计算第二门限η₂：

η₂＝{η₂:F(η₂|N,1/(K-L+D+1))＝P_fa}，

其中，

t表示积分因子，！表示阶乘运算；

(6)对传输到融合中心数据B中的元素进行求和，得到全局检验统计量E，将全局检验统计量E与第二门限η₂比较：若E≥η₂，则判决目标存在，若E＜η₂，则判决目标不存在。

2.根据权利要求1所述的方法，其中步骤(1)中第n个局部雷达站中的原始数据X_n为L×(K+1)维的矩阵，其表示如下：

其中，x_nls表示第n个局部雷达站中第l个脉冲的第s+1个数据，l＝1,2,...,L，s＝0,1,...,K，L为一个非相干积累周期中所包含的脉冲个数，s＝0时，x_nls代表检测单元信息，s＝1,2,...,K时，x_nls代表检测单元附近参考单元信息，K代表检测单元附近参考单元个数。

3.根据权利要求1所述的方法，其中步骤(2)中对数据矩阵X中的元素进行广义似然比检测，按如下步骤进行：

2.1)根据已知的L×D维结构化干扰子空间矩阵H_n，计算L×(L-D)维干扰补空间U_n：

2.1.1)根据已知的L×D维结构化干扰子空间H_n计算L×L维补空间映射E_n：

E_n＝I-H_nH_n⁺，

其中，D表示局部雷达站的干扰源个数，I表示L×L维单位矩阵，(·)⁺表示共轭伪矩阵；

2.1.2)取L×L维补空间映射E_n的前(L-D)列元素，构成L×(L-D)维干扰补空间矩阵U_n；

2.2)根据L×(L-D)维干扰补空间U_n，计算局部检验统计量Z_n：

其中，x_n0表示原始数据X_n中第一列数据元素，x_ni表示原始数据X_n中第i+1列数据元素，表示求共轭矩阵，(·)^-1表示求矩阵逆运算，(·)^Τ表示求矩阵转置，|·|²表示求数据模平方。