CN102981144B - 空中运动平台对目标的三维无源定位方法 - Google Patents
空中运动平台对目标的三维无源定位方法 Download PDFInfo
- Publication number
- CN102981144B CN102981144B CN201210476641.5A CN201210476641A CN102981144B CN 102981144 B CN102981144 B CN 102981144B CN 201210476641 A CN201210476641 A CN 201210476641A CN 102981144 B CN102981144 B CN 102981144B
- Authority
- CN
- China
- Prior art keywords
- target
- main website
- extension station
- location
- station
- 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.)
- Expired - Fee Related
Links
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种空中运动平台对目标的三维无源定位方法,主要解决现有技术中空中运动平台不能同时对自身和目标进行无源定位的问题,其实现步骤为:1)对接收目标散射回波的主站和各个辅站进行时间同步;2)雷达发射站辐射目标检测所需的电磁波;3)主站记录发射直达波的到达时刻并记录目标散射回波的到达时刻和各个辅站记录主站或目标的散射回波的到达时刻;4)辅站将自身位置信息连同主站或目标的散射回波的到达时刻信息上传到主站;5)主站对自身和目标进行定位;6)将目标在全局坐标系下的位置转换为局部坐标系下的位置。本发明具有定位精度高,隐蔽性强的优点,可用于空中运动平台对目标的三维无源定位。
Description
技术领域
本发明属于目标定位技术领域,特别涉及目标的三维无源定位方法,可用于空中运动平台对目标的三维无源定位。
背景技术
目前空中运动平台普遍采用有源的工作方式对目标进行定位,由于其自身有辐射源,因此很容易被对方发现并实施有效打击,另外空中运动平台自身的位置往往需要借助其他定位系统如GPS获得,这使得当其他定位系统失灵时空中运动平台将失去对自身定位的能力从而无法获得目标相对于自身的位置信息。因此同时对空中运动平台和目标进行无源定位具有非常重要的现实意义和应用前景,但同时对空中运动平台和目标进行定位是多年来国内外的技术难题,尤其是当有源定向干扰存在时,常规有源雷达和无源雷达往往得到错误的目标波达方向估计,这为空中运动平台对目标进行准确定位带来了更大的难度。
在现有的无源定位技术中,国内外的研究主要是针对空中运动平台自身或者仅仅针对目标进行定位,其定位方法主要有以下几种:
(1)测向目标定位方法,二维平面测向定位通过测量两个不同位置上雷达辐射电磁波的方向DOA,各站测得的雷达方向数据结合两个雷达站之间的基线距离按三角测量法交会计算出雷达的位置;三维空间测向定位则是利用飞机或卫星上的侦察设备,测出地面雷达信号的俯仰角和方向角,并利用导航数据提供的飞机或卫星高度来确定雷达的地理位置。测向目标定位的缺点是定位精度较低,而且易受定向干扰源的影响。
(2)到达时间定位方法,利用运动平台上的接收机通过观测位置已知的四颗卫星信号到接收机的到达时间TOA,利用空间球面交会的方法确定接收机所在点位置的三维坐标信息,也即空中运动平台当前的位置。这种方法主要是根据各个发射站的信号到达接收机的时间算出接收机相对于发射站的位置信息,但系统比较复杂,对各个发射站的时间同步精度要求较高,相应的信号处理系统复杂,不能对没有部署接收机的目标进行定位。应用此技术进行目标定位的一个例子是GPS全球定位系统。
(3)时差定位方法,由设在地面的1个主台与3个副台合成的台链和运动平台上的接收机组成。接收机测定主、副台发射的两个脉冲信号的时间差TDOA即可获得运动平台到主、副台的距离差。距离差保持不变的航迹是一个双曲面,根据多个双曲面的交点可以定出运动平台的位置,但该系统同样仅仅是针对装有接收机的空中运动平台进行定位,不能对没有安装接收机的目标进行定位,罗兰C远程定位系统就是采用时差定位系统实现的。
(4)频差定位方法,该方法是利用来波信号到达各传感器之间的多普勒频差信息FDOA,并以此来确定目标的位置,此方法的缺点是需要各个站相对于目标运动,一般频差定位难以单独使用达到满意的定位效果,需要配合其他无源定位方法联合使用。
(5)外辐射源的无源定位技术,该技术是利用已知地理位置的民用调频广播FM、电视TV、卫星通信、手机基站等外辐射源信号来实现无源定位。此类定位系统通过设有多个接收站分别接收这些直射信号和目标反射信号,测量它们的TDOA、DOA和多普勒频移,FDOA实现目标的定位。由于这些外辐射源信号的基带调制信号是时变的,不能在接收端进行相干处理,因此此类定位技术的定位区域较小,定位精度较低。
上述各方法主要是针对安装接收机的运动平台自身进行无源定位,或者是多个接收站对空中运动平台和目标进行无源定位,不能实现空中运动平台对空中目标进行无源定位,对于空中运动平台对目标进行无源定位的方法目前还未见报道。
发明内容
本发明的目的在于克服传统时差定位中存在的定位模糊缺点,提供一种空中运动平台对目标的三维无源定位方法,以提高定位精度,增强隐蔽性。
为实现上述目的,本发明包括如下步骤:
1)对接收目标散射回波的主站和各个辅站进行时间同步,给目标回波到达时刻的测量提供统一的时间基准,所述主站为空中运动平台,所述辅站为地面上的三个接收站;
2)雷达发射站辐射目标检测所需的电磁波;
3)主站测量并记录雷达发射站直达波的到达时刻TD0;
4)主站记录目标散射回波的到达时刻t0,各个辅站记录主站或目标的散射回波的到达时刻si,i=1,2,3;
5)辅站将主站或目标的散射回波的到达时刻si,i=1,2,3连同辅站在全局坐标系下的位置(xi,yi,zi)T进行编码和调制,并通过无线信道上传到主站,i=1,2,3,其中()T表示转置;
6)主站对辅站上传的数据解调和解码后,剔除主站自身产生的回波到达时刻tDi,i=1,2,3,并对主站自身进行定位;
7)主站结合目标散射回波到达主站的时刻t0和辅站上传的目标散射回波到达辅站的时刻ti,i=1,2,3对目标进行定位,即将求解目标位置的非线性方程组转化为与其对应的最优化问题,通过求取该最优化问题的最优解得到目标的定位结果
8)将目标在全局坐标系下的定位结果转换为以主站为坐标原点的局部坐标系下的定位结果
本发明与现有技术相比,具有如下优点:
(1)本发明利用时差定位算法的特性,将目标回波信号到一个主站以及三个辅站之间的到达时间差转变为距离差,从而对目标位置进行求解,此方法中主站不辐射电磁信号,不易被检测,隐蔽性强。
(2)本发明定位方法通过合理布站采用三个双曲面的单页的交点对目标进行定位,采用三个半球面的交点对主站自身进行定位,采用直接求解非线性方程组的方式无精度损失的对目标位置进行求解,通过合理的选择求解的初值可以使传统时差定位算法中存在的定位模糊问题得到消除,不需要引入其他解定位模糊的手段,提高了目标真实位置的估计精度。
(4)本发明中目标到辅站的到达时间以及全局坐标系下辅站的站点位置信息和站号等信息采用编码和调制的方式发送给主站,主站在已知编码方案以及调制方式的具体参数时才能得到辅站提供的定位信息,辅站到主站的通信信息不易被侦察。
(5)本发明采用的定位技术可以获得主站自身位置,从而可以以主站为局部坐标系原点建立局部坐标系,得到目标在局部坐标系下的高精度估计值。
附图说明
图1是本发明的应用场景图;
图2是本发明的实现总流程图;
图3是应用本发明在三个等时差双曲面只有一个交点时对目标三维无源定位的仿真结果;
图4是应用本发明在三个等时差双曲面有两个交点时对目标三维无源定位的仿真结果;
图5是应用本发明在图3模型下的定位精度几何稀释GDOP;
图6是应用本发明在图4模型下的定位精度几何稀释GDOP;
图7是用不同方法在图3模型下目标定位精度随时差测量误差变化的定位精度仿真图。
图8是用不同方法在图4模型下目标定位精度随时差测量误差变化的定位精度仿真图。
具体实施方式
参照图1,本发明实施例所用的主站3为空中运动平台,发射站1辐射目标检测所需的电磁波,2是待定位目标,主站3负责接收发射站辐射的直达波以及待定位目标2的散射回波,辅站4、辅站5、辅站6分别是地面上的三个接收站。
参照图2,本实施例的具体实现步骤如下:
步骤1,对主站3、辅站4、辅站5、辅站6进行时间同步。
如图1所示,本发明是通过三个辅站上传到主站数据的到达时间对主站3自身进行定位,通过待定位目标2散射回波到达各个站的时间差实现主站3对待定位目标2的定位。因此本发明首先利用GPS时钟信号对主站3、辅站4、辅站5、辅站6进行时间同步,目的是给整个系统提供一个统一的时间基准,以确保空中运动平台自身定位精度以及目标定位精度。
步骤2,雷达发射站辐射目标检测所需的电磁波。
通过发射站1向待定位目标2所在区域辐射相位编码电磁波信号,该电磁波应具有足够大的功率,以满足待定位目标2检测信噪比的要求。
步骤3,主站3测量并记录发射直达波的到达时刻TD0。
本实例以直达波脉冲包络前沿的到达主站的时刻作为直达波的到达时刻TD0,其表达式为:
其中m是脉冲数,T是脉冲重复周期,β=v/c,α=(1+β)/(1-β),c是光速,RD0是发射站1到主站3的距离,v是主站3与发射站1之间的径向速度,eD0是零均值的高斯白噪声。
步骤4,主站3记录待定位目标2散射回波到达时刻t0,各个辅站记录主站3或待定位目标2的散射回波到达时刻si,i=1,2,3。
发射的电磁波碰到待定位目标2以及主站3后分别发生散射,主站3记录待定位目标2散射回波的到达时刻t0。如果主站3的散射回波到达辅站4、辅站5、辅站6的时刻分别为si,i=1,2,3,则辅站4、辅站5、辅站6各自记录主站3散射回波到各个辅站的到达时刻tDi,i=1,2,3,此时令si=tDi;如果待定位目标2的散射回波在si,i=1,2,3时刻分别到达辅站4、辅站5、辅站6,则辅站4、辅站5、辅站6各自记录待定位目标2散射回波到各个辅站的到达时刻ti,i=1,2,3,此时令si=ti。
步骤5,各个辅站将主站3或待定位目标2的散射回波到达时刻si,i=1,2,3连同各个辅站自身在全局坐标系下的位置(xi,yi,zi)T进行编码和调制,并通过无线信道上传到主站3,i=1,2,3。
将主站3或待定位目标2的散射回波到达时刻si,i=1,2,3连同辅站4、辅站5、辅站6各自在全局坐标系下的坐标(xi,yi,zi)T,i=1,2,3依次经过Gold Code编码和二相相移键控BPSK调制后在si+Δτ,i=1,2,3时刻分别通过无线信道发送给主站(x0,y0,z0)T,其中Δτ是编码和调制占用的信号处理时间。
步骤6,主站3对辅站上传的定位数据进行解调、解码。
主站3在时刻Ti,i=1,2,3接收到辅站4、辅站5、辅站6各自上传的数据并对其依次进行二相相移键控BPSK解调和Gold Code解码,获得此数据中包含的主站3或目标2的散射回波到达时刻si,i=1,2,3,以及辅站4、辅站5、辅站6各自在全局坐标系下的坐标(xi,yi,zi)T,i=1,2,3。
步骤7,主站3对辅站上传的主站3自身产生的目标回波到达时刻tDi,i=1,2,3进行剔除。
辅站上传的定位数据中的回波到达时刻si,i=1,2,3有可能是主站3散射回波到达时刻tDi,i=1,2,3,也有可能是待定位目标2散射回波到达时刻ti,i=1,2,3。给定一个接近于零的正数γ,计算表达式Ei=|si-Ti+Δτ|-|si-TD0|,i=1,2,3的数值,如果|Ei|≤γ说明si是主站3自身散射回波到达辅站的时刻tDi,i=1,2,3,将其剔除;如果|Ei|>γ则说明si是目标2散射回波到达辅站的时刻ti,将其保存,其中||表示取绝对值,Ti,i=1,2,3分别是辅站4、辅站5、辅站6上传的数据到主站3的到达时刻。
步骤8,主站3对自身进行无源定位。
8.1)将辅站4、辅站5、辅站6各自上传到主站3数据的到达时间Ti,i=1,2,3代入如下所示的三个半球面方程组:
其中c代表光速,T1是辅站4上传数据到达主站3的时刻,T2是辅站5上传数据到达主站3的时刻,T3是辅站6上传数据到达主站3的时刻,s1是主站3或待定位目标2的散射回波到达辅站4的时刻,s2是主站3或待定位目标2的散射回波到达辅站5的时刻,s3是主站3或待定位目标2的散射回波到达辅站6的时刻,Δτ是辅站编码和调制占用的信号处理时间,ρ1是主站3到辅站4的距离,ρ2是主站3到辅站5的距离,ρ3是主站3到辅站6的距离,w1是与ρ1对应的零均值噪声,w2是与ρ2对应的零均值噪声,w3是与ρ3对应的零均值噪声。
8.2)对以上方程组进行求解,该方程组的解(x0,y0,z0)T就是主站3自身的定位结果,其中()T表示转置。
步骤9,主站3对待定位目标2进行无源定位。
9.1建立等效的目标位置求解优化模型。
将待定位目标2到各个辅站的距离Ri,i=1,2,3,待定位目标2到主站3的距离R0,待定位目标2散射回波到主站3的到达时刻t0及待定位目标2散射回波到各个辅站的到达时刻ti,i=1,2,3代入如下三个等时差双曲面方程:
式中ΔRi是待定位目标2到各个辅站的距离Ri与待定位目标2到主站3的距离R0之差,ei是零均值的量测噪声,i=1,2,3,根据最小均方误差准则,将所述三个双叶双曲面方程组的求解转化为如下最优化问题的求解:
式(4)中J(u)是代价函数,g(u)=[g1,g2,g3]T为辅助函数,u=(x,y,z)T为待定位目标2的真实位置,为待定位目标2真实位置u的估计值,为待定位目标2位置x分量的估计值,为待定位目标2位置y分量的估计值,为待定位目标2位置z分量的估计值,gi=Ri-R0-c(ti-t0)-ei,i=1,2,3分别为辅助函数g(u)包含的3个分量, 分别为待定位目标2到各个辅站的距离,为待定位目标2到主站3的距离。
9.2对优化问题求解。
Step1:选取迭代初值,即当待定位目标2的位置高于主站3时,迭代初值u0=[x0,y0,z0]T应该满足初值的z轴分量z0大于主站3的z轴分量z0;当待定位目标2低于主站3时,初值u0应满足初值的z轴分量z0小于主站的z轴分量z0;
Step2:将初值u0带入式(4),计算代价函数值J(u0);
Step3:对于精度ε>0,若存在J(u0)<ε,则将u0作为待定位目标2位置的估计值否则令迭代序号k=0继续执行下面的步骤;
Step4:根据式(5)计算辅助函数g(u)在第k步目标位置估计值uk处的雅可比矩阵F(uk)
其中uk=[xk,yk,zk]T,xk是第k步待定位目标2位置估计值uk的x轴分量,yk是uk的y轴分量,zk是uk的z轴分量;
Step5:根据式(6)计算代价函数J在第k步待定位目标2位置估计值uk处的梯度
Step6:根据式(7)迭代计算第k+1步待定位目标2位置估计值uk+1
其中为代价函数Jk的负梯度,参数μk+1为自适应时变步长;
其中α是一个加速因子,其取值范围为1≤α≤2,||||表示向量的2范数;
Step7:置k←k+1,即将迭代序号k加1,重复执行Step4到Step7直到代价函数Jk满足精度要求J(uk)<ε,将该uk作为待定位目标2位置的估计值
步骤10,将待定位目标2在全局坐标系下的定位结果转换为以主站3建立的局部坐标系下的定位结果
获得待定位目标2在全局直角坐标系下的定位结果后,通过式(9)将其转换为待定位目标2在局部直角坐标系下的定位结果其中[]T表示转置。
其中ux=(ux1,ux2,ux3)T是局部直角坐标系X'轴在全局坐标系下的单位向量,uy=(uy1,uy2,uy3)T是局部直角坐标系Y'轴在全局坐标系下的单位向量,uz=(uz1,uz2,uz3)T是局部直角坐标系Z'轴在全局坐标系下的单位向量,局部坐标系的坐标轴指向矩阵U=[ux,uy,uz]是单位正交矩阵。
本发明的效果通过以下仿真试验进一步说明:
1.仿真内容
仿真1,将本发明方法应用于三个等时差双曲面只有一个交点对应的定位模型,对该方法的定位性能进行仿真。
设辅站4在全局坐标系下的坐标位置为(x1,y1,z1)T=(-30,42,0)Tkm,辅站5在全局坐标系下的坐标位置为(x2,y2,z2)T=(-30,-42,0)Tkm,辅站6在全局坐标系下的坐标位置为(x3,y3,z3)T=(48,0,0)Tkm;各个辅站上传数据到主站经历的时间分别为0.2187845211ms,0.2579405616ms以及0.2380476143ms;假设目标位置为(x,y,z)T=(30,40,20)Tkm,迭代初值取u0=(x0,y0,z0)T=(0,0,7)Tkm,时差测量误差为10ns,迭代收敛精度为ε=0.01。将这些参数代入式(4)对目标位置进行求解,经过44次迭代得到目标在全局坐标系下的坐标位置 结果如图3所示,其对应的定位精度的几何稀释GDOP如图5所示。在不同时差测量误差下,分别对Chan算法和本发明做1000次蒙特卡罗试验,目标定位精度随时差测量误差变化的定位精度结果如图7所示。以空中运动平台当前位置(x0,y0,z0)T为原点建立局部坐标系,局部坐标系的坐标轴指向矩阵U为单位阵I,得到目标在局部坐标系下的坐标位置为
仿真2,将本发明方法应用于三个等时差双曲面有两个交点对应的定位模型,对该方法的定位性能进行仿真。
设辅站4在全局坐标系下的坐标位置为(x1,y1,z1)T=(-30,42,0)Tkm,辅站5在全局坐标系下的坐标位置为(x2,y2,z2)T=(-30,-42,0)Tkm,辅站6在全局坐标系下的坐标位置为(x3,y3,z3)T=(48,0,0)Tkm;各个辅站上传数据到主站的经历的时间分别为0.1427118309ms,0.1975685535ms,0.1707825128ms;假设目标位置为(x,y,z)T=(30,40,20)Tkm,迭代初值u0=(x0,y0,z0)T=(0,0,7)Tkm,时差测量误差为10ns。将这些参数代入式(4)对目标位置进行求解,经过129次迭代得到正确的目标位置的估计值为结果如图4所示,其对应的定位精度的几何稀释GDOP,如图6所示。在不同时差测量误差下,分别对Chan算法和本方法做1000次蒙特卡罗试验,目标定位精度随时差测量误差变化的定位精度结果如图8所示。在不同时差测量误差下,分别对Chan算法和本方法做1000次蒙特卡罗试验,目标定位精度随时差测量误差变化的定位精度结果如图7所示。以空中运动平台当前位置(x0,y0,z0)T为原点建立局部坐标系,局部坐标系的坐标轴指向矩阵U为单位阵I,得到目标在局部坐标系下的坐标位置为
2.仿真结果分析
从图3可以看出,对于三个等时差双曲面只有一个交点对应的定位模型,三个等时差双叶双曲面仅形成一个交点,此交点就是目标位置,此时在定位区域内任选一初值均可通过该定位方法得到目标位置。
从图5可以看出,此种情况下时差测量误差为10纳秒时,整个定位范围中的定位误差不超过20米,最高定位精度在定位范围中心,可达2.4米。从图7可以看出在三个等时差双曲面只有一个交点时,本发明在目标位置的x,y,z分量上估计均方根误差都比Chan算法的小,并且在时差测量误差较小时能达到克拉美罗下界CRLB,这说明当时差测量误差较大时,本发明优于Chan算法,可以达到更小的定位误差。
从图4可以看出,对于三个等时差双曲面有两个交点对应的定位模型,三个等时差双叶双曲面形成两个交点,此时采用本发明的初值选取方法选取初值然后通过该方法可获得目标位置。
从图6可以看出,此种情况下时差测量误差为10纳秒时,整个定位范围中的定位误差不超过160米,各个辅站附近的定位误差在30米左右,最高定位精度在定位范围中心附近,可达8.2米。
从图8可以看出,在三个等时差双曲面有两个交点时,本方法在目标位置的x,y,z分量上估计均方根误差都比Chan算法的小,并且在时差测量误差较小时能达到克拉美罗下界CRLB,这说明当时差测量误差较大时,本发明的目标定位算法优于Chan算法,可以达到更小的定位误差。
综合以上分析过程可以得出如下结论:应用本发明可以实现空中运动平台对目标进行高精度的三维无源定位,该方法不存在定位模糊,具有较好的实时性和鲁棒性。
Claims (3)
1.一种空中运动平台对目标的三维无源定位方法,其特征在于:包括如下过程:
1)对接收目标散射回波的主站和各个辅站进行时间同步,给目标散射回波到达时刻的测量提供统一的时间基准,所述主站为空中运动平台,所述辅站为地面上的三个接收站;
2)雷达发射站辐射目标检测所需的电磁波;
3)主站测量并记录雷达发射站直达波的到达时刻TD0;
4)主站记录目标散射回波的到达时刻t0,各个辅站记录主站或目标的散射回波的到达时刻si,i=1,2,3;
5)辅站将主站或目标的散射回波的到达时刻si,i=1,2,3连同辅站在全局坐标系下的位置(xi,yi,zi)T进行编码和调制,并通过无线信道上传到主站,i=1,2,3,其中()T表示转置;
6)主站对辅站上传的数据解调和解码后,剔除主站自身产生的散射回波到达时刻tDi,i=1,2,3,并对主站自身进行定位;
7)主站结合目标散射回波到达主站的时刻t0和辅站上传的目标散射回波到达辅站的时刻ti,i=1,2,3对目标进行定位,即将求解目标位置的非线性方程组转化为与其对应的最优化问题;
8)通过合理的选取求解初值对上述最优化问题进行求解,其最优解即为目标的定位结果
9)将目标在全局坐标系下的定位结果转换为以主站为坐标原点的局部坐标系下的定位结果
2.根据权利要求1所述的空中运动平台对目标的三维无源定位方法,其中步骤7)所述的将求解目标位置的非线性方程组转化为与其对应的最优化问题,按如下步骤进行:
7a)将目标到各个辅站的距离Ri,i=1,2,3,目标到主站的距离R0,目标散射回波到主站的到达时刻t0及目标散射回波到各个辅站的到达时刻ti,i=1,2,3代入如下三个双叶双曲面方程组:
式中ei是零均值的量测噪声,i=1,2,3,(x0,y0,z0)T是主站的坐标位置,c是光速;
7b)根据最小均方误差准则,将所述三个双叶双曲面方程组的求解转化为如下最优化问题的求解:
式(2)中J(u)是代价函数,g(u)=[g1,g2,g3]T为辅助函数,u=(x,y,z)T为目标真实位置,为目标真实位置u的估计值,为目标位置x分量的估计值,为目标位置y分量的估计值,为目标位置z分量的估计值,辅助函数g(u)的各个分量为gi=Ri-R0-c(ti-t0)-ei,分别为目标到各个辅站的距离,为目标到主站的距离。
3.根据权利要求2所述空中运动平台对目标的三维无源定位方法,其中步骤8)所述的通过合理的选取求解初值对上述最优化问题进行求解,按如下步骤进行:
8a)选取迭代初值,即当目标位置高于主站时,迭代初值u0=[x0,y0,z0]T应该满足初值的z轴分量z0大于主站的z轴分量z0;当目标低于主站时,初值u0应满足初值的z轴分量z0小于主站的z轴分量z0;
8b)将初值u0带入式(2),计算代价函数值J(u0);
8c)对于精度ε>0,若存在J(u0)<ε,则将u0作为目标位置的估计值否则令迭代序号k=0继续执行下面的步骤;
8d)根据式(3)计算辅助函数g(u)在第k步目标位置估计值uk处的雅可比矩阵F(uk)
其中uk=[xk,yk,zk]T,xk是第k步目标位置估计值uk的x轴分量,yk是uk的y轴分量,zk是uk的z轴分量;
8e)根据式(4)计算代价函数J在第k步目标位置估计值uk处的梯度▽J(uk);
8f)根据式(5)迭代计算第k+1步目标位置估计值uk+1
其中为代价函数J在第k步目标位置估计值uk处的梯度,参数μk+1为自适应时变步长;
其中α是一个加速因子,其取值范围为1≤α≤2,||||表示向量的2范数;
8g)置k←k+1,即将迭代序号k加1,重复执行8d)到8g)直到代价函数J满足精度要求J(uk)<ε,将该uk作为目标位置的估计值
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210476641.5A CN102981144B (zh) | 2012-11-21 | 2012-11-21 | 空中运动平台对目标的三维无源定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210476641.5A CN102981144B (zh) | 2012-11-21 | 2012-11-21 | 空中运动平台对目标的三维无源定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102981144A CN102981144A (zh) | 2013-03-20 |
CN102981144B true CN102981144B (zh) | 2015-03-04 |
Family
ID=47855349
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210476641.5A Expired - Fee Related CN102981144B (zh) | 2012-11-21 | 2012-11-21 | 空中运动平台对目标的三维无源定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102981144B (zh) |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106019348A (zh) * | 2015-08-25 | 2016-10-12 | 中国运载火箭技术研究院 | 一种基于多滞空平台组网的高精度定位装置及方法 |
CN105509665B (zh) * | 2015-11-30 | 2018-01-12 | 公安部天津消防研究所 | 一种基于超声波原理的测点空间位移测量装置及方法 |
CN106371078B (zh) * | 2016-08-16 | 2018-06-26 | 电子科技大学 | 基于发射波形和位置联合估计的多发射源被动定位方法 |
CN106501767A (zh) * | 2016-10-13 | 2017-03-15 | 哈尔滨工程大学 | 一种运动多站无源时差定位方法 |
CN106772237B (zh) * | 2016-11-18 | 2019-06-21 | 中国船舶重工集团公司第七二四研究所 | 长基线时差定位干涉仪测向辅助解模糊方法 |
CN106535217B (zh) * | 2016-12-15 | 2020-02-28 | 中铁二院工程集团有限责任公司 | 基于gdop分析的gsm-r干扰源融合定位方法 |
CN107294548B (zh) * | 2017-06-06 | 2019-09-17 | 武汉船舶通信研究所(中国船舶重工集团公司第七二二研究所) | 一种罗兰c发射机主电路及其发射波形生成方法 |
CN107247279B (zh) * | 2017-07-12 | 2019-03-29 | 中国航天科工集团八五一一研究所 | 存在站址误差下的时差体制定位校正方法 |
CN108375766A (zh) * | 2018-02-08 | 2018-08-07 | 四川九洲电器集团有限责任公司 | 一种基于协同识别高精度定位融合方法 |
CN108717184B (zh) * | 2018-04-27 | 2020-09-01 | 杭州电子科技大学 | 基于误差校正的联合doa与toa单站无源定位方法 |
CN108919185A (zh) * | 2018-07-12 | 2018-11-30 | 中国船舶重工集团公司第七二四研究所 | 非对称布站多元微分时差定位误差分析方法 |
CN111505573B (zh) * | 2020-03-18 | 2022-09-13 | 中国民用航空总局第二研究所 | 一种分布式定位系统的航迹生成方法及装置 |
CN111352444A (zh) * | 2020-04-23 | 2020-06-30 | 上海交通大学 | 基于无线导航的无人机室外移动平台自主降落方法及系统 |
CN112114341B (zh) * | 2020-08-13 | 2022-01-25 | 中国人民解放军军事科学院国防科技创新研究院 | 低轨卫星协同测频无源定位方法 |
CN112083406B (zh) * | 2020-09-17 | 2022-05-17 | 电子科技大学 | 一种对目标高度稳健的外辐射源二维定位方法 |
CN112083407B (zh) * | 2020-09-17 | 2022-03-29 | 电子科技大学 | 一种利用时差与一维方位测量的外辐射源三维定位方法 |
CN112327248B (zh) * | 2020-11-05 | 2021-07-30 | 中国人民解放军32802部队 | 一种窄带连续波信号多站定位与跟踪方法 |
CN112327251B (zh) * | 2020-11-05 | 2021-07-30 | 中国人民解放军32802部队 | 低速率通信条件下的无人值守式信号时差定位方法 |
CN112954591B (zh) * | 2021-02-10 | 2022-04-15 | 北京理工大学 | 一种协同分布式定位方法及系统 |
CN113176606B (zh) * | 2021-06-02 | 2023-09-26 | 中国恩菲工程技术有限公司 | 微震震源定位方法、系统、设备及存储介质 |
CN114994603A (zh) * | 2022-03-17 | 2022-09-02 | 西安电子科技大学 | 一种长波定位方法、系统、存储介质、计算机设备及终端 |
CN114779168A (zh) * | 2022-04-21 | 2022-07-22 | 上海炉石信息科技有限公司 | 一种设备阵列的自动定位方法、装置、存储介质 |
CN115436905B (zh) * | 2022-09-30 | 2023-07-18 | 中国科学院空天信息创新研究院 | 基于无源接收系统的基线估计方法、系统及电子设备 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102608569A (zh) * | 2012-03-29 | 2012-07-25 | 清华大学 | 基于双观测点的空间匹配无源定位方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
IT1398976B1 (it) * | 2009-06-19 | 2013-03-28 | Sie Soc It Elettronica | Metodo e apparecchiatura per la localizzazione passiva di emittenti di segnali radio. |
-
2012
- 2012-11-21 CN CN201210476641.5A patent/CN102981144B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102608569A (zh) * | 2012-03-29 | 2012-07-25 | 清华大学 | 基于双观测点的空间匹配无源定位方法 |
Non-Patent Citations (2)
Title |
---|
"Accurate Passive Location Estimation Using TOA Measurements";Junyang Shen et al.;《IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS》;20120630;第11卷(第6期);第2182-2192页 * |
"布站形式对TDOA无源区域定位系统定位精度的影响";孙宝国 等;《火力与指挥控制》;20110930;第36卷(第9期);第129-132页 * |
Also Published As
Publication number | Publication date |
---|---|
CN102981144A (zh) | 2013-03-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102981144B (zh) | 空中运动平台对目标的三维无源定位方法 | |
US7170441B2 (en) | Target localization using TDOA distributed antenna | |
Johnson et al. | Ultra-wideband aiding of GPS for quick deployment of anchors in a GPS-denied ad-hoc sensor tracking and communication system | |
CN103308934B (zh) | 一种利用wifi反射信号实现室内移动人员定位的方法 | |
CN108351400A (zh) | 用于超宽带位置定位的系统和方法 | |
CN204166130U (zh) | 射频定位装置和系统 | |
CN103995250B (zh) | 射频标签轨迹追踪方法 | |
Kloeden et al. | Vehicle localization using cooperative RF-based landmarks | |
CN103097906A (zh) | 用于确定移动式发射器位置和方向的方法和装置 | |
CN105572635B (zh) | 基于最小二乘法的单站无源快速定位方法 | |
CN103412286B (zh) | 基于mimo雷达的发射极化优化doa估计方法 | |
CN104849737A (zh) | 一种卫星定位系统和定位方法 | |
US20210011109A1 (en) | Method and apparatus for determining the direction of arrival of radio or acoustic signals, and for transmitting directional radio or acoustic signals | |
US6762721B2 (en) | Urban terrain geolocation system | |
CN101446634A (zh) | 一种高精度位置、方位角和俯仰角的组合测量方法及装置 | |
CN104537202B (zh) | 基于卫星编队协作的空间天线阵列合成方法 | |
Fokin | Passive geolocation with unmanned aerial vehicles using TDOA-AOA measurement processing | |
CN105487072A (zh) | 一种基于t2/r的时差与多普勒频移联合定位方法及系统 | |
CN102879762A (zh) | 基于射频接收信号强度值的隧道内车辆的动态定位方法 | |
CN110187333A (zh) | 一种基于合成孔径雷达技术的rfid标签定位方法 | |
CN104459675A (zh) | 一种基于测距的物体定位和跟踪方法及使用该方法的定位设备 | |
CN105487094A (zh) | 一种数据链与卫星导航协同定位方法及定位系统 | |
CN108732534A (zh) | 一种基于加权mds的多标签协作定位方法 | |
CN101858977A (zh) | 基于双红外系统的室内空间定位方法与系统 | |
Lee et al. | Simulation of UWB radar-based positioning performance for a UAV in an urban area |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150304 Termination date: 20211121 |
|
CF01 | Termination of patent right due to non-payment of annual fee |