发明名称 双站合成孔径雷达海浪方向谱的仿真方法
摘要 本发明公布了一种双站合成孔径雷达海浪方向谱的仿真方法。包括以下步骤:根据海浪参数,计算JONSWAP海浪谱和海浪铅直位移;依据雷达波调制理论和BiSAR系统参数,得到含有海浪轨道速度和加速度的BiSAR回波信号;然后根据BiSAR距离-多普勒成像算法,对回波信号进行匹配滤波处理,得到BiSAR海浪图像的解析式;并采用数值积分的方法来求得BiSAR海浪图像。最后,对BiSAR海浪图像进行二维FFT运算可得到BiSAR海浪方向谱。本发明不需要在时域内产生回波信号和成像处理而是采用数值积分的方法来求得BiSAR海浪图像,具有较高的运算效率。
申请公布号 CN101697011B 申请公布日期 2011.12.07
申请号 CN200910233593.5 申请日期 2009.10.29
申请人 江苏科技大学 发明人 杨永红;张尤赛;林明;陆南
分类号 G01S13/90(2006.01)I;G01S7/41(2006.01)I;G06F17/50(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 楼高潮
主权项 1.一种双站合成孔径雷达海浪方向谱的仿真方法,其特征在于包括以下步骤:步骤一,根据海浪参数,计算JONSWAP海浪谱和海面铅直位移:JONSWAP海浪谱的表达式为<img file="FSB00000583960100011.GIF" wi="625" he="155" />其中,<img file="FSB00000583960100012.GIF" wi="484" he="161" /><img file="FSB00000583960100013.GIF" wi="439" he="194" />ω<sub>m</sub>为谱峰值角频率,ω为海浪的角频率,γ为峰升因子,δ为峰形参量,g为引力常数,α为尺度系数;根据Longuet-Higgins模型,海面铅直位移ζ为<img file="FSB00000583960100014.GIF" wi="1232" he="166" />其中,c.c.表示复共轭算子,i为虚数,t为时间,<img file="FSB00000583960100015.GIF" wi="32" he="46" />为海面空间向量,N为海浪波数的个数,j表示索引数,<img file="FSB00000583960100016.GIF" wi="65" he="70" />为第j个海浪的波数向量,φ<sub>j</sub>为第j个海浪的相位噪声,S(ω<sub>j</sub>)表示角频率为ω<sub>j</sub>的JONSWAP海浪谱的数值,Δω<sub>j</sub>=ω<sub>j+1</sub>-ω<sub>j</sub>表示为第j+1与第j海浪角频率的差;步骤二,根据双站SAR几何参数和雷达波参数,则含有轨道速度和加速度扰动的海面的回波信号可表示为:<img file="FSB00000583960100017.GIF" wi="1872" he="180" />其中,x<sub>0</sub>为海面的某一位置,σ(x<sub>0</sub>)表示海面在x<sub>0</sub>处的雷达波散射系数,<img file="FSB00000583960100018.GIF" wi="128" he="70" />为海浪对雷达波的调制系数,<img file="FSB00000583960100019.GIF" wi="162" he="77" />为发射机所产生的相位,<img file="FSB000005839601000110.GIF" wi="164" he="70" />为接收机所产生的相位;当收发天线为正侧视,则发射信号或接收信号的相位可表示为<img file="FSB000005839601000111.GIF" wi="1290" he="139" />k为雷达波的波数;下标n={1,2},n=1表示的是发射机,n=2表示的是接收机,V<sub>n</sub>为发射机/接收机的速度,R<sub>n</sub>为发射机/接收机与海面中心的最近距离,<img file="FSB000005839601000112.GIF" wi="60" he="91" />为发射机/接收机的常数相位项,U<sub>n</sub>为发射机/接收机的海浪轨道速度,A<sub>n</sub>为发射机/接收机的海浪轨道的加速度;步骤三,根据双站SAR距离-多普勒成像算法,对回波信号进行匹配滤波处理,得到双站SAR的海浪图像:当发射天线的方向图为<img file="FSB000005839601000113.GIF" wi="565" he="142" />和接收天线的方向图为<img file="FSB000005839601000114.GIF" wi="657" he="227" />其中,T表示为合成孔径时间,方位向匹配滤波器为<img file="FSB00000583960100021.GIF" wi="841" he="127" />双站SAR输出功率密度函数为r(t,x<sub>0</sub>)=G<sub>2</sub>(t,x<sub>0</sub>)G<sub>1</sub>(t,x<sub>0</sub>)e(t,x<sub>0</sub>)*h(t,x<sub>0</sub>),其中,*表示卷积运算,双站SAR海浪图像I(x′)可表式为,<img file="FSB00000583960100022.GIF" wi="2127" he="448" />其中,<img file="FSB00000583960100023.GIF" wi="199" he="144" />为发射机或接收机所产生的方位分辨率,λ为雷达波波长,x′为双站SAR图像位置,<img file="FSB00000583960100024.GIF" wi="723" he="150" />为由轨道加速度导致图像分辨率下降因子;步骤四,将双站SAR海浪图像进行二维FFT运算后,得到双站SAR海浪方向谱<img file="FSB00000583960100025.GIF" wi="22" he="27" />
地址 212003 江苏省镇江市梦溪路2号江苏科技大学