CN114020003A - 用于测控天线海上轴系参数标校的无人机航路规划方法 - Google Patents

用于测控天线海上轴系参数标校的无人机航路规划方法 Download PDF

Info

Publication number
CN114020003A
CN114020003A CN202111595157.XA CN202111595157A CN114020003A CN 114020003 A CN114020003 A CN 114020003A CN 202111595157 A CN202111595157 A CN 202111595157A CN 114020003 A CN114020003 A CN 114020003A
Authority
CN
China
Prior art keywords
coordinate system
waypoint
pitch angle
coordinates
waypoints
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
Application number
CN202111595157.XA
Other languages
English (en)
Other versions
CN114020003B (zh
Inventor
顾新锋
刘童岭
徐荣
田英国
刘洋
黄坤
周新舟
韩子龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Satelite Maritime Measurement And Control Bureau
Original Assignee
China Satelite Maritime Measurement And Control Bureau
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Satelite Maritime Measurement And Control Bureau filed Critical China Satelite Maritime Measurement And Control Bureau
Priority to CN202111595157.XA priority Critical patent/CN114020003B/zh
Publication of CN114020003A publication Critical patent/CN114020003A/zh
Application granted granted Critical
Publication of CN114020003B publication Critical patent/CN114020003B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/08Control of attitude, i.e. control of roll, pitch, or yaw
    • G05D1/0808Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/10Simultaneous control of position or course in three dimensions
    • G05D1/101Simultaneous control of position or course in three dimensions specially adapted for aircraft
    • G05D1/106Change initiated in response to external conditions, e.g. avoidance of elevated terrain or of no-fly zones

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Navigation (AREA)

Abstract

本申请涉及一种用于测控天线海上轴系参数标校的无人机航路规划方法。所述方法根据轴系参数标校需要满足不同方位俯仰角度的跟踪要求,根据标校参数采用保持相对距离不变的绕圈升高和绕圈下降的方式设计航路,得到航点的相对位置关系,并通过航点的相对位置关系和测量船当前地理位置,计算地理坐标系下的航点位置坐标,再将计算得到的地理坐标系下的航点位置坐标转换成地面站航路文件,上传到无人机,完成航路规划。利用本发明给出的计算模型,可以通过编程的方式实现快速高精度的跟踪无人机轴系参数标校应用的航路设计,提高低成本无人机用于轴系参数标校的效率。

Description

用于测控天线海上轴系参数标校的无人机航路规划方法
技术领域
本发明涉及航天测控技术领域,特别涉及一种用于测控天线海上轴系参数标校的无人机航路规划方法。
背景技术
为了实现海上高精度的角度测量,船载测控***每次在进行海上测控之前需要通对设备进行标校。对于无线电设备,标校内容为相位标定和轴系参数的标定。其中,海上轴系参数的标定,目前主要有两种方法:一种是通过跟踪过境卫星的方式直接得到电轴修正参数;另一种方法是先以经纬仪为基准,通过同步测恒星法得到各测控天线标校电视光轴修正参数,然后再以光轴为基准,通过放球的方式进行光电偏差的标定,得到电轴的参数。第一种方法虽然可以直接得到电轴参数,但受过境星频点的限制,如果其频点与本次测控频点不一致,还需要通过放球的方式得到不同频点之间光电参数的差异进行等效计算。轴系参数标校涉及时间长,效率低。
随着近几年无人机技术的发展,其价格不断降低、应用领域不断拓展,同样也可以应用于海上标校。通过无人机携带信标进行海上标校,不仅可以利用无人机可回收的优势降低成本,还可以利用无人机地面站和飞机之间的通信链路进行信标频率的远程切换,实现一次飞行完成多个频点的设备标校,提高标校效率;并且无人机飞行平稳,跟踪信号稳定,标校随机误差小,并且可以通过无人机的位置信息计算无人机相对天线的引导角度和距离。
但是利用无人机进行轴系参数标校,需要覆盖不同象限方位和俯仰角度,每次标校,船舶的位置也不同,需要根据实时的船位进行航路设计,传统的航路规划方法很难满足海上轴系参数标校应用的要求。
发明内容
基于此,有必要针对上述技术问题,提供一种低成本、高效率的用于测控天线海上轴系参数标校的无人机航路规划方法。
一种用于测控天线海上轴系参数标校的无人机航路规划方法,所述方法包括:
获取测量船当前地理位置,设置标校距离、最低俯仰角、最高俯仰角、上升过程飞行圈数、下降过程飞行圈数、上升过程航点数以及下降过程航点数。
根据所述最低俯仰角、所述最高俯仰角、所述标校距离、所述上升过程飞行圈数、所述上升过程航点数、所述下降过程飞行圈数以及所述下降过程航点数,采用保持各航点与测量船的相对距离不变的绕圈升高和绕圈下降的方式设计航路,得到所有航点的相对位置关系。
根据所述测量船当前地理位置和所述相对位置关系,得到每个航点的地理坐标位置;
设置返航点地理坐标位置,并将所述返航点地理坐标位置和所有航点的地理坐标位置转换成地面站航路文件,将所述地面站航路文件导入到地面站,并通过地面站软件上传至无人机,完成无人机航路规划。
在其中一个实施例中,航点的相对位置关系包括:航点与测量船的相对距离、航点的方位角和俯仰角。
根据所述最低俯仰角、所述最高俯仰角、所述标校距离、所述上升过程飞行圈数、所述上升过程航点数、所述下降过程飞行圈数以及所述下降过程航点数,采用保持各航点与测量船的相对距离不变的绕圈升高和绕圈下降的方式设计航路,得到所有航点的相对位置关系,包括:
设置初始方位角,设置各航点与测量船的相对距离为标校距离。
当为上升过程时,根据初始方位角、所述最低俯仰角、所述最高俯仰角、所述上升过程飞行圈数以及所述上升过程航点数,得到上升航点的方位角和俯仰角;上升航点的方位角和俯仰角的计算表达式为:
Ai=A0+360×i×(N1/M1)
Ei=E0+i×(Emax-E0)/M1
其中,A0为初始方位角;Ai为第i个上升航点的方位角,i为大于0且小于等于上升过程航点数的整数;N1为上升过程飞行圈数;M1为上升过程航点数;Emax为最高俯仰角;E0为最低俯仰角;Ei为第i个上升航点的俯仰角。
当为下降过程时,根据上升过程中第M1个上升航点的方位角、所述最低俯仰角、所述最高俯仰角、所述下降过程飞行圈数以及所述下降过程航点数,得到下降航点的方位角和俯仰角;下降航点的方位角和俯仰角的计算表达式为:
Aj=AM1+360×j×(N2/M2)
Ej=Emax-i×(Emax-E0)/M2
其中,AM1为第M1个上升航点的方位角;Aj为第j个下降航点的方位角,j为大于0且小于等于下降过程航点数的整数;N2为下降过程飞行圈数;M2为下降过程航点数;Ej为第j个下降航点的俯仰角。
在其中一个实施例中,根据所述测量船当前地理位置和所述相对位置关系,得到每个航点的地理坐标位置,包括:
将航点的相对位置关系转换成航点的直角坐标系坐标,并将测量船当前地理位置转换成当前位置的地心固连坐标系坐标。
根据航点的直角坐标系坐标与测量船当前位置的地心固连坐标系坐标,计算航点的地心固连坐标系坐标。
将航点的地心固连坐标系坐标转换成航点的地理坐标位置。
在其中一个实施例中,将航点的相对位置关系转换成航点的直角坐标系坐标,并将测量船当前地理位置转换成当前位置的地心固连坐标系坐标,包括:
将上升过程航点的相对位置关系按照俯仰角增大的顺序进行排序,下降过程航点的相对位置关系按照俯仰角减小的顺序进行排序,并将得到的两个序列组合,得到航点的相对位置关系序列。
根据航点的相对位置关系序列和相对位置关系到直角坐标系坐标的转换公式,得到航点的直角坐标系坐标;所述相对位置关系到直角坐标系坐标的转换公式为:
Figure BDA0003430296460000031
其中,(R,Ak,Ek)为航点的相对位置关系序列中第k个航点的相对位置关系,
其中,R为航点与测量船的相对距离,Ak为第k个航点的方位角,Ek为第k个航点的俯仰角,k为航点的序号,k为大于等于1小于等于上升过程航点数与下降过程航点数之和的整数;(Xkc,Ykc,Zkc)为第k个航点的直角坐标系坐标。
根据测量船当前地理位置和从地理坐标系到地心固连坐标系的转换公式,得到测量船当前地理位置的地心固连坐标系坐标;从地理坐标系到地心固连坐标系的转换公式为:
Figure BDA0003430296460000041
其中,(L0,B0,H0)为测量船当前地理位置,(Xoc,Yoc,Zoc)为测量船当前地理位置的地心固连坐标系坐标,e为地球椭球偏心率,e2=f(2-f),f为地球椭球的极扁率,N为测量船当前地理位置的卯酉圈曲率。
在其中一个实施例中,根据航点的直角坐标系坐标与测量船当前位置的地心固连坐标系坐标,计算航点的地心固连坐标系坐标,步骤中航点的地心固连坐标系坐标的计算表达式为:
Figure BDA0003430296460000042
其中,(Xkt,Ykt,Zkt)为第k个航点的地心固连坐标系坐标。
在其中一个实施例中,将航点的地心固连坐标系坐标转换成航点的地理坐标位置,包括:
将航点位置的地心固连坐标系坐标转换成航点的经度、纬度和高度。
航点的经度为:
Figure BDA0003430296460000043
航点的纬度和高度是通过迭代计算公式得到的,所述迭代计算公式为:
Figure BDA0003430296460000051
Figure BDA0003430296460000052
Figure BDA0003430296460000053
其中,u为大于等于1的整数,初始值为
Figure BDA0003430296460000054
N0=a,
Figure BDA0003430296460000055
a、b分别为地球椭球的半长轴和半短轴;迭代收敛条件为:|Hu-Hu-1|<ε1,|Bu-Bu-1|<ε2,其中ε1和ε2为迭代收敛精度。
在其中一个实施例中,设置返航点地理坐标位置,并将所述返航点地理坐标位置和所有航点的地理坐标位置转换成地面站航路文件,将所述地面站航路文件导入到地面站,并通过地面站软件上传至无人机,完成无人机航路规划,包括:
根据测量船的位置设置返航点地理坐标位置。
根据地面站软件的航路文件的接口要求,将所有航点的地理坐标位置转化成航路文件。
将所述地面站航路文件导入到地面站软件中,并通过地面站软件与无人机的通信链路,将所述地面站航路文件加载到无人机,完成无人机航路规划。
上述一种用于测控天线海上轴系参数标校的无人机航路规划方法,所述方法根据轴系参数标校需要满足不同方位俯仰角度的跟踪要求,根据标校参数采用保持相对距离不变的绕圈升高和绕圈下降的方式设计航路,得到航点的相对位置关系,并通过航点的相对位置关系和测量船当前地理位置,计算地理坐标系下的航点位置坐标,再将计算得到的地理坐标系下的航点位置坐标转换成地面站航路文件,上传到无人机,完成航路规划。利用本发明给出的计算模型,可以通过编程的方式实现快速高精度的跟踪无人机轴系参数标校应用的航路设计,提高低成本无人机用于轴系参数标校的效率
附图说明
图1为一个实施例中一种用于测控天线海上轴系参数标校的无人机航路规划方法的流程示意图;
图2为另一个实施例中跟踪无人机轴系参数标校应用航路设计流程;
图3为另一个实施例中轴系标校应用的等距绕圈航路规划结果。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
在一个实施例中,如图1所示,提供了一种用于测控天线海上轴系参数标校的无人机航路规划方法,该方法包括以下步骤:
步骤100:获取测量船当前地理位置,设置标校距离、最低俯仰角、最高俯仰角、上升过程飞行圈数、下降过程飞行圈数、上升过程航点数以及下降过程航点数。
具体的,测量船当前地理位置P0(L0,B0,H0)直接通过读取无人机下传并在地面站上显示为当前位置得到。
标校参数包括:标校距离R、最低俯仰角E0以及最高俯仰角Emax。标校参数根据具体设备标校的要求设置。
上升过程飞行圈数N1和下降过程飞行圈数N2根据无人机的飞行性能设置,具体值可以是整数,也可取实数。
在本实施例中在方位上每10度设置1个航点,则上升过程航点数=上升过程飞行圈数×(360/10)+1),下降过程航点数=下降过程飞行圈数×(360/10),最高点包含在上升过程中,不包含在下降过程中。
步骤102:根据最低俯仰角、最高俯仰角、标校距离、上升过程飞行圈数、上升过程航点数、下降过程飞行圈数以及下降过程航点数,采用保持各航点与测量船的相对距离不变的绕圈升高和绕圈下降的方式设计航路,得到所有航点的相对位置关系。
步骤104:根据测量船当前地理位置和相对位置关系,得到每个航点的地理坐标位置。
步骤106:设置返航点地理坐标位置,并将返航点地理坐标位置和所有航点的地理坐标位置转换成地面站航路文件,将地面站航路文件导入到地面站,并通过地面站软件上传至无人机,完成无人机航路规划。
上述一种用于测控天线海上轴系参数标校的无人机航路规划方法中,所述方法根据轴系参数标校需要满足不同方位俯仰角度的跟踪要求,根据标校参数采用保持相对距离不变的绕圈升高和绕圈下降的方式设计航路,得到航点的相对位置关系,并通过航点的相对位置关系和测量船当前地理位置,计算地理坐标系下的航点位置坐标,再将计算得到的地理坐标系下的航点位置坐标转换成地面站航路文件,上传到无人机,完成航路规划。利用本发明给出的计算模型,可以通过编程的方式实现快速高精度的跟踪无人机轴系参数标校应用的航路设计,提高低成本无人机用于轴系参数标校的效率
在其中一个实施例中,航点的相对位置关系包括:航点与测量船的相对距离、航点的方位角和俯仰角;步骤102包括:设置初始方位角,设置各航点与测量船的相对距离为标校距离;当为上升过程时,根据初始方位角、最低俯仰角、最高俯仰角、上升过程飞行圈数以及上升过程航点数,得到上升航点的方位角和俯仰角;上升航点的方位角和俯仰角的计算表达式为:
Ai=A0+360×i×(N1/M1)
Ei=E0+i×(Emax-E0)/M1
其中,A0为初始方位角;Ai为第i个上升航点的方位角,i为大于0且小于等于上升过程航点数的整数;N1为上升过程飞行圈数;M1为上升过程航点数;Emax为最高俯仰角;E0为最低俯仰角;Ei为第i个上升航点的俯仰角。
当为下降过程时,根据上升过程中第M1个上升航点的方位角、最低俯仰角、最高俯仰角、下降过程飞行圈数以及下降过程航点数,得到下降航点的方位角和俯仰角;下降航点的方位角和俯仰角的计算表达式为:
Aj=AM1+360×j×(N2/M2)
Ej=Emax-i×(Emax-E0)/M2
其中,AM1为第M1个上升航点的方位角;Aj为第j个下降航点的方位角,j为大于0且小于等于下降过程航点数的整数;N2为下降过程飞行圈数;M2为下降过程航点数;Ej为第j个下降航点的俯仰角。
在其中一个实施例中,步骤104包括:将航点的相对位置关系转换成航点的直角坐标系坐标,并将测量船当前地理位置转换成当前位置的地心固连坐标系坐标;根据航点的直角坐标系坐标与测量船当前位置的地心固连坐标系坐标,计算航点的地心固连坐标系坐标;将航点的地心固连坐标系坐标转换成航点的地理坐标位置。
在其中一个实施例中,步骤104包括:将上升过程航点的相对位置关系按照俯仰角增大的顺序进行排序,下降过程航点的相对位置关系按照俯仰角减小的顺序进行排序,并将得到的两个序列组合,得到航点的相对位置关系序列;根据航点的相对位置关系序列和相对位置关系到直角坐标系坐标的转换公式,得到航点的直角坐标系坐标;相对位置关系到直角坐标系坐标的转换公式为:
Figure BDA0003430296460000081
其中,(R,Ak,Ek)为航点的相对位置关系序列中第k个航点的相对位置关系,
其中,R为航点与测量船的相对距离,Ak为第k个航点的方位角,Ek为第k个航点的俯仰角,k为航点的序号,k为大于等于1小于等于上升过程航点数与下降过程航点数之和的整数;(Xkc,Ykc,Zkc)为第k个航点的直角坐标系坐标;
根据测量船当前地理位置和从地理坐标系到地心固连坐标系的转换公式,得到测量船当前地理位置的地心固连坐标系坐标;从地理坐标系到地心固连坐标系的转换公式为:
Figure BDA0003430296460000091
其中,(L0,B0,H0)为测量船当前地理位置,(Xoc,Yoc,Zoc)为测量船当前地理位置的地心固连坐标系坐标,e为地球椭球偏心率,e2=f(2-f),f为地球椭球的极扁率,f=1/298.257223565,N为测量船当前地理位置的卯酉圈曲率
Figure BDA0003430296460000092
在其中一个实施例中,步骤104中航点的地心固连坐标系坐标的计算表达式为:
Figure BDA0003430296460000093
其中,(Xkt,Ykt,Zkt)为第k个航点的地心固连坐标系坐标。
在其中一个实施例中,将航点的地心固连坐标系坐标转换成航点的地理坐标位置,包括:
将航点位置的地心固连坐标系坐标转换成航点的经度、纬度和高度;航点的经度为:
Figure BDA0003430296460000094
航点的纬度和高度是通过迭代计算公式得到的,迭代计算公式为:
Figure BDA0003430296460000095
Figure BDA0003430296460000096
Figure BDA0003430296460000097
其中,u为大于等于1的整数,初始值为
Figure BDA0003430296460000101
N0=a,
Figure BDA0003430296460000102
a、b分别为地球椭球的半长轴和半短轴;迭代收敛条件为:|Hu-Hu-1|<ε1,|Bu-Bu-1|<ε2,其中ε1和ε2为迭代收敛精度。作为优选ε1=10-3,ε2=10-4
在其中一个实施例中,步骤106包括:根据测量船的位置设置返航点地理坐标位置;根据地面站软件的航路文件的接口要求,将所有航点的地理坐标位置转化成航路文件;将地面站航路文件导入到地面站软件中,并通过地面站软件与无人机的通信链路,将地面站航路文件加载到无人机,完成无人机航路规划。
应该理解的是,虽然图1的流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,图1中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
在一个具体的实施例中,如图2所示,提供了一种跟踪无人机进行轴系参数标校的航路规划方法,包括以下步骤:
S1、获取测量船当前地理位置P0(L0,B0,H0),设置标校距离为R(本实例中设置为5000m),最低俯仰角为E0(本实例中设置为10°),最高俯仰角为Emax(本实例中设置为60°),上升过程飞行圈数为N1(本实例中设置为1),下降过程飞行圈数为N2(本实例中设置为1),上升过程航点数为M1(本实例中设置M1为36),下降过程航点数为M2(本实例中设置M2为36)。
S2、计算第1个航点的地理位置P1(L1,B1,H1),该位置相对测量船距离为R,俯仰角为E0,方位角A0
具体的,将相对位置关系转换成航点的具体的地理坐标位置,先将相对位置关系(R,A0,E0)转换成直角坐标系坐标(Xc,Yc,Zc),将测量船当前地理位置P0(L0,B0,H0)转换成地心固连坐标系坐标(Xoc,Yoc,Zoc),再计算第一个航点对应的地心固连坐标系坐标Pt(Xt,Yt,Zt),将Pt转换地理坐标位置P1(L1,B1,H1);坐标转换的数学计算表达式如下:
1)相对位置关系(R,A0,E0)转换成直角坐标系坐标(Xc,Yc,Zc)
Figure BDA0003430296460000111
2)测量船当前地理位置P0(L0,B0,H0)转换成地心固连坐标系坐标(Xoc,Yoc,Zoc)
Figure BDA0003430296460000112
3)根据1)和2)的结果计算第一个航点对应的地心固连坐标系坐标Pt(Xt,Yt,Zt)
Figure BDA0003430296460000113
4)将航点位置的地心固连坐标系Pt(Xt,Yt,Zt)坐标转换成航点的经度、纬度和高度。经度的计算公式为,
Figure BDA0003430296460000114
纬度和高度采用迭代法计算,初始值设置为N0=a,
Figure BDA0003430296460000115
Figure BDA0003430296460000116
随后,每次迭代按下述公式进行:
Figure BDA0003430296460000117
Figure BDA0003430296460000121
Figure BDA0003430296460000122
迭代收敛条件为:|Hu-Hu-1|<ε1,|Bu-Bu-1|<ε2,其中ε1和ε2为迭代收敛精度。这里取ε1=10-3,ε2=10-4
S3、计算上升过程中的第i(i=1,2,...,M1-1)个上升航点Pi(Li,Bi,Hi),该位置相对测量船距离为R,俯仰角为Ei,方位角为Ai
S4、计算俯仰角最高点对应的航点PM(LM,BM,HM),该位置相对测量船距离为R,俯仰角为Emax,方位角为AM
S5、计算下降升过程中的第j(j=1,2,...,M2)个下降航点Pj(Lj,Bj,Hj),该位置相对测量船距离为R,俯仰角为Ej,方位角为Aj
S6、设置反航点PL(LL,BL,HL)。
S7、将计算得到的航点位置信息转换成地面站可识别的航路文件,导入到地面站,并通过地面站软件上传至无人机完成规划。
通过上述步骤可完成低成本无人机1个架次的航路飞行。第2个架次飞行时,采用同样的方法设计航路,所不同的是,第一个航点的方位角相比第1个架次,增加180度,这样,通过两个架次飞行,每个象限的跟踪俯仰角均匀覆盖E0到Emax。图3给出了测量船位置为(东经100°,南纬12°)时,按照本实例设计的航路,无人机以顺时针方向飞行,其中三角形的点表示无人机上升过程中的航点,圆形的点表示无人机下降过程中的航点。
本发明通过设计满足轴系标校的应用的等距离变高度和方向的相对位置关系,再将相对位置关系转换成航点坐标,进一步转换成航点文件,上传至无人机实现航路规划,每个架次飞行1圈上升,1圈下降,通过2个架次飞行就可以实现不同象限均匀覆盖E0到Emax俯仰角。可有效提高航路设计的精度和效率,实现跟踪无人机进行轴系参数标校的应用,有助于提高标校效率和标校精度,降低标校成本。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。

