CN106289259A - 一种基于地理坐标的自适应分布式系统定位方法 - Google Patents
一种基于地理坐标的自适应分布式系统定位方法 Download PDFInfo
- Publication number
- CN106289259A CN106289259A CN201610604334.9A CN201610604334A CN106289259A CN 106289259 A CN106289259 A CN 106289259A CN 201610604334 A CN201610604334 A CN 201610604334A CN 106289259 A CN106289259 A CN 106289259A
- Authority
- CN
- China
- Prior art keywords
- aircraft
- max
- fitness
- airport
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000005457 optimization Methods 0.000 claims abstract description 17
- 239000002245 particle Substances 0.000 claims description 16
- 238000005070 sampling Methods 0.000 claims description 7
- 230000003044 adaptive effect Effects 0.000 claims description 5
- 238000003032 molecular docking Methods 0.000 claims 1
- 239000004744 fabric Substances 0.000 abstract 1
- 230000008685 targeting Effects 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000012163 sequencing technique Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明涉及民航空中交通管制技术领域,尤其是一种基于地理坐标的自适应分布式系统定位方法。本发明针对现有技术存在的问题,提供一种基于地理坐标的自适应分布式系统定位方法。本发明针对现有技术存在的问题,提供一种基于地理坐标的自适应分布式系统布站方法。在GDOP图的基础上,从google地图上引入机场地理坐标,计算遮挡因子,加入遮挡因子后,根据飞机停靠点覆盖区域Rj中存在的障碍遮挡,对目标函数进行加权得到当前飞机停靠点位置的适应度f(xi,yi,zi):对适应度f(xi,yi,zi)通过基于粒子群工程优化方法或凸优化优化方法进行优化,得到不同布站方式下的适应度fg(xi,yi,zi);根据fg(xi,yi,zi)对各个站点xi,yi,zi的位置进行排序,并给出飞机停靠与机场内任意位置的定位方案。
Description
技术领域
本发明涉及民航空中交通管制技术领域,尤其是一种基于地理坐标的自适应分布式系统定位方法。
背景技术
航空运输市场迅猛发展,航空运输量的快速增长对空域和机场的使用提出更高的要求。高密度交通流量对空域的需求与可供使用空域资源不足的矛盾日益突出。多点定位(MLAT)系统可以提供高精度、高覆盖、高数据率的监视服务,满足空中交通服务对监视的需求,确保航空运输安全,提高空域容量和航空器运行效率。
MLAT系统利用多个接收站,接收民航飞行器发射的1090MHz应答信号,通过测量信号到达各接收站的时间,由到达时间差(TDOA)定位体制完成定位。对于该体制而言,其接收站点的位置选择尤为重要。现有的选站技术,主要是基于TDOA定位的几何精度因子(GDOP)的选站,不同的站点数和站点位置,对应于不同的GDOP分布图,根据机场对于不同区域的定位精度要求即可完成选站。
然而,该方法存在以下缺点:
1)没有将地形遮挡因素考虑进来;
2)对于某一特定的区域,无法给出各站参与定位时对于最终结果贡献程度的多少,从而无法获知某个站主要对哪些区域起作用;
3)无法对各站进行综合评分,从而对站点位置进行排序、优选;
4)无法给出目标处于任意位置时的最优定位方案(以谁为主站,哪些站 参与定位,定位结果最优)。
发明内容
本发明所要解决的技术问题是:针对现有技术存在的问题,提供一种基于地理坐标的自适应分布式系统定位方法。在GDOP图的基础上,从google地图上引入机场地理坐标,计算遮挡因子,加入遮挡因子后,根据飞机停靠点覆盖区域Rj中存在的障碍遮挡,对目标函数进行加权得到当前飞机停靠点位置的适应度f(xi,yi,zi):对适应度f(xi,yi,zi)通过基于粒子群工程优化方法或凸优化优化方法进行优化,得到不同布站方式下的适应度fg(xi,yi,zi);根据fg(xi,yi,zi)对各个站点xi,yi,zi的位置进行排序,并给出飞机停靠与机场内任意位置的定位方案。
本发明采用的技术方案如下:
一种基于地理坐标的自适应分布式系统布站方法包括:
步骤1:在GDOP图的基础上,从google地图上引入机场地理坐标中站点位置坐标(xi,yi,zi),计算遮挡物遮住飞机和建筑物视线的遮挡因子,并将其引入目标函数,因此目标函数写为如下形式:
其中,σx,σy,σz分别表示飞机停靠位置(x,y,z)到第i站址坐标(xi,yi,zi)的距离在视线方向上的定位误差;其中i=1,...,n,xmin<x<xmax,ymin<y<ymax,zmin<z<zmax,xmin,ymin,zmin是目标分别处于机场地理坐标系(笛卡尔坐标系)X轴、Y轴、Z轴位置分量的下限,xmax,ymax,zmax是目标分别处于X轴、Y轴、Z轴位置分量的上限;
步骤2:加入遮挡因子后,根据飞机停靠点覆盖区域Rj中存在的障碍遮挡,对目标函数进行加权得到当前飞机停靠点位置的适应度f(xi,yi,zi):
其中,Wj∈[0,为遮挡因子;飞机停靠点覆盖区域Rj指的是x∈(xmin,xmax),y∈(ymin,ymax),z∈(zmin,zmax);
步骤3:对适应度f(xi,yi,zi)通过基于粒子群工程优化方法或凸优化优化方法进行优化,得到不同布站方式下的适应度fg(xi,yi,zi);
步骤4:根据fg(xi,yi,zi)对各个站点xi,yi,zi的位置进行排序,并给出飞机停靠与机场内任意位置的定位方案;其中fg(xi,yi,zi)较大值对应的是飞机优先停靠的位置,fg(xi,yi,zi)较小值对应的是飞机次于优先停靠位置的位置。
进一步的,所述计算遮挡因子具体过程是:在GDOP图的基础上,从google地图上引入机场地理坐标,根据机场建筑物的地理坐标,计算其与飞机与站点连线之间的方位角与俯仰角;计算飞机与站点连线被机场建筑物遮挡的程度,如果不存在遮挡,则Wj=1,完全被遮挡,则Wj=0,如果存在部分遮挡,则根据被遮挡程度,Wj在[0,1]之间进行取值。
进一步的,所述步骤3通过基于粒子群工程优化方法优化具体步骤包括:
步骤31:在机场观测区域内飞机停靠点覆盖区域Rj,对飞机停靠位置(x,y,z)进行随机采样,执行步34;
步骤32:计算m个飞机停靠位置(x,y,z)对应的代价函数,令代价函数值最小时,飞机停靠位置(x,y,z)为最优位置pi,设置pi为其当前位置,并计算其适应度f(xi,yi,zi);
步骤33:对步骤32进行n次迭代,得到m个飞机停靠位置最优位置pi,计算这m个最优位置pi的代价函数,令代价函数最小值对应的飞机停靠位置为全局最优位置pg;
步骤34:根据粒子群工程优化方法中的速度和位置更新公式更新所有粒子 的位置,并计算其适应度f(xi,yi,zi);
步骤35:基于步:32及步骤34进行判断,如果当前第i个飞机停靠点位置对应的适应度值大于pi对应的适应度值,则设置pi为第i个飞机的当前停靠位置;否则,不更新pi;
步骤36:基于步骤34及步骤35进行判断,如果pi对应的适应度值大于pg对应的适应度值,则用pg更新pi的值,执行步骤37;否则,不更新pg
步骤37:输出pg为第i个飞机当前停靠位置,同时计算其对应的适应度值fg(xi,yi,zi),否则转步骤34。
进一步的,所述GDOP图,通常可通过Google Earth进行截图;机场轨道采样点,可到现场实际测量,也可在Google Earth上点选;遮挡物采样点,可到现场实际测量,也可在Google Earth上点选;站点位置坐标,可到现场实际测量,也可在Google Earth上点选。
综上所述,由于采用了上述技术方案,本发明的有益效果是:
1)将地形遮挡因素考虑进来了;
2)对于某一特定的区域,可以给出各站参与定位时对于最终结果贡献程度的多少,从而获知某个站主要对哪些区域起作用;
3)可以对各站进行综合评分,从而对站点位置进行排序、优选;
4)可以给出目标处于任意位置时的最优定位方案(以谁为主站,哪些站参与定位,定位结果最优)。
5)本发明考虑的因素更全面,呈现的要素更丰富,特别是细节呈现能力远高于传统技术,能够有力地支撑选站和机场建设,进而提高民航系统的空中交通管制能力。
具体实施方式
本说明书中公开的所有特征,或公开的所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以以任何方式组合。
本说明书中公开的任一特征,除非特别叙述,均可被其他等效或具有类似目的的替代特征加以替换。即,除非特别叙述,每个特征只是一系列等效或类似特征中的一个例子而已。
本发明相关说明:
1、机场地理坐标(包括建筑物坐标,飞机位置坐标,站点位置坐标(xi,yi,zi)。
本发明实施过程是:
步骤1:在GDOP图的基础上,从google地图上引入机场地理坐标中站点位置坐标(xi,yi,zi),计算遮挡物遮住飞机和建筑物视线的遮挡因子,并将其引入目标函数,因此目标函数写为如下形式:
其中,σx,σy,σz分别表示飞机停靠位置(x,y,z)到第i站址坐标(xi,yi,zi)的距离在视线方向上的定位误差;其中i=1,...,n,xmin<x<xmax,ymin<y<ymax,zmin<z<zmax,xmin,ymin,zmin是目标分别处于笛卡尔坐标系中X轴、Y轴、Z轴位置分量的下限,xmax,ymax,zmax是目标分别处于笛卡尔坐标系中X轴、Y轴、Z轴位置分量的上限;
步骤2:加入遮挡因子后,根据飞机停靠点覆盖区域Rj中存在的障碍遮挡,对目标函数进行加权得到当前飞机停靠点位置的适应度f(xi,yi,zi):
其中,Wj∈[0,为遮挡因子;飞机停靠点覆盖区域Rj指的是x∈(xmin,xmax),y∈(ymin,ymax),z∈(zmin,zmax);
步骤3:对适应度f(xi,yi,zi)通过基于粒子群工程优化方法或凸优化优化方 法进行优化,得到不同布站方式下的适应度fg(xi,yi,zi);
步骤4:根据fg(xi,yi,zi)对各个站点xi,yi,zi的位置进行排序,并给出飞机停靠与机场内任意位置的定位方案;其中fg(xi,yi,zi)较大值对应的是飞机优先停靠的位置,fg(xi,yi,zi)较小值对应的是飞机次于优先停靠位置的位置。
所述步骤3通过基于粒子群工程优化方法优化具体步骤包括:
步骤31:在机场观测区域内飞机停靠点覆盖区域Rj,对飞机停靠位置(x,y,z)进行随机采样,执行步34;
步骤32:计算m个飞机停靠位置(x,y,z)对应的代价函数,令代价函数值最小时,飞机停靠位置(x,y,z)为最优位置pi,设置pi为其当前位置,并计算其适应度f(xi,yi,zi);
步骤33:对步骤32进行n次迭代,得到m个飞机停靠位置最优位置pi,计算这m个最优位置pi的代价函数,令代价函数最小值对应的飞机停靠位置为全局最优位置pg;
步骤34:根据粒子群工程优化方法中的速度和位置更新公式更新所有粒子的位置,并计算其适应度f(xi,yi,zi);
步骤35:基于步:32及步骤34进行判断,如果当前第i个飞机停靠点位置对应的适应度值大于pi对应的适应度值,则设置pi为第i个飞机的当前停靠位置;否则,不更新pi;
步骤36:基于步骤34及步骤35进行判断,如果pi对应的适应度值大于pg对应的适应度值,则用pg更新pi的值,执行步骤37;否则,不更新pg
步骤37:输出pg为第i个飞机当前停靠位置,同时计算其对应的适应度值fg(xi,yi,zi),否则转步骤34。
本发明并不局限于前述的具体实施方式。本发明扩展到任何在本说明书中 披露的新特征或任何新的组合,以及披露的任一新的方法或过程的步骤或任何新的组合。
Claims (4)
1.一种基于地理坐标的自适应分布式系统定位方法,其特征在于包括:
步骤1:在GDOP图的基础上,从google地图上引入机场地理坐标中站点位置坐标(xi,yi,zi),计算遮挡物遮住飞机和建筑物视线的遮挡因子,并将其引入目标函数,因此目标函数写为如下形式:
其中,σx,σy,σz分别表示飞机停靠位置(x,y,z)到第i站址坐标(xi,yi,zi)的距离在视线方向上的定位误差;其中i=1,...,n,xmin<x<xmax,ymin<y<ymax,zmin<z<zmax,xmin,ymin,zmin是目标分别处于笛卡尔坐标系中X轴、Y轴、Z轴位置分量的下限,xmax,ymax,zmax是目标分别处于笛卡尔坐标系中X轴、Y轴、Z轴位置分量的上限;
步骤2:加入遮挡因子后,根据飞机停靠点覆盖区域Rj中存在的障碍遮挡,对目标函数进行加权得到当前飞机停靠点位置的适应度f(xi,yi,zi):
其中,Wj∈[0,1为]遮挡因子;飞机停靠点覆盖区域Rj指的是x∈(xmin,xmax),y∈(ymin,ymax),z∈(zmin,zmax);
步骤3:对适应度f(xi,yi,zi)通过基于粒子群工程优化方法或凸优化优化方法进行优化,得到不同布站方式下的适应度fg(xi,yi,zi);
步骤4:根据fg(xi,yi,zi)对各个站点xi,yi,zi的位置进行排序,并给出飞机停靠于机场内任意位置的定位方案;其中fg(xi,yi,zi)较大值对应的是飞机优先停靠的位置,fg(xi,yi,zi)较小值对应的是飞机次于优先停靠位置的位置。
2.根据权利要求1所述的一种基于地理坐标的自适应分布式系统定位方法,其特征在于所述计算遮挡因子具体过程是:在GDOP图的基础上,从 google地图上引入机场地理坐标,根据机场建筑物的地理坐标,计算其与飞机与站点连线之间的方位角与俯仰角;计算飞机与站点连线被机场建筑物遮挡的程度,如果不存在遮挡,则Wj=1,完全被遮挡,则Wj=0,如果存在部分遮挡,则根据被遮挡程度,Wj在[0,1]之间进行取值。
3.根据权利要求1所述的一种基于地理坐标的自适应分布式系统定位方法,其特征在于所述步骤3通过基于粒子群工程优化方法优化具体步骤包括:
步骤31:在机场观测区域内飞机停靠点覆盖区域Rj,对飞机停靠位置(x,y,z)进行随机采样,执行步34;
步骤32:计算m个飞机停靠位置(x,y,z)对应的代价函数,令代价函数值最小时,飞机停靠位置(x,y,z)为最优位置pi,设置pi为其当前位置,并计算其适应度f(xi,yi,zi);
步骤33:对步骤32进行n次迭代,得到m个飞机停靠位置最优位置pi,计算这m个最优位置pi的代价函数,令代价函数最小值对应的飞机停靠位置为全局最优位置pg;
步骤34:根据粒子群工程优化方法中的速度和位置更新公式更新所有粒子的位置,并计算其适应度f(xi,yi,zi);
步骤35:基于步:32及步骤34进行判断,如果当前第i个飞机停靠点位置对应的适应度值大于pi对应的适应度值,则设置pi为第i个飞机的当前停靠位置;否则,不更新pi;
步骤36:基于步骤34及步骤35进行判断,如果pi对应的适应度值大于pg对应的适应度值,则用pg更新pi的值,执行步骤37;否则,不更新pg
步骤37:输出pg为第i个飞机当前停靠位置,同时计算其对应的适应度值fg(xi,yi,zi),否则转步骤34。
4.根据权利要求1所述的一种基于地理坐标的自适应分布式系统定位方法,其特征在于GDOP图,通常可通过Google Earth进行截图;机场轨道采样点,可到现场实际测量,也可在Google Earth上点选;遮挡物采样点,可到现场实际测量,也可在Google Earth上点选;站点位置坐标,可到现场实际测量,也可在Google Earth上点选。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610604334.9A CN106289259B (zh) | 2016-07-28 | 2016-07-28 | 一种基于地理坐标的自适应分布式系统定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610604334.9A CN106289259B (zh) | 2016-07-28 | 2016-07-28 | 一种基于地理坐标的自适应分布式系统定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106289259A true CN106289259A (zh) | 2017-01-04 |
CN106289259B CN106289259B (zh) | 2019-03-08 |
Family
ID=57662756
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610604334.9A Active CN106289259B (zh) | 2016-07-28 | 2016-07-28 | 一种基于地理坐标的自适应分布式系统定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106289259B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110493853A (zh) * | 2019-08-15 | 2019-11-22 | 卓望信息技术(北京)有限公司 | 一种高铁专网基站的排序方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6353744B1 (en) * | 1999-08-16 | 2002-03-05 | Nortel Networks Limited | Method for positioning with low geometrical dilution of position (GDOP) |
CN101706569A (zh) * | 2009-11-05 | 2010-05-12 | 民航数据通信有限责任公司 | 多点定位系统的覆盖精度评估方法和装置 |
CN104023394A (zh) * | 2014-06-17 | 2014-09-03 | 河海大学 | 基于自适应惯性权重的wsn定位方法 |
CN104640206A (zh) * | 2015-02-15 | 2015-05-20 | 中国民航大学 | 一种偏心式星形布站体制的多点定位方法 |
-
2016
- 2016-07-28 CN CN201610604334.9A patent/CN106289259B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6353744B1 (en) * | 1999-08-16 | 2002-03-05 | Nortel Networks Limited | Method for positioning with low geometrical dilution of position (GDOP) |
CN101706569A (zh) * | 2009-11-05 | 2010-05-12 | 民航数据通信有限责任公司 | 多点定位系统的覆盖精度评估方法和装置 |
CN104023394A (zh) * | 2014-06-17 | 2014-09-03 | 河海大学 | 基于自适应惯性权重的wsn定位方法 |
CN104640206A (zh) * | 2015-02-15 | 2015-05-20 | 中国民航大学 | 一种偏心式星形布站体制的多点定位方法 |
Non-Patent Citations (2)
Title |
---|
李东新等: "基于粒子群算法的伪卫星布站优化设计", 《导航定位与授时》 * |
李家蓬等: "多点定位系统的布站设计", 《电光与控制》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110493853A (zh) * | 2019-08-15 | 2019-11-22 | 卓望信息技术(北京)有限公司 | 一种高铁专网基站的排序方法 |
CN110493853B (zh) * | 2019-08-15 | 2021-07-27 | 卓望信息技术(北京)有限公司 | 一种高铁专网基站的排序方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106289259B (zh) | 2019-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP7190329B2 (ja) | 航空機の飛行計画において微細風況を考慮するシステム | |
CN110243359B (zh) | 基于低空风预测模型的安全航迹规划方法 | |
CN110111608B (zh) | 基于雷达轨迹构建机坪场面运动目标运行意图识别的方法 | |
CN104834769B (zh) | 基于航空噪声限制区域拟合的消噪飞行程序设计方法 | |
CN101201870A (zh) | 空中交通飞行态势动态仿真的方法 | |
CN108716919A (zh) | 基于任意多边形无障碍区域的植保无人机航迹规划方法 | |
CN103093649A (zh) | 用于推断飞行器参数的方法和系统 | |
CN106872942A (zh) | 用于分布式多点定位监视系统的定位精度解算方法 | |
CN109141427A (zh) | 在非视距环境下基于距离和角度概率模型的ekf定位方法 | |
CN109814405B (zh) | 一种测控布站方案综合量化评估方法 | |
CN109918818A (zh) | 一种基于性能导航的pbn单发程序保护区划设方法 | |
CN112146650A (zh) | 一种用于无人蜂群协同导航的构型优化方法 | |
CN108860631A (zh) | 一种基于固定翼飞机的性能管理系统 | |
CN108596397B (zh) | 一种基于水文模型多机航线规划方法 | |
CN114860863B (zh) | 一种用于座位选择的位置确定方法、存储介质及电子设备 | |
CN104833343B (zh) | 基于多旋翼飞行器的复杂地形边界与面积估计系统与方法 | |
CN111948675A (zh) | 一种地基增强卫星导航系统的伪卫星布设方法 | |
CN110175793B (zh) | 一种基于地面需求分析的人工增雨飞机航线设计方法 | |
CN107784866A (zh) | 一种飞行管理系统横向导航精度试飞空域规划方法 | |
CN106289259B (zh) | 一种基于地理坐标的自适应分布式系统定位方法 | |
CN112146615B (zh) | 基于多架无人机的边坡监测方法 | |
CN105066986B (zh) | 一种多模式上传体制的航空地面着陆系统及其控制方法 | |
KR20230078097A (ko) | 도심항공교통 실증 관리를 위한 디지털 트윈 기술에 기반한 3d 가시화 방법 | |
CN112002012B (zh) | 城市区域的可视性分析方法 | |
CN117433527A (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 |