当前位置: 首页 > 专利查询>武汉大学专利>正文

一种高频天地波雷达海洋表面流反演方法技术

技术编号:17096407 阅读:69 留言:0更新日期:2018-01-21 07:39
本发明专利技术公开一种基于傅里叶级数展开的高频天地波雷达海洋表面流反演算法,该方法为天地波超视距雷达面临的电离层倾斜、反射点高度难以确定、多普勒污染严重等问题提供了一个有效的解决方案。考虑到不同散射元对应的反射高度、电离层倾斜角的不同,提出傅里叶级数展开算法用以提取天波模式下各散射元的一阶Bragg频率,再利用多站组网接收模式遍历搜索公共散射元,从而求解出各个散射元对应的俯仰角和双基地角,实现精准定位。本发明专利技术摆脱了无倾斜平面电离层的假设,去除了由于电离层倾斜和反射点高度变化引起的探测误差,为天地波模式下海洋动力学参数的反演提供了一个新颖的反演算法,并提高了探测精度。

【技术实现步骤摘要】
一种高频天地波雷达海洋表面流反演方法
本专利技术属于雷达
,特别涉及一种高频天地波雷达海洋表面流反演方法,用于解决天地波超视距雷达面临的电离层状态未知以至于难以实现精准定位的问题。
技术介绍
高频天地波超视距雷达是一种适用于大范围海态监测的遥感新体制雷达设备,因其探测距离远、信噪比高等特点,近年来已经被广泛应用于目标探测、海洋动力学参数反演等领域。天地波超视距雷达工作于天波反射/地波绕射的新传播模式,由内陆放置的发射机发射高频无线电波并斜入射到电离层,经过电离层反射到达地(海)面,再经地(海)面绕射传播到地波接收站。地波接收站可以放置在远离发射站的海岸或舰船上,相比传统的高频地波雷达具有很大的灵活性。然而,由于电离层具有分层、非平稳、时变等特性,导致无线电波的反射点位置难以确定,从而给探测与反演工作带来了一定的困难。目前国内已有关于高频天地波雷达提取海流方面的研究,而这些研究都是建立于平面无倾斜电离层反射模型的基础之上,假设无线电波的波束照射在同一反射高度的平面电离层,通过解析几何来确定海面散射元的位置。这种简易模型忽略了由于电离层倾斜和高度变化引起的探测误差,会使散射元定位不准确,严重影响探测精度。这将使得高频天地波雷达的应用受到了极大的限制,即只有在电离层非常平稳无倾斜的状态下该简易模型才是有效的,而在实际情况中难以实现。高频段(3~30MHz)的无线波通常会在E/F层(200~300km)的高度发生反射,电离层的小角度倾斜和反射高度的微小变化都将引起无线电波在地(海)投射距离的巨大变化。因此,基于平面无倾斜电离层反射模型的天地波海态探测研究结果目前仍存在较大误差。
技术实现思路
本专利技术针对
技术介绍
存在的问题,提出了一种高频天地波雷达海洋表面流反演新方法,利用傅里叶级数展开法重新提取天地波模式下的真实一阶Bragg频率,并利用最小二乘约束求解出天地波超视距雷达在发生海面散射时难以确定的俯仰角和双基地角,从而实现精准定位。为达到上述目的,本专利技术采用如下技术方案:一种高频天地波雷达海洋表面流反演方法,包括如下步骤,步骤1:海洋回波数据预处理。步骤2:将一阶Bragg散射频率fB和由径向流产生的多普勒偏移量Δf分别用傅里叶级数展开为关于群距离r和波达方向θ的函数。步骤3:根据海洋回波的散射传播机理和多普勒效应建立了线性方程组,利用最小二乘拟合求解超定方程,提取一阶Bragg频率。步骤4:双站遍历搜索公共散射元,求解双基地角和俯仰角。步骤5:联立双站的径向流结果以及投影关系,合成天地波模式下矢量流结果。上述技术方案的特点在于:步骤1中对实验数据作前期预处理,包括距离匹配滤波、通道校准、谱区划分及DOA(DirectionofArrival)估计等关键技术,得到群距离r、波达方向θ和海洋回波的谱点信息f+/f-。步骤2中将一阶Bragg散射频率fB和由径向流产生的多普勒偏移量Δf看作在海洋表面上关于r和θ的分布函数,由若干不同频率的正弦波叠加组合而成。根据二维傅里叶级数展开公式可以得到:其中,P1,P2是Δf傅里叶展开的阶数,Q1,Q2是fB傅里叶展开的阶数,ε1和ε2分别为Δf和fB傅里叶级数展开的残余量,可忽略;其中,k=1,2,3,4x1,mn=amn,x2,mn=bmn,x3,mn=cmn,x4,mn=dmn;y1,mn=αmn,y2,mn=βmn,y3,mn=γmn,y4,mn=ηmn;A1,mn=cos(mω1r)cos(nω2θ),A2,mn=cos(mω1r)sin(nω2θ),。A3,mn=sin(mω1r)cos(nω2θ),A4,mn=sin(mω1r)sin(nω2θ)步骤3中根据海洋回波的散射传播机理和多普勒效应建立线性方程组。在实际的海洋探测过程中,海洋回波谱分布在零频两侧,且满足如下关系式:即其中,f+和f-分别是正负一阶海洋回波的谱点,我们可以通过高频超视距雷达获得。将步骤2中fB和Δf的傅里叶级数表达式带入上式方程组,可得到如下超定方程组:其中,A和Β分别是Δf和fB的傅里叶级数展开系数矩阵,下标的正负号分别对应正谱点和负谱点,为保证线性方程组有解,方程组的个数不能少于未知变量的个数。然而,在长期的雷达观测试验中发现,探测结果会存在一定的随机误差,因此我们需要使用尽可能多的谱点将随机误差平均减弱。此时,问题就变成了超定方程组求解,这里采用最小二乘的方法求解:minα||Λα-F||2其中,得到[XY]T=(ΛTΛ)-1ΛTF为获得方程组的最优解,傅立叶展开阶数的选择也至关重要。高阶最小二乘解经常会带来数值计算的不稳定性,阶数太小又会导致收敛误差过大。在实际的探测试验中,Δf比fB变化更剧烈,所以我们考虑适当减小fB的展开阶数,增大Δf的展开阶数。综合以上几点,加上仿真测试,我们得到fB和Δf的展开阶数经验值为:P1=3,P2=7,Q1=Q2=3于是,我们可以提取真实的一阶Bragg频率fB为:fBinversion=[B+YB-Y]T步骤4中假设照射在同一散射元上的电磁波发生各个方向的散射,两个雷达接收站分别接收到该散射元散射的回波,两个回波在空中的传播路径认为是相同的,沿海面散射的路径分别指向雷达接收站,如附图2所示。利用两站雷达数据群距离的差值进行匹配。假设对应散射元S,雷达接收站A的群距离为G1,雷达接收站B的群距离为G2。散射元S到雷达接收站A的实际距离为R1,到雷达接收站B的实际距离为R2。在假设空中传播路径相同的情况下,群距离与散射元实际距离应满足下式:R1-R2=G1-G2由此可以确定两个天波雷达接收站的公共散射元S。傅里叶级数展开法可以确定两个雷达接收站的一阶Bragg频率fB1和fB2。假设照射公共散射元的是同一列波束,则双站在公共散射元的地方具有相同的俯仰角γ,且双基地角满足:β2=β1+θ3其中θ3为散射元到两个雷达接收站连线的夹角,如附图2。天地波模式下的一阶Bragg频率表达式如下:其中,f0为雷达工作频率,c为电磁波传播速度,γ为掠射角,β为双基地角。基于天地波模式下的一阶Bragg频率表达式,可以联立求解俯仰角γ和双基地角β1和β2。进而通过解析几何求解散射元距离雷达接收站的实际距离R1或R2,实现坐标配准。步骤5中首先计算出单站Bragg散射方向的矢量流分量,立体几何关系如附图3。其中为天波照射海面的入射波,为反射波,则Bragg散射方向的海流分量流vb可用下式计算:其中vr为单站雷达接收到的径向流速,vB为vb在入射波与反射波角平分线上的分量,α为vB与vb的夹角,通过解析立体几何,可得到双站接收时,将两个分量进行合成即可得到矢量流,如附图4所示。其中标号1和2分别代表两个接收站的标识,通过解如下方程组,即可求得矢量流的大小和方向。其中κ1和κ2分别为两个接收站Bragg散射方向的分量流vb1和vb2与矢量流方向的夹角,v是矢量流,κ是分量流vb1和vb2的夹角。与现有技术相比,本专利技术的优势在于:本专利技术提出的高频天地波雷达海洋表面流反演方法,通过傅里叶级数展开模型和双站搜索公共散射元的方法,实现精准定位,提高了天地波模式下的海流探测精度。该方法摆脱了平面无倾斜电离层的假设,求解出受电离层影响后的真实一阶Bragg频率,进而求解海本文档来自技高网
...
一种高频天地波雷达海洋表面流反演方法

