发明名称 提高MIMO-STAP检测性能的稳健波形优化方法
摘要 本发明公开了一种提高MIMO-STAP检测性能的稳健波形优化方法,步骤包括:(1)建立MIMO-STAP模型(2)针对空时频导向矢量存在估计误差的场景,建立实际空时频导向矢量模型,据此表述稳健波形优化问题;(2)对<img file="2014100397853100004DEST_PATH_IMAGE001.GIF" wi="17" he="15" />应用对角加载方法,将稳健波形优化问题松弛为SDP问题,从而可以利用诸多成熟优化工具箱获得高效求解。本发明所提出的稳健波形优化方法考虑了求解波形优化问题所需参数存在估计误差时的稳健设计问题,更贴近于实际场景,因而更具有较高的实用价值,并且利用对角加载方法将此稳健优化问题松弛为SDP问题,从而可以获得高效求解。
申请公布号 CN103852749A 申请公布日期 2014.06.11
申请号 CN201410039785.3 申请日期 2014.01.28
申请人 大连大学 发明人 王洪雁;裴炳南;汪祖民;刘艳;裴腾达
分类号 G01S7/02(2006.01)I 主分类号 G01S7/02(2006.01)I
代理机构 大连八方知识产权代理有限公司 21226 代理人 任洪成
主权项 1.提高MIMO-STAP检测性能的稳健波形优化方法,包括如下步骤:(1)表述输出SINR1a) MIMO-STAP接收信号描述MIMO雷达的发射阵列和接收阵列均为均匀线阵,且平行放置,阵元数分别为<img file="348620DEST_PATH_IMAGE001.GIF" wi="22" he="18" />和<img file="423237DEST_PATH_IMAGE002.GIF" wi="20" he="20" />,阵元间距分别为<img file="157975DEST_PATH_IMAGE003.GIF" wi="21" he="25" />和<img file="953761DEST_PATH_IMAGE004.GIF" wi="22" he="25" />,雷达平台沿平行于发射、接收阵列的方向匀速直线飞行,飞行高度和速度分别为<img file="64937DEST_PATH_IMAGE005.GIF" wi="14" he="20" />和<img file="262569DEST_PATH_IMAGE006.GIF" wi="13" he="16" />,目标沿与发射、接收阵列法线夹角为<img file="168208DEST_PATH_IMAGE007.GIF" wi="17" he="25" />的直线匀速运动,速度为<img file="388974DEST_PATH_IMAGE008.GIF" wi="16" he="25" />,且与雷达平台处于同一平面,在一个相干处理间隔(CPI)内,各发射阵元同时辐射由<img file="100578DEST_PATH_IMAGE009.GIF" wi="16" he="18" />个脉冲组成的脉冲串波形,且脉冲重复间隔(PRI)为<img file="155646DEST_PATH_IMAGE010.GIF" wi="16" he="18" />,将距离环离散化为<img file="232186DEST_PATH_IMAGE011.GIF" wi="116" he="24" />个小单元,则第<img file="2565DEST_PATH_IMAGE012.GIF" wi="9" he="20" />个PRI内的接收数据可表示为:<img file="189964DEST_PATH_IMAGE013.GIF" wi="286" he="42" />式中,<img file="362188DEST_PATH_IMAGE014.GIF" wi="194" he="34" />和<img file="609630DEST_PATH_IMAGE015.GIF" wi="208" he="36" />分别表示目标以及位于<img file="867305DEST_PATH_IMAGE016.GIF" wi="17" he="22" />的杂波发射导向矢量;<img file="858394DEST_PATH_IMAGE017.GIF" wi="182" he="34" />和<img file="882195DEST_PATH_IMAGE018.GIF" wi="196" he="36" />分别表示目标以及位于<img file="362855DEST_PATH_IMAGE016.GIF" wi="17" he="22" />的杂波的接收导向矢量;<img file="592979DEST_PATH_IMAGE019.GIF" wi="89" he="40" />,<img file="574711DEST_PATH_IMAGE020.GIF" wi="137" he="40" />,<img file="455948DEST_PATH_IMAGE021.GIF" wi="94" he="40" />,<img file="45192DEST_PATH_IMAGE022.GIF" wi="49" he="44" />,以及<img file="277459DEST_PATH_IMAGE023.GIF" wi="61" he="44" />;<img file="610352DEST_PATH_IMAGE024.GIF" wi="169" he="25" />表示每个PRI内发射波形矩阵;<img file="349025DEST_PATH_IMAGE025.GIF" wi="64" he="24" />为第<img file="374750DEST_PATH_IMAGE026.GIF" wi="18" he="16" />个发射阵元在每个PRI内发射的复基带信号的离散形式,<img file="641783DEST_PATH_IMAGE027.GIF" wi="18" he="18" />为波形采样数;<img file="27634DEST_PATH_IMAGE028.GIF" wi="20" he="22" />和<img file="368617DEST_PATH_IMAGE029.GIF" wi="40" he="22" />分别为所考虑的距离环内目标的复幅度以及位于<img file="814510DEST_PATH_IMAGE030.GIF" wi="17" he="25" />的杂波反射系数;<img file="506523DEST_PATH_IMAGE031.GIF" wi="16" he="20" />为波形中心波长;<img file="430485DEST_PATH_IMAGE032.GIF" wi="201" he="33" />,<img file="891554DEST_PATH_IMAGE033.GIF" wi="70" he="26" />表示第<img file="528857DEST_PATH_IMAGE034.GIF" wi="14" he="16" />个接收阵元在第<img file="442586DEST_PATH_IMAGE012.GIF" wi="9" he="20" />个PRI内接收的干扰以及噪声,假设<img file="107923DEST_PATH_IMAGE035.GIF" wi="21" he="24" />的每一列独立同分布,且服从于均值为0,方差为<img file="751393DEST_PATH_IMAGE036.GIF" wi="18" he="21" />的循环对称复高斯分布;1b)感兴趣距离环内空时快拍表述为得到用于目标检测的统计量,利用<img file="273511DEST_PATH_IMAGE037.GIF" wi="82" he="29" />作为匹配滤波器,则相应矢量化的匹配输出可表示如下:<img file="127066DEST_PATH_IMAGE038.GIF" wi="504" he="42" />式中,<img file="333444DEST_PATH_IMAGE039.GIF" wi="129" he="30" />;<img file="769104DEST_PATH_IMAGE040.GIF" wi="121" he="25" />即为波形协方差矩阵(WCM),由上式可得距离环内总的空时快拍为:<img file="462123DEST_PATH_IMAGE041.GIF" wi="509" he="69" />式中,<img file="616024DEST_PATH_IMAGE042.GIF" wi="196" he="36" />和<img file="685480DEST_PATH_IMAGE043.GIF" wi="210" he="37" />分别为目标以及位于<img file="975647DEST_PATH_IMAGE016.GIF" wi="17" he="22" />的杂波的多普勒导向矢量;<img file="105146DEST_PATH_IMAGE044.GIF" wi="18" he="22" />表示元素全部为1的<img file="480763DEST_PATH_IMAGE045.GIF" wi="33" he="18" />维矢量;且<img file="85401DEST_PATH_IMAGE046.GIF" wi="74" he="25" />1c)最优MIMO-STAP处理器条件下输出SINR表述对于最优的MIMO-STAP处理器,输出SINR可表示为:<img file="495654DEST_PATH_IMAGE047.GIF" wi="517" he="32" />式中,<img file="530475DEST_PATH_IMAGE048.GIF" wi="394" he="128" />1d)杂波高斯分布,且与干扰不相关条件下输出SINR表述简化在杂波与干扰加噪声项不相关以及<img file="393389DEST_PATH_IMAGE049.GIF" wi="18" he="22" />独立同分布,且服从均值为0,方差为<img file="804647DEST_PATH_IMAGE050.GIF" wi="22" he="24" />的高斯分布条件下,输出SINR可简化为如下表达式:<img file="69407DEST_PATH_IMAGE051.GIF" wi="233" he="30" />式中,<img file="540708DEST_PATH_IMAGE052.GIF" wi="114" he="22" />;<img file="890918DEST_PATH_IMAGE053.GIF" wi="156" he="30" />;<img file="843218DEST_PATH_IMAGE054.GIF" wi="86" he="25" />;<img file="962484DEST_PATH_IMAGE055.GIF" wi="122" he="24" />;<img file="604687DEST_PATH_IMAGE056.GIF" wi="132" he="32" />;<img file="442193DEST_PATH_IMAGE057.GIF" wi="117" he="25" />;且<img file="195254DEST_PATH_IMAGE058.GIF" wi="164" he="26" />;<b> </b>(2) 建立MIMO-STAP的稳健波形优化模型                          由上式可得,SINR的计算需要以下参量的确切值:波形相关条件下的空时导向矢量<img file="231343DEST_PATH_IMAGE059.GIF" wi="18" he="22" />(即<img file="795179DEST_PATH_IMAGE007.GIF" wi="17" he="25" />,<img file="369249DEST_PATH_IMAGE060.GIF" wi="22" he="25" />),杂波协方差矩阵<img file="595176DEST_PATH_IMAGE061.GIF" wi="25" he="24" />,以及干扰加噪声的协方差矩阵<img file="875984DEST_PATH_IMAGE036.GIF" wi="18" he="21" />;2a)真实空时频导向矢量建模假设<img file="797673DEST_PATH_IMAGE059.GIF" wi="18" he="22" />不确知,但存在于一确知的凸紧支集内,即<img file="344192DEST_PATH_IMAGE062.GIF" wi="117" he="28" />式中,<img file="704635DEST_PATH_IMAGE063.GIF" wi="74" he="25" />;<img file="715316DEST_PATH_IMAGE064.GIF" wi="18" he="25" />和<img file="748519DEST_PATH_IMAGE065.GIF" wi="18" he="25" />分别为真实的以及对应假设的空时导向矢量;<img file="844651DEST_PATH_IMAGE066.GIF" wi="16" he="16" />为刻画<img file="477626DEST_PATH_IMAGE064.GIF" wi="18" he="25" />的误差的未知复矢量;而<img file="546076DEST_PATH_IMAGE067.GIF" wi="14" he="18" />则为<img file="871884DEST_PATH_IMAGE064.GIF" wi="18" he="25" />模的下界;2b)稳健波形优化问题表述基于以上讨论,提高MIMO-STAP检测性能的稳健优化波形问题可表述如下:基于发射总功率约束,在凸集<img file="392995DEST_PATH_IMAGE068.GIF" wi="18" he="20" />上优化WCM以最大化最差情况下的输出SINR,即<img file="829662DEST_PATH_IMAGE069.GIF" wi="257" he="108" />式中,<img file="752618DEST_PATH_IMAGE070.GIF" wi="17" he="18" />为总发射功率,       由上式可知,此稳健波形优化问题是关于发射波形以及空时频导向矢量误差的最大最小联合优化问题;(3) 基于对角加载求解稳健波形优化问题3a)求解内层优化问题观察上式可知,当真实空时导向矢量<img file="246398DEST_PATH_IMAGE064.GIF" wi="18" he="25" />位于对应于矩阵<img file="441756DEST_PATH_IMAGE071.GIF" wi="141" he="26" />最小特征值的特征矢量所指方向时,就会出现MIMO-STAP检测性能最差的情况,因此,2b)中的优化波形问题可重写为:<img file="432846DEST_PATH_IMAGE072.GIF" wi="236" he="82" />式中,<img file="459577DEST_PATH_IMAGE073.GIF" wi="45" he="25" />表示矩阵的最小特征值;3b)基于DL方法的正定化<img file="877920DEST_PATH_IMAGE074.GIF" wi="24" he="25" />对<img file="357311DEST_PATH_IMAGE074.GIF" wi="24" he="25" />应用对角加载方法,使得<img file="152092DEST_PATH_IMAGE075.GIF" wi="117" he="26" />式中,<img file="846379DEST_PATH_IMAGE076.GIF" wi="94" he="25" />,<img file="687820DEST_PATH_IMAGE077.GIF" wi="46" he="25" />表示矩阵的最大特征值;3c)基于正定化<img file="795453DEST_PATH_IMAGE074.GIF" wi="24" he="25" />简化稳健波形优化问题由于<img file="377613DEST_PATH_IMAGE078.GIF" wi="50" he="26" />,可得<img file="926406DEST_PATH_IMAGE079.GIF" wi="165" he="26" />,利用矩阵求逆定理,稳健波形优化问题可简化为:<img file="139082DEST_PATH_IMAGE080.GIF" wi="162" he="84" />3d)基于凸优化求解稳健波形优化问题基于上述讨论,稳健波形优化问题可等价表示为如下SDP问题:<img file="343798DEST_PATH_IMAGE081.GIF" wi="173" he="134" />式中,<img file="484578DEST_PATH_IMAGE082.GIF" wi="9" he="17" />为辅助优化变量,从而,此波形优化问题可利用诸多成熟优化工具箱获得高效求解。
地址 116622 辽宁省大连市金州新区学府大街10号