发明名称 警戒雷达抗同型干扰方法
摘要 本发明公布了一种基于克隆选择的警戒雷达抗同型干扰方法。该方法基于克隆选择框架构造启发式搜索,求解一个有约束的整数二次规划问题,从任意的初始解出发,通过克隆、变异、选择逐代进化,直至达到最大进化代数或最优解无改善,得到各部雷达选用各套参数配置的最优频度,设计最优工作序列表,实现抗同型干扰同时不降低每部雷达性能的目标。具有良好的普适性、稳健性和较高的效率,是一种实际可行、可靠、高效的方法,并且程序简单,容易在嵌入式数字处理平台上实现,具有较高的实用价值。
申请公布号 CN101839975A 申请公布日期 2010.09.22
申请号 CN201010157716.4 申请日期 2010.04.21
申请人 清华大学 发明人 陈思;张卫杰;杨健;周广益
分类号 G01S7/36(2006.01)I 主分类号 G01S7/36(2006.01)I
代理机构 北京众合诚成知识产权代理有限公司 11246 代理人 黄家俊
主权项 1.一种基于克隆选择的警戒雷达抗同型干扰方法,其特征在于,包括以下步骤:步骤1:将归一化的平均同型干扰式(1)作为优化目标式,根据实际情况确定各个参量,<maths num="0001"><![CDATA[<math><mrow><mover><mi>I</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mrow><mi>P</mi><mrow><mo>(</mo><mi>P</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><msub><mi>p</mi><mn>1</mn></msub><mo>=</mo><mn>1</mn></mrow><mi>P</mi></munderover><mrow><mmultiscripts><munderover><munder><mi>&Sigma;</mi><mrow><msub><mi>p</mi><mn>2</mn></msub><mo>=</mo><mn>1</mn></mrow></munder><mrow><msub><mi>p</mi><mn>2</mn></msub><mo>&NotEqual;</mo><msub><mi>p</mi><mn>1</mn></msub></mrow><mi>P</mi></munderover></mmultiscripts><mfrac><mn>1</mn><msup><mi>N</mi><mn>2</mn></msup></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><msub><mi>p</mi><mn>1</mn></msub></msub></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><msub><mi>p</mi><mn>2</mn></msub></msub></munderover><msubsup><mi>&sigma;</mi><mi>ij</mi><mrow><mrow><mo>(</mo><msub><mi>p</mi><mn>1</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>p</mi><mn>2</mn></msub><mo>)</mo></mrow></mrow></msubsup></mrow><msubsup><mi>n</mi><mi>i</mi><mrow><mo>(</mo><msub><mi>p</mi><mn>1</mn></msub><mo>)</mo></mrow></msubsup><msubsup><mi>n</mi><mi>j</mi><mrow><mo>(</mo><msub><mi>p</mi><mn>2</mn></msub><mo>)</mo></mrow></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,P为系统中雷达数量,各部雷达可选参数配置数目为M<sub>p</sub>(p=1,2,...,P),N为频度量化数,<img file="FSA00000113896700012.GIF" wi="230" he="104" />为归一化同型干扰系数,该系数需经过实际测试确定并归一化至[0,1]区间;步骤2:列出解<img file="FSA00000113896700013.GIF" wi="95" he="83" />(p=1,2,...,P,i=1,2,...,M<sub>p</sub>)应满足的约束条件,至少包括<img file="FSA00000113896700014.GIF" wi="92" he="81" />为非负整数以及<maths num="0002"><![CDATA[<math><mrow><mmultiscripts><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mi>p</mi></msub></munderover></mmultiscripts><msubsup><mi>n</mi><mi>i</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mi>N</mi><mo>,</mo><mo>&ForAll;</mo><mi>p</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>P</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>步骤3:任意给定一组满足所述约束条件的初始解,作为第1代抗体群,记为<img file="FSA00000113896700016.GIF" wi="64" he="62" />(n=1,2,...,N<sub>A</sub>),常数N<sub>A</sub>为抗体群规模;步骤4:用式(1)计算所述第1代抗体群中每个抗体对应的平均归一化同型干扰<img file="FSA00000113896700017.GIF" wi="343" he="68" />计算抗体适应度<maths num="0003"><![CDATA[<math><mrow><mi>f</mi><mrow><mo>(</mo><msubsup><mi>A</mi><mi>n</mi><mn>1</mn></msubsup><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>-</mo><mover><mi>I</mi><mo>-</mo></mover><mrow><mo>(</mo><msubsup><mi>A</mi><mi>n</mi><mn>1</mn></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>并记录其最大值;步骤5:对于第k代抗体群<img file="FSA00000113896700019.GIF" wi="66" he="65" />中的每个抗体,用式(4)计算其克隆数<img file="FSA00000113896700021.GIF" wi="1441" he="272" />并构造克隆抗体群<img file="FSA00000113896700022.GIF" wi="1050" he="66" />其中,常数N<sub>C</sub>为克隆抗体群规模;步骤6:对克隆抗体群中的每个抗体<maths num="0004"><![CDATA[<math><mo>[</mo><mrow><msubsup><mi>n</mi><mn>1</mn><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>n</mi><msub><mi>M</mi><mn>1</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>n</mi><mn>1</mn><mrow><mo>(</mo><mi>P</mi><mo>)</mo></mrow></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>n</mi><msub><mi>M</mi><mi>P</mi></msub><mrow><mo>(</mo><mi>P</mi><mo>)</mo></mrow></msubsup><mo>]</mo><mo>---</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>随机选取p<sub>0</sub>∈{1,2,...,P},并随机选取<img file="FSA00000113896700024.GIF" wi="479" he="109" />构造变异抗体<maths num="0005"><![CDATA[<math><mrow><mo>[</mo><msubsup><mi>n</mi><mn>1</mn><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>n</mi><msub><mi>M</mi><mn>1</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>n</mi><mn>1</mn><mrow><mo>(</mo><msub><mi>p</mi><mn>0</mn></msub><mo>)</mo></mrow></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mover><mi>n</mi><mo>^</mo></mover><mi>i</mi><mrow><mo>(</mo><msub><mi>p</mi><mn>0</mn></msub><mo>)</mo></mrow></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mover><mi>n</mi><mo>^</mo></mover><mi>j</mi><mrow><mo>(</mo><msub><mi>p</mi><mn>0</mn></msub><mo>)</mo></mrow></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>n</mi><msub><mi>M</mi><msub><mi>p</mi><mn>0</mn></msub></msub><mrow><mo>(</mo><msub><mi>p</mi><mn>0</mn></msub><mo>)</mo></mrow></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>n</mi><mn>1</mn><mrow><mo>(</mo><mi>P</mi><mo>)</mo></mrow></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>n</mi><msub><mi>M</mi><mi>P</mi></msub><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow></msubsup><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>满足<maths num="0006"><![CDATA[<math><mrow><msubsup><mover><mi>n</mi><mo>^</mo></mover><mi>i</mi><mrow><mo>(</mo><msub><mi>p</mi><mn>0</mn></msub><mo>)</mo></mrow></msubsup><mo>+</mo><msubsup><mover><mi>n</mi><mo>^</mo></mover><mi>j</mi><mrow><mo>(</mo><msub><mi>p</mi><mn>0</mn></msub><mo>)</mo></mrow></msubsup><mo>=</mo><msubsup><mi>n</mi><mi>i</mi><mrow><mo>(</mo><msub><mi>p</mi><mn>0</mn></msub><mo>)</mo></mrow></msubsup><mo>+</mo><msubsup><mi>n</mi><mi>j</mi><mrow><mo>(</mo><msub><mi>p</mi><mn>0</mn></msub><mo>)</mo></mrow></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>在和不变的条件下两个分量的重新分配也是随机的;步骤7:计算所有所述变异抗体的适应度,从所有由同一个抗体<img file="FSA00000113896700027.GIF" wi="62" he="64" />经过克隆和变异得到的变异抗体中选出适应度最大的,记为<img file="FSA00000113896700028.GIF" wi="83" he="64" />步骤8:依概率从<img file="FSA00000113896700029.GIF" wi="65" he="65" />知<img file="FSA000001138967000210.GIF" wi="64" he="65" />二者中选择其一进入下一代,选择<img file="FSA000001138967000211.GIF" wi="63" he="64" />的概率为<maths num="0007"><![CDATA[<math><mrow><mi>Pr</mi><mrow><mo>(</mo><msubsup><mi>A</mi><mi>n</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mi>B</mi><mi>n</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>=</mo></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn><mo>,</mo></mtd><mtd><mi>f</mi><mrow><mo>(</mo><msubsup><mi>A</mi><mi>n</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>&lt;</mo><mi>f</mi><mrow><mo>(</mo><msubsup><mi>B</mi><mi>n</mi><mi>k</mi></msubsup><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mrow><mi>f</mi><mrow><mo>(</mo><msubsup><mi>A</mi><mi>n</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mi>f</mi><mrow><mo>(</mo><msubsup><mi>B</mi><mi>n</mi><mi>k</mi></msubsup><mo>)</mo></mrow></mrow><mi>&alpha;</mi></mfrac><mo>)</mo></mrow><mo>,</mo></mtd><mtd><mi>f</mi><mrow><mo>(</mo><msubsup><mi>A</mi><mi>n</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>&GreaterEqual;</mo><mi>f</mi><mrow><mo>(</mo><msubsup><mi>B</mi><mi>n</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>,</mo><mi>f</mi><mrow><mo>(</mo><msubsup><mi>A</mi><mi>n</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>&NotEqual;</mo><munder><mi>max</mi><mi>m</mi></munder><mi>f</mi><mrow><mo>(</mo><msubsup><mi>A</mi><mi>m</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>,</mo><mi>m</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>N</mi><mi>A</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn><mo>,</mo></mtd><mtd><mi>f</mi><mrow><mo>(</mo><msubsup><mi>A</mi><mi>n</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>&GreaterEqual;</mo><mrow><mo>(</mo><msubsup><mi>B</mi><mi>n</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>,</mo><mi>f</mi><mrow><mo>(</mo><msubsup><mi>A</mi><mi>n</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>=</mo><munder><mi>max</mi><mi>m</mi></munder><mi>f</mi><mrow><mo>(</mo><msubsup><mi>A</mi><mi>m</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>,</mo><mi>m</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>N</mi><mi>A</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,常数α>0控制选择的随机性;步骤9:记录第k+1代抗体群<img file="FSA000001138967000214.GIF" wi="114" he="67" />(n=1,2,...,N<sub>A</sub>)的适应度最大值;步骤10:k:=k+1,若k=K<sub>max</sub>,或连续K<sub>stop</sub>代抗体群适应度最大值没有改善,则终止计算,当前代抗体群中适应度最大的抗体即为求得的最优解;否则返回步骤5继续进化。
地址 100084 北京市100084-82信箱