【技术保护点】
一种高频天地波雷达海洋表面流反演方法,其特征在于,包括如下步骤:步骤1:海洋回波数据预处理;步骤2:将一阶Bragg散射频率fB和由径向流产生的多普勒偏移量Δf分别用傅里叶级数展开为关于群距离r和波达方向θ的函数;步骤3:根据海洋回波的散射传播机理和多普勒效应建立线性方程组,最小二乘拟合求解超定方程,提取一阶Bragg频率;步骤4:双站遍历搜索公共散射元,求解双基地角和俯仰角;步骤5:联立双站的径向流结果以及投影关系,合成天地波模式下的矢量流。

【技术特征摘要】
1.一种高频天地波雷达海洋表面流反演方法,其特征在于,包括如下步骤:步骤1:海洋回波数据预处理;步骤2:将一阶Bragg散射频率fB和由径向流产生的多普勒偏移量Δf分别用傅里叶级数展开为关于群距离r和波达方向θ的函数;步骤3:根据海洋回波的散射传播机理和多普勒效应建立线性方程组,最小二乘拟合求解超定方程,提取一阶Bragg频率;步骤4:双站遍历搜索公共散射元,求解双基地角和俯仰角;步骤5:联立双站的径向流结果以及投影关系,合成天地波模式下的矢量流。2.根据权利要求1所述的一种高频天地波雷达海洋表面流反演方法,其特征在于:所述步骤1,高频天地波雷达接收到的海洋回波经过匹配滤波和傅里叶变换得到距离多普勒谱,数据预处理包括:通道校准、谱区划分及DOA估计,得到群距离r、波达方向θ和海洋回波的谱点信息f+/f-。3.根据权利要求2所述的一种高频天地波雷达海洋表面流反演方法,其特征在于:所述步骤2中将fB和Δf用傅里叶级数展开是基于海洋状态的不确定性可以用正弦函数的叠加来拟合的特性,fB和Δf可看作是关于r和θ的非周期性函数,其展开形式为:

【专利技术属性】
技术研发人员:吴雄斌李苗张兰岳显昌
申请(专利权)人:武汉大学
类型:发明
国别省市:湖北,42

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1