CN105487048B - 基于置信椭圆的二站纯方位定位模糊区方法 - Google Patents
基于置信椭圆的二站纯方位定位模糊区方法 Download PDFInfo
- Publication number
- CN105487048B CN105487048B CN201510731399.5A CN201510731399A CN105487048B CN 105487048 B CN105487048 B CN 105487048B CN 201510731399 A CN201510731399 A CN 201510731399A CN 105487048 B CN105487048 B CN 105487048B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- location
- direction finding
- confidence
- 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.)
- Active
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)
- Analysing Materials By The Use Of Radiation (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明公开了一种基于置信椭圆的二站纯方位定位模糊区方法，要解决的技术问题是提高定位模糊区的定位精度和定位置信度。技术方案是先建立基于置信椭圆定位模糊区模型，然后实验求解，得到置信椭圆定位模糊区模型中测向误差系数k的精确取值范围为1.75～2.5，确定置信椭圆定位模糊区，同时确定目标辐射源O3位于置信椭圆定位模糊区中的置信度。采用本发明一方面可以达到工程实践对定位置信度的要求，具有较好的定位置信水平，另一方面，相对于现有的定位模糊区方法，本发明具有较高的定位精度，同时可以体现定位误差的分布情况。
Description
技术领域
本发明涉及无源定位领域，特指在二站纯方位背景下，采用基于置信椭圆的定位模糊区方法对目标点进行定位。
背景技术
在现代化战争中，争夺电磁领域的主动权是争夺战场主动权的核心，对辐射源进行定位有助于提高精确打击武器的使用性能，为最终摧毁敌方提供有力保障。因此，如何提高对威胁目标的定位精度以及如何提高定位效率具有重要的意义。
无源探测系统自身不发射电磁波，主要结合时差定位和天线方向来实现测角进而得到目标地理位置，其优点是有效范围广、隐蔽性好而且抗干扰能力强且节能对环境友好。无源定位技术分为纯方位定位技术和联合定位技术。纯方位定位技术主要利用测向角度参数进行目标点定位，联合定位技术在纯方位思想的基础上，结合多种信息进行联合探测定位，其定位原理基本一样。纯方位定位是无源定位技术的基础，具有最广泛深入的研究。二站纯方位定位作为一种最基本的多站纯方位定位方式，在目标领域是研究的热点和重点。
作为最重要的一种纯方位定位方法，二站测向交叉定位受到国内外许多学者的重视并在工程和理论研究上取得了许多成果。但是许多学者的研究重点是解决如何排除虚假定位点以及如何提高优化布站效率问题，针对定位模糊区方法的研究较少而且主要集中在定位模糊区面积研究上。针对二站纯方位定位模糊区的研究主要有基于交叉定位方法的四边形定位模糊区方法和基于GDOP方法的圆形定位模糊区方法。基于GDOP的圆形定位模糊区方法由于其定位面积较大，定位精度不高，所以在纯方位定位领域不常采用。基于交叉定位的四边形定位模糊区方法参见修建娟，王国宏，何友，修建华在《系统工程与电子技术学报》2005年第8 卷发表的《纯方位系统中的定位模糊区分析》。基于交叉定位的四边形定位模糊区方法中的四边形定位模糊区几何示意图如图1所示：由测向站O1和测向站O2构成二站纯方位定位系统，以测向站O1和测向站O2连线的中点O为坐标原点，以两测向站连线所在直线为x轴，右向为 x轴正方向，以两测向站连线垂直方向为y轴建立平面直角坐标系。R为两测向站之间的距离，则两个测向站O1(－R/2,0)和O2(R/2,0)的位置已知,O3(x,y)为一个需要定位的目标辐射源。考虑由于测向不准所引入的定位误差，两个测向站的方位测量均服从正态分布，Δθ1为测向站O1方位测量误差的标准差，Δθ2为测向站O2方位测量误差的标准差，Δθ1和Δθ2通过阅读测向站中传感器的使用规格说明可以得到，若两测向站使用的传感器相同，则Δθ1＝Δθ2。θ1为测向站O1对目标辐射源O3实施测向的均方位角；θ2为测向站O2对目标O3辐射源实施测向的均方位角(其中θ2'为θ2的补角)。
取测向站O1测向误差为±kΔθ1，测向站O2测向误差为±kΔθ2，k为测向误差系数。则以 tan(θ1+kΔθ1)为斜率，过O1形成的直线l1为O1的一条测向边界线；以tan(θ1-kΔθ1)为斜率，过 O1点形成的直线l2为O1的另一条测向边界线。以tan(θ2+kΔθ2)为斜率，过O2点形成的直线l3为 O2的一条测向边界线；以tan(θ2-kΔθ2)为斜率，过O2点形成的直线l4为O2的另一条测向边界线。l1与l3、l4形成的交点分别为点B和点A；l2与l3、l4形成的交点分别为点C和点D。由 ABCD形成的四边形区域ABCD称为定位模糊区ABCD(图1 所示阴影区)。为了方便描述，用LAB表示l1，LCD表示l2，LBC表示l3，LAD表示l4。当目标辐射源位于x轴上时不能确定目标的位置。但在实际工程应用中，目标辐射源一般都处于不断运动的状态，其运动方向会随时间不断改变，目标辐射源位于x轴上的情况只会存在于某个时刻，因此从时间角度来说，这种情况并不会对定位产生影响。另一方面，实际工程应用中多采用多站定位，当目标辐射源位于某两个测向站连线上时，可以调用其它测向站进行定位，因此从空间角度来说，这种问题也可以得到解决。综上，目标辐射源位于测向站连线上这种情况并不会对基于交叉定位的四边形定位模糊区方法在应用上造成影响。
但是，基于交叉定位的四边形定位模糊区方法只能反映模糊区的面积大小，不能反映模糊区误差分布等其它特征；由于存在定位误差，在四边形定位模糊区以外的目标辐射源容易被遗漏丢失导致该方法定位置信水平不高；由于四边形定位模糊区表示参数较为复杂，导致其不易于在工程实践中应用。
如何提高定位模糊区的置信水平和定位精度是本领域技术人员正在努力解决的重要技术问题。
发明内容
本发明要解决的技术问题在于：针对目前现有基于交叉定位的四边形定位模糊区方法定位精度不高，定位置信度较低，无法体现定位误差分布特性等问题，在二站定位背景下，提出一种新的定位方法——基于置信椭圆的二站纯方位定位模糊区方法，有效提高定位模糊区的定位精度和定位置信度。
本发明包括以下两个步骤：
第一步、建立基于置信椭圆定位模糊区模型，方法是：
1.1 对于由测向站O1和测向站O2构成的二站纯方位定位系统，如图4所示，测向站O1和测向站O2的测角标准差Δθ1和Δθ2为已知量，通过测量获取测向站O1和测向站O2对目标定位点的测量参数，即测角θ1和θ2以及测向站O1和测向站O2之间的距离R；
1.2 建立四边形定位模糊区ABCD的测向边界线LAB、LCD、LBC、LAD方程：
其中θ2'为θ2的补角，k为测向误差系数。
1.3 建立四边形定位模糊区ABCD的外接椭圆簇方程，得到基于置信椭圆的定位模糊区模型：
四边形定位模糊区ABCD的外接椭圆簇方程为:
LABLCD+λLBCLAD＝0 (2)
将式(1)中LAB、LCD、LBC、LAD的直线方程带入式(2)，得到：
将式(3)转化为椭圆簇方程形式，得到基于置信椭圆的定位模糊区模型为：
其中:a＝tan(θ1+kΔθ1),b＝tan(θ1-kΔθ1),c＝tan(θ2'-kΔθ2),d＝tan(θ2'+kΔθ2)
第二步、实验求解，得到置信椭圆定位模糊区模型中测向误差系数k的精确取值范围，确定置信椭圆定位模糊区，同时确定目标辐射源O3位于置信椭圆定位模糊区中的置信度：
2.1 获得测向误差系数k的粗略取值范围：
2.1.1 取测向站O1和测向站O2之间的距离R＝10km，取两测向站的测角标准差Δθ1＝Δθ2＝0.5°，取测角θ1＝30°，θ'2＝30°，测向误差系数k依次取值为1、1.5、2、2.5、3，代入基于置信椭圆的定位模糊区模型(即公式(4))，共进行5次实验，采用matlab数学软件进行模型求解，每次实验得到相应的置信椭圆定位模糊区，用蒙特卡罗方法(参见裴鹿成在第五届全国蒙特卡罗学术交流会上发表的《蒙特卡罗方法及其应用》第3页第2段)对置信椭圆定位模糊区进行仿真实验，根据5次仿真实验结果求解出该定位模糊区的置信度和面积并记录，得到5个置信度和5个面积。
2.1.2 设置R＝20km，按步骤2.1.1的方法再进行5次实验，得到5个置信度和5个面积。
2.1.3 设置R＝30km，按步骤2.1.1的方法再进行5次实验，得到5个置信度和5个面积。
2.1.4 观察上述15个置信度，发现这些置信度均处于0.5～1这个范围，且置信度随着k的增加而增加，当k小于1时，置信度小于0.5，当k接近3时，置信度逐渐趋近于1；观察15个面积，发现面积随着k的增加而增加，且增加速度逐渐变大。因此得出结论：当k 取值在1～3之间置信度可满足0.5～1。故为了同时保证定位置信度和定位精度，k的粗略取值范围为1～3。
2.2 获得测向误差系数k的精确取值范围：
2.2.1 初始化变量：设置测向站O1和测向站O2之间的距离变量R＝10km；设置测向站O1方位测量误差的标准差变量Δθ1＝0.1°，测向站O2方位测量误差的标准差变量Δθ2＝0.1°；设置第一测角变量θ1＝30°，第二测角补角变量θ'2＝30°；设置测向误差系数变量k＝1。
2.2.2 将R、Δθ1、Δθ2、θ1、θ'2、k带入基于置信椭圆的定位模糊区模型(即公式(4))，采用matlab数学软件进行模型求解，得到相应的置信椭圆定位模糊区，用蒙特卡罗方法对该置信椭圆定位模糊区进行仿真实验，根据仿真结果求解出该定位模糊区的置信度并记录，得到置信度和面积。
2.2.3 令k＝k+0.1。
2.2.4 判断k是否小于等于3，若是，则转入步骤2.2.2；若否，则令k＝1，θ'2＝θ'2+5°，转步骤2.2.5。
2.2.5 判断θ'2是否小于等于140°，若是，则转步骤2.2.6；若否，则令θ'2＝30°，θ1＝θ1+5°，转步骤2.2.7。
2.2.6 判断θ'2+θ1＝180°是否成立，若成立，则θ'2＝θ'2+5°，转步骤2.2.5；若否，则转步骤2.2.2。
2.2.7 判断θ1是否小于等于140°，若是，则转步骤2.2.2；若否，则θ1＝30°，Δθ1＝Δθ1+0.05°，Δθ2＝Δθ2+0.05°，转步骤2.2.8。
2.2.8 判断Δθ1是否小于等于1°且Δθ2是否小于等于1°，若是，则转步骤2.2.2；若否，则Δθ1＝Δθ2＝0.1°，R＝R+5，转步骤2.2.9。
2.2.9 判断R是否小于等于100，若是，则转入步骤2.2.2；若否，则实验结束，执行2.2.10。
2.2.10 对比实验记录的结果，发现当测向误差系数k在1.75～2.5之间取值可以满足工程实践中定位模糊区置信度在0.950～0.999的要求，同时具有较高定位精度。得到了k，也即确定了基于置信椭圆的定位模糊区模型的各参数，图4 的四边形定位模糊区ABCD的外接椭圆也就确定了，也即得到了置信椭圆定位模糊区。
采用本发明可以达到以下效果；
1、本发明一方面可以达到工程实践对定位置信度的要求，具有较好的定位置信水平，另一方面，通过仿真实验的结果分析和对比，如图6所示，相对于现有的定位模糊区方法，本发明具有较高的定位精度，同时可以体现定位误差的分布情况。在现有的无源定位中，可以有效提高纯方位定位系统的目标定位效率和准确性，同时为无源定位提供更多的定位方向和定位分布上的信息，有助于对辐射源进行定位，从而提高精确打击武器的使用性能，为最终摧毁敌方提供有力保障。
2、本发明第一步只利用测向站的测角标准差Δθ1和Δθ2和测向角度θ1和θ2，以及测向站之间的距离R进行建模，该模型不依赖于定位分布情况的变化而变化，具有较好的稳健性。因此可使用的范围较为广泛，而且定位所需参数较少，定位速度快，操作简单，便于工程实践中的实现。
3、本发明第二步通过多次实验给出给出了测向误差系数k的取值范围在1.75—2.5之间。当k在该区间取值时，置信椭圆定位模糊区可以满足工程实践中定位模糊区置信度达到0.95 以上的要求，同时具有较高定位精度，加强了本发明在工程实践中的可应用性。
附图说明
图1是背景技术基于交叉定位的四边形定位模糊区方法示意图。
图2是本发明基于置信椭圆定位模糊区的定位方法流程图。
图3是本发明第二步中步骤2.2的实验方法流程图。
图4是本发明第一步建立的置信椭圆定位模糊区模型示意图。
图5是本发明第二步中当测向误差系数k＝1.75，k＝2，k＝2.5时置信椭圆定位模糊区示意图。
图6是三种定位模糊区方法仿真结果对比示意图。
图7是本发明第二步中步骤2.2的部分实验数据。
具体实施方式
图1是背景技术基于交叉定位的四边形定位模糊区方法示意图。以测向站O1和测向站O2连线的中点为坐标原点，以两测向站连线所在直线为x轴，右向为x轴正方向，以两测向站连线垂直方向为y轴建立平面直角坐标系。R为两测向站之间的距离，则两个测向站 O1(－R/2,0)和O2(R/2,0)的位置已知,O3(x,y)为一个需要定位的目标辐射源。考虑由于测向不准所引入的定位误差，两个测向站的方位测量均服从正态分布，Δθ1为测向站O1方位测量误差的标准差；Δθ2为O2测向站方位测量误差的标准差，Δθ1和Δθ2通过阅读测向站中传感器的使用规格说明可以得到，若两测向站使用的传感器相同，则Δθ1＝Δθ2。θ1为测向站O1对目标O3辐射源实施测向的均方位角；θ2为测向站O2对目标O3辐射源实施测向的均方位角(其中θ2'为θ2的补角)。
取测向站O1测向误差为±kΔθ1，k为测向误差系数；测向站O2测向误差为±kΔθ2。则以 tan(θ1+kΔθ1)为斜率，过O1点形成的直线l1为测向站O1的一条测向边界线；以tan(θ1-kΔθ1)为斜率，过O1点形成的直线l2为测向站O1的另一条测向边界线。以tan(θ2+kΔθ2)为斜率，过O2点形成的直线l3为测向站O2的一条测向边界线；以tan(θ2-kΔθ2)为斜率，过O2点形成的直线l4为测向站O2的另一条测向边界线。l1与l3、l4形成的交点分别为点B和点A；l2与l3、l4形成的交点分别为点D和点C。由ABCD形成的四边形区域ABCD称为定位模糊区ABCD(图1所示阴影区)。为了方便描述，用LAB表示l1，LCD表示l2，LBC表示l3，LAD表示l4。
图2是本发明基于置信椭圆定位模糊区的定位方法流程图。本发明包括两个步骤：第一步确定交叉定位的四边形定位模糊区的最小外接椭圆，建立置信椭圆定位模糊区模型。首先获取两个测量基站O1和O2对目标定位点的测量参数，即测角θ1和θ2以及两测量基站的测角标准差Δθ1和Δθ2，两基站之间的距离R。建立四边形定位模糊区ABCD的测向边界线LAB、LCD、LBC、LAD方程，然后建立四边形定位模糊区ABCD的外接椭圆簇方程，得到基于置信椭圆的定位模糊区模型。第二步是实验求解，得到置信椭圆定位模糊区模型中测向误差系数k的取值范围为1.75～2.5，此时置信椭圆定位模糊区可以满足工程实践中定位模糊区置信度达到0.95以上的要求，同时具有较高定位精度。
图3是本发明第二步中步骤2.2的实验方法流程图。
2.2 获得测向误差系数k的精确取值范围：
2.2.1 初始化变量：设置测向站O1和测向站O2之间的距离变量R＝10km；设置测向站O1方位测量误差的标准差变量Δθ1＝0.1°，测向站O2方位测量误差的标准差变量Δθ2＝0.1°；设置第一测角变量θ1＝30°，第二测角补角变量θ'2＝30°；设置测向误差系数变量k＝1。
2.2.2 将R、Δθ1、Δθ2、θ1、θ'2、k带入基于置信椭圆的定位模糊区模型(即公式(4))，采用matlab数学软件进行模型求解，得到相应的置信椭圆定位模糊区，用蒙特卡罗方法对该置信椭圆定位模糊区进行仿真实验，根据仿真结果求解出该定位模糊区的置信度并记录，得到置信度和面积。
2.2.3 令k＝k+0.1。
2.2.4 判断k是否小于等于3，若是，则转入步骤2.2.2；若否，则令k＝1，θ'2＝θ'2+5°，转步骤2.2.5。
2.2.5判断θ'2是否小于等于140°，若是，则转步骤2.2.6；若否，则令θ'2＝30°，θ1＝θ1+5°，转步骤2.2.7。
2.2.6 判断θ'2+θ1＝180°是否成立，若成立，则θ'2＝θ'2+5°，转步骤2.2.5；若否，则转步骤2.2.2。
2.2.7 判断θ1是否小于等于140°，若是，则转步骤2.2.2；若否，则θ1＝30°，Δθ1＝Δθ1+0.05°，Δθ2＝Δθ2+0.05°，转步骤2.2.8。
2.2.8 判断Δθ1是否小于等于1°且Δθ2是否小于等于1°，若是，则转步骤2.2.2；若否，则Δθ1＝Δθ2＝0.1°，R＝R+5，转步骤2.2.9。
2.2.9 判断R是否小于等于100，若是，则转入步骤2.2.2；若否，则实验结束，执行2.2.10。
2.2.10 对比实验记录的结果，得出测向误差系数k为1.75～2.5的结论。
图4是本发明第一步建立置信椭圆定位模糊区模型示意图。
与图1的区别是增加了四边形定位模糊区ABCD的最小外界椭圆，即置信椭圆定位模糊区。
图5是本发明第二步中当测向误差系数k＝1.75，k＝2，k＝2.5时置信椭圆定位模糊区示意图。取两个测向站之间的距离即R＝30km，Δθ1＝Δθ2＝0.5°，θ1＝30°，θ2＝150°。采用蒙特卡洛方法进行仿真实验。图5 (a)为k＝1.75，即测向误差系数为1.75时，进行10000次仿真实验得到的置信椭圆定位模糊区示意图。其中椭圆为置信椭圆定位模糊区，散布的点为实验仿真得到的定位点的分布。实验结果表明其置信度为0.9518，置信椭圆定位模糊区面积为0.5079km2。图5 (b)为k＝2，即测向误差系数为2时，进行10000次仿真实验得到的置信椭圆定位模糊区示意图。其中椭圆为置信椭圆定位模糊区，散布的点为实验仿真得到的定位点的分布。实验结果表明其置信度为0.9834，置信椭圆定位模糊区面积为0.6636km2。图5 (c)为k＝2.5，即测向误差系数为2.5时，进行10000次仿真实验得到的置信椭圆定位模糊区示意图。其中椭圆为置信椭圆定位模糊区，散布的点为实验仿真得到的定位点的分布。实验结果表明其置信度为0.9981，置信椭圆定位模糊区面积为1.0373km2。三种情况均能满足定位置信度0.95的要求，且定位面积较小，定位精度高。定位置信度随着k的增加而增加，最后逐渐趋近于1，定位面积也随着k的增加而增加，而且增加速度逐渐变大。
图6是三种定位模糊区方法(即基于交叉定位的四边形定位模糊区方法、基于GDOP的圆形定位模糊区方法和本发明基于置信椭圆的定位模糊区方法)仿真结果对比示意图。取两个测向站之间的距离即R＝30km，Δθ1＝Δθ2＝0.5°，θ1＝30°，θ2＝150°，k＝2。采用蒙特卡洛方法进行10000次仿真实验，得到定位点的分布即图中散布的点。然后根据基于交叉定位的四边形定位模糊区方法做出四边形定位模糊区即图6 中四边形ABCD；根据基于GDOP的圆形定位模糊区方法做出圆形定位模糊区即图6 中的圆；根据本发明提出的基于置信椭圆的定位模糊区方法做出相应的置信椭圆定位模糊区即图6 中的椭圆。实验结果为：四边形ABCD定位模糊区置信度为0.9072,定位模糊区面积为0.4221km2；圆形定位模糊区置信度为0.9981,定位模糊区面积为1.4972km2；置信椭圆定位模糊区置信度为0.9855,定位模糊区面积为0.6636km2。综上可见，置信椭圆定位模糊区方法相较于现有的两种方法，不仅可以保证较高的定位置信水平，同时定位面积较小，定位精度高。
图7是本发明第二步中步骤2.2的部分实验数据。其中R为两测向站之间的距离，Δθ1和Δθ2为测角标准差，θ1为测向站O1对目标O3辐射源实施测向的均方位角；θ2为测向站O2对目标O3辐射源实施测向的均方位角(其中θ2'为θ2的补角)。k为测向误差系数。S为相应的置信椭圆定位模糊区面积，P为相应的置信椭圆定位模糊区的置信度。其中，R在取值范围分别取为10Km，50km，100km，Δθ1和Δθ2在取值范围内分别取为0.5°，0.1°，1°，θ1为30°，θ2'在取值范围内取分别取为30°，60°，120°，k在取值范围内分别取为1.75,2,2.5。分别得到81组相应置信椭圆定位模糊区的面积S和置信度P。观察数据可发现，这些置信椭圆定位模糊区的置信度P均满足0.95～0.99且定位面积较小，定位精度高。
Claims (3)
1.一种基于置信椭圆的二站纯方位定位模糊区方法，其特征在于包括以下步骤：
第一步、建立基于置信椭圆定位模糊区模型，方法是：
1.1对于由测向站O1和测向站O2构成的二站纯方位定位系统，通过测量获取测向站O1和测向站O2对目标定位点的测量参数，即测角θ1和θ2以及测向站O1和测向站O2之间的距离R，Δθ1为测向站O1方位测量误差的标准差，Δθ2为测向站O2方位测量误差的标准差，Δθ1和Δθ2通过阅读测向站中传感器的使用规格说明得到，若两测向站使用的传感器相同，则Δθ1＝Δθ2；
1.2建立四边形定位模糊区ABCD的测向边界线LAB、LCD、LBC、LAD方程：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>A</mi>
<mi>B</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>k&Delta;&theta;</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<mfrac>
<mi>R</mi>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>C</mi>
<mi>D</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>k&Delta;&theta;</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<mfrac>
<mi>R</mi>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>B</mi>
<mi>C</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msub>
<mi>k&Delta;&theta;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<mfrac>
<mi>R</mi>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>A</mi>
<mi>D</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>+</mo>
<msub>
<mi>k&Delta;&theta;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<mfrac>
<mi>R</mi>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中θ2'为θ2的补角，k为测向误差系数。
1.3建立四边形定位模糊区ABCD的外接椭圆簇方程，得到基于置信椭圆的定位模糊区模型：
四边形定位模糊区ABCD的外接椭圆簇方程为:
LABLCD+λLBCLAD＝0 (2)
将式(1)中LAB、LCD、LBC、LAD的直线方程带入式(2)，得到：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mo>&lsqb;</mo>
<mi>y</mi>
<mo>-</mo>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>k&Delta;&theta;</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<mfrac>
<mi>R</mi>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>&lsqb;</mo>
<mi>y</mi>
<mo>-</mo>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>k&Delta;&theta;</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<mfrac>
<mi>R</mi>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mi>&lambda;</mi>
<mo>&lsqb;</mo>
<mi>y</mi>
<mo>+</mo>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msub>
<mi>k&Delta;&theta;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<mfrac>
<mi>R</mi>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>&lsqb;</mo>
<mi>y</mi>
<mo>+</mo>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>+</mo>
<msub>
<mi>k&Delta;&theta;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<mfrac>
<mi>R</mi>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
将式(3)转化为椭圆簇方程形式，得到基于置信椭圆的定位模糊区模型为：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mi>b</mi>
<mo>+</mo>
<mi>&lambda;</mi>
<mi>c</mi>
<mi>d</mi>
<mo>)</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mo>(</mo>
<mi>&lambda;</mi>
<mi>d</mi>
<mo>+</mo>
<mi>&lambda;</mi>
<mi>c</mi>
<mo>-</mo>
<mi>b</mi>
<mo>-</mo>
<mi>a</mi>
<mo>)</mo>
<mi>x</mi>
<mi>y</mi>
<mo>+</mo>
<mo>(</mo>
<mi>&lambda;</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mi>b</mi>
<mi>R</mi>
<mo>-</mo>
<mi>&lambda;</mi>
<mi>c</mi>
<mi>d</mi>
<mi>R</mi>
<mo>)</mo>
<mi>x</mi>
<mo>+</mo>
<mo>(</mo>
<mo>-</mo>
<mi>b</mi>
<mi>R</mi>
<mo>/</mo>
<mn>2</mn>
<mo>-</mo>
<mi>a</mi>
<mi>R</mi>
<mo>/</mo>
<mn>2</mn>
<mo>-</mo>
<mi>&lambda;</mi>
<mi>d</mi>
<mi>R</mi>
<mo>/</mo>
<mn>2</mn>
<mo>-</mo>
<mi>&lambda;</mi>
<mi>c</mi>
<mi>R</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
<mi>y</mi>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>abR</mi>
<mn>2</mn>
</msup>
<mo>/</mo>
<mn>4</mn>
<mo>+</mo>
<msup>
<mi>&lambda;cdR</mi>
<mn>2</mn>
</msup>
<mo>/</mo>
<mn>4</mn>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:a＝tan(θ1+kΔθ1),b＝tan(θ1-kΔθ1),c＝tan(θ2'-kΔθ2),d＝tan(θ2'+kΔθ2)
第二步、实验求解，得到置信椭圆定位模糊区模型中测向误差系数k的精确取值范围，确定置信椭圆定位模糊区，同时确定目标辐射源O3位于置信椭圆定位模糊区中的置信度：
2.1获得测向误差系数k的粗略取值范围为1～3；
2.2获得测向误差系数k的精确取值范围为1.75～2.5，得到了k，也即确定了基于置信椭圆的定位模糊区模型的各参数，四边形定位模糊区ABCD的外接椭圆也就确定了，也即得到了置信椭圆定位模糊区。
2.如权利要求1所述的基于置信椭圆的二站纯方位定位模糊区方法，其特征在于2.1步获得测向误差系数k的粗略取值范围的方法是：
2.1.1取测向站O1和测向站O2之间的距离R＝10km，取两测向站的测角标准差Δθ1＝Δθ2＝0.5°，取测角θ1＝30°，θ'2＝30°，测向误差系数k依次取值为1、1.5、2、2.5、3，代入公式(4)，共进行5次实验，采用matlab数学软件进行模型求解，每次实验得到相应的置信椭圆定位模糊区，用蒙特卡罗方法对置信椭圆定位模糊区进行仿真实验，根据5次仿真实验结果求解出该定位模糊区的置信度和面积并记录，得到5个置信度和5个面积；
2.1.2设置R＝20km，按步骤2.1.1的方法再进行5次实验，得到5个置信度和5个面积；
2.1.3设置R＝30km，按步骤2.1.1的方法再进行5次实验，得到5个置信度和5个面积；
2.1.4观察上述15个置信度，发现这些置信度均处于0.5～1这个范围，且置信度随着k的增加而增加，当k小于1时，置信度小于0.5，当k接近3时，置信度逐渐趋近于1；得出结论：k的粗略取值范围为1～3。
3.如权利要求1所述的基于置信椭圆的二站纯方位定位模糊区方法，其特征在于2.2步获得测向误差系数k的精确取值范围的方法是：
2.2.1初始化变量：设置测向站O1和测向站O2之间的距离变量R＝10km；设置测向站O1方位测量误差的标准差变量Δθ1＝0.1°，测向站O2方位测量误差的标准差变量Δθ2＝0.1°；设置第一测角变量θ1＝30°，第二测角补角变量θ'2＝30°；设置测向误差系数变量k＝1；
2.2.2将R、Δθ1、Δθ2、θ1、θ'2、k带入公式(4)，采用matlab数学软件进行模型求解，得到相应的置信椭圆定位模糊区，用蒙特卡罗方法对该置信椭圆定位模糊区进行仿真实验，根据仿真结果求解出该定位模糊区的置信度并记录，得到置信度和面积；
2.2.3令k＝k+0.1；
2.2.4判断k是否小于等于3，若是，则转入步骤2.2.2；若否，则令k＝1，θ'2＝θ'2+5°，转步骤2.2.5；
2.2.5判断θ'2是否小于等于140°，若是，则转步骤2.2.6；若否，则令θ'2＝30°，θ1＝θ1+5°，转步骤2.2.7；
2.2.6判断θ'2+θ1＝180°是否成立，若成立，则θ'2＝θ'2+5°，转步骤2.2.5；若否，则转步骤2.2.2；
2.2.7判断θ1是否小于等于140°，若是，则转步骤2.2.2；若否，则θ1＝30°，Δθ1＝Δθ1+0.05°，Δθ2＝Δθ2+0.05°，转步骤2.2.8；
2.2.8判断Δθ1是否小于等于1°且Δθ2是否小于等于1°，若是，则转步骤2.2.2；若否，则Δθ1＝Δθ2＝0.1°，R＝R+5，转步骤2.2.9；
2.2.9判断R是否小于等于100，若是，则转入步骤2.2.2；若否，则实验结束，执行2.2.10；
2.2.10对比实验记录的结果，发现当测向误差系数k在1.75～2.5之间取值可以满足工程实践中定位模糊区置信度在0.950～0.999的要求，同时具有较高定位精度，得出结论：k的精确取值范围为1.75～2.5。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510731399.5A CN105487048B (zh) | 2015-11-02 | 2015-11-02 | 基于置信椭圆的二站纯方位定位模糊区方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510731399.5A CN105487048B (zh) | 2015-11-02 | 2015-11-02 | 基于置信椭圆的二站纯方位定位模糊区方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105487048A CN105487048A (zh) | 2016-04-13 |
| CN105487048B true CN105487048B (zh) | 2018-02-02 |
Family
ID=55674157
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510731399.5A Active CN105487048B (zh) | 2015-11-02 | 2015-11-02 | 基于置信椭圆的二站纯方位定位模糊区方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105487048B (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102175991A (zh) * | 2011-01-16 | 2011-09-07 | 哈尔滨工程大学 | 基于最大定位似然传感器配置的目标定位方法 |
| CN102508197A (zh) * | 2011-09-29 | 2012-06-20 | 哈尔滨工程大学 | 基于信道容量的无源目标定位方法 |
| CN103440353A (zh) * | 2013-03-05 | 2013-12-11 | 上海大学 | 一种基于vc++软件的针对双基站对无线电源三角定位及其误差分析方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4266669B2 (ja) * | 2003-02-28 | 2009-05-20 | 日本電気株式会社 | バイスタティック方位検出システム及び検出方法 |
-
2015
- 2015-11-02 CN CN201510731399.5A patent/CN105487048B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102175991A (zh) * | 2011-01-16 | 2011-09-07 | 哈尔滨工程大学 | 基于最大定位似然传感器配置的目标定位方法 |
| CN102508197A (zh) * | 2011-09-29 | 2012-06-20 | 哈尔滨工程大学 | 基于信道容量的无源目标定位方法 |
| CN103440353A (zh) * | 2013-03-05 | 2013-12-11 | 上海大学 | 一种基于vc++软件的针对双基站对无线电源三角定位及其误差分析方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105487048A (zh) | 2016-04-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106126794B (zh) | 一种三角网格曲面下切割面动态调整的射线寻迹方法 | |
| CN104318622B (zh) | 一种室内场景非均匀三维点云数据的三角网格建模方法 | |
| CN106446422B (zh) | 一种基于对数似然估计的无源定位跟踪新方法 | |
| CN103558602B (zh) | 一种用于多基地声纳配置方式的模拟退火定位方法 | |
| CN106019217A (zh) | 基于aoa的二维无线传感器网络半定规划定位方法 | |
| CN105242252B (zh) | 基于图像匹配的下降轨聚束sar雷达定位方法 | |
| CN107690184A (zh) | 联合tdoa‑aoa无线传感器网络半定规划定位方法 | |
| CN104848833A (zh) | 一种电子经纬仪和激光跟踪仪联合测量系统的建立方法 | |
| CN108413988A (zh) | 机器人末端经纬仪坐标系快速标定方法 | |
| CN103292784B (zh) | 一种基于声线轨迹的水下地表地形可视性分析方法 | |
| CN109917333A (zh) | 融合aoa观测量与tdoa观测量的无源定位方法 | |
| CN103217688A (zh) | 一种基于不规则三角网机载激光雷达点云平差计算方法 | |
| CN107271957A (zh) | 基于tdoa和toa的室内三维定位方法 | |
| CN108168559A (zh) | 一种基于分布式天线的室内定位系统及方法 | |
| CN103529428A (zh) | 一种基于空间十元阵的点声源被动声定位方法 | |
| CN104569915A (zh) | 多光电探测系统中基于目标运动模型的定位方法 | |
| CN105372676A (zh) | 一种三维场景导航信号多径预测方法 | |
| CN105627916A (zh) | 一种建立跟踪仪地理坐标系及六自由度测量的方法 | |
| Gao et al. | TDOA collaborative localization algorithm based on PSO and Newton iteration in WGS-84 coordinate system | |
| CN105487048B (zh) | 基于置信椭圆的二站纯方位定位模糊区方法 | |
| CN110087307A (zh) | 基于测距修正的水下传感器网络定位方法 | |
| CN107884745A (zh) | 一种针对uwb定位的大气误差处理方法 | |
| CN107504862A (zh) | 一种全方位高精度激光定位方法 | |
| CN110501672A (zh) | 稀疏los/nlos场景下基于距离观测值定位的约束l1范数方法 | |
| CN102564383B (zh) | 基于轴面角确定三轴磁传感器轴间角的方法 |
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 |