CN109916406A - 一种基于无人机群的环绕目标定位方法 - Google Patents
一种基于无人机群的环绕目标定位方法 Download PDFInfo
- Publication number
- CN109916406A CN109916406A CN201910024119.5A CN201910024119A CN109916406A CN 109916406 A CN109916406 A CN 109916406A CN 201910024119 A CN201910024119 A CN 201910024119A CN 109916406 A CN109916406 A CN 109916406A
- Authority
- CN
- China
- Prior art keywords
- target
- pitch angle
- azimuth
- matrix equation
- matrix
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种基于无人机群环绕目标的定位方法,通过无人机对目标观测所得的方位角和俯仰角,构建以方位角和俯仰角的正余弦值为参数的线性矩阵方程,利用最小二乘法求解获得目标位置。与过去以测量角度的正切值为参数的侧向交叉定位算法相比,本发明提出的方法有效减小了矩阵方程系数矩阵的条件数,不仅防止了无人机环绕目标定位时可能出现的病态问题,还进一步提高了环绕区域内目标的定位精度。
Description
技术领域
本发明涉及无人机目标定位领域,具体涉及一种基于无人机群的环绕目标定位方法。
背景技术
无人机是近年来发展十分迅速的一项技术,已经在控制精度、带载能力、抗风能力、续航能力上都有了很大的进步。长航时高空无人机可在空中停留几天甚至几个月以获取目标区完整的信息。无人机具有受目标上空云雾的影响较小、可获取高分辨率的图像、较卫星应用灵活、成本低等优点,可以随时选定目标区域,能接近目标区域,提供更加精确的目标位置。
利用无人机进行目标定位在军事和民用方面均有广泛的应用。在民用方面,可以应用于地质勘测、海洋测绘、海上石油作业、飞机失事救援、消防灭火等方面。在军事方面,可以应用于导弹定位、海上军事演练等方面。
目前的目标定位算法多基于雷达或红外,大多需要俯仰角,方位角以及距离信息。在通过纯角度信息进行定位的算法中,侧向交叉定位是一种典型的方法,主要用于二维定位。在观测点分布在目标的同一侧时,侧向交叉定位的定位精度较好。但当观测点环绕目标时,该算法有可能遇到严重的病态问题,其定位精度也有提高空间。对目标进行环绕观测是基于无人机群对目标进行定位的常见形式,需要给出更为合适的定位方法。本次发明主要针对环绕观测环境,基于条件数扰动分析得到传统侧向定位方法出现病态问题的原因,并给出防止病态问题的方法,从而达到了对目标进行环绕定位的精度。
发明内容
本发明的目的是提供一种基于无人机群的环绕目标定位算法,该算法可以实现三维空间的定位,提高基于角度传感器的定位精度,尤其是能够克服算法在特殊角度定位精度较差的缺点。
为了实现上述目的,本发明的技术实现方案如下:
有n(n≥2)架无人机对目标进行定位,每架无人机自带GPS定位系统,可以较为精准的进行自定位,同时携带摄像头对目标进行拍摄,通过对拍摄图片的处理,每个摄像机都可以测得一个方位角和一个俯仰角,所以至少有两架无人机才可以对目标进行三维定位。
进一步的,假设第i架无人机得到的方位角和俯仰角信息为(αi,βi),其中,方位角,俯仰角分别表示为:
进一步的,由于难以将x,y拆分,为了便于求解,将更改为那么俯仰角变为:
进一步的,根据方位角和俯仰角的定义(式(1)和(3)),可得:
-xtanαi+y=yi-xi tanαi (4)
ytanβi-zsinαi=yi tanβi-zi sinαi (5)
进一步的,将上述式(4)和式(5)表示为矩阵形式:
进一步的,将式(6)扩展为i=1,2…n,得到的矩阵方程符合HX=Z形式,可以利用最小二乘法求解落点位置X,从而得到落点位置坐标。
进一步的,当方位角和俯仰角位于特殊位置(90°)时,tan90°无穷大,导致利用最小二乘法得到的解在该位置误差较大,所以将tanα,tanβ分别改为 重新代入到式(1)和式(3),重新进行推导,可以得到新的H矩阵,在利用最小二乘法进行求解后,可大大提升其在特殊位置的定位精度。
本发明与现有技术相比的优点(有益效果)是:将测向定位交叉算法从两维平面扩展为三维,将交叉定位算法中的tan参数改为sin和cos参数。更改后,有效解决了以tan为参数的交叉定位算法中位于特殊位置时误差过大的问题,同时,当目标处于环绕区域内,较原有算法,本发明的算法精度更高。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为两架无人机定位示意图。
图2为相机针孔模型示意图。
图3为本目标定位方法的流程图。
图4为以tan参数的目标定位算法的仿真图。
图5为本发明提出的目标定位算法的仿真图。
具体实施方式
下面根据具体实施实例,并结合附图,进一步说明本发明,使本发明的实现步骤更加易于明白了解。
假设在三维空间中,有i(i=1,2…n)架无人机对目标进行定位,坐标系为OXYZ,目标C的坐标为(x,y,z),各无人机的三维坐标为(xi,yi,zi),本文定义方位角为无人机和目标连线在OXY平面的投影和X正半轴的夹角,俯仰角为无人机和目标连线和OXY平面的夹角。第i架无人机测得的方位角和俯仰角为(αi,βi)。为了便于理解原理,使用两架无人机A和B对目标C进行定位,如图1所示。
如图2所示,每架无人机对目标进行拍摄,从机载摄像机所拍摄的图像中,利用图像处理以及坐标转换,得到方位角和俯仰角(αi,βi)。
如图3所示流程,根据得到的俯仰角和方位角进行建模,进行目标定位。将俯仰角中的更改为并将式(6)中的i扩展为i=1,2…n,得到的矩阵方程为:
令:
则式(7)可以表示为:
HX=Z (11)
考虑到H为奇异矩阵,为了提高定位精度,因此可利用最小二乘法求解落点位置坐标X,从而得到落点位置坐标:
X=(HTH)-1HTZ (12)
当H矩阵形式为式(7)所示形式,由于为无穷大这一特殊性,H矩阵中很多元素均由tan组成,那么在某些特殊位置(即方位角或俯仰角为)时,可能会造成该位置的求解精确度较差的问题,具体推导(为方便推导,取n为2)如下所示:
由式(7)的H矩阵可以得到:
可以看到,tanαi越大,或tanαi值大的观测数目越多,HTH中的前两行越相关,其最小特征值也越趋于0,意味着H的条件数越大,误差传播的上限也越大。所以,当α1,α2,β1,β2中有一个或多个越趋进时,该位置的误差会越大,即该位置的落点位置坐标精确度较低。
因此,本发明提出将tanα(tanβ)替换为形式,方位角和俯仰角的表示形式变为:
重新进行整理,得到矩阵方程如下所示:
新的H矩阵变为:
那么重新计算HTH(同样取n=2)得到:
从式(18)可以看出,此方法可以有效防止式(13)出现的问题,说明采用此方法,可以有效提高无人机处于特殊位置时对目标的估计精度。
仿真实验过程和结果分析:
为更好的展示本发明的优点,在matlab平台上对本发明中提出的算法进行仿真。仿真实验环境设置为由6架无人机组成正六边形编队对航天测量船进行船姿估计,六边形平行于海平面,中心位于测量船惯导甲板坐标系原点的正上方。为了简化分析,我们假设无人机编队与测量船之间的相对位置关系不变,同时编队的所有无人机均可通过光电经纬仪观测到。综合考虑无人机观测、通信能力和安全性的需求,假设无人机编队与测量船相对位置固定,中心处于航天测量船正上方,每架无人机直线距离航天测量船1km,光电经纬仪对每个无人机的观测仰角为45°,无人机的定位精度为0.5m,测量的角度精度为1°。
由图4所示,500m是目标与六边形编队中两架无人机处于一条直线时的位置,可以看到500m处的误差突然变高,充分展示出原目标定位算法在特殊位置时误差较大的问题。从图5可以看出,本发明提出的基于无人机群的环绕目标定位方法有效的解决了上述问题,不管在特殊位置还是其他位置,误差都处于一个稳定的区间,同时通过对图4和图5的比较可以看出,原算法在非特殊位置的精度为40m左右,而本发明提出的算法的精度为25m左右。充分证明了本发明的优点,即不仅解决了特殊位置时误差变高的问题,同时提高了总体的目标定位精度。
Claims (3)
1.一种基于无人机群的环绕目标定位方法,其特征在于,包括如下步骤:
步骤1,多架无人机对目标进行拍摄,通过相机成像模型以及坐标转换,获得多组方位角和俯仰角;
步骤2,n架无人机获得n组方位角和俯仰角(αi,βi),通过对俯仰角βi的定义公式进行变形,获得HX=Z形式的矩阵方程求解问题;
所述的步骤2具体为:
步骤2.1,假设第i架无人机得到的方位角和俯仰角信息为(αi,βi),
其中,
将俯仰角中的替换成为得到方位角,俯仰角的表示形式:
步骤2.2,将方位角和俯仰角的形式进行整理得到:
-x tanαi+y=yi-xi tanαi
y tanβi-z sinαi=yi tanβi-zi sinαi
步骤2.3,将步骤2.2的方程组表示为矩阵形式,i扩展为i=1,2…n,得到的矩阵方程为:
即得到了符合HX=Z形式的矩阵方程;
步骤3,对步骤2得到的矩阵方程中的正切值进行变换,保证H矩阵在方位角或俯仰角为90°时不再是病态方程;
步骤4,根据最终得到的矩阵方程利用最小二乘法进行求解,得到目标的位置坐标。
2.根据权利要求1所述的环绕目标定位方法,其特征在于所述的步骤1中,采用n架无人机同一时刻对同一目标进行拍摄,得到该时刻目标所在的n组方位角和俯仰角,采用图像处理中的目标定位的方法得到目标在图片中的位置,利用相机的成像原理得到角度。
3.根据权利要求1所述的环绕目标定位方法,其特征在于所述的步骤3具体为,将tanα和tanβ分别替换为和对矩阵方程进行处理,得到新的矩阵方程:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910024119.5A CN109916406B (zh) | 2019-01-10 | 2019-01-10 | 一种基于无人机群的环绕目标定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910024119.5A CN109916406B (zh) | 2019-01-10 | 2019-01-10 | 一种基于无人机群的环绕目标定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109916406A true CN109916406A (zh) | 2019-06-21 |
CN109916406B CN109916406B (zh) | 2020-10-13 |
Family
ID=66960232
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910024119.5A Active CN109916406B (zh) | 2019-01-10 | 2019-01-10 | 一种基于无人机群的环绕目标定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109916406B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111123341A (zh) * | 2019-11-15 | 2020-05-08 | 西安电子科技大学 | 无人机群三维协同定位方法 |
CN112229380A (zh) * | 2020-10-15 | 2021-01-15 | 西北工业大学 | 一种基于多旋翼无人机协同的无源目标实时定位方法 |
CN112880685A (zh) * | 2021-01-19 | 2021-06-01 | 中国人民解放军陆军边海防学院 | 无人机群对探测目标的联合定位方法 |
CN113960581A (zh) * | 2021-10-26 | 2022-01-21 | 众芯汉创(北京)科技有限公司 | 一种应用于变电站并结合雷达的无人机目标探测系统 |
CN116164754A (zh) * | 2023-04-25 | 2023-05-26 | 中国电子科技集团公司第五十四研究所 | 一种基于无人机运动状态的云台虚拟化方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102589549A (zh) * | 2011-12-22 | 2012-07-18 | 中国人民解放军海军航空工程学院 | 地球坐标系空间中三站测向交叉定位跟踪算法 |
CN102997911A (zh) * | 2012-12-13 | 2013-03-27 | 中国航空无线电电子研究所 | 一种被动传感器组网探测多目标方法 |
CN107065895A (zh) * | 2017-01-05 | 2017-08-18 | 南京航空航天大学 | 一种植保无人机定高技术 |
CN108571981A (zh) * | 2018-03-28 | 2018-09-25 | 西安爱生技术集团公司 | 一种侦察无人机目标定位集成检校方法 |
-
2019
- 2019-01-10 CN CN201910024119.5A patent/CN109916406B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102589549A (zh) * | 2011-12-22 | 2012-07-18 | 中国人民解放军海军航空工程学院 | 地球坐标系空间中三站测向交叉定位跟踪算法 |
CN102997911A (zh) * | 2012-12-13 | 2013-03-27 | 中国航空无线电电子研究所 | 一种被动传感器组网探测多目标方法 |
CN107065895A (zh) * | 2017-01-05 | 2017-08-18 | 南京航空航天大学 | 一种植保无人机定高技术 |
CN108571981A (zh) * | 2018-03-28 | 2018-09-25 | 西安爱生技术集团公司 | 一种侦察无人机目标定位集成检校方法 |
Non-Patent Citations (6)
Title |
---|
YANG MENG等: "A UAV Aided Method for Ship Attitude Estimation", 《PROCEEDINGS OF THE 37TH CHINESE CONTROL CONFERENCE》 * |
ZONG JUNJUN等: "Algorithm and Accuracy Analysis of Weighted Maximum Likelihood Estimation in Multi-station DF Crossing Localization", 《4TH INTERNATIONAL CONFERENCE ON COMPUTER,MECHATRONICS,CONTROL AND ELECTRONIC ENGINEERING(ICCMCEE 2015)》 * |
杨洪康: "海上编队辐射源融合定位与关联模型研究及仿真实现", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
沈晓峰等: "基于正性约束的测向交叉定位方法", 《电子科技大学学报》 * |
谢鑫: "测向交叉定位最优布站方案分析", 《电子科技》 * |
黄剑伟等: "一种改进的测向交叉定位方法", 《航天电子对抗》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111123341A (zh) * | 2019-11-15 | 2020-05-08 | 西安电子科技大学 | 无人机群三维协同定位方法 |
CN111123341B (zh) * | 2019-11-15 | 2023-03-10 | 西安电子科技大学 | 无人机群三维协同定位方法 |
CN112229380A (zh) * | 2020-10-15 | 2021-01-15 | 西北工业大学 | 一种基于多旋翼无人机协同的无源目标实时定位方法 |
CN112880685A (zh) * | 2021-01-19 | 2021-06-01 | 中国人民解放军陆军边海防学院 | 无人机群对探测目标的联合定位方法 |
CN113960581A (zh) * | 2021-10-26 | 2022-01-21 | 众芯汉创(北京)科技有限公司 | 一种应用于变电站并结合雷达的无人机目标探测系统 |
CN113960581B (zh) * | 2021-10-26 | 2024-06-04 | 众芯汉创(北京)科技有限公司 | 一种应用于变电站并结合雷达的无人机目标探测系统 |
CN116164754A (zh) * | 2023-04-25 | 2023-05-26 | 中国电子科技集团公司第五十四研究所 | 一种基于无人机运动状态的云台虚拟化方法 |
CN116164754B (zh) * | 2023-04-25 | 2023-06-23 | 中国电子科技集团公司第五十四研究所 | 一种基于无人机运动状态的云台虚拟化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109916406B (zh) | 2020-10-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109916406A (zh) | 一种基于无人机群的环绕目标定位方法 | |
CN105353772B (zh) | 一种无人机机动目标定位跟踪中的视觉伺服控制方法 | |
CN106871927B (zh) | 一种无人机光电吊舱安装误差标校方法 | |
CN107247458A (zh) | 无人机视频图像目标定位系统、定位方法及云台控制方法 | |
CN105698762A (zh) | 一种单机航迹上基于不同时刻观测点的目标快速定位方法 | |
CN109597427A (zh) | 一种基于无人机的炸弹随遇攻击规划方法及系统 | |
CN103674029B (zh) | 一种基于水声通信的多艇协同导航编队构型的方法 | |
CN106155081B (zh) | 一种旋翼无人机大范围目标监测和精确定位方法 | |
CN107490364A (zh) | 一种大角度倾斜成像航空相机对地目标定位方法 | |
CN102928861A (zh) | 机载设备用目标定位方法及装置 | |
CN110542407A (zh) | 一种航拍图像任意像素点定位信息获取方法 | |
CN107783544A (zh) | 一种控制单旋翼植保无人机避障飞行的方法 | |
CN117455960B (zh) | 时变观测噪声条件下机载光电系统对地无源定位滤波方法 | |
CN106885573A (zh) | 面向四旋翼飞行器的运动捕捉系统实时测姿方法 | |
CN115876197A (zh) | 一种系留升空光电成像目标定位方法 | |
CN106403934A (zh) | 一种弹载地磁姿态测量处理算法 | |
CN105424060B (zh) | 一种飞行器星敏感器与捷联惯组安装误差的测量方法 | |
CN105261023B (zh) | 一种大气折射修正状态下目标重建方法 | |
CN108106597B (zh) | 全捷联激光导引头在目标出线性视场情况下角度测量方法 | |
CN107884781B (zh) | 一种双无人机跟踪测距方法 | |
CN110579744B (zh) | 一种基于dom卫片的场景建模仿真方法 | |
CN115493598B (zh) | 运动过程中的目标定位方法、装置及存储介质 | |
CN110989655A (zh) | 一种舰载侦察较射无人机起降阶段目标跟踪方法 | |
CN110794391A (zh) | 一种基于无人机集群组网平台的无源定位优化布站方法 | |
CN105243639B (zh) | 卧装的光电吊舱的图像调整方法、装置以及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |