CN104316903B - 一种三站时差定位性能试验评估方法 - Google Patents
一种三站时差定位性能试验评估方法 Download PDFInfo
- Publication number
- CN104316903B CN104316903B CN201410380137.4A CN201410380137A CN104316903B CN 104316903 B CN104316903 B CN 104316903B CN 201410380137 A CN201410380137 A CN 201410380137A CN 104316903 B CN104316903 B CN 104316903B
- Authority
- CN
- China
- Prior art keywords
- positioning
- station
- radiation source
- error
- time difference
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明属于无源探测技术领域,公开一种三站时差定位性能试验评估方法,是建立三站时差定位性能试验评估模型,设计时差定位***的部署和机载辐射源或机载雷达信号模拟器的飞行航线,然后基于整条航线的定位误差试验结果,解算综合时间测量误差模型中未知参数,最后利用综合时间测量误差模型和定位误差GDOP模型,得到任意位置的定位误差或任意平面的定位误差,本发明给出了综合时间测量误差模型中相关参数的求解方法,及任意布站模式和任意辐射源位置的定位误差。解决了背景技术中的三站时差定位性能试验评估问题。弥补了理想或固定的综合时间测量误差模型的不足。还能够推广到不需辐射源挂飞的地面试验模式。
Description
技术领域
本发明属于无源探测技术领域,尤其适用于一种三站时差定位性能试验评估方法。
背景技术
时差定位***能够对空间辐射源进行侦察定位,具有侦察距离远、抗干扰能力强、定位精度高和反应速度快等特点,因此被广泛应用到导航、航空、航天和电子战等领域。时差定位***由多个基站组成,包括一个主站和若干副站,主站是时差定位***信息处理中心,在无约束情况下,要实现任意空间辐射源的三维无源定位至少需要四站。辐射源信号同时被两个不同位置的接收机接收,时差相同的点在平面上为双曲线,在三维空间为双曲面,多站时差定位通过求解曲线或曲面的交点获取辐射源位置信息。三站时差定位主站和副站侦察接收到空间辐射源信号,并获取侦察信号时差数据,然后主站通过求解两条曲线的交点获取辐射源位置信息。由于三站时差定位***结构简单,因此得到广泛应用。
影响时差定位性能的因素包括基站空间位置、站址测量误差、定位模型误差、时间测量误差等,其中时间测量误差是影响时差定位***定位性能的关键因素之一。如何根据试验数据评估多站时差定位性能是一个紧迫的课题。由于试验时间和经费的限制,试验不可能直接获取多站时差定位***对空间任意一点的定位性能,时差定位性能试验评估目的就是根据试验航线的外场时差定位数据,推算辐射源任意位置的定位精度。为了评估时差定位***的定位性能,通常采用时差定位***地面布站方式,机载辐射源按照预定航线飞行,通过试验数据处理得到任意位置的定位误差或任意平面的定位误差几何稀释GDOP,下面用定位误差GDOP表示定位误差几何稀释。
现有技术一的技术方案
文献[1]:北京国防工业出版社,1996.孙仲康,周一宇,何黎星著的单多基地有源无源定位技术[M].给出了三站时差定位性能试验评估方法。
现有技术一的缺点
常规时差定位误差模型通常将时差测量误差假定为常数[1][2],文献[2]给出了三站时差定位直线布站情况下的定位精度推算模型,但没有考虑时差测量误差随侦察信号信噪比的变化,因此评估方法存在模型误差缺陷。时差测量误差模型是时差定位***GDOP计算的基础,文献[3]给出了时间/频率测量误差和信噪比模型、文献[4]对时间测量误差模型进行了完善,文献[5][6][7][8]也给出了SNR模型,以上文献认为采样时间段内信号是连续的,没有考虑信号占空比对SNR的影响,例如对于雷达信号,相干处理时间内可能存在一个或多个雷达脉冲情况。实际时间测量误差方差和辐射源***参数及信号参数、侦察设备***参数、两者距离等因素有关。***参数确定后,距离远近决定了时间测量能力和测量误差。无限远处侦察设备不能侦察到信号,时差测量无从谈起,随着距离的减小,逐渐能侦察到信号且时差测量误差逐渐减小。
文献[2]:电子学报,2004,32(9):1452-1455.陈永光,李昌锦,李修和著的三站时差定位的精度分析与推算模型[J]。给出了三站时差定位直线布站情况下的定位精度推算模型,采用了固定时差测量误差模型或称为固定系数模型,没有考虑时差测量误差随侦察信噪比的变化,即随着基站与信号源距离变化,信噪比也变化,最终导致时间测量误差方差的变化,因此评估方法存在模型误差缺陷。
发明内容
为了克服背景技术中的不足,本发明提供一种三站时差定位性能试验评估方法,是针对平面三站协同时差定位性能试验评估需求,综合考虑时间测量误差模型,给出了多站时差定位性能试验评估方法,将解决试验过程中时差定位性能评估问题,同时给出了三站时差定位性能试验评估方法。
为了实现上述发明目的,本发明采用技术方案如下:
一种三站时差定位性能试验评估方法,是建立三站时差定位性能试验评估模型,并且通过三站时差定位性能试验评估模型,设计外场时差定位***的部署方法和机载辐射源的飞行航线,然后基于整条航线的定位误差试验结果,解算综合时间测量误差模型中未知参数,最后利用综合时间测量误差模型和定位误差模型,得到任意位置的定位误差或任意平面的定位误差GDOP,其具体步骤如下:
步骤一试验布站设计
按照直线、三角形、T形或任意空间位置部署三站时差定位***的定位站空间位置,根据定位主站和副站间通讯能力和定位***确定最小基线长度,基线长度要满足***的站间通信距离以及发射站的信号辐射功率、波束宽度条件,同时满足最小基线长度和站间的通讯;另外定位站的选择也要满足实际空域试验的飞行航线;
步骤二定位误差的归一化GDOP
根据步骤一的定位站部署位置,按照三站时差定位***的定位方法,计算机载辐射源飞行高度H或飞行器载辐射源的巡航高度的平面的归一化定位误差GDOP;
设置三个时差定位站和辐射源位于同一平面,基站和辐射源坐标用二维坐标表示,基站位置分别为R0(x0,y0),R1(x1,y1)和R2(x2,y2),其中R0(x0,y0)为主站坐标,辐射源位置为T(x,y);
辐射源信号到各基站的时刻分别为t0,t1和t2,时差测量表达式为:
式中:c为光速,辐射源信号到达第i副站与到达主站的距离差di=c(ti-t0),i=1,2;通过求解该方程组可得到辐射源位置;
(1)式两边求微分,得到
其中
写成矩阵形式
cdT=AdX+dXs (3)
其中dT=[dt1-dt0 dt2-dt0]T,dX=[dx dy]T,dXs=[k1+k0 k2+k0]T;
利用(3)式,解得定位误差为:
dX=A-1(cdT-dXs) (4)
进而得到定位误差方差为:
PdX=E[dX(dX)T]=A-1Pε[A-1]T (5)
Pε=E[(cdT-dXs)(cdT-dXs)T] (6)
=c2E[dTdTT]+E[dXs(dXs)T]
其中Pε为综合时间测量误差方差矩阵;E[]为求均值函数(即求数学期望);
设置各测量误差是零均值且不相关的高斯白噪声,基站位置测量和时间测量之间相互独立,站址测量误差方差时间测量误差方差为:
其中i=0,1,2,由于受基站和辐射源天线动态扫描、两者之间的距离、辐射源发射信号样式和功率资源调度等因素的影响,侦察信号强度是时变,侦察信号的互相关输出也是时变的,时差测量精度是一个动态过程,因此是不同的,而且均随侦察距离变化而变化。常规定位误差模型将时间测量误差方差假定为常数[2],即称之为固定时差测量模型。
通常站址测量误差值满足三维正态分布,即利用得到
其中In为n阶单位矩阵,En为n阶方阵,n=2。对于固定时差测量模型,
综合时间测量误差方差矩阵的迹为:
三站时差定位***中,参考附件A(A8式)多个基站的综合时间测量误差方差可近似表示为综合时间测量误差方差矩阵迹的形式,即:
其中E[]为求均值函数,下标m表示均值,σPt=σs/c为基站的站址测量误差(已经转换为时间);
平面(XY)三站时差定位误差GDOP为:
用单个基站综合时间测量误差对GDOP归一化,得到归一化GDOP为:
其中归一化GDOP为时的定位误差。
假设可以得到从而得到归一化GDOP为:
证明归一化GDOP与基站和辐射源位置有关,与综合时间测量误差方差矩阵Pε无关;
步骤三机载辐射源的直线航线设计
利用步骤1,2获取巡航高度H平面的归一化定位误差几何稀释设置直线飞行航线,选取的航线变化趋势平缓,用于避免通过定位误差梯度变换比较快的位置,以便分段统计定位误差;同时当目标飞机沿设计的航线飞行时,要满足在有效航线上辐射源信号始终能被定位基站接收到;包括:航路最远点、航线长度、航线高度以及有效飞行航次;
所述航线最远点能够满足飞行空域限制、飞机有效航程、最远侦察距离因素的限制;所述航线高度按照飞机巡航高度,试验航线最近点位置相对最近基站的俯仰角小于三站时差定位***俯仰覆盖的范围;
飞行架次的确定:飞行架次根据距离单元的大小进行推算,首先计算所需的飞行航次数
式中:Fn为试验航次,△R为航线距离取样区间大小(m),N为航线距离取样区间△R内需统计的定位数据总量,通常为60,V为目标机飞行速度(m/s),T为定位***数据录取周期(s)。
步骤四试验数据的获取与处理
获取试验数据的数据量,根据统计数据量设置目标飞机一架沿航线水平、直线、往返飞行的航线往返次数Fn;
当机载辐射源开机,被试装备主站、副站天线对准预定航线进行扇扫,分别对机载辐射源的目标信号进行侦收;副站将数据传送到主站,主站进行数据处理,用相关处理法得到主站和所有副站的侦察信号时差数据,然后主站按照***定位算法联合处理时差信息,获取空间辐射源的位置;
标准位置测量设备为机载GPS***或精密测量雷达,全航路跟踪机载辐射源目标,获取辐射源真实位置;
当机载辐射源不在定位平面时,三站时差定位存在***模型误差,由于时间测量误差和***模型误差是相互独立的,因此实际三站时差定位***的理论定位误差模型为:
其中为辐射源高度引入的定位模型误差方差;
以三个基站为平面建立定位平面(XY),其中主站坐标为(x0,y0,0),两副站坐标分别为(x1,y1,0)和(x2,y2,0),辐射源坐标为(x,y,z),解算出实际辐射源到各定位基站的距离或时延,为:
然后根据三站定位算法解算出目标位置测量值(x′,y′,0);
由目标的真实位置和测量位置得到辐射源高度引起的定位模型误差方差为:
由(13,16)式,得到综合时间测量误差方差的试验统计值为:
其中c为光速,σX 2由三维定位误差试验结果统计得到,不同航线位置σX 2不同的,为基站和辐射源坐标位置确定的定位模型误差方差,GDOPe为基站和辐射源坐标位置确定的归一化误差。
用相关处理方法,计算综合时间测量误差方差
其中σPt=σs/c为站址测量误差,B为辐射信号带宽,τ=1/B,TC为相参处理时间长度,fs为采样率,Ri为辐射源距离第i个基站的距离,和K0为固定常数,同等信号参数情况下,K0越小,距离产生的影响越大;
当辐射源***参数和各基站***参数确定后,K0为常数;另外由于三站时差定位各主副站间距一般比较小,而侦察天线波束比较宽,在远距离各站的参数Ki基本相同,即:
其中Pt为辐射源峰值功率,Gt(φi)为辐射源在第i个侦察接收机方向的增益,Gi(θi)为第i个侦察接收机的接收增益,λ为雷达工作波长,Ljr为侦察接收机功率接收损耗,Lt为发射机***损耗,为第i个接收机噪声功率,NP为侦察脉冲数目,Tr为辐射信号脉冲宽度,fs为采样率;
三站时差定位性能的评估,就是通过试验结果求取参数和K0的值,通过(19,20)式,建立方程:
其中c为光速,σX 2为试验获取的定位误差方差,GDOPe为由辐射源位置和定位站位置得到的归一化误差,和K0为固定常数,和由辐射源位置和定位站位置得到,TC为相关处理长度,fs为采样率,B为信号带宽,为的估计数值,通过试验能够获取综合时间测量误差方差的估计值f;
综合时间测量误差的方差f分布是非平稳的,采用等距离段分段统计,判断距离段内是平稳随机信号;
设置第m个距离段得到L个定位实测结果Xi,对应真实位置为X0i,以及L个定位点对应的辐射源和定位站位置;
通过(23)式、辐射源和定位站位置直接计算gm和qm参数均值;对分段数据f用3σ准则剔除异常误差,统计得到各距离段内fm;第m个距离段的综合时间测量误差方差fm为:
其中,L为距离段内定位点数,||||为向量范数或取向量长度,为第i个辐射源位置对应的定位模型误差方差,GDOPei为第i个辐射源位置归一化定位误差;
由(22)式得到第m个距离段的综合时间测量误差方差为:
其中m为1~M;得到M个位置的综合时间测量误差方差fm,及参数gm和qm;
步骤五参数化求解过程
将方程(25)简化为:
fm=a+bxm+b2ym (26)
其中
按照最小二乘法,定义:
判断a,b最优估计的问题,就是求解Q(a,b)最小值的问题,Q(a,b)分别对a,b求偏导,并令它们等于零,解方程组就能够得到a,b的估计数值,即
解(28)式可得到b的方程为:
U(b)=Sxf+Syfxxb+Sxyb2+Syyb3=0 (29)
根据三次方程求解方法得到该方程的解析解,或采用MATLAB中roots函数多项式求根方法;再根据辐射源和时差定位***参数,估计出K0的数量级或近似数值范围;
根据b的取值范围,可得到b的一个有效解,解得b,代入下式得到a为:
按照最小二乘法得到待估计参数和K0;
或通过曲线插值方法获取两位置的综合时间测量误差方差参数fC,fD,由基站和辐射源位置,按照(23)式得到C和D两个位置参数gC,qC的gD,qD,通过求二元方程组得到估计综合时间测量误差参数和K0;
步骤六任意布站模式下定位性能
三站时差定位存在***模型误差,由于时间测量噪声误差和***模型误差是相互独立的,则三站时差定位***定位精度的任意航线或任意X位置的定位误差,即综合时间测量误差模型定位误差为:
其中(GDOPe)X由步骤二(14)式获取,由步骤五(22)式获取,由(18)式获取。
一种三站时差定位性能试验评估方法,所述三站时差定位性能试验评估模型,包括:建立综合时间测量误差模型和定位误差模型,所述综合时间测量误差模型由信噪比引起的时差测量误差和综合时间测量误差组成,综合时间测量误差包括基站站址测量误差和***固有时间测量误差。
一种三站时差定位性能试验评估方法,所述飞行辐射源设置为在两固定位置上的固定辐射源;固定辐射源在地面上支撑两个杆的固定位置T0、T1上,通过两个杆的固定位置上辐射源在两固定位置T0、T1的多次测试结果,计算辐射源任意位置的定位精度;通过试验获取两固定位置的综合时间测量误差方差参数为fC,fD,由基站和辐射源位置,按照(23)式得到C和D两个位置参数gC,qC的gD,qD,通过求二元方程组得到估计综合时间测量误差参数和K0,进而按照步骤六得到任意布站模式和任意辐射源位置的定位误差。
由于采用如上所述的技术方案,本发明技术方案具有如下优越性:
一种三站时差定位性能试验评估方法,是针对时差定位性能试验评估需求,建立了综合时间测量误差模型和定位误差模型,给出了多站时差定位性能试验评估方法。即基于外场固定布站模式情况下单条飞行航线的多站时差定位试验数据,给出了综合时间测量误差模型中相关参数的求解方法,推算定位***的由信噪比引起的时差测量误差和其他综合时间测量误差,并给出了任意布站模式和任意辐射源位置的定位误差。
本发明给出了三站时差定位性能试验评估方法,是根据单条航线的外场时差定位数据,推算辐射源任意位置的定位精度。采用三站时差定位试验数据处理和时差定位性能评估,解决了背景技术中的三站时差定位性能试验评估问题。
本发明还能够推广到不需辐射源挂飞的地面试验模式,通过辐射源在两固定位置的多次测试结果,推算辐射源任意位置的定位精度。
本发明的创新点基于信噪比的综合时间测量误差模型和时差定位性能试验评估模型,弥补了理想或固定的综合时间测量误差模型的不足。本发明给出时差定位***对试验辐射源的定位性能,相关研究可以应用到辐射源替代推算理论中,为拓展等效替代理论奠定基础。
附图说明
图1三站时差定位性能试验评估模型流程图;
图2三站时差定位接收机与辐射源集合关系示意图;
图3归一化GDOP与航线设计示意图;
图4a定位误差测量结果图;图4b定位误差分段统计结果图;图4c综合时间测量误差方差分段统计结果及其统计结果图;
图5最小二乘综合时间测量误差方差估计曲线与理论曲线对比图;
图6a模型误差分布图;图6b***误差分布图;图6c理论综合定位误差图;6d综合时差测量误差模型法定位误差/试验结果图;单位km;
图7固定综合时间测量误差情况下的GDOP图;
图8不同时差测量误差模型对定位误差的影响对比图;
图9三站时差定位接收机与两个固定辐射源集合关系示意图。
具体实施方案
由于试验时间和经费的限制,试验不可能直接获取多站时差定位***对空间任意一点的定位性能。为了评估时差定位***的定位性能,通常采用被试时差定位***地面布站方式,机载辐射源按照预定航线飞行,通过试验数据处理得到任意位置的定位误差或任意平面的定位误差GDOP。
时差定位性能试验评估技术的关键是建立综合时间测量误差模型和定位误差模型,综合时间测量误差模型包括由信噪比引起的时差测量误差和其他综合时间测量误差。时差定位性能试验评估技术的过程是首先设计外场时差定位***部署方案和机载辐射源飞行航线,然后基于整条航线的定位误差试验结果,解算综合时间测量误差模型中未知参数,最后利用综合时间测量误差模型和定位误差GDOP模型,得到任意位置的定位误差或任意平面的定位误差GDOP。时差定位性能试验评估模型及流程图如图1,步骤如下:
步骤一试验布站设计
根据定位主站和副站间通讯能力和定位***确定最小基线长度,多站空间位置部署可以按照直线、三角形、T形或任意空间位置布站。基线长度要考虑***的站间通信距离以及发射站的信号辐射功率、波束宽度等条件,同时满足最小基线长度和站间通讯能力要求。另外阵地的选择也要兼顾试验飞行航线的设计,满足实际空域要求。
步骤二定位误差的归一化GDOP
根据步骤一的定位站部署位置,按照三站时差定位***的定位方法计算机载辐射源飞行高度H,即飞行器巡航高度,平面的归一化定位误差GDOP。
假设三个时差定位站的位置和辐射源位于同一平面,基站和辐射源坐标可以用二维坐标表示,基站位置分别为R0(x0,y0),R1(x1,y1)和R2(x2,y2),其中R0(x0,y0)为主站坐标,辐射源位置为T(x,y),如图2所示。
辐射源信号到各基站的时间分别为t0,t1和t2,时差测量表达式为:
式中:c为光速,辐射源信号到达第i副站与到达主站的距离差di=c(ti-t0),i=1,2。通过求解该方程组可得到辐射源位置。
(1)式两边求微分,得到
其中
写成矩阵形式
cdT=AdX+dXs (3)
其中dT=[dt1-dt0 dt2-dt0]T,dX=[dx dy]T,dXs=[k1+k0 k2+k0]T。
利用(3)式,解得定位误差为:
dX=A-1(cdT-dXs) (4)
进而得到定位误差方差为:
PdX=E[dX(dX)T]=A-1Pε[A-1]T (5)
Pε=E[(cdT-dXs)(cdT-dXs)T] (6)
=c2E[dTdTT]+E[dXs(dXs)T]
其中Pε为综合时间测量误差方差矩阵。
假设各测量误差是零均值且不相关的高斯白噪声,基站位置测量和时间测量之间相互独立,站址测量误差方差时间测量误差方差为:
其中i=0,1,2,由于受基站和辐射源天线动态扫描、两者之间的距离、辐射源发射信号样式和功率资源调度等因素的影响,侦察信号强度是时变,侦察信号的互相关输出也是时变的,因此时差测量精度是一个动态过程,是不同的,而且均随侦察距离变化而变化。常规定位误差将时间测量误差方差假定为常数[2],即称之为固定时差测量模型。
通常站址测量误差值满足三维正态分布,即利用得到
其中In为n阶单位矩阵,En为n阶方阵,n=2。对于固定时差测量模型,
综合时间测量误差方差矩阵的迹为:
三站时差定位***中,参考附件A(A8式)多个基站的综合时间测量误差方差可近似表示为综合时间测量误差方差矩阵迹的形式,即:
其中σPt=σs/c为基站的站址测量误差(已经转换为时间)。
平面(XY)三站时差定位误差GDOP为:
用单个基站综合时间测量误差对GDOP归一化,得到归一化GDOP为:
其中归一化GDOP为时的定位误差。
假设可以得到从而得到归一化GDOP为:
可以证明归一化GDOP和基站、辐射源位置有关,与综合时间测量误差方差矩阵Pε无关。
举例:假设三站时差定位***直线布站,主站位置为(0,0,0),2个副站位置分别为(-10km,0,0),(10km,0,0)。试验通常设计直线航线,辐射源载机航线高度为H=6km,按照该布站方式,得到6km高度三站时差定位的归一化GDOP。归一化GDOP与航线设计示意图如图3所示,航线起始点为A点,终止点为B点,为了便于分段统计定位误差,因此航线设计时避免通过定位误差梯度变换比较快的位置。
步骤三机载辐射源直线航线设计
利用步骤1,2,获取巡航高度H平面的归一化定位误差GDOP。试验通常设计直线飞行航线,选取的航线变化趋势平缓,航线应该避免通过定位误差梯度变换比较快的位置,以便分段统计定位误差;同时当目标飞机沿设计航线飞行时,要满足在有效航线上辐射源信号始终能被定位基站接收到。
设计的主要内容包括航路最远点、航线长度、航线高度以及有效飞行航次等。航线最远点可以满足飞行空域限制、飞机有效航程、最远侦察距离等因素的限制。航线高度按照飞机巡航高度,试验航线最近点位置相对最近基站的俯仰角应该小于多站时差定位***俯仰覆盖范围要求。
飞行架次的确定:对任一固定航路而言,其任一点的多次试验样本的定位误差分布函数上是平稳随机函数,但全航路的定位误差分布函数是非平稳随机函数。在非平稳信号处理的理论分析和工程实践过程中,通常采用划分间隔区间的方式,确保定位误差在一定范围内满足平稳随机过程的要求。然后对区间内的数据进行数值估计,在满足一定置信度和置信区间的估值要求,每个范围内的试验数据量应达到一定的数量要求,要确保精度统计的正确性和合理性,因此需要计算合适的飞行架次,飞行架次需根据距离单元的大小进行推算,首先计算所需的飞行航次数:
式中:Fn为试验航次,△R为航线距离取样区间大小(m),N为航线距离取样区间△R内需统计的定位数据总量,通常为60,V为目标机飞行速度(m/s),T为定位***数据录取周期(s)。
步骤四试验数据的获取与处理
具体试验方法为目标飞机一架沿预定设计航线水平、直线、往返飞行,根据统计数据量需求,航线往返次数Fn,以满足数据统计需要。机载辐射源开机,被试装备主站、副站天线对准预定航线进行扇扫,分别对目标信号进行侦收。副站将数据传送到主站,主站进行数据处理,用相关处理法得到主站和所有副站的侦察信号时差数据,然后主站按照***定位算法联合处理时差信息,获取空间辐射源的位置。标准位置测量设备为机载GPS***或精密测量雷达,全航路跟踪目标获取辐射源真实位置[10]。
当辐射源不在定位平面时,三站时差定位存在***模型误差。由于时间测量误差和***模型误差是相互独立的,因此实际三站时差定位***的理论定位误差模型为:
其中为辐射源高度引入的定位模型误差方差。
以三个基站为平面建立定位平面(XY),其中主站坐标为(x0,y0,0),两副站坐标分别为(x1,y1,0)和(x2,y2,0),辐射源坐标为(x,y,z),解算出实际辐射源到各定位基站的距离(或时延)为:
然后根据三站定位算法解算出目标位置测量值(x′,y′,0)。
由目标的真实位置和测量位置就可以得到辐射源高度引起的定位模型误差方差为:
由(13,16)式,得到综合时间测量误差方差的试验统计值为:
其中c为光速,σX 2由三维定位误差试验结果统计得到,不同航线位置σX 2不同的,为基站和辐射源坐标位置确定的定位模型误差方差,GDOPe为基站和辐射源坐标位置确定的归一化误差。
相关处理方法,理论上综合时间测量误差方差为:
其中σPt=σs/c为站址测量误差,B为辐射信号带宽,τ=1/B,TC为相参处理时间长度,fs为采样率,Ri为辐射源距离第i个基站的距离,和K0为固定常数,为未知参数。同等信号参数情况下,K0越小,距离产生的影响越大。
辐射源***参数和各基站***参数确定后,统计意义上K0为常数。另外由于多站时差定位各主副站间距一般比较小,而侦察天线波束比较宽,在远距离各站的参数Ki基本相同,即:
其中Pt为辐射源峰值功率,Gt(φi)为辐射源在第i个侦察接收机方向的增益,Gi(θi)为第i个侦察接收机的接收增益,λ为雷达工作波长,Ljr为侦察接收机功率接收损耗,Lt为发射机***损耗,为第i个接收机噪声功率,NP为侦察脉冲数目,Tr为辐射信号脉冲宽度,fs为采样率。
时差定位性能评估的问题,就是通过试验结果求取参数和K0的问题,参考(19,20)式,可以建立方程,即综合时间测量误差模型为:
其中c为光速,σX 2为试验获取的定位误差方差,GDOPe为由辐射源位置和定位站位置得到的归一化误差,和K0为固定常数,和可由辐射源位置和定位站位置得到,TC为相关处理长度,fs为采样率,B为信号带宽,为的估计数值,通过试验可以获取综合时间测量误差方差的估计值f。
综合时间测量误差方差f分布是非平稳的,采用等距离段(例如间隔6km)分段统计可以认为距离段内是平稳随机信号。假设第m个距离段得到L个定位实测结果Xi,对应真实位置为X0i,以及L个定位点对应的辐射源和定位站位置,参考(23)式,利用辐射源和定位站位置可以直接计算gm和qm参数均值;对分段数据f用3σ准则剔除异常误差[11],统计得到各距离段内fm,第m个距离段的综合时间测量误差方差fm为:
其中,L为距离段内定位点数,||||为向量范数或取向量长度,为第i个辐射源位置对应的定位模型误差方差,GDOPei为第i个辐射源位置归一化定位误差。
由(22)式可得到第m个距离段的综合时间测量误差方差为:
其中m为1~M。
仿真表明在小于20km的距离近端,综合时间测量误差的估计数值逐渐增大,这和归一化的定位误差变化增大有关,另外距离远端(超过侦察能力)存在无法定位情况比较多,为了提高统计段误差稳定性,需要剔除距离近端和距离远端定位误差测量数据。
假设通过步骤4,得到M个位置的综合时间测量误差方差fm,参数gm和qm。
举例:布站方式同上,辐射源载机航线高度为6km,按照设计航线沿着Y轴移动,假设时间测量误差为10ns,站址测量误差方差辐射源最远点坐标(0,150km),最近点坐标(0,0),时差测量***对运动辐射源连续时差定位,辐射源每隔100m定位一次,采用6km分段统计,段内测量点为60个。通过仿真得到定位误差及其统计结果,即三站时差定位误差试验结果和综合时间测量误差方差统计结果如图中:4a)为定位误差测量结果;4b)为定位误差分段统计结果,表明定位误差是平面模型误差和***误差的合成,定位误差分段统计结果和理论定位误差很好地吻合,证明了分段统计的可行性;4c)综合时间测量误差方差分段统计结果。表明综合时间测量误差方差分段统计结果和理论值很好地吻合,同样可以证明分段统计方法的可行性。
步骤五参数化求解过程
方程(25)简化为:
fm=a+bxm+b2ym (26)
其中
按照最小二乘法,定义:
a,b最优估计的问题,就是求解Q(a,b)最小值的问题。Q(a,b)分别对a,b求偏导,并令它们等于零,解方程组就可以得到a,b的估计数值,即
解(28)式可得到b的方程为:
U(b)=Sxf+Syfxxb+Sxyb2+Syyb3=0 (29)
参考数学手册,文献高等教育出版社,2010:88-89,王连祥,方德植,张鸣镛著的数学手册[M].根据三次方程求解方法可得到该方程的解析解,或采用MATLAB中roots函数多项式求根方法。由于K0>>0,因此0<b<<1,另外根据辐射源和时差定位***参数,也可估计出K0的数量级或近似数值范围。根据b的取值范围,可得到b的一个有效解。解得b后,代入下式得到a为:
按照最小二乘法可得到待估计参数和K0。
当然,也可以通过曲线插值方法获取两位置的综合时间测量误差方差参数fC,fD,由基站和辐射源位置,按照(23)式得到C和D两个位置参数gC,qC的gD,qD,通过求二元方程组得到参数和K0。
举例:有效航线试验数据分段统计得到M个位置的参数fm测量值,利用已知参数可以得到按照最小二乘法可得到待估计参数和K0。最小二乘法参数估计结果与实际仿真参数对比列表如表1。图5为最小二乘综合时间测量误差方差估计曲线与理论曲线对比图,可见最小二乘参数估计法可估计出综合时差测量误差参数,且估计精度高。实际最小二乘需要分段统计平均,相当于平滑处理,否则误差很大。
仿真表明在小于20km的距离近端,综合时间测量误差的估计数值逐渐增大,这和归一化的定位误差变化增大有关,因此应该利用远距离段进行参数估计。另外距离远端(超过侦察能力)存在无法定位情况比较多,为了提高统计段误差稳定性,需要剔除距离远端定位误差测量数据。
表1 最小二乘法参数估计结果与实际仿真参数对比列表
通过步骤5,得到估计综合时间测量误差参数和K0。
步骤六任意布站模式下定位性能
三站时差定位存在***模型误差,时间测量噪声误差和***模型误差是相互独立的。三站时差定位***定位精度的任意航线或任意X位置的定位误差,即综合时间测量误差模型定位误差为:
其中(GDOPe)X由步骤二获取,由步骤五获取,由(18)式获取。
举例:基于时差定位模型,由站址测量误差、时差测量误差以及各定位主副站相对位置关系,就可以得到任意空间辐射源位置所对应的定位精度。利用最小二乘法得到和K0=4.0336×1013,利用(31)式就可以得到同样直线布站情况下的H=6km平面的GDOP,即理论定位误差和综合时差测量误差模型法的定位误差GDOP,如图中:6a)为模型误差分布,6b)为***误差分布,6c)为理论定位误差GDOP,理论定位误差是模型误差和***误差的合成,6d)为最小二乘法得到定位精度GDOP的估计。理论定位误差GDOP和本发明的综合时差测量误差模型定位误差GDOP重合性很好,证明该方法的有效性。
常规方法认为f=σPt 2+σt0 2+σt1 2为常数,即全程时差测量为常数,利用定位误差仿真试验结果,可以得到利用(31)式得到固定时差定位误差情况下三站时差定位性能。固定综合时间测量误差情况下的GDOP如图7所示,可见理想固定系数模型得到的三站时差定位误差偏离了理论定位误差曲线(图6(c)),存在***误差。
为了对比基于时间测量误差固定系数模型定位误差、本发明综合时间测量误差模型定位误差和理论定位误差性能,仿真了单条航线情况下三者的定位误差GDOP对比结果如图8所示,其中航线为本发明直线试验航线。可见本发明综合时间测量误差模型定位误差和理论定位误差非常接近,而固定系数模型定位误差与理论定位误差相差很大,这体现了用本发明时差定位性能试验评估方法进行定位性能研究的必要性。
本发明还能够推广到不需辐射源挂飞的地面试验模式,通过辐射源在两固定位置的多次测试结果,推算辐射源任意位置的定位精度,如图9所示,所述飞行辐射源设置为在两固定位置上的固定辐射源;固定辐射源在地面上支撑两个杆的固定位置T0、T1上,通过两个杆的固定位置上辐射源在两固定位置T0、T1的多次测试结果,计算辐射源任意位置的定位精度。例如通过试验获取两固定位置的综合时间测量误差方差参数为fC,fD,由基站和辐射源位置,按照(23)式得到C和D两个位置参数gC,qC的gD,qD,通过求二元方程组得到估计综合时间测量误差参数和K0,进而按照步骤六得到任意布站模式和任意辐射源位置的定位误差。
附件A基于信噪比的综合时间测量误差模型
根据多站时差定位理论和模型研究,可知时差定位性能最终影响因素包括布站方式、站址测量误差、定位模型误差、时差测量误差、电波传播扰动误差等。其中时差测量误差和主副站时统误差、***信号采样率、时差测量方法、两站接收信号多普勒和侦察信号强度等因素有关。时差测量误差包括固有时间测量误差和信噪比引起的时间测量误差,其中固有时间测量误差和主副站时统误差、***信号采样率、时差测量方法、两站接收信号多普勒等因素有关,忽略两站接收信号多普勒处理误差,固有时间测量误差近似为常数。信噪比引起的时间测量误差和辐射源***参数、侦察设备***参数、侦察距离等因素有关。
不考虑多路径影响,忽略电波传播扰动误差,两站接收信号多普勒处理误差,将站址测量误差和时间测量误差统称为综合时间测量误差。第i个基站的综合时间测量误差表示为:
其中σPt=σs/c为基站的站址测量误差(已经转换为时间),σ0为***固有时差测量误差(或时统误差),该数值是常数,包括主副站时统误差、***信号采样率和时差测量方法引起的时差测量误差等,σti为第i个基站侦察信号相关信噪比引起的时差测量误差,该数值是动态量。
第i个基站侦察信号相关信噪比引起的时间测量误差为:
式中:τ=1/B为相关输出脉宽,B为信号带宽,SNRC为相关处理信噪比。时差测量条件是相关信噪比SNRc≥13dB。
由侦察方程可知第i个基站接收信号信噪比
其中Pt为辐射源峰值功率,Gt(φi)为辐射源在第i个侦察接收机方向的增益,Gi(θi)为第i个侦察接收机的接收增益,Lt为发射机***损耗,Ljr为侦察接收机功率接收损耗(包括侦察接收机接收损耗Lr,单程大气损耗Ldq及波束损耗LBS等),Ri为侦察接收机和雷达的距离,λ为雷达工作波长;为第i个接收机噪声功率。
参考文献Stein S.Differential delay/Doppler ML estimation with unknownsignals[J].IEEE Trans.on Signal Process,1993,41(8):2717-2719,考虑了对时间测量误差的影响,经过理论分析和仿真验证可知第i和j基站的相关法输出信噪比为:
其中Tr为辐射信号脉冲宽度,NP为侦察脉冲数目,fs为采样率,Ri为辐射源距离第i个定位基站的距离,相参处理时间长度TC,一般两侦察接收机接收信号延迟Td远小于采样信号长度TL,即实际相关处理长度TC=TL-Td≈TL。相关输出信噪比和***参数Ki和Kj,距离Ri和Rj,相关处理长度TC,采样率fs有关,而参数Ki和接收机带宽(采样率)fs、采样脉冲数NP、脉冲宽度Tr、侦察接收机接收信号信噪比有关。
第i个基站的***参数Ki为:
其中Pt为辐射源峰值功率,Gt(φi)为辐射源在第i个侦察接收机方向的增益,Gi(θi)为第i个侦察接收机的接收增益,λ为雷达工作波长,Ljr为侦察接收机功率接收损耗,Lt为发射机***损耗,为第i个接收机噪声功率,NP为侦察脉冲数目,Tr为辐射信号脉冲宽度,fs为采样率。
辐射源***参数和各基站***参数确定后,统计意义上K0为常数。另外由于多站时差定位各主副站间距一般比较小,而侦察天线波束比较宽,在远距离各站的参数Ki基本相同,即:
其中,E[]为求均值函数。通常***参数K0为对辐射源副瓣侦察情况。
各站相关处理均以主站j=0为参考,得到相关处理信噪比为:
参考A1式第i个基站的综合时间测量误差方差为:
其中σPt=σs/c为站址测量误差,c为光速,σ0 2为***固时差测量误差,参数K0=E[Ki],fs为采样率,B为辐射信号带宽,τ=1/B,相参处理时间长度TC,Ri为辐射源距离第i个定位基站的距离。
Claims (3)
1.一种三站时差定位性能试验评估方法,其特征在于:建立三站时差定位性能试验评估模型,并且通过三站时差定位性能试验评估模型,设计时差定位***的部署和机载辐射源或机载雷达信号模拟器的飞行航线,然后基于整条航线的定位误差试验结果,解算综合时间测量误差模型中未知参数,最后利用综合时间测量误差模型和定位误差GDOP模型,得到任意位置的定位误差或任意平面的定位误差GDOP,其具体步骤如下:
步骤一 试验布站设计
按照直线、三角形、T形或任意空间位置部署三站时差定位***的定位站空间位置,根据定位主站和副站间通讯能力和定位***确定最小基线长度,基线长度要满足***的站间通信距离以及发射站的信号辐射功率、波束宽度条件,同时满足最小基线长度和站间的通讯;另外定位站的选择也要满足实际空域试验的飞行航线;
步骤二 定位误差的归一化GDOP
根据步骤一的定位站部署位置,按照三站时差定位***的定位方法,计算机载辐射源飞行高度H或飞行器载辐射源的巡航高度平面的归一化定位误差GDOP;
设置三个时差定位站的位置和辐射源位于同一平面,基站和辐射源坐标用二维坐标表示,基站位置分别为R0(x0,y0),R1(x1,y1)和R2(x2,y2),其中R0(x0,y0)为主站坐标,辐射源位置为T(x,y);
辐射源信号到各基站的时间分别为t0,t1和t2,时差测量表达式为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mi>c</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mi>c</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:c为光速,辐射源信号到达第i副站与到达主站的距离差di=c(ti-t0),i=1,2;通过求解该方程组可得到辐射源位置;
(1)式两边求微分,得到
<mrow>
<mi>c</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>dt</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>dt</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>x</mi>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>y</mi>
<mo>+</mo>
<msub>
<mi>k</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中
写成矩阵形式
cdT=AdX+dXs (3)
其中dT=[dt1-dt0 dt2-dt0]T,dX=[dx dy]T,
dXs=[k1+k0 k2+k0]T;
利用(3)式,解得定位误差为:
dX=A-1(cdT-dXs) (4)
进而得到定位误差方差为:
PdX=E[dX(dX)T]=A-1Pε[A-1]T (5)
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mi>&epsiv;</mi>
</msub>
<mo>=</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>c</mi>
<mi>d</mi>
<mi>T</mi>
<mo>-</mo>
<msub>
<mi>dX</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>c</mi>
<mi>d</mi>
<mi>T</mi>
<mo>-</mo>
<msub>
<mi>dX</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<mi>E</mi>
<mo>&lsqb;</mo>
<msup>
<mi>dTdT</mi>
<mi>T</mi>
</msup>
<mo>&rsqb;</mo>
<mo>+</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<msub>
<mi>dX</mi>
<mi>s</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>dX</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
其中Pε为综合时间测量误差方差矩阵;E[]为求均值函数;
设置各测量误差是零均值且不相关的高斯白噪声,基站位置测量和时间测量之间相互独立,站址测量误差时间测量误差方差为:
<mrow>
<mi>E</mi>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>dt</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>=</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>T</mi>
<mi>i</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中i=0,1,2,由于受基站和辐射源天线动态扫描、两者之间的距离、辐射源发射信号样式和功率资源调度因素的影响,侦察信号强度是时变的,侦察信号的互相关输出也是时变的,因此时差测量精度是一个动态过程,因此是不同的,而且均随侦察距离变化而变化;常规定位误差模型将时间测量误差方差假定为常数,即称之为固定时差测量模型;
通常站址测量误差值满足三维正态分布,即利用得到
<mrow>
<mi>E</mi>
<mo>&lsqb;</mo>
<msub>
<mi>dX</mi>
<mi>s</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>dX</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>&rsqb;</mo>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>2</mn>
<msubsup>
<mi>&sigma;</mi>
<mi>s</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
<mtd>
<msubsup>
<mi>&sigma;</mi>
<mi>s</mi>
<mn>2</mn>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>&sigma;</mi>
<mi>s</mi>
<mn>2</mn>
</msubsup>
</mtd>
<mtd>
<mrow>
<mn>2</mn>
<msubsup>
<mi>&sigma;</mi>
<mi>s</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>s</mi>
<mn>2</mn>
</msubsup>
<msub>
<mi>I</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>s</mi>
<mn>2</mn>
</msubsup>
<msub>
<mi>E</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>E</mi>
<mo>&lsqb;</mo>
<msup>
<mi>dTdT</mi>
<mi>T</mi>
</msup>
<mo>&rsqb;</mo>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>T</mi>
<mn>0</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>T</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
<mtd>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>T</mi>
<mn>0</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>T</mi>
<mn>0</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>T</mi>
<mn>0</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>T</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
其中In为n阶单位矩阵,En为n阶方阵,n=2;对于固定时差测量模型,
综合时间测量误差方差矩阵的迹为:
<mrow>
<mi>t</mi>
<mi>r</mi>
<mi>a</mi>
<mi>c</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>&epsiv;</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>T</mi>
<mn>0</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>T</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>T</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>4</mn>
<msubsup>
<mi>&sigma;</mi>
<mi>s</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
三站时差定位***中,多个基站的综合时间测量误差方差表示为综合时间测量误差方差矩阵迹的形式为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&sigma;</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>&rsqb;</mo>
<mo>=</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>P</mi>
<mi>t</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>T</mi>
<mi>i</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>&rsqb;</mo>
<mo>=</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<mfrac>
<msubsup>
<mi>&sigma;</mi>
<mi>s</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>T</mi>
<mi>i</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>T</mi>
<mn>0</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>T</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>T</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<msubsup>
<mi>&sigma;</mi>
<mi>s</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>4</mn>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mi>t</mi>
<mi>r</mi>
<mi>a</mi>
<mi>c</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>&epsiv;</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
其中σPt=σs/c为基站的站址测量误差,即已经转换为时间;
XY平面三站时差定位误差GDOP为:
<mrow>
<msub>
<mi>GDOP</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msubsup>
<mi>&sigma;</mi>
<mi>x</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>y</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>=</mo>
<msqrt>
<mrow>
<mi>t</mi>
<mi>r</mi>
<mi>a</mi>
<mi>c</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>d</mi>
<mi>X</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
用单个基站综合时间测量误差对GDOP归一化,得到归一化GDOP为:
<mrow>
<msub>
<mi>GDOP</mi>
<mi>e</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>GDOP</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
</msub>
</mrow>
<msqrt>
<mrow>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msubsup>
<mi>&sigma;</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mfrac>
<mo>=</mo>
<mfrac>
<msqrt>
<mrow>
<mi>t</mi>
<mi>r</mi>
<mi>a</mi>
<mi>c</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>d</mi>
<mi>X</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
<msqrt>
<mrow>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<mi>t</mi>
<mi>r</mi>
<mi>a</mi>
<mi>c</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>&epsiv;</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
其中归一化GDOP为时的定位误差;
假设可以得到从而得到归一化GDOP为:
<mrow>
<msub>
<mi>GDOP</mi>
<mi>e</mi>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<mi>t</mi>
<mi>r</mi>
<mi>a</mi>
<mi>c</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>d</mi>
<mi>X</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
<mo>=</mo>
<msqrt>
<mrow>
<mi>t</mi>
<mi>r</mi>
<mi>a</mi>
<mi>c</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>A</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msub>
<mi>P</mi>
<mi>&epsiv;</mi>
</msub>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msup>
<mi>A</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
证明归一化GDOP与基站、辐射源位置有关,与综合时间测量误差方差矩阵Pε无关;
步骤三 机载辐射源的直线航线设计
利用步骤一和步骤二获取巡航高度H平面的归一化定位误差几何稀释,设置直线飞行航线,选取的航线变化趋势平缓,用于避免通过定位误差梯度变换比较快的位置,以便分段统计定位误差;同时当目标飞机沿设计的航线飞行时,要满足在有效航线上辐射源信号始终能被定位基站接收到;包括:航线最远点、航线长度、航线高度以及有效飞行航次;
所述航线最远点能够满足飞行空域限制、飞机有效航程、最远侦察距离因素的限制;所述航线高度按照飞机巡航高度,试验航线最近点位置相对最近基站的俯仰角设置小于三站时差定位***俯仰覆盖的范围;
飞行架次的确定:飞行架次根据距离单元的大小进行推算,首先计算所需的航线往返次数
<mrow>
<msub>
<mi>F</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>V</mi>
<mo>&times;</mo>
<mi>T</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>R</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:Fn为航线往返次数,ΔR为航线距离取样区间大小,单位为m,N为航线距离取样区间ΔR内需统计的定位数据总量,通常为60,V为目标机飞行速度,单位为m/s,T为定位***数据录取周期,单位为s;
步骤四试验数据的获取与处理
获取试验数据的数据量,根据统计数据量设置目标飞机一架沿航线水平、直线、往返飞行的航线往返次数Fn;
当机载辐射源开机,被试装备主站、副站天线对准预定航线进行扇扫,分别对机载辐射源的目标信号进行侦收;副站将数据传送到主站,主站进行数据处理,用相关处理法得到主站和所有副站的侦察信号时差数据,然后主站按照***定位算法联合处理时差信息,获取空间辐射源的位置;
标准位置测量设备为机载GPS***或精密测量雷达,全航路跟踪目标获取辐射源真实位置;
当机载辐射源不在定位平面时,三站时差定位存在***模型误差,由于时间测量误差和***模型误差是相互独立的,因此实际三站时差定位***的理论定位误差模型为:
<mrow>
<msub>
<mi>&sigma;</mi>
<mi>X</mi>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msubsup>
<mi>&sigma;</mi>
<mi>x</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>y</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>z</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>=</mo>
<msqrt>
<mrow>
<msub>
<mi>GDOP</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>z</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msubsup>
<mi>&sigma;</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>GDOP</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>z</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
其中为辐射源高度引入的定位模型误差方差;
以三个基站为平面建立XY定位平面,其中主站坐标为(x0,y0,0),两副站坐标分别为(x1,y1,0)和(x2,y2,0),辐射源坐标为(x,y,z),解算出实际辐射源到各定位基站的距离或时延,为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mi>c</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mi>c</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
然后根据三站定位算法解算出目标位置测量值(x′,y′,0);
由目标的真实位置和测量位置就可以得到辐射源高度引起的定位模型误差方差为:
<mrow>
<msubsup>
<mi>&sigma;</mi>
<mi>z</mi>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msup>
<mi>y</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
</mrow>
由(13)、(16)式,得到综合时间测量误差方差的试验统计值为:
<mrow>
<msubsup>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<msub>
<mi>&sigma;</mi>
<mi>X</mi>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>z</mi>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msubsup>
<mi>GDOP</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>19</mn>
<mo>)</mo>
</mrow>
</mrow>
其中c为光速,σX 2由三维定位误差试验结果统计得到,不同航线位置σX 2不同,为基站和辐射源坐标位置确定的定位模型误差方差,GDOPe为基站和辐射源坐标位置确定的归一化误差;
用相关处理方法,计算多个基站的综合时间测量误差方差
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&sigma;</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>4</mn>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mi>t</mi>
<mi>r</mi>
<mi>a</mi>
<mi>c</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>&epsiv;</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>T</mi>
<mn>0</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>T</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>T</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msup>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>P</mi>
<mi>t</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>t</mi>
<mn>0</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>t</mi>
<mn>1</mn>
</mrow>
<mn>4</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>t</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msup>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>P</mi>
<mi>t</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msubsup>
<mi>&sigma;</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msup>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>P</mi>
<mi>t</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msup>
<mi>B</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>4</mn>
<msub>
<mi>K</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mn>6</mn>
<msubsup>
<mi>R</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>4</mn>
<msup>
<msub>
<mi>K</mi>
<mn>0</mn>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<msub>
<mi>T</mi>
<mi>C</mi>
</msub>
<msub>
<mi>f</mi>
<mi>s</mi>
</msub>
<msup>
<msub>
<mi>R</mi>
<mn>0</mn>
</msub>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>R</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>20</mn>
<mo>)</mo>
</mrow>
</mrow>
其中σPt=σs/c为站址测量误差,B为辐射信号带宽,τ=1/B,TC为相参处理时间长度,fs为采样率,Ri为辐射源距离第i个基站的距离,和K0为固定常数,同等信号参数情况下,K0越小,距离产生的影响越大;
当辐射源***参数和各基站***参数确定后,K0为常数;另外由于三站时差定位各主副站间距一般比较小,而侦察天线波束比较宽,在远距离各站的参数Ki基本相同,即:
<mrow>
<msub>
<mi>K</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<msub>
<mi>K</mi>
<mi>i</mi>
</msub>
<mo>&rsqb;</mo>
<mo>=</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mi>t</mi>
</msub>
<msub>
<mi>G</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&phi;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>G</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mi>&pi;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msub>
<mi>L</mi>
<mrow>
<mi>j</mi>
<mi>r</mi>
</mrow>
</msub>
<msub>
<mi>L</mi>
<mi>t</mi>
</msub>
<msubsup>
<mi>&sigma;</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<msub>
<mi>N</mi>
<mi>p</mi>
</msub>
<msub>
<mi>T</mi>
<mi>r</mi>
</msub>
<msub>
<mi>f</mi>
<mi>s</mi>
</msub>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>21</mn>
<mo>)</mo>
</mrow>
</mrow>
其中Pt为辐射源峰值功率,Gt(φi)为辐射源在第i个侦察接收机方向的增益,Gi(θi)为第i个侦察接收机的接收增益,λ为雷达工作波长,Ljr为侦察接收机功率接收损耗,Lt为发射机***损耗,为第i个接收机噪声功率,NP为侦察脉冲数目,Tr为辐射信号脉冲宽度,fs为采样率;
三站时差定位性能的评估,就是通过试验结果求取参数和K0的值,通过(19)、(20)式,建立方程,即综合时间测量误差模型为:
<mrow>
<mi>f</mi>
<mo>=</mo>
<msubsup>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>&sigma;</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msup>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>P</mi>
<mi>t</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msup>
<mi>B</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>K</mi>
<mn>0</mn>
</msub>
</mfrac>
<mi>g</mi>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msup>
<msub>
<mi>K</mi>
<mn>0</mn>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<msub>
<mi>T</mi>
<mi>C</mi>
</msub>
<msub>
<mi>f</mi>
<mi>s</mi>
</msub>
<mi>q</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>22</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>f</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
</mfrac>
<mfrac>
<mrow>
<msup>
<msub>
<mi>&sigma;</mi>
<mi>X</mi>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>z</mi>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>GDOP</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>g</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mn>6</mn>
<msubsup>
<mi>R</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>q</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<msup>
<msub>
<mi>R</mi>
<mn>0</mn>
</msub>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>R</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>23</mn>
<mo>)</mo>
</mrow>
</mrow>
其中c为光速,σX 2为试验获取的定位误差方差,GDOPe为由辐射源位置和定位站位置得到的归一化误差,和K0为固定常数,和由辐射源位置和定位站位置得到,TC为相参处理时间长度,fs为采样率,B为信号带宽,为的估计数值,通过试验能够获取综合时间测量误差方差的估计值f;
综合时间测量误差的方差估计值f分布是非平稳的,采用等距离段分段统计,判断距离段内是平稳随机信号;
设置第m个距离段得到L个定位实测结果Xi,对应真实位置为X0i,以及L个定位点对应的辐射源和定位站位置;
通过(23)式、辐射源和定位站位置直接计算gm和qm参数均值;对分段数据综合时间测量误差的方差估计值f用3σ准则剔除异常误差,统计得到各距离段内的估计值fm;第m个距离段的综合时间测量误差方差的估计值fm为:
<mrow>
<msub>
<mi>f</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>L</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>L</mi>
</munderover>
<mfrac>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>X</mi>
<mrow>
<mn>0</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>-</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>z</mi>
<mi>i</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msup>
<msub>
<mi>GDOP</mi>
<mrow>
<mi>e</mi>
<mi>i</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>24</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,L为距离段内定位点数,|| ||为向量范数或取向量长度,为第i个辐射源位置对应的定位模型误差方差,GDOPei为第i个辐射源位置归一化定位误差;
由(22)式得到第m个距离段的综合时间测量误差方差的估计值为:
<mrow>
<msub>
<mi>f</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<msubsup>
<mi>&sigma;</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msup>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>P</mi>
<mi>t</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msup>
<mi>B</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>K</mi>
<mn>0</mn>
</msub>
</mfrac>
<msub>
<mi>g</mi>
<mi>m</mi>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msup>
<msub>
<mi>K</mi>
<mn>0</mn>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<msub>
<mi>T</mi>
<mi>C</mi>
</msub>
<msub>
<mi>f</mi>
<mi>s</mi>
</msub>
<msub>
<mi>q</mi>
<mi>m</mi>
</msub>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>25</mn>
<mo>)</mo>
</mrow>
</mrow>
其中m的取值范围为1~M;得到M个位置的综合时间测量误差方差的估计值fm,及参数gm和qm;
步骤五参数化求解过程
将方程(25)简化为:
fm=a+bxm+b2ym (26)
其中
按照最小二乘法,定义:
<mrow>
<mi>Q</mi>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>+</mo>
<msub>
<mi>bx</mi>
<mi>m</mi>
</msub>
<mo>+</mo>
<msup>
<mi>b</mi>
<mn>2</mn>
</msup>
<msub>
<mi>y</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>27</mn>
<mo>)</mo>
</mrow>
</mrow>
判断a,b最优估计的问题,就是求解Q(a,b)最小值的问题,Q(a,b)分别对a,b求偏导,并令它们等于零,解方程组就能够得到a,b的估计数值,即
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>Q</mi>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>a</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mn>2</mn>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>+</mo>
<msub>
<mi>bx</mi>
<mi>m</mi>
</msub>
<mo>+</mo>
<msup>
<mi>b</mi>
<mn>2</mn>
</msup>
<msub>
<mi>y</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>Q</mi>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>b</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mn>2</mn>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>+</mo>
<msub>
<mi>bx</mi>
<mi>m</mi>
</msub>
<mo>+</mo>
<msup>
<mi>b</mi>
<mn>2</mn>
</msup>
<msub>
<mi>y</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>m</mi>
</msub>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>by</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>28</mn>
<mo>)</mo>
</mrow>
</mrow>
解(28)式可得到b的方程为:
U(b)=Sxf+Syfxxb+Sxyb2+Syyb3=0 (29)
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>x</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>x</mi>
<mi>m</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>f</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>f</mi>
<mi>m</mi>
</msub>
<msub>
<mi>x</mi>
<mi>m</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>y</mi>
<mi>f</mi>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>f</mi>
<mi>m</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>y</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>f</mi>
<mi>m</mi>
</msub>
<msub>
<mi>y</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>x</mi>
<mi>m</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>x</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>x</mi>
<mi>m</mi>
</msub>
<msub>
<mi>x</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<mn>3</mn>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>x</mi>
<mi>m</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>y</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>x</mi>
<mi>m</mi>
</msub>
<msub>
<mi>y</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>y</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>y</mi>
<mi>m</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>y</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>y</mi>
<mi>m</mi>
</msub>
<msub>
<mi>y</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>30</mn>
<mo>)</mo>
</mrow>
</mrow>
根据三次方程求解方法得到该方程的解析解,或采用MATLAB中roots函数多项式求根方法;再根据辐射源和时差定位***参数,估计出K0的数量级或近似数值范围;
根据b的取值范围,可得到b的一个有效解,解得b,代入(31)式得到a为:
<mrow>
<mi>a</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>f</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<mi>b</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>x</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<msup>
<mi>b</mi>
<mn>2</mn>
</msup>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>y</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>31</mn>
<mo>)</mo>
</mrow>
</mrow>
按照最小二乘法得到待估计参数和K0;
或通过曲线插值方法获取两位置的综合时间测量误差方差参数fC,fD,由基站和辐射源位置,按照(23)式得到C和D两个位置参数gC,qC的gD,qD,通过求二元方程组得到估计综合时间测量误差参数和K0;
步骤六任意布站模式下定位性能
三站时差定位存在***模型误差,由于时间测量噪声误差和***模型误差是相互独立的,则三站时差定位***定位的任意航线误差或任意X位置的定位误差,即综合时间测量误差模型定位误差,为:
<mrow>
<msub>
<mi>&sigma;</mi>
<mi>X</mi>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msubsup>
<mi>&sigma;</mi>
<mi>x</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>y</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>z</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mi>X</mi>
</msub>
<msup>
<msub>
<mrow>
<mo>(</mo>
<msub>
<mi>GDOP</mi>
<mi>e</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>X</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>z</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>32</mn>
<mo>)</mo>
</mrow>
</mrow>
其中(GDOPe)X由步骤二获取,由步骤五获取,由(18)式获取。
2.根据权利要求1所述的一种三站时差定位性能试验评估方法,其特征在于:所述三站时差定位性能试验评估模型包括:建立综合时间测量误差模型和定位误差模型,所述综合时间测量误差模型由信噪比引起的时间测量误差和综合时间测量误差组成,综合时间测量误差包括基站站址测量误差和***固有时间测量误差。
3.根据权利要求1所述的一种三站时差定位性能试验评估方法,其特征在于:所述辐射源设置为在两固定位置上的固定辐射源;固定辐射源在地面上支撑两个杆的固定位置T0、T1上,通过两个杆的固定位置上辐射源在两固定位置T0、T1的多次测试结果,计算辐射源任意位置的定位精度;通过试验获取两固定位置的综合时间测量误差方差参数为fC,fD,由基站和辐射源位置,按照(23)式得到C和D两个位置参数gC,qC的gD,qD,通过求二元方程组得到估计综合时间测量误差参数和K0,进而按照步骤六得到任意布站模式和任意辐射源位置的定位误差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410380137.4A CN104316903B (zh) | 2014-08-01 | 2014-08-01 | 一种三站时差定位性能试验评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410380137.4A CN104316903B (zh) | 2014-08-01 | 2014-08-01 | 一种三站时差定位性能试验评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104316903A CN104316903A (zh) | 2015-01-28 |
CN104316903B true CN104316903B (zh) | 2018-03-09 |
Family
ID=52372156
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410380137.4A Expired - Fee Related CN104316903B (zh) | 2014-08-01 | 2014-08-01 | 一种三站时差定位性能试验评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104316903B (zh) |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105044669A (zh) * | 2015-05-22 | 2015-11-11 | 中国电子科技集团公司第十研究所 | 三站测时差立体定位方法 |
CN105093174B (zh) * | 2015-08-31 | 2018-02-06 | 成都金本华电子有限公司 | 一种基于2.5g无线网络信号损益的定位算法 |
CN105403854B (zh) * | 2015-11-04 | 2017-08-25 | 哈尔滨工业大学 | 基于天线周期扫描时间差无源定位方法 |
CN105242262B (zh) * | 2015-11-04 | 2017-11-03 | 哈尔滨工业大学 | 一种基于天线周期扫描时间差无源定位方法 |
CN106842184B (zh) * | 2015-12-03 | 2019-12-17 | 中国航空工业集团公司雷华电子技术研究所 | 一种基于波束调度的多目标探测与跟踪方法 |
CN106249197B (zh) * | 2016-05-03 | 2019-11-29 | 电子科技大学 | 一种多点定位***中接收机位置误差的自校准方法 |
CN106855631B (zh) * | 2016-12-02 | 2019-06-04 | 中国人民解放军63891部队 | 一种转发式gnss动态测量精度测试评估方法 |
CN106872942A (zh) * | 2016-12-23 | 2017-06-20 | 安徽四创电子股份有限公司 | 用于分布式多点定位监视***的定位精度解算方法 |
CN107357758B (zh) * | 2017-06-29 | 2021-04-13 | 中国人民解放军63796部队 | 定位信息的多项式最小二乘回归记忆性求速方法 |
CN107315891B (zh) * | 2017-07-17 | 2019-09-13 | 电子科技大学 | 基于辐射源的机载通信干扰航线建模与干扰效能评估方法 |
CN107907860A (zh) * | 2017-11-15 | 2018-04-13 | 电子科技大学 | 一种利用频控阵技术抗时差测量的定位欺骗方法 |
EP3588135B1 (en) | 2018-06-28 | 2021-05-05 | Aptiv Technologies Limited | Method of determining an alignment error of an antenna and vehicle with an antenna and a detection device |
CN108919211B (zh) * | 2018-07-12 | 2022-03-11 | 中国船舶重工集团公司第七二四研究所 | 一种面向三站协同定位的定位精度评估方法 |
CN110611878B (zh) * | 2019-08-23 | 2021-10-19 | 华清科盛(北京)信息技术有限公司 | 一种基于基站一维部署的高精度定位***及方法 |
CN110441779B (zh) * | 2019-09-18 | 2023-05-23 | 哈尔滨工程大学 | 一种多声纳浮标分布式协同定位方法 |
CN111273223A (zh) * | 2020-01-20 | 2020-06-12 | 西安交通大学 | 一种基于双天线时延无源定位对抗方法 |
CN111371516B (zh) * | 2020-02-28 | 2021-11-12 | 中国民用航空总局第二研究所 | 一种分布式定位***远端接收站性能监测方法与装置 |
CN111273274B (zh) * | 2020-03-12 | 2022-03-18 | 四川九洲电器集团有限责任公司 | 多基协同定位方法、存储介质、雷达及雷达定位*** |
CN111815031A (zh) * | 2020-06-16 | 2020-10-23 | 智慧航海(青岛)智能***工程有限公司 | 一种基于岸基雷达***探测能力的雷达优化布站方法 |
CN111880145B (zh) * | 2020-08-03 | 2023-07-07 | 中国电子科技集团公司第三十六研究所 | 辐射源时差定位方法、装置及电子设备 |
CN112147574B (zh) * | 2020-09-25 | 2022-02-01 | 四川九洲电器集团有限责任公司 | 一种多站定位算法验证方法和一种模拟*** |
-
2014
- 2014-08-01 CN CN201410380137.4A patent/CN104316903B/zh not_active Expired - Fee Related
Non-Patent Citations (5)
Title |
---|
Differential Delay/Doppler ML Estimation with Unknown Signals;Seymour Stein;《IEEE TRANSACTIONS ON SIGNAL PROCESSING》;19930831;第41卷(第8期);2717-2719 * |
三站时差定位的精度分析与推算模型;陈永光 等;《电子学报》;20040930;第32卷(第9期);1452-1455 * |
四站时差定位精度分析;俞志强;《空军雷达学院学报》;20101231;第24卷(第6期);400-402 * |
时差无源定位精度分析及其等效推算试验方法研究;张政超 等;《中国电子科学研究院学报》;20130430;第8卷(第2期);196-200 * |
有源无源一体模式下目标定位及其精度分析;张政超 等;《电讯技术》;20121130;第52卷(第11期);1758-1762 * |
Also Published As
Publication number | Publication date |
---|---|
CN104316903A (zh) | 2015-01-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104316903B (zh) | 一种三站时差定位性能试验评估方法 | |
US11246010B2 (en) | Method and system for positioning and tracking using a plurality of sensors | |
KR101240629B1 (ko) | Ads-b 시스템이 탑재된 항공기를 이용한 미지신호 검출 및 발생원 위치 추정방법 | |
CN112073895A (zh) | 一种基于csi高精度人员定位跟踪方法 | |
Strohmeier et al. | A k-NN-based localization approach for crowdsourced air traffic communication networks | |
CN105759241B (zh) | 基于时差频差的直接定位方法 | |
KR101234177B1 (ko) | 사용자 단말의 위치 측정 방법 | |
CN110954865A (zh) | 一种基于电离层信息的短波时差定位方法 | |
Zwirello et al. | Localization in industrial halls via ultra-wideband signals | |
Sinha et al. | Impact of 3D antenna radiation patterns on TDOA-based wireless localization of UAVs | |
Keunecke et al. | Deriving 2D TOA/TDOA IEEE 802.11 g/n/ac location accuracy from an experimentally verified fading channel model | |
CN106569180A (zh) | 一种基于Prony方法的方位估计算法 | |
Galati et al. | Multilateration (local and wide area) as a distributed sensor system: lower bounds of accuracy | |
Stefanski | Asynchronous wide area multilateration system | |
Strohmeier et al. | A localization approach for crowdsourced air traffic communication networks | |
Yaro et al. | Position estimation error performance model for a minimum configuration 3-D multilateration | |
CN102638810A (zh) | 基于多维信道分量功率谱密度的信道建模仿真平台 | |
RU2657237C1 (ru) | Однопозиционный способ определения координат местоположения источников радиоизлучения | |
Kelner et al. | Influence of the frequency stability on the emitter position in SDF method | |
CN102998674A (zh) | 多通道sar地面慢动目标检测方法及装置 | |
CN104698444A (zh) | 确定阵元位置误差对阵列天线三维成像质量影响的方法 | |
Kim et al. | A study on the implementation of a 3-dimensional positioning system on indoor environments | |
Marinov et al. | An algorithm optimizing the geometric distribution of the stations in aeronautical multilateration systems. | |
RU2625094C1 (ru) | Способ определения пеленга и дальности до источника сигналов | |
Yaro et al. | Position Estimation Bias Analysis of a Multilateration System with a Reference Station Selection Technique |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180309 Termination date: 20180801 |
|
CF01 | Termination of patent right due to non-payment of annual fee |