CN107607936A - 一种高频天地波雷达海洋表面流反演方法 - Google Patents
一种高频天地波雷达海洋表面流反演方法 Download PDFInfo
- Publication number
- CN107607936A CN107607936A CN201710775319.5A CN201710775319A CN107607936A CN 107607936 A CN107607936 A CN 107607936A CN 201710775319 A CN201710775319 A CN 201710775319A CN 107607936 A CN107607936 A CN 107607936A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- mtd
- mtr
- theta
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开一种基于傅里叶级数展开的高频天地波雷达海洋表面流反演算法,该方法为天地波超视距雷达面临的电离层倾斜、反射点高度难以确定、多普勒污染严重等问题提供了一个有效的解决方案。考虑到不同散射元对应的反射高度、电离层倾斜角的不同,提出傅里叶级数展开算法用以提取天波模式下各散射元的一阶Bragg频率,再利用多站组网接收模式遍历搜索公共散射元,从而求解出各个散射元对应的俯仰角和双基地角,实现精准定位。本发明摆脱了无倾斜平面电离层的假设,去除了由于电离层倾斜和反射点高度变化引起的探测误差,为天地波模式下海洋动力学参数的反演提供了一个新颖的反演算法,并提高了探测精度。
Description
技术领域
本发明属于雷达技术领域,特别涉及一种高频天地波雷达海洋表面流反演方法,用于解决天地波超视距雷达面临的电离层状态未知以至于难以实现精准定位的问题。
背景技术
高频天地波超视距雷达是一种适用于大范围海态监测的遥感新体制雷达设备,因其探测距离远、信噪比高等特点,近年来已经被广泛应用于目标探测、海洋动力学参数反演等领域。
天地波超视距雷达工作于天波反射/地波绕射的新传播模式,由内陆放置的发射机发射高频无线电波并斜入射到电离层,经过电离层反射到达地(海)面,再经地(海)面绕射传播到地波接收站。地波接收站可以放置在远离发射站的海岸或舰船上,相比传统的高频地波雷达具有很大的灵活性。然而,由于电离层具有分层、非平稳、时变等特性,导致无线电波的反射点位置难以确定,从而给探测与反演工作带来了一定的困难。
目前国内已有关于高频天地波雷达提取海流方面的研究,而这些研究都是建立于平面无倾斜电离层反射模型的基础之上,假设无线电波的波束照射在同一反射高度的平面电离层,通过解析几何来确定海面散射元的位置。这种简易模型忽略了由于电离层倾斜和高度变化引起的探测误差,会使散射元定位不准确,严重影响探测精度。这将使得高频天地波雷达的应用受到了极大的限制,即只有在电离层非常平稳无倾斜的状态下该简易模型才是有效的,而在实际情况中难以实现。高频段(3~30MHz)的无线波通常会在E/F层(200~300km)的高度发生反射,电离层的小角度倾斜和反射高度的微小变化都将引起无线电波在地(海)投射距离的巨大变化。因此,基于平面无倾斜电离层反射模型的天地波海态探测研究结果目前仍存在较大误差。
发明内容
本发明针对背景技术存在的问题,提出了一种高频天地波雷达海洋表面流反演新方法,利用傅里叶级数展开法重新提取天地波模式下的真实一阶Bragg频率,并利用最小二乘约束求解出天地波超视距雷达在发生海面散射时难以确定的俯仰角和双基地角,从而实现精准定位。
为达到上述目的,本发明采用如下技术方案:
一种高频天地波雷达海洋表面流反演方法,包括如下步骤,
步骤1:海洋回波数据预处理。
步骤2:将一阶Bragg散射频率fB和由径向流产生的多普勒偏移量Δf分别用傅里叶级数展开为关于群距离r和波达方向θ的函数。
步骤3:根据海洋回波的散射传播机理和多普勒效应建立了线性方程组,利用最小二乘拟合求解超定方程,提取一阶Bragg频率。
步骤4:双站遍历搜索公共散射元,求解双基地角和俯仰角。
步骤5:联立双站的径向流结果以及投影关系,合成天地波模式下矢量流结果。
上述技术方案的特点在于:
步骤1中对实验数据作前期预处理,包括距离匹配滤波、通道校准、谱区划分及DOA(Direction of Arrival)估计等关键技术,得到群距离r、波达方向θ和海洋回波的谱点信息f+/f-。
步骤2中将一阶Bragg散射频率fB和由径向流产生的多普勒偏移量Δf看作在海洋表面上关于r和θ的分布函数,由若干不同频率的正弦波叠加组合而成。根据二维傅里叶级数展开公式可以得到:
其中,P1,P2是Δf傅里叶展开的阶数,Q1,Q2是fB傅里叶展开的阶数,ε1和ε2分别为Δf和fB傅里叶级数展开的残余量,可忽略;
其中,k=1,2,3,4
x1,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
其中,
得到
[X Y]T=(ΛTΛ)-1ΛTF
为获得方程组的最优解,傅立叶展开阶数的选择也至关重要。高阶最小二乘解经常会带来数值计算的不稳定性,阶数太小又会导致收敛误差过大。在实际的探测试验中,Δf比fB变化更剧烈,所以我们考虑适当减小fB的展开阶数,增大Δf的展开阶数。综合以上几点,加上仿真测试,我们得到fB和Δf的展开阶数经验值为:
P1=3,P2=7,Q1=Q2=3
于是,我们可以提取真实的一阶Bragg频率fB为:
fBinversion=[B+Y B-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是双站接收天地波模式的信道传播示意图;
图3是电磁波在海面发生Bragg散射时海流分量关系示意图;
图4是矢量流合成示意图;
图5是东山站天波海洋回波距离多普勒谱图;
图6是天地波雷达与地波雷达矢量流结果对比图;图6(a)是2017年4月4日10:04海流对比图;图6(b)是2017年4月4日10:24海流对比图。
具体实施方式
下面结合实例对本发明作进一步的详细描述,此处所描述的实施示例仅用于说明和解释本发明,并不用于限定本发明。
一种高频天地波雷达海洋表面流反演方法,该实现方法主要包括雷达回波信号预处理和算法的具体实现,具体包括如下步骤,
步骤1:海洋回波数据预处理。武汉大学于2017年4月4进行了高频天地波雷达组网实验。天波发射站位于武汉,地波接收站分别位于福建东山和赤湖附近的一艘船只上。
对两个接收站的回波数据分别作处理。包括回波数据的匹配滤波、通道校准、谱区划分、DOA(Direction of Arrival)估计等,从而得到群距离r、波达方向θ和海洋回波的谱点信息f+/f-。附图5为东山站接收到的2017年4月4日10:04的天地波海洋回波距离多普勒谱。船载回波数据还包括晃动补偿与由船体金属引起的天线方向图校正的前期处理。
步骤2:将一阶Bragg散射频率fB和由径向流产生的多普勒偏移量Δf分别用傅里叶级数展开为关于群距离r和波达方向θ的函数。根据二维傅里叶级数展开公式可以得到:
其中,P1,P2是Δf傅里叶展开的阶数,Q1,Q2是fB傅里叶展开的阶数,ε1和ε2分别为Δf和fB傅里叶级数展开的残余量,可忽略;
其中,k=1,2,3,4
x1,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θ)
将步骤1中得到的回波数据的群距离r、频点f+/f-和方位信息θ带入以上展开式。
步骤3:根据海洋回波的散射传播机理和多普勒效应建立线性方程组,最小二乘拟合求解超定方程,提取一阶Bragg频率。在实际的海洋探测过程中,海洋回波谱分布在零频两侧,且满足如下关系式:
其中,f+和f-分别是正负一阶海洋回波的谱点,我们可以通过高频超视距雷达获得。
将步骤2中fB和Δf的傅里叶级数表达式带入步骤上式方程组,可得到如下超定方程组:
其中,A和Β分别是Δf和fB的傅里叶级数展开系数矩阵,下标的正负号分别对应正谱点和负谱点,
利用最小二乘求解超定方程,得到:
[X Y]T=(ΛTΛ)-1ΛTF
其中,
为获得方程组的最优解,傅立叶展开阶数的选择也至关重要。高阶最小二乘解经常会带来数值计算的不稳定性,阶数太小又会导致收敛误差过大。在实际的探测试验中,Δf比fB变化更剧烈,所以我们考虑适当减小fB的展开阶数,增大Δf的展开阶数。综合以上几点,加上仿真测试,我们得到fB和Δf的展开阶数经验值为:
P1=3,P2=7,Q1=Q2=3
于是,我们可以提取真实的一阶Bragg频率fB为:
fBinversion=[B+Y B-Y]T
步骤4:双站遍历搜索公共散射元,求解双基地角和俯仰角。利用两站雷达数据群距离的差值进行匹配。假设对应散射元S,雷达接收站A的群距离为G1,雷达接收站B的群距离为G2。散射元S到雷达接收站A的实际距离为R1,到雷达接收站B的实际距离为R2。在假设空中传播路径相同的情况下,群距离与散射元实际距离应满足下式:
R1-R2=G1-G2
由于天波群距离的分辨率ΔD=10km,要达到上式完全相等会出现散射元堆叠与模糊的现象。于是我们修正约束条件为:
其中ΔR=R1-R2,ΔG=G1-G2。
通过遍历谱点,初步确定两个天波雷达接收站的公共散射元S。傅里叶级数展开法可以确定两个雷达接收站的一阶Bragg频率fB1和fB2。假设照射公共散射元的是同一列波束,则双站在公共散射元的地方具有相同的俯仰角γ,且双基地角满足:
β2=β1+θ3
其中θ3为散射元到两个雷达接收站连线的夹角,如附图2。
天地波模式下的一阶Bragg频率表达式如下:
其中,f0为雷达工作频率,c为电磁波传播速度,γ为掠射角,β为双基地角。
联立fB1和fB2,求解俯仰角γ、双基地角β1和β2,以及实际散射元距离R。利用直达波反射高度h和R计算无倾斜时俯仰角γ0,约束条件返回修正为γ0-3<γ<γ0+3。进而通过解析几何求解散射元距离雷达接收站的实际距离R1或R2,实现坐标配准。
步骤5:联立双站的径向流结果以及投影关系,合成天地波模式下的矢量流。首先由径向流计算出单站Bragg散射方向的矢量流分量vb,由下式计算:
其中vr为单站雷达接收到的径向流速,vB为vb在入射波与反射波角平分线上的分量,α为vB与vb的夹角,通过解析立体几何,可得到
双站接收时,将两个分量进行合成即可得到矢量流。其中标号1和2分别代表两个接收站的标识,通过解如下方程组,即可求得矢量流的大小和方向。
其中κ1和κ2分别为两个接收站Bragg散射方向的分量流vb1和vb2与矢量流方向的夹角,v是矢量流,κ是分量流vb1和vb2的夹角。
附图6为2017年4月4日10:04和10:24天波矢量流与地波矢量流在空间海平面上的对比结果。其公共覆盖区域的流速大小均方根误差分别为16.4cm/s和12.8cm/s,流向均方根误差分别为15.7°和38.8°。
以上实验结果表明,本发明所提的一种高频天地波雷达海洋表面流反演方法对高频天地波雷达海流探测是适用的,并提高探测精度,具备与地波雷达相当的探测水平。
Claims (7)
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和θ的非周期性函数,其展开形式为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mo>&lsqb;</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>m&omega;</mi>
<mn>1</mn>
</msub>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>n&omega;</mi>
<mn>2</mn>
</msub>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>b</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>m&omega;</mi>
<mn>1</mn>
</msub>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>n&omega;</mi>
<mn>2</mn>
</msub>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>m&omega;</mi>
<mn>1</mn>
</msub>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>n&omega;</mi>
<mn>2</mn>
</msub>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>d</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>m&omega;</mi>
<mn>1</mn>
</msub>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>n&omega;</mi>
<mn>2</mn>
</msub>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>X</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>X</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<msub>
<mi>X</mi>
<mn>3</mn>
</msub>
<mo>,</mo>
<msub>
<mi>X</mi>
<mn>4</mn>
</msub>
<mo>&rsqb;</mo>
<mo>&CenterDot;</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>A</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>A</mi>
<mn>4</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>A</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>A</mi>
<mn>4</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>&CenterDot;</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>X</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>X</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<msub>
<mi>X</mi>
<mn>3</mn>
</msub>
<mo>,</mo>
<msub>
<mi>X</mi>
<mn>4</mn>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mi>A</mi>
<mo>&CenterDot;</mo>
<mi>X</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mi>B</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<msub>
<mi>Q</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<msub>
<mi>Q</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mo>&lsqb;</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>m&omega;</mi>
<mn>1</mn>
</msub>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>n&omega;</mi>
<mn>2</mn>
</msub>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>m&omega;</mi>
<mn>1</mn>
</msub>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>n&omega;</mi>
<mn>2</mn>
</msub>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msub>
<mi>&gamma;</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>m&omega;</mi>
<mn>1</mn>
</msub>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>n&omega;</mi>
<mn>2</mn>
</msub>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>m&omega;</mi>
<mn>1</mn>
</msub>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>n&omega;</mi>
<mn>2</mn>
</msub>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>B</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>B</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>B</mi>
<mn>4</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>&CenterDot;</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>Y</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>Y</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<msub>
<mi>Y</mi>
<mn>3</mn>
</msub>
<mo>,</mo>
<msub>
<mi>Y</mi>
<mn>4</mn>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mi>B</mi>
<mo>&CenterDot;</mo>
<mi>Y</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,P1,P2是Δf傅里叶展开的阶数,Q1,Q2是fB傅里叶展开的阶数,ε1和ε2分别为Δf和fB傅里叶级数展开的残余量;
<mrow>
<msub>
<mi>X</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mn>11</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mn>12</mn>
</mrow>
</msub>
<mo>,</mo>
<mo>...</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<mo>&rsqb;</mo>
<mo>,</mo>
<msub>
<mi>Y</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mn>11</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mn>12</mn>
</mrow>
</msub>
<mo>,</mo>
<mo>...</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<msub>
<mi>Q</mi>
<mn>1</mn>
</msub>
<msub>
<mi>Q</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<mo>&rsqb;</mo>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>A</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mn>11</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mn>12</mn>
</mrow>
</msub>
<mo>,</mo>
<mo>...</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<mo>&rsqb;</mo>
<mo>,</mo>
<msub>
<mi>B</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mn>11</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mn>12</mn>
</mrow>
</msub>
<mo>,</mo>
<mo>...</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<msub>
<mi>Q</mi>
<mn>1</mn>
</msub>
<msub>
<mi>Q</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<mo>&rsqb;</mo>
<mo>;</mo>
</mrow>
其中,k=1,2,3,4
x1,mn=amn,x2,mn=bmn,x3,mn=cmn,x4,mn=dmn;
y1,mn=αmn,y2,mn=βmn,y3,mn=γmn,y4,mn=ηmn;
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>m&omega;</mi>
<mn>1</mn>
</msub>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>n&omega;</mi>
<mn>2</mn>
</msub>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>A</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>m&omega;</mi>
<mn>1</mn>
</msub>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>n&omega;</mi>
<mn>2</mn>
</msub>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mn>3</mn>
<mo>,</mo>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>m&omega;</mi>
<mn>1</mn>
</msub>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>n&omega;</mi>
<mn>2</mn>
</msub>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>A</mi>
<mrow>
<mn>4</mn>
<mo>,</mo>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>m&omega;</mi>
<mn>1</mn>
</msub>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>n&omega;</mi>
<mn>2</mn>
</msub>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>.</mo>
</mrow>
4.根据权利要求3所述的一种高频天地波雷达海洋表面流反演方法,其特征在于:
所述步骤3中海洋回波谱分布在零频两侧,且满足如下关系式:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>f</mi>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>B</mi>
</msub>
<mo>=</mo>
<msup>
<mi>f</mi>
<mo>+</mo>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>f</mi>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>B</mi>
</msub>
<mo>=</mo>
<msup>
<mi>f</mi>
<mo>-</mo>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,f+和f-分别是正负一阶海洋回波的谱点;
将步骤2中所述fB和Δf的傅里叶级数展开式带入上式,得到如下超定方程组:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mo>+</mo>
</msub>
</mtd>
<mtd>
<msub>
<mi>B</mi>
<mo>+</mo>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mo>-</mo>
</msub>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>B</mi>
<mo>-</mo>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>Y</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msup>
<mi>f</mi>
<mo>+</mo>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mi>f</mi>
<mo>-</mo>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,A和Β分别是Δf和fB的傅里叶级数展开系数矩阵,下标的正负号分别对应正谱点和负谱点,
最小二乘求解超定方程组,提取fB的结果为:
fBinversion=[B+Y B-Y]T。
5.根据权利要求4所述的一种高频天地波雷达海洋表面流反演方法,其特征在于:所述最小二乘求解过程包括:
减小fB的展开阶数,增大Δf的展开阶数;得到fB和Δf的展开阶数经验值为:
P1=3,P2=7,Q1=Q2=3。
6.根据权利要求4所述的一种高频天地波雷达海洋表面流反演方法,其特征在于:
所述步骤4中遍历搜索公共散射元是利用两站雷达数据群距离的差值进行匹配;
假设对应散射元S,雷达接收站A的群距离为G1,雷达接收站B的群距离为G2;散射元S到雷达接收站A的实际距离为R1,到雷达接收站B的实际距离为R2;在假设空中传播路径相同的情况下,群距离与散射元实际距离应满足下式:
R1-R2=G1-G2
傅里叶级数展开法确定两个雷达接收站的一阶Bragg频率fB1和fB2,联立求解俯仰角γ和双基地角β1和β2,并计算出散射元的实际距离R。
7.根据权利要求6所述的一种高频天地波雷达海洋表面流反演方法,其特征在于:
所述步骤5中解如下方程组,合成天地波模式下的矢量流v:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>v</mi>
<mi> </mi>
<msub>
<mi>cos&kappa;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>b</mi>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>v</mi>
<mi> </mi>
<msub>
<mi>cos&kappa;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>b</mi>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&kappa;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&kappa;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mi>&kappa;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中κ1和κ2分别为两个接收站Bragg散射方向的分量流vb1和vb2与矢量流方向的夹角,κ是分量流vb1和vb2的夹角。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710775319.5A CN107607936B (zh) | 2017-08-31 | 2017-08-31 | 一种高频天地波雷达海洋表面流反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710775319.5A CN107607936B (zh) | 2017-08-31 | 2017-08-31 | 一种高频天地波雷达海洋表面流反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107607936A true CN107607936A (zh) | 2018-01-19 |
CN107607936B CN107607936B (zh) | 2019-12-24 |
Family
ID=61056527
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710775319.5A Active CN107607936B (zh) | 2017-08-31 | 2017-08-31 | 一种高频天地波雷达海洋表面流反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107607936B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108562896A (zh) * | 2018-04-23 | 2018-09-21 | 广西师范大学 | 一种基于三维正压浅海大陆架模型的深层海流反演方法 |
CN111610573A (zh) * | 2020-05-20 | 2020-09-01 | 北京遥测技术研究所 | 一种带有运动补偿的安检成像方法 |
CN111650665A (zh) * | 2020-05-20 | 2020-09-11 | 北京遥测技术研究所 | 一种带有运动补偿的安检成像系统及其使用方法 |
CN113901741A (zh) * | 2021-10-20 | 2022-01-07 | 哈尔滨工业大学(威海) | 海洋-电离层耗散动力学仿真方法、仿真系统、设备及应用 |
CN114330163A (zh) * | 2021-12-27 | 2022-04-12 | 哈尔滨工业大学(威海) | 高频地波超视距雷达台风-电离层扰动动力学模型建模方法 |
CN114325625A (zh) * | 2021-11-23 | 2022-04-12 | 哈尔滨工业大学(威海) | 海洋-电离层混沌动力学特征提取方法、装置、介质、应用 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9423495B1 (en) * | 2014-02-20 | 2016-08-23 | Raytheon Company | Ship-based over-the-horizon radar |
CN104749563B (zh) * | 2015-03-26 | 2017-02-22 | 武汉大学 | 从高频地波雷达海洋回波一阶Bragg谐振峰中提取浪高的方法 |
CN106226762B (zh) * | 2016-09-29 | 2018-07-27 | 武汉大学 | 一种用于确定高频天地波超视距雷达探测区域空间分布的方法 |
-
2017
- 2017-08-31 CN CN201710775319.5A patent/CN107607936B/zh active Active
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108562896A (zh) * | 2018-04-23 | 2018-09-21 | 广西师范大学 | 一种基于三维正压浅海大陆架模型的深层海流反演方法 |
CN111610573A (zh) * | 2020-05-20 | 2020-09-01 | 北京遥测技术研究所 | 一种带有运动补偿的安检成像方法 |
CN111650665A (zh) * | 2020-05-20 | 2020-09-11 | 北京遥测技术研究所 | 一种带有运动补偿的安检成像系统及其使用方法 |
CN113901741A (zh) * | 2021-10-20 | 2022-01-07 | 哈尔滨工业大学(威海) | 海洋-电离层耗散动力学仿真方法、仿真系统、设备及应用 |
CN113901741B (zh) * | 2021-10-20 | 2023-08-18 | 哈尔滨工业大学(威海) | 海洋-电离层耗散动力学仿真方法、仿真系统、设备 |
CN114325625A (zh) * | 2021-11-23 | 2022-04-12 | 哈尔滨工业大学(威海) | 海洋-电离层混沌动力学特征提取方法、装置、介质、应用 |
CN114325625B (zh) * | 2021-11-23 | 2024-03-12 | 哈尔滨工业大学(威海) | 海洋-电离层混沌动力学特征提取方法、装置、介质、应用 |
CN114330163A (zh) * | 2021-12-27 | 2022-04-12 | 哈尔滨工业大学(威海) | 高频地波超视距雷达台风-电离层扰动动力学模型建模方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107607936B (zh) | 2019-12-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107607936A (zh) | 一种高频天地波雷达海洋表面流反演方法 | |
CN104391281A (zh) | 提高天波雷达海面船舶目标跟踪定位精度的方法 | |
CN105445730B (zh) | 一种基于角度分集的海洋流场反演星载sar系统及其方法 | |
CN104035095B (zh) | 基于空时最优处理器的低空风切变风速估计方法 | |
CN103869311B (zh) | 实波束扫描雷达超分辨成像方法 | |
CN108828550A (zh) | 一种海缆区界船舶目标位置检测方法 | |
CN102288964A (zh) | 一种星载高分辨率合成孔径雷达的成像处理方法 | |
CN105487074B (zh) | 一种双基地合成孔径雷达数值距离多普勒成像方法 | |
CN106093935B (zh) | 俯仰宽波束雷达地面和低空目标分类方法 | |
CN114942428B (zh) | 多目标跟踪的混合体制雷达数据关联算法 | |
CN102749621A (zh) | 一种双基地合成孔径雷达频域成像方法 | |
Inggs et al. | Planning and design phases of a commensal radar system in the FM broadcast band | |
CN104076338A (zh) | 基于数字高程和数字地表覆盖的机载雷达杂波仿真方法 | |
CN106526583A (zh) | 一种基于天线方向图信息的地面运动目标定位方法 | |
CN103744081A (zh) | 一种机载圆迹合成孔径雷达高精度三维成像补偿方法 | |
Sun et al. | Track-to-track association based on maximum likelihood estimation for T/RR composite compact HFSWR | |
CN104101869B (zh) | 一种极坐标下的地波雷达运动目标仿真建模方法 | |
CN106291547B (zh) | 基于天线方向图辅助的多普勒模糊分量自适应抑制方法 | |
CN103091682A (zh) | 基于时频分析InISAR多动目标成像和运动轨迹重建法 | |
Maresca et al. | Theoretical limitation of the sea on the detection of low Doppler targets by over-the-horizon radar | |
CN103245950A (zh) | 配接于多功能雷达的气象终端 | |
CN104360335B (zh) | 基于ais目标指示的电离层参数反演方法 | |
CN104391288B (zh) | 拟合干涉相位的虚警剔除方法 | |
Buck et al. | Steps towards the preparation of a Wavemill mission | |
CN107329117A (zh) | 一种基于改进omp的双基机载雷达自适应杂波谱补偿方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |