发明名称 一种基于集中式MIMO雷达搜索模式下的资源管理方法
摘要 本发明属于通信雷达技术领域,特别涉及基于新体制MIMO雷达搜索模式下的资源管理方法。本发明通过对影响MIMO雷达搜索性能的各个参数进行分析,提出了一种基于集中式MIMO雷达搜索模式下的资源管理方法,即选取子阵划分个数、搜索帧周期、信号占空比和波束驻留时间作为可控参数,在满足雷达搜索性能要求的前提下,最小化雷达资源消耗量。其中雷达搜索性能不仅考虑了目标径向的累积检测概率,并首次考虑了目标切向的累积检测概率。最后采用遗传算法求解上述优化问题,得到具体场景中搜索参数配置结果。
申请公布号 CN105182317A 申请公布日期 2015.12.23
申请号 CN201510512693.7 申请日期 2015.08.20
申请人 电子科技大学 发明人 程婷;武俊青;杨少委;张洁;张宇轩
分类号 G01S7/523(2006.01)I;H04L7/06(2006.01)I;H04B7/08(2006.01)I 主分类号 G01S7/523(2006.01)I
代理机构 电子科技大学专利中心 51203 代理人 张杨
主权项 一种基于集中式MIMO雷达搜索模式下的资源管理方法,包括以下步骤:步骤1、建立MIMO信号处理模型:针对收发阵列共址的MIMO雷达,将阵列分为K个子阵,各个子阵中包含L个阵元,则阵列共包含阵元个数M=KL,推导出雷达接收信号的信噪比为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>S</mi><mi>N</mi><mi>R</mi><mo>=</mo><mfrac><mrow><msub><mi>G</mi><mi>t</mi></msub><msub><mi>G</mi><mi>r</mi></msub><msub><mi>Mp</mi><mi>t</mi></msub><msub><mi>&eta;t</mi><mi>B</mi></msub><msup><mi>&sigma;&lambda;</mi><mn>2</mn></msup></mrow><mrow><msup><mrow><mo>(</mo><mn>4</mn><mi>&pi;</mi><mo>)</mo></mrow><mn>3</mn></msup><msub><mi>N</mi><mn>0</mn></msub><msup><mi>R</mi><mn>4</mn></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000785083090000011.GIF" wi="805" he="209" /></maths>其中雷达单个天线的峰值功率为p<sub>t</sub>,发射天线增益为G<sub>t</sub>,天线接收增益为G<sub>r</sub>,目标的雷达截面积(RCS)为σ,λ为信号波长,N<sub>0</sub>为噪声功率谱密度,R为目标到雷达的径向距离,η为信号的占空比、t<sub>B</sub>为波束驻留时间;步骤2、计算切向累积检测概率:令MIMO雷达的搜索范围为[θ<sub>1</sub>,θ<sub>2</sub>],依据跟踪起始距离的概念,定义跟踪起始角度θ<sub>max</sub>,θ<sub>max</sub>=(θ<sub>2</sub>‑θ<sub>1</sub>)Q,Q∈(0,1)  (2)即对于给定目标,当雷达的切向累积检测概率达到累计检测概率门限时目标飞过的角度;角度上的飞过θ<sub>max</sub>对应的时间段为[t<sub>0</sub>,t<sub>max</sub>],其中,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>t</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>=</mo><msub><mi>t</mi><mn>0</mn></msub><mo>+</mo><mfrac><mrow><msub><mi>&theta;</mi><mi>max</mi></msub><msub><mi>R</mi><mi>x</mi></msub></mrow><msub><mi>v</mi><mi>x</mi></msub></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000785083090000012.GIF" wi="492" he="182" /></maths>切向积累检测概率的计算步骤如下:A、计算[t<sub>0</sub>,t<sub>max</sub>]时间段内,雷达发射搜索波束个数:目标沿切向飞入搜索区域到达到累积检测概率所用时间为[t<sub>0</sub>,t<sub>max</sub>],此期间雷达总共发射搜索波束个数I=I<sub>1</sub>+I<sub>2</sub>,其中,<img file="FDA0000785083090000021.GIF" wi="1153" he="228" /><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>I</mi><mn>2</mn></msub><mo>=</mo><mi>m</mi><mi>i</mi><mi>n</mi><mo>&lsqb;</mo><mrow><mfrac><mrow><mi>mod</mi><mrow><mo>(</mo><mfrac><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><msub><mi>T</mi><mi>f</mi></msub></mfrac><mo>)</mo></mrow><mo>*</mo><msub><mi>T</mi><mi>f</mi></msub></mrow><msub><mi>t</mi><mi>B</mi></msub></mfrac><mo>,</mo><msub><mi>N</mi><mi>B</mi></msub></mrow><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000785083090000022.GIF" wi="1184" he="451" /></maths>t<sub>0</sub>时刻波束计数器值n<sub>count</sub>,n<sub>count</sub>∈[0,I‑1]开始计数,直到n<sub>count</sub>=I‑1循环以下计算;B、计算每个搜索波束在目标上的驻留时间以及单次检测概率,目标飞入搜索区域t<sub>0</sub>时刻,角度位置为θ<sub>1</sub>,此时搜索波束编号为n<sub>0</sub>,<img file="FDA0000785083090000023.GIF" wi="799" he="413" />在发射第n<sub>count</sub>个雷达搜索波束时,波束指向范围为[θ<sub>s1</sub>(t<sub>0</sub>),θ<sub>s2</sub>(t<sub>0</sub>)]=[θ<sub>1</sub>+(n<sub>search</sub>(t<sub>0</sub>)‑1)B<sub>w</sub>,θ<sub>1</sub>+n<sub>search</sub>(t<sub>0</sub>)B<sub>w</sub>]  (6)对应的搜索波束执行时间为<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mo>&lsqb;</mo><msub><mi>t</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>t</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>&rsqb;</mo><mo>=</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>&lsqb;</mo><mrow><mo>(</mo><mrow><msub><mi>n</mi><msub><mi>T</mi><mi>f</mi></msub></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><msub><mi>T</mi><mi>f</mi></msub><mo>+</mo><mrow><mo>(</mo><mrow><msub><mi>n</mi><mrow><mi>s</mi><mi>e</mi><mi>a</mi><mi>r</mi><mi>c</mi><mi>h</mi></mrow></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><msub><mi>t</mi><mi>B</mi></msub><mo>,</mo><mrow><mo>(</mo><mrow><msub><mi>n</mi><msub><mi>T</mi><mi>f</mi></msub></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><msub><mi>T</mi><mi>f</mi></msub><mo>+</mo><msub><mi>n</mi><mrow><mi>s</mi><mi>e</mi><mi>a</mi><mi>r</mi><mi>c</mi><mi>h</mi></mrow></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><msub><mi>t</mi><mi>B</mi></msub><mo>&rsqb;</mo></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000785083090000024.GIF" wi="1808" he="269" /></maths>其中,搜索帧周期计数器<img file="FDA0000785083090000025.GIF" wi="90" he="87" />为<img file="FDA0000785083090000026.GIF" wi="263" he="109" /><img file="FDA0000785083090000031.GIF" wi="810" he="209" />搜索波位编号n<sub>search</sub>为(n<sub>search</sub>≥1)<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>n</mi><mrow><mi>s</mi><mi>e</mi><mi>a</mi><mi>r</mi><mi>c</mi><mi>h</mi></mrow></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>mod</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>n</mi><mn>0</mn></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>n</mi><mrow><mi>c</mi><mi>o</mi><mi>u</mi><mi>n</mi><mi>t</mi></mrow></msub></mrow><msub><mi>N</mi><mi>B</mi></msub></mfrac><mo>)</mo></mrow><mo>,</mo><msub><mi>n</mi><mn>0</mn></msub><mo>+</mo><msub><mi>n</mi><mrow><mi>c</mi><mi>o</mi><mi>u</mi><mi>n</mi><mi>t</mi></mrow></msub><mo>&NotEqual;</mo><msub><mi>kN</mi><mi>B</mi></msub><mo>,</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>...</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>N</mi><mi>B</mi></msub><mo>,</mo><msub><mi>n</mi><mn>0</mn></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>n</mi><mrow><mi>c</mi><mi>o</mi><mi>u</mi><mi>n</mi><mi>t</mi></mrow></msub><mo>=</mo><msub><mi>kN</mi><mi>B</mi></msub><mo>,</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>..</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000785083090000032.GIF" wi="1919" he="330" /></maths>运动目标所在角度与时间的函数关系式为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>&theta;</mi><mi>t</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>v</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow><msub><mi>r</mi><mi>x</mi></msub></mfrac><mo>+</mo><msub><mi>&theta;</mi><mn>1</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000785083090000033.GIF" wi="785" he="194" /></maths>然后,求解该搜索波束在目标上的驻留时间;令<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>&theta;</mi><msub><mi>s</mi><mn>1</mn></msub></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>v</mi><mrow><mo>(</mo><msub><mover><mi>t</mi><mo>~</mo></mover><mn>1</mn></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow><mi>r</mi></mfrac><mo>+</mo><msub><mi>&theta;</mi><mn>1</mn></msub><mo>,</mo><msub><mi>&theta;</mi><msub><mi>s</mi><mn>2</mn></msub></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>v</mi><mrow><mo>(</mo><msub><mover><mi>t</mi><mo>~</mo></mover><mn>2</mn></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow><mi>r</mi></mfrac><mo>+</mo><msub><mi>&theta;</mi><mn>1</mn></msub><mo>,</mo></mrow>]]></math><img file="FDA0000785083090000034.GIF" wi="1261" he="179" /></maths>求解可以分别求得<img file="FDA0000785083090000035.GIF" wi="57" he="101" />和<img file="FDA0000785083090000036.GIF" wi="93" he="101" />判断第n<sub>count</sub>个波束的照射时间<img file="FDA00007850830900000312.GIF" wi="213" he="118" />和<img file="FDA0000785083090000037.GIF" wi="176" he="107" />是否存在交集,若存在交集则代表目标被该波束照射到,交集的长度t<sub>B</sub>′则代表被照射的时间长度,若t<sub>B</sub>′≠0,P<sub>d</sub>(t<sub>B</sub>′)=0。其中,<img file="FDA0000785083090000038.GIF" wi="386" he="158" />C、将I个搜索波束的目标检测概率进行非相干积累,得到目标切向累计检测概率,<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><munderover><mo>&Pi;</mo><mrow><msub><mi>n</mi><mrow><mi>c</mi><mi>o</mi><mi>u</mi><mi>n</mi><mi>t</mi></mrow></msub><mo>=</mo><mn>1</mn></mrow><mi>I</mi></munderover><mo>&lsqb;</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>d</mi><mi>c</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>P</mi><mi>d</mi></msub><mrow><mo>(</mo><mrow><msup><msub><mi>t</mi><mi>B</mi></msub><mo>&prime;</mo></msup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>n</mi><mrow><mi>c</mi><mi>o</mi><mi>u</mi><mi>n</mi><mi>t</mi></mrow></msub></mrow><mo>)</mo></mrow><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000785083090000039.GIF" wi="1046" he="190" /></maths>D、令t<sub>0</sub>在[0,T<sub>f</sub>]内均匀分布,于是,以<img file="FDA00007850830900000310.GIF" wi="400" he="168" />作为跟踪起始时刻(t<sub>max</sub>&gt;T<sub>f</sub>),雷达对目标的平均切向累计检测概率为:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>P</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><msub><mi>T</mi><mi>f</mi></msub></mfrac><msubsup><mo>&Integral;</mo><mn>0</mn><msub><mi>T</mi><mi>f</mi></msub></msubsup><mrow><mrow><mo>{</mo><mrow><mn>1</mn><mo>-</mo><munderover><mo>&Pi;</mo><mrow><msub><mi>n</mi><mrow><mi>c</mi><mi>o</mi><mi>u</mi><mi>n</mi><mi>t</mi></mrow></msub><mo>=</mo><mn>1</mn></mrow><mi>I</mi></munderover><mo>&lsqb;</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>d</mi><mi>c</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>P</mi><mi>d</mi></msub><mrow><mo>(</mo><mrow><msup><msub><mi>t</mi><mi>B</mi></msub><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>n</mi><mrow><mi>c</mi><mi>o</mi><mi>u</mi><mi>n</mi><mi>t</mi></mrow></msub></mrow><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mo>}</mo></mrow><mi>d</mi><mi>t</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow><mo>;</mo></mrow>]]></math><img file="FDA00007850830900000311.GIF" wi="1656" he="220" /></maths>步骤3、构建MIMO雷达搜索模式下的优化模型:一帧内搜索波束所用时间与搜索帧周期的比值描述了系统在时域上的资源分布,信号占空比则描述了每次搜索能量消耗,子阵划分个数则描述了雷达硬件资源的消耗,三者按照各自权重相加作为目标函数来描述MIMO雷达资源消耗,同时,MIMO雷达要保持搜索性能,则以目标径向累积检测概率,目标切向累积检测概率大于累计检测概率门限值,时间资源消耗合理性为约束条件,优化模型为<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><munder><mi>min</mi><mrow><mo>(</mo><mrow><mi>K</mi><mo>,</mo><mi>&eta;</mi><mo>,</mo><msub><mi>t</mi><mi>B</mi></msub><mo>,</mo><msub><mi>T</mi><mi>f</mi></msub></mrow><mo>)</mo></mrow></munder><msub><mi>c</mi><mn>1</mn></msub><mfrac><mrow><msub><mi>t</mi><mi>B</mi></msub><msub><mi>N</mi><mi>B</mi></msub></mrow><msub><mi>T</mi><mi>f</mi></msub></mfrac><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><mi>&eta;</mi><mo>+</mo><msub><mi>c</mi><mn>3</mn></msub><mfrac><mi>K</mi><mi>M</mi></mfrac></mrow>]]></math><img file="FDA0000785083090000041.GIF" wi="759" he="187" /></maths><maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>P</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>R</mi><mi>t</mi></msub><mo>)</mo></mrow><mo>&GreaterEqual;</mo><msub><mi>P</mi><mi>D</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>P</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>max</mi></msub><mo>)</mo></mrow><mo>&GreaterEqual;</mo><msub><mi>P</mi><mi>D</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mrow><msub><mi>N</mi><mi>B</mi></msub><msub><mi>t</mi><mi>B</mi></msub></mrow><msub><mi>T</mi><mi>f</mi></msub></mfrac><mo>&le;</mo><mn>1</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000785083090000042.GIF" wi="787" he="516" /></maths>步骤4、采用遗传算法求解上述优化问题:通过求解该优化模型的优化过程,得到一组参数[K<sub>opt</sub>,η<sub>opt</sub>,t<sub>Bopt</sub>,T<sub>fopt</sub>],在满足搜索性能要求的前提下,消耗资源最少。
地址 611731 四川省成都市高新区(西区)西源大道2006号