发明名称 干扰机与接收机联合组网的SAR欺骗干扰方法
摘要 干扰机和接收机联合组网的SAR欺骗干扰方法,步骤如下:设置若干台接收机与干扰机组网,通过干扰机和接收机截获雷达信号,分析SAR信号参数;利用截获的雷达时域信号计算各接收机到干扰机的瞬时斜距差;计算各接收机对应的瞬时位置距离参数;构造系数矩阵,计算系数矩阵的逆;根据系数矩阵及瞬时位置距离参数计算调制参数;计算干扰场景中的点相对于干扰机的瞬时斜距差;对干扰机截获的信号进行干扰调制生成干扰信号;转发欺骗干扰信号,对方雷达得到真实场景雷达信号与虚假场景欺骗干扰信号叠加的欺骗干扰信号。本发明方法避免了近似误差对SAR欺骗干扰效果的影响,可获得逼真的干扰图像,实现良好的干扰效果,具有简单,高效,精确的优点。
申请公布号 CN103760532A 申请公布日期 2014.04.30
申请号 CN201410015031.4 申请日期 2014.01.13
申请人 西安电子科技大学 发明人 周峰;赵博;崔乐园;王海兵
分类号 G01S7/38(2006.01)I;G01S13/90(2006.01)I 主分类号 G01S7/38(2006.01)I
代理机构 广东秉德律师事务所 44291 代理人 杨焕军
主权项 1.干扰机和接收机联合组网的SAR欺骗干扰方法,其特征在于,包括以下步骤:步骤1、设置若干台接收机,将接收机与干扰机组网,通过干扰机和接收机分别截获雷达信号,分析SAR信号参数;通过干扰机以及p台接收机A<sub>l</sub>截获雷达时域信号,l=1,…,p,p≥3,分别得到以距离为行向量、以方位为列向量的雷达信号数据:包括干扰机雷达信号数据s<sub>O</sub>(t<sub>r</sub>,t<sub>a</sub>)和接收机雷达信号数据s<sub>Al</sub>(t<sub>r</sub>,t<sub>a</sub>),对雷达信号数据进行分析,得到SAR信号参数,其中,t<sub>r</sub>为距离向快时间,t<sub>a</sub>为方位向慢时间;步骤2、利用截获的雷达时域信号,分别计算各台接收机到干扰机的瞬时斜距差w<sub>l</sub>(t<sub>a</sub>);<maths num="0001"><![CDATA[<math><mrow><msub><mi>w</mi><mi>l</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mi>c</mi><mrow><mn>2</mn><mo>&CenterDot;</mo><msub><mi>F</mi><mi>s</mi></msub></mrow></mfrac><mo>&CenterDot;</mo><mrow><mo>(</mo><mfrac><msub><mi>N</mi><mi>r</mi></msub><mn>2</mn></mfrac><mo>-</mo><msub><mi>max</mi><mi>ID</mi></msub><mo>{</mo><mi>IFFT</mi><mo>[</mo><mi>FFT</mi><mrow><mo>(</mo><msub><mi>s</mi><mi>o</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>r</mi></msub><mo>,</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>conj</mi><mrow><mo>(</mo><mi>FFT</mi><mrow><mo>(</mo><msub><mi>s</mi><mi>Al</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>r</mi></msub><mo>,</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>)</mo></mrow><mo>]</mo><mo>}</mo><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>其中,c表示光速,N<sub>r</sub>为距离采样点数,F<sub>s</sub>为距离向采样频率,s<sub>O</sub>(t<sub>r</sub>,t<sub>a</sub>)为干扰机的雷达信号数据,s<sub>Al</sub>(t<sub>r</sub>,t<sub>a</sub>)为第l台接收机的雷达信号数据,max<sub>ID</sub>(·)表示返回最大值位置的运算,IFFT[·]表示逆傅里叶变换,FFT(·)表示傅里叶变换,conj(·)表示共轭;步骤3、计算各台接收机对应的瞬时位置距离参数e<sub>l</sub>(t<sub>a</sub>);根据接收机的位置和对应接收机的瞬时斜距差,分别计算各接收机的瞬时位置距离参数e<sub>l</sub>(t<sub>a</sub>):e<sub>l</sub>(t<sub>a</sub>)=-x<sub>l</sub><sup>2</sup>-y<sub>l</sub><sup>2</sup>+(w<sub>l</sub>(t<sub>a</sub>))<sup>2</sup>,式中的x<sub>l</sub>,y<sub>l</sub>为第l台接收机A<sub>l</sub>的位置坐标,w<sub>l</sub>(t<sub>a</sub>)为第l台接收机到干扰机的瞬时斜距差;步骤4、构造系数矩阵P(t<sub>a</sub>),并计算系数矩阵的逆P<sup>-1</sup>(t<sub>a</sub>);<maths num="0002"><![CDATA[<math><mrow><mi>P</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mn>1</mn></msub></mtd><mtd><msub><mi>y</mi><mn>1</mn></msub></mtd><mtd><msub><mi>w</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msub><mi>x</mi><mn>1</mn></msub></mtd><mtd><msub><mi>y</mi><mi>l</mi></msub></mtd><mtd><msub><mi>w</mi><mi>l</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msub><mi>x</mi><mi>p</mi></msub></mtd><mtd><msub><mi>y</mi><mi>p</mi></msub></mtd><mtd><msub><mi>w</mi><mi>p</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>式中的x<sub>l</sub>,y<sub>l</sub>为第l台接收机A<sub>l</sub>的位置坐标,w<sub>l</sub>(t<sub>a</sub>)为第l台接收机到干扰机的瞬时斜距差;步骤5、根据系数矩阵以及接收机的瞬时位置距离参数计算调制参数Q;<maths num="0003"><![CDATA[<math><mrow><mi>Q</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>a</mi></mtd></mtr><mtr><mtd><mi>b</mi></mtd></mtr><mtr><mtd><mi>c</mi></mtd></mtr></mtable></mfenced><mo>=</mo><msup><mi>P</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow><mi>e</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>,</mo><mi>e</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>e</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msub><mi>e</mi><mi>l</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msub><mi>e</mi><mi>p</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>其中,e<sub>l</sub>(t<sub>a</sub>)为第l台接收机的瞬时位置距离参数;步骤6、计算干扰场景中的点相对于干扰机的瞬时斜距差ΔR<sub>Jk</sub>(t<sub>a</sub>);设计干扰场景中点的分布位置,干扰场景中任意一点相对于干扰机的瞬时斜距差<maths num="0004"><![CDATA[<math><mrow><mi>&Delta;</mi><msub><mi>R</mi><mi>Jk</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>=</mo><msqrt><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>Jk</mi></msub><mo>+</mo><mfrac><mi>a</mi><mn>2</mn></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>Jk</mi></msub><mo>+</mo><mfrac><mi>b</mi><mn>2</mn></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mrow><mo>(</mo><mo>-</mo><mfrac><msup><mi>a</mi><mn>2</mn></msup><mn>4</mn></mfrac><mo>-</mo><mfrac><msup><mi>b</mi><mn>2</mn></msup><mn>4</mn></mfrac><mo>+</mo><mfrac><msup><mi>c</mi><mn>2</mn></msup><mn>4</mn></mfrac><mo>)</mo></mrow></msqrt><mo>+</mo><mfrac><mi>c</mi><mn>2</mn></mfrac><mo>,</mo></mrow></math>]]></maths>其中,x<sub>Jk</sub>、y<sub>Jk</sub>为干扰场景中第k个点的位置坐标,k=1,2,…,N,N表示干扰场景中设置的点个数,a、b、c为调制参数Q中的元素;步骤7、对干扰机截获的信号进行干扰调制,生成干扰信号;根据干扰场景中各点相对于干扰机的瞬时斜距差,计算延时调制系数<img file="FDA0000456259550000024.GIF" wi="557" he="171" />和相位调制系数<img file="FDA0000456259550000025.GIF" wi="588" he="168" />对干扰机截获的信号进行干扰调制,得到虚假目标干扰信号s<sub>J</sub>(t<sub>r</sub>,t<sub>a</sub>):<maths num="0005"><![CDATA[<math><mrow><msub><mi>s</mi><mi>J</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>r</mi></msub><mo>,</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>IFFT</mi><mo>{</mo><mi>FFT</mi><mo>[</mo><msub><mi>s</mi><mi>o</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>r</mi></msub><mo>,</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>]</mo><mo>&CenterDot;</mo><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>r</mi></msub><mfrac><mrow><mn>2</mn><mo>&CenterDot;</mo><mi>&Delta;</mi><msub><mi>R</mi><mi>Jk</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow></mrow><mi>c</mi></mfrac><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>c</mi></msub><mfrac><mrow><mn>2</mn><mo>&CenterDot;</mo><mi>&Delta;</mi><msub><mi>R</mi><mi>Jk</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow></mrow><mi>c</mi></mfrac><mo>)</mo></mrow><mo>]</mo><mo>}</mo><mo>;</mo></mrow></math>]]></maths>其中,s<sub>O</sub>(t<sub>r</sub>,t<sub>a</sub>)为干扰机雷达信号数据,j表示虚数单位,f<sub>r</sub>为距离向频率采样序列,c表示光速,f<sub>c</sub>为雷达载频,ΔR<sub>Jk</sub>(t<sub>a</sub>)为干扰场景中第k个点相对于干扰机的瞬时斜距差;步骤8、转发虚假目标干扰信号;将虚假场景的欺骗干扰信号转发,对方雷达得到真实场景雷达信号与虚假场景欺骗干扰信号叠加的欺骗干扰信号。
地址 710068 陕西省西安市太白南路2号