发明名称 基于分数阶傅里叶域信道化的宽带线性调频信号多假目标干扰方法
摘要 本发明涉及基于分数阶傅里叶域信道化的宽带线性调频信号多假目标干扰方法,一种基于分数阶傅里叶域信道化,在脉内对宽带线性调频信号进行多假目标欺骗的干扰方法,属于电子干扰技术领域。本发明提出的基于分数阶傅里叶域信道化的宽带LFM信号多假目标干扰方法利用分数阶傅里叶变换对非平稳信号的能量聚焦性,可以将宽带线性调频信号均聚焦到一个信道内输出,保证了信号的完整性,很好地解决了宽带LFM信号在传统傅里叶域信道化中能量会溢出到两个或更多的信道中,造成信号的能量削弱和波形的失真;同时降低了DRFM系统的瞬时带宽要求。
申请公布号 CN103532656B 申请公布日期 2016.11.09
申请号 CN201310344201.9 申请日期 2013.08.08
申请人 北京理工大学 发明人 陶然;赵兴浩;戚士斌
分类号 H04K3/00(2006.01)I 主分类号 H04K3/00(2006.01)I
代理机构 代理人
主权项 基于分数阶傅里叶域信道化的宽带线性调频信号多假目标干扰方法,其特征在于:步骤一、根据雷达发射的线性调频信号x(n)的调频率μ选定分数阶傅里叶域的变换阶次p及分数阶傅里叶域相对于傅里叶域的逆时针旋转角度α,α=pπ/2,构建分数阶傅里叶域滤波器组{h<sub>k,p</sub>(n)}<sub>k=0,1,…K‑1</sub>;<maths num="0001"><math><![CDATA[<mrow><msub><mi>h</mi><mrow><mi>k</mi><mo>,</mo><mi>p</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>h</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>cot</mi><mrow><mo>(</mo><mi>&alpha;</mi><mo>)</mo></mrow><msup><mrow><mo>(</mo><mi>n</mi><mi>&Delta;</mi><mi>t</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></msup><mo>=</mo><msub><mi>h</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mi>j</mi><mi>k</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>K</mi></mfrac><mi>n</mi></mrow></msup><mo>&CenterDot;</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>cot</mi><mrow><mo>(</mo><mi>&alpha;</mi><mo>)</mo></mrow><msup><mrow><mo>(</mo><mi>n</mi><mi>&Delta;</mi><mi>t</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000984984420000011.GIF" wi="1617" he="134" /></maths>其中<img file="FDA0000984984420000012.GIF" wi="549" he="108" />为传统傅里叶域的带通滤波器组,<img file="FDA0000984984420000013.GIF" wi="485" he="103" />为传统傅里叶域的高通滤波器;h<sub>k</sub>(n)=h<sub>0</sub>(n)为传统傅里叶域的低通滤波器;K为信道数,K为2的幂次方;n为滤波器的点数;△t为采样间隔,△t=1/f<sub>s</sub>,f<sub>s</sub>为采样频率;p=‑2arccot(2πμ)/π;步骤二、开启干扰机的接收模块,对截获到的雷达信号x(n)进行间歇采样,采样波门函数为p(n),其中:<maths num="0002"><math><![CDATA[<mrow><mi>p</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>r</mi><mo>=</mo><mn>0</mn></mrow><mrow><mo>+</mo><mi>&infin;</mi></mrow></munderover><munderover><mo>&Sigma;</mo><mrow><msup><mi>k</mi><mo>&prime;</mo></msup><mo>=</mo><mn>0</mn></mrow><msub><mi>N</mi><mn>1</mn></msub></munderover><mi>&delta;</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><msup><mi>k</mi><mo>&prime;</mo></msup><mo>+</mo><msub><mi>rN</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000984984420000014.GIF" wi="1429" he="158" /></maths>采样完成后关闭接收模块;采样得到的信号为s(n)<maths num="0003"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mi>s</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mi>x</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>&times;</mo><mi>p</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mi>x</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>&times;</mo><mo>&lsqb;</mo><munderover><mo>&Sigma;</mo><mrow><mi>r</mi><mo>=</mo><mn>0</mn></mrow><mrow><mo>+</mo><mi>&infin;</mi></mrow></munderover><munderover><mo>&Sigma;</mo><mrow><msup><mi>k</mi><mo>&prime;</mo></msup><mo>=</mo><mn>0</mn></mrow><msub><mi>N</mi><mn>1</mn></msub></munderover><mi>&delta;</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><msup><mi>k</mi><mo>&prime;</mo></msup><mo>+</mo><msub><mi>rN</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>r</mi><mo>=</mo><mn>0</mn></mrow><mrow><mo>+</mo><mi>&infin;</mi></mrow></munderover><munderover><mo>&Sigma;</mo><mrow><msup><mi>k</mi><mo>&prime;</mo></msup><mo>=</mo><mn>0</mn></mrow><msub><mi>N</mi><mn>1</mn></msub></munderover><mi>x</mi><mrow><mo>(</mo><msup><mi>k</mi><mo>&prime;</mo></msup><mo>-</mo><msub><mi>rN</mi><mn>2</mn></msub><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000984984420000015.GIF" wi="1710" he="350" /></maths>其中,k'为离散点数,取值为0到N<sub>1</sub>中的自然数;p(n)是宽度为N<sub>1</sub>、重复间隔为N<sub>2</sub>的周期方波序列,N2的取值范围为(1/B~T/2)f<sub>s</sub>,B为雷达脉冲信号的带宽,T为雷达脉冲信号的重复周期,周期序列占空比为<img file="FDA0000984984420000021.GIF" wi="107" he="126" />r为p(n)的循环周期数目;步骤三、用步骤一中得到的分数阶傅里叶域滤波器组{h<sub>k,p</sub>(n)}<sub>k=0,1,…K‑1</sub>对步骤二中间歇采样得到的信号s(n)进行分数阶傅里叶域滤波并对每个子信道的输出进行M倍下抽取,M为抽取因子,得到分数阶傅里叶域信道化各个子信道的输出信号<img file="FDA0000984984420000022.GIF" wi="535" he="110" /><img file="FDA0000984984420000023.GIF" wi="1838" he="335" />其中,<img file="FDA0000984984420000024.GIF" wi="47" he="79" />表示p阶分数阶傅里叶域上的卷积,<img file="FDA0000984984420000025.GIF" wi="44" he="44" />表示傅里叶域上的卷积;K=F·M,F为比例因子,F取值为1或2;通常取F=1,则K=M;步骤四、利用自相关积累算法,对步骤三中的K个通道的输出信号进行信道检测,方法为:对每个信道的输出信号<img file="FDA0000984984420000026.GIF" wi="392" he="104" />分别选取N个点做自相关积累,得到信号z<sub>k,p</sub>(m)|<sub>k=0,1...K‑1</sub><img file="FDA0000984984420000027.GIF" wi="1469" he="151" />其中,i为0至N‑1中的自然数;然后根据选定的门限TH,对K个信道进行门限判决,如果第k个信道的自相关积累z<sub>k,p</sub>(m)的模值有连续D个点大于门限TH,则判决该信道有信号输出,并记下信道号k;否则,认为该信道没有信号输出;门限TH是根据TH=μ<sub>1</sub>+ξσ<sub>1</sub>计算得到,其中ξ是由给定的恒虚警概率P<sub>fa</sub>确定的门限系数,其中<img file="FDA0000984984420000031.GIF" wi="1078" he="157" />N为自相关积累的点数,σ<sup>2</sup>为噪声功率,由雷达干扰机接收端的带宽确定;步骤五、根据步骤四确定的信道号k,对分数阶傅里叶域信道化输出的信号<img file="FDA0000984984420000032.GIF" wi="157" he="87" />进行M倍零值内插,得到信号y<sub>k,p</sub>(m):<img file="FDA0000984984420000033.GIF" wi="1494" he="262" />步骤六、对步骤五内插完成的信号y<sub>k,p</sub>(m)进行幅度调制,调制系数为C,得到干扰信号J(m):J(m)=C·y<sub>k,p</sub>(m)   (7)C大于或等于M;步骤七、将步骤六调制完成的信号J(m)经过高速数模转换器后转换成模拟干扰信号,开启干扰机发送模块,发送模块将干扰信号通过发射天线转发出去;步骤八、转发完成后关闭发送模块并同时开启接收模块,转到步骤二,对下一个雷达发射脉冲信号进行后续步骤处理,直到不再接收到雷达脉冲信号。
地址 100081 北京市海淀区中关村南大街5号