[发明专利]基于北斗作为外辐射源的无源雷达直达波干扰抑制方法

权利要求书

1.一种基于北斗作为外辐射源的无源雷达直达波干扰抑制方法，其特征是：

(1)首先对M个阵元接收的信号分别进行下变频；

(2)将M个接收阵元组成互相重叠的K个子线阵；

(3)第1个子线阵在频域内的输出为X_l(ω)；

(4)将每个子线阵的输出移相并取模的平方；

(5)对此时的K个方程依次递减；

(6)将递减后的方程组转化为矩阵形式；

(7)利用最小二乘法求解均方意义下误差最小的目标回波估计值R＝B⁺X，再转换到时域即可得到估计的信号。

2.根据权利要求1所述的基于北斗作为外辐射源的无源雷达直达波干扰抑制方法，其特征是其具体步骤为：

雷达接收系统采用M个接收阵元组成互相重叠的K个子线阵，其中M为接收系统天线阵元个数，K为阵元组成的子线阵个数；

子线阵1：阵元1到阵元(M-K+1)

子线阵2：阵元2到阵元(M-K+2)

……        ……

子线阵K：阵元K到阵元M

第l(1≤l≤K)个子线阵的输出信号为：

Xl(ω)=R(ω)+Dl(ω)+Σi=1qMli(ω)---(1)]]>

式中，X_l(ω)为第l个子线阵在频域内输出的信号，D_l(ω)为第l个子线阵接收到的频域内直达波干扰，R(ω)为频域内雷达目标回波信号，M_li(ω)为频域内第l个子线阵接收到的第i路多径干扰信号，1≤i≤q，q为多径干扰总的个数；

对上式中第1个子线阵输出的信号移相并取模平方可得：

|Xl|2+|R|2-2Re(Xl)Re(R)-2Im(Xl)Im(R)=|Dl+Σi=1qMli|2---(2)]]>

式中，l＝1，2，·，K为子线阵序号，X_l为第l个子线阵在频域内输出的信号，R为频域内雷达目标回波信号，D_l为第l个子线阵接收到的频域内直达波干扰，M_li为频域内第l个子线阵接收到的第i路多径干扰信号，1≤i≤q，q为多径干扰总的个数，Re()表示一个信号的实部，Im()表示一个信号的虚部，||表示对信号取模；

将式(2)得到的K个方程依次相减可得：

2Re(X_p+1-X_p)Re(R)+2×Im(X_p+1-X_p)Im(R)＝|X_p+1|²-X_p|²-ε_p  (3)

式中，p＝1，2，·，K-1为子线阵的序号。为扰动项。

将式(3)转化为矩阵形式如下所示：

BR＝X+ε                            (4)

式中，R＝[Re(R)，Im(R)]^T为雷达目标回波信号矩阵，其中Re(R)和Im(R)分别为频域内目标回波信号的实部和虚部；

ε＝[ε₁，ε₂，ε₃，·，ε_p]                            (5)

X＝[|X₂|²-|X₁|²，|X₃|²-|X₂|²，·，|X_p+1|²-|X_p|²]        (6)

B=2Re(X2-X1)Im(X2-X1)Re(X3-X2)Im(X3-X2)··Re(Xp+1-Xp)Im(Xp+1-Xp),]]>p＝1，2，3，·，K    (7)

利用最小二乘法求解得到均方意义下误差最小的目标回波估计值：

R＝B⁺X                                                  (8)

式中，B⁺表示矩阵B的Moore-Penrose广义逆，R为目标回波信号估计值的频域表示形式，再将其转换到时域就得到目标回波信号基于相位匹配抑制直达波后的估计信号。