[发明专利]干扰机与接收机联合组网的SAR欺骗干扰方法

权利要求书

1.干扰机和接收机联合组网的SAR欺骗干扰方法，其特征在于，包括以下步骤：

步骤1、设置若干台接收机，将接收机与干扰机组网，通过干扰机和接收机分别截获雷达信号，分析SAR信号参数；

通过干扰机以及p台接收机A_l截获雷达时域信号，l=1,…,p，p≥3，分别得到以距离为行向量、以方位为列向量的雷达信号数据：包括干扰机雷达信号数据s_O(t_r,t_a)和接收机雷达信号数据s_Al(t_r,t_a)，对雷达信号数据进行分析，得到SAR信号参数，其中，t_r为距离向快时间，t_a为方位向慢时间；

步骤2、利用截获的雷达时域信号，分别计算各台接收机到干扰机的瞬时斜距差w_l(t_a)；

wl(ta)=c2·Fs·(Nr2-maxID{IFFT[FFT(so(tr,ta))·conj(FFT(sAl(tr,ta)))]});]]>

其中，c表示光速，N_r为距离采样点数，F_s为距离向采样频率，s_O(t_r,t_a)为干扰机的雷达信号数据，s_Al(t_r,t_a)为第l台接收机的雷达信号数据，max_ID(·)表示返回最大值位置的运算，IFFT[·]表示逆傅里叶变换，FFT(·)表示傅里叶变换，conj(·)表示共轭；

步骤3、计算各台接收机对应的瞬时位置距离参数e_l(t_a)；

根据接收机的位置和对应接收机的瞬时斜距差，分别计算各接收机的瞬时位置距离参数e_l(t_a)：

e_l(t_a)=-x_l²-y_l²+(w_l(t_a))²，

式中的x_l，y_l为第l台接收机A_l的位置坐标，w_l(t_a)为第l台接收机到干扰机的瞬时斜距差；

步骤4、构造系数矩阵P(t_a)，并计算系数矩阵的逆P^-1(t_a)；

P(ta)=x1y1w1(ta)·········x1ylwl(ta)·········xpypwp(ta);]]>

式中的x_l，y_l为第l台接收机A_l的位置坐标，w_l(t_a)为第l台接收机到干扰机的瞬时斜距差；

步骤5、根据系数矩阵以及接收机的瞬时位置距离参数计算调制参数Q；

Q=abc=P-1(ta)e(ta),e(ta)=e1(ta)···el(ta)···ep(ta);]]>

其中，e_l(t_a)为第l台接收机的瞬时位置距离参数；

步骤6、计算干扰场景中的点相对于干扰机的瞬时斜距差ΔR_Jk(t_a)；

设计干扰场景中点的分布位置，干扰场景中任意一点相对于干扰机的瞬时斜距差ΔRJk(ta)=(xJk+a2)2+(yJk+b2)2+(-a24-b24+c24)+c2,]]>其中，x_Jk、y_Jk为干扰场景中第k个点的位置坐标，k=1,2,…,N，N表示干扰场景中设置的点个数，a、b、c为调制参数Q中的元素；

步骤7、对干扰机截获的信号进行干扰调制，生成干扰信号；

根据干扰场景中各点相对于干扰机的瞬时斜距差，计算延时调制系数和相位调制系数对干扰机截获的信号进行干扰调制，得到虚假目标干扰信号s_J(t_r,t_a)：

sJ(tr,ta)=IFFT{FFT[so(tr,ta)]·[Σk=1Nexp(-j2πfr2·ΔRJk(ta)c)·exp(-j2πfc2·ΔRJk(ta)c)]};]]>

其中，s_O(t_r,t_a)为干扰机雷达信号数据，j表示虚数单位，f_r为距离向频率采样序列，c表示光速，f_c为雷达载频，ΔR_Jk(t_a)为干扰场景中第k个点相对于干扰机的瞬时斜距差；

步骤8、转发虚假目标干扰信号；将虚假场景的欺骗干扰信号转发，对方雷达得到真实场景雷达信号与虚假场景欺骗干扰信号叠加的欺骗干扰信号。