发明名称 一种海洋杂波环境中电大尺寸目标雷达成像方法
摘要 本发明涉及一种海洋杂波环境中电大尺寸目标雷达成像方法,具体为:本发明的内容是通过建立多干扰环境中面—体复合目标散射模型,并采用数值与解析混合算法的电磁仿真计算和快速后向投影成像方式对海洋环境中的电大尺寸目标进行雷达成像。本发明中电大尺寸是指目标的物理尺寸跟波长的比值大于10。雷达工作的频段一般是在高频区域,例如L波段频率范围是1-2GHz,S波段在2-4GHz。一些军事目标,包括战斗机、隐身飞机、航空母舰等在上述波段都属于电大尺寸。本发明方法中,通过复杂目标几何建模,电磁仿真计算和快速后向投影成像算法以及Kaiser窗边缘滤波,即可实现对海洋环境中的电大尺寸目标快速有效的雷达成像。
申请公布号 CN102901965B 申请公布日期 2014.09.17
申请号 CN201210361270.6 申请日期 2012.09.26
申请人 同济大学 发明人 吴霞
分类号 G01S13/90(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 上海正旦专利代理有限公司 31200 代理人 张磊
主权项 一种海洋杂波环境中电大尺寸目标雷达成像方法,其特征在于具体步骤如下:(1),在仿真软件主界面菜单中选取目标模型,在选取的目标模型基础上建立海洋上空电大尺寸目标几何模型,采用GMSH方法对目标几何模型中的目标进行初始网格剖分,设定剖分类型、面元剖分尺寸以及计算频率参数;(2),根据步骤(1)得到的初始网格剖分设定的剖分类型、面元剖分尺寸以及计算频率参数通过海洋谱建模生成三维动态起伏海洋表面,采用多层快速多极子和物理光学数值混合算法进行电磁仿真计算,同时设置各种雷达工作参数和海洋表面生成,得到散射场计算结果,所述雷达工作参数包括雷达距离分辨率,入射角、方位角、场景尺寸以及方位参数;散射场计算过程为:设入射角为θi 、 φi,散射角为 θs 、φs,引入Fresnel反射系数<i>R<sub>h</sub></i>与<i>R<sub>v</sub></i>分别表示平行于入射面与垂直于入射面,因此任意极化方向的入射电场表示为平行极化与垂直极化的叠加;散射场可近似表达为:<img file="905004dest_path_image001.GIF" wi="517" he="45" />(1)其中:<img file="113262dest_path_image002.GIF" wi="202" he="28" />(2a)<img file="301536dest_path_image003.GIF" wi="169" he="26" />(2b)<img file="186315dest_path_image004.GIF" wi="176" he="28" />(2c)<img file="266398dest_path_image005.GIF" wi="218" he="28" />(2d)<img file="527615dest_path_image006.GIF" wi="200" he="57" />(2e)<img file="570395dest_path_image007.GIF" wi="202" he="57" />(2f)其中:<img file="376808dest_path_image008.GIF" wi="15" he="20" />是面元法向矢量,<img file="193455dest_path_image009.GIF" wi="18" he="25" />和<img file="298508dest_path_image010.GIF" wi="17" he="25" />为面元局部入射坐标系,<img file="634943dest_path_image011.GIF" wi="15" he="25" />为入射场极化矢量,分为水平极化入射<img file="127104dest_path_image012.GIF" wi="30" he="26" />和垂直极化入射<img file="742631dest_path_image013.GIF" wi="34" he="26" />, 并且(<img file="830804dest_path_image012.GIF" wi="30" he="26" />,<img file="848176dest_path_image013.GIF" wi="34" he="26" />,<img file="58709dest_path_image014.GIF" wi="46" he="39" />)构成一组正交坐标系;<img file="849947dest_path_image015.GIF" wi="37" he="32" />和<img file="240346dest_path_image016.GIF" wi="41" he="34" />分别是相对介电常数和相对磁导率;所述海洋谱建模表达式为:<img file="613690dest_path_image017.GIF" wi="291" he="48" />(3)其中,<i>k</i> 是电磁波数,<i>g</i>是重力加速度,<i>U_wind</i> 是海面风速,<img file="765097dest_path_image018.GIF" wi="253" he="56" />, s是常数,取 1.5×10<sup>‑4</sup> m<sup>2</sup>,<b><i>φ</i></b>为海洋表面风向角;海面风速<i>U_wind</i> 的计算公式如下:<img file="778053dest_path_image019.GIF" wi="171" he="52" />(4)其中,<img file="473608dest_path_image020.GIF" wi="287" he="46" />,<img file="199993dest_path_image021.GIF" wi="26" he="20" />是摩擦速度,z是海面高度;(3),将步骤(2)得到的散射计算结果采用快速后向投影成像算法结合Kaiser窗函数进行边缘滤波选取合适的参数做目标特征提取,建立目标散射特性与几何结构特征的对应关系, 为目标识别提供图像数据验证,经过极化处理,获得全极化雷达成像结果;根据快速后向投影成像算法分析不同入射、雷达的频率参数下的目标特征;成像借助信号处理,从回波信号中还原出目标函数:<img file="17907dest_path_image022.GIF" wi="193" he="41" />(5)<img file="518159dest_path_image023.GIF" wi="197" he="26" />(6)<img file="312677dest_path_image024.GIF" wi="241" he="47" />(7)其中,<img file="660613dest_path_image025.GIF" wi="28" he="20" />是发射信号, <i>s</i>(<i>t</i>,<i>u</i>) 是t时刻、u处的回波信号,<img file="85647dest_path_image026.GIF" wi="16" he="16" />代表卷积,星号代表复共轭;对成像区域的每一个像素<img file="886244dest_path_image027.GIF" wi="52" he="28" />,计算出在对应的一个合成孔径时间内,它与每个位置之间的时延<img file="172869dest_path_image028.GIF" wi="42" he="28" />后,在距离压缩后的数据矩阵中找出相应的累积曲线,并沿此曲线将所有信号进行相干迭加,得到这个像素的能量和, 即为该点的图像。
地址 200092 上海市杨浦区四平路1239号