Claims (7)

1.一种用于测控天线海上轴系参数标校的无人机航路规划方法,其特征在于,所述方法包括:
获取测量船当前地理位置,设置标校距离、最低俯仰角、最高俯仰角、上升过程飞行圈数、下降过程飞行圈数、上升过程航点数以及下降过程航点数;
根据所述最低俯仰角、所述最高俯仰角、所述标校距离、所述上升过程飞行圈数、所述上升过程航点数、所述下降过程飞行圈数以及所述下降过程航点数,采用保持各航点与测量船的相对距离不变的绕圈升高和绕圈下降的方式设计航路,得到所有航点的相对位置关系;
根据所述测量船当前地理位置和所述相对位置关系,得到每个航点的地理坐标位置;
设置返航点地理坐标位置,并将所述返航点地理坐标位置和所有航点的地理坐标位置转换成地面站航路文件,将所述地面站航路文件导入到地面站,并通过地面站软件上传至无人机,完成无人机航路规划。
2.根据权利要求1所述的方法,其特征在于,航点的相对位置关系包括:航点与测量船的相对距离、航点的方位角和俯仰角;
根据所述最低俯仰角、所述最高俯仰角、所述标校距离、所述上升过程飞行圈数、所述上升过程航点数、所述下降过程飞行圈数以及所述下降过程航点数,采用保持各航点与测量船的相对距离不变的绕圈升高和绕圈下降的方式设计航路,得到所有航点的相对位置关系,包括:
设置初始方位角,设置各航点与测量船的相对距离为标校距离;
当为上升过程时,根据初始方位角、所述最低俯仰角、所述最高俯仰角、所述上升过程飞行圈数以及所述上升过程航点数,得到上升航点的方位角和俯仰角;上升航点的方位角和俯仰角的计算表达式为:
Ai=A0+360×i×(N1/M1)
Ei=E0+i×(Emax-E0)/M1
其中,A0为初始方位角;Ai为第i个上升航点的方位角,i为大于0且小于等于上升过程航点数的整数;N1为上升过程飞行圈数;M1为上升过程航点数;Emax为最高俯仰角;E0为最低俯仰角;Ei为第i个上升航点的俯仰角;
当为下降过程时,根据上升过程中第M1个上升航点的方位角、所述最低俯仰角、所述最高俯仰角、所述下降过程飞行圈数以及所述下降过程航点数,得到下降航点的方位角和俯仰角;下降航点的方位角和俯仰角的计算表达式为:
Figure FDA0003430296450000022
Ej=Emax-i×(Emax-E0)/M2
其中,
Figure FDA0003430296450000023
为第M1个上升航点的方位角;Aj为第j个下降航点的方位角,j为大于0且小于等于下降过程航点数的整数;N2为下降过程飞行圈数;M2为下降过程航点数;Ej为第j个下降航点的俯仰角。
3.根据权利要求2所述的方法,其特征在于,根据所述测量船当前地理位置和所述相对位置关系,得到每个航点的地理坐标位置,包括:
将航点的相对位置关系转换成航点的直角坐标系坐标,并将测量船当前地理位置转换成当前位置的地心固连坐标系坐标;
根据航点的直角坐标系坐标与测量船当前位置的地心固连坐标系坐标,计算航点的地心固连坐标系坐标;
将航点的地心固连坐标系坐标转换成航点的地理坐标位置。
4.根据权利要求3所述的方法,其特征在于,将航点的相对位置关系转换成航点的直角坐标系坐标,并将测量船当前地理位置转换成当前位置的地心固连坐标系坐标,包括:
将上升过程航点的相对位置关系按照俯仰角增大的顺序进行排序,下降过程航点的相对位置关系按照俯仰角减小的顺序进行排序,并将得到的两个序列组合,得到航点的相对位置关系序列;
根据航点的相对位置关系序列和相对位置关系到直角坐标系坐标的转换公式,得到航点的直角坐标系坐标;所述相对位置关系到直角坐标系坐标的转换公式为:
Figure FDA0003430296450000021
其中,(R,Ak,Ek)为航点的相对位置关系序列中第k个航点的相对位置关系,其中,R为航点与测量船的相对距离,Ak为第k个航点的方位角,Ek为第k个航点的俯仰角,k为航点的序号,k为大于等于1小于等于上升过程航点数与下降过程航点数之和的整数;(Xkc,Ykc,Zkc)为第k个航点的直角坐标系坐标;
根据测量船当前地理位置和从地理坐标系到地心固连坐标系的转换公式,得到测量船当前地理位置的地心固连坐标系坐标;从地理坐标系到地心固连坐标系的转换公式为:
Figure FDA0003430296450000031
其中,(L0,B0,H0)为测量船当前地理位置,(Xoc,Yoc,Zoc)为测量船当前地理位置的地心固连坐标系坐标,e为地球椭球偏心率,e2=f(2-f),f为地球椭球的极扁率,N为测量船当前地理位置的卯酉圈曲率。
5.根据权利要求4所述的方法,其特征在于,根据航点的直角坐标系坐标与测量船当前位置的地心固连坐标系坐标,计算航点的地心固连坐标系坐标,步骤中航点的地心固连坐标系坐标的计算表达式为:
Figure FDA0003430296450000032
其中,(Xkt,Ykt,Zkt)为第k个航点的地心固连坐标系坐标。
6.根据权利要求5所述的方法,其特征在于,将航点的地心固连坐标系坐标转换成航点的地理坐标位置,包括:
将航点位置的地心固连坐标系坐标转换成航点的经度、纬度和高度;
航点的经度为:
Figure FDA0003430296450000033
航点的纬度和高度是通过迭代计算公式得到的,所述迭代计算公式为:
Figure FDA0003430296450000041
Figure FDA0003430296450000042
Figure FDA0003430296450000043
其中,u为大于等于1的整数,初始值为
Figure FDA0003430296450000044
N0=a,
Figure FDA0003430296450000045
a、b分别为地球椭球的半长轴和半短轴;迭代收敛条件为:|Hu-Hu-1|<ε1,|Bu-Bu-1|<ε2,其中ε1和ε2为迭代收敛精度。
7.根据权利要求1所述的方法,其特征在于,设置返航点地理坐标位置,并将所述返航点地理坐标位置和所有航点的地理坐标位置转换成地面站航路文件,将所述地面站航路文件导入到地面站,并通过地面站软件上传至无人机,完成无人机航路规划,包括:
根据测量船的位置设置返航点地理坐标位置;
根据地面站软件的航路文件的接口要求,将所有航点的地理坐标位置转化成航路文件;
将所述地面站航路文件导入到地面站软件中,并通过地面站软件与无人机的通信链路,将所述地面站航路文件加载到无人机,完成无人机航路规划。
CN202111595157.XA 2021-12-23 2021-12-23 用于测控天线海上轴系参数标校的无人机航路规划方法 Active CN114020003B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111595157.XA CN114020003B (zh) 2021-12-23 2021-12-23 用于测控天线海上轴系参数标校的无人机航路规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111595157.XA CN114020003B (zh) 2021-12-23 2021-12-23 用于测控天线海上轴系参数标校的无人机航路规划方法

Publications (2)

Publication Number Publication Date
CN114020003A true CN114020003A (zh) 2022-02-08
CN114020003B CN114020003B (zh) 2024-04-19

Family

ID=80069012

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111595157.XA Active CN114020003B (zh) 2021-12-23 2021-12-23 用于测控天线海上轴系参数标校的无人机航路规划方法

Country Status (1)

Country Link
CN (1) CN114020003B (zh)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20130068399A (ko) * 2011-12-15 2013-06-26 한국전자통신연구원 레이더 시스템의 오차 보정 장치 및 방법
CN103207622A (zh) * 2012-01-11 2013-07-17 波音公司 飞行器的安静着陆姿态调节装置
CN106054920A (zh) * 2016-06-07 2016-10-26 南方科技大学 一种无人机飞行路径规划方法和装置
CN107515003A (zh) * 2017-07-19 2017-12-26 中国南方电网有限责任公司超高压输电公司检修试验中心 一种规划飞机巡视输电线路飞行航线的方法
CN109062233A (zh) * 2018-08-15 2018-12-21 广东电网有限责任公司 一种输电线路无人机自动驾驶巡检方法
RU2692350C1 (ru) * 2018-07-19 2019-06-24 Федеральное государственное бюджетное учреждение науки Институт астрономии Российской академии наук Способ высокоточного позиционирования аппарата на поверхности Луны и устройство для его осуществления
CN109947123A (zh) * 2019-02-27 2019-06-28 南京航空航天大学 一种基于视线导引律的无人机路径跟踪与自主避障方法
CN110806210A (zh) * 2019-11-18 2020-02-18 中国航空工业集团公司沈阳飞机设计研究所 一种无人机校准测向设备航线规划方法及其飞行控制方法
CN110908405A (zh) * 2019-12-18 2020-03-24 中国人民解放军总参谋部第六十研究所 一种固定翼无人机进行同心圆飞行时的控制方法
RU2749194C1 (ru) * 2020-12-15 2021-06-07 Общество с ограниченной ответственностью "Опытно-конструкторское бюро УЗГА" (ООО "ОКБ УЗГА") Способ дистанционного определения координат местоположения наземного (надводного) объекта
CN113093801A (zh) * 2021-04-02 2021-07-09 南京航空航天大学 一种多无人机协同编队的航路文件制作方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20130068399A (ko) * 2011-12-15 2013-06-26 한국전자통신연구원 레이더 시스템의 오차 보정 장치 및 방법
CN103207622A (zh) * 2012-01-11 2013-07-17 波音公司 飞行器的安静着陆姿态调节装置
CN106054920A (zh) * 2016-06-07 2016-10-26 南方科技大学 一种无人机飞行路径规划方法和装置
CN107515003A (zh) * 2017-07-19 2017-12-26 中国南方电网有限责任公司超高压输电公司检修试验中心 一种规划飞机巡视输电线路飞行航线的方法
RU2692350C1 (ru) * 2018-07-19 2019-06-24 Федеральное государственное бюджетное учреждение науки Институт астрономии Российской академии наук Способ высокоточного позиционирования аппарата на поверхности Луны и устройство для его осуществления
CN109062233A (zh) * 2018-08-15 2018-12-21 广东电网有限责任公司 一种输电线路无人机自动驾驶巡检方法
CN109947123A (zh) * 2019-02-27 2019-06-28 南京航空航天大学 一种基于视线导引律的无人机路径跟踪与自主避障方法
CN110806210A (zh) * 2019-11-18 2020-02-18 中国航空工业集团公司沈阳飞机设计研究所 一种无人机校准测向设备航线规划方法及其飞行控制方法
CN110908405A (zh) * 2019-12-18 2020-03-24 中国人民解放军总参谋部第六十研究所 一种固定翼无人机进行同心圆飞行时的控制方法
RU2749194C1 (ru) * 2020-12-15 2021-06-07 Общество с ограниченной ответственностью "Опытно-конструкторское бюро УЗГА" (ООО "ОКБ УЗГА") Способ дистанционного определения координат местоположения наземного (надводного) объекта
CN113093801A (zh) * 2021-04-02 2021-07-09 南京航空航天大学 一种多无人机协同编队的航路文件制作方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
刘童岭 等: "船载雷达无塔动态相位标校方法研究", 科学技术与工程, 31 December 2011 (2011-12-31), pages 5899 *
王益忠;叶玲玲;: "船载测控设备测角精度标定方法", 舰船科学技术, no. 11, 15 November 2012 (2012-11-15), pages 86 *
郭倩倩;周晓慧;: "基于GPS的无人机天线自动跟踪***", 计算机应用与软件, no. 06, 15 June 2017 (2017-06-15), pages 91 *

Also Published As

Publication number Publication date
CN114020003B (zh) 2024-04-19

Similar Documents

Publication Publication Date Title
CN106871927B (zh) 一种无人机光电吊舱安装误差标校方法
US20160187140A1 (en) Coordinated route distribution systems and methods
US6784840B2 (en) Method for determining azimuth and elevation angles using a single axis direction finding system
CN110487301A (zh) 一种雷达辅助机载捷联惯性导航***初始对准方法
CN109478930B (zh) 光束跟踪***及其操作方法
CN110398242A (zh) 一种高旋高过载条件飞行器的姿态角确定方法
CN108489466B (zh) 一种山区航测无人机像控点坐标测量方法
CN105929400A (zh) 基于精确测量的恒定prf下变航速机载sar数据补偿方法
CN113624197A (zh) 一种基于无人机平台的测控天线大盘不水平度测量方法
CN113048939B (zh) 跟踪无人机的光电偏差标定方法、装置和计算机设备
CN114838740A (zh) 一种考虑不同经纬度区域的卫星图像几何定标方法
RU2513900C1 (ru) Способ и устройство определения координат объектов
CN110579744B (zh) 一种基于dom卫片的场景建模仿真方法
CN114001756A (zh) 一种小视场星敏感器外场地面寻星方法
CN111721282B (zh) 一种基于天文导航原理的捷联惯导坐标系动态取齐方法
CN114020003A (zh) 用于测控天线海上轴系参数标校的无人机航路规划方法
CN106643726B (zh) 一种统一惯性导航解算方法
JP2005249629A (ja) 電波到来方向特定システム
CN109471102B (zh) 一种惯组误差修正方法
CN114353831B (zh) 一种船载usb测角误差修正模型及其实现方法
CN108489483B (zh) 一种船载星光定向仪单星次优修正算法
CN114355396A (zh) 一种基于北斗***的usb轴系参数标定方法
CN113721188A (zh) 拒止环境下的多无人机自身定位与目标定位方法
CN112859053A (zh) 一种标定激光雷达时变参数的方法及***
CN111027646A (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