CN111221347B - 垂直起降固定翼无人机姿态估计中加速度补偿方法及系统 - Google Patents

垂直起降固定翼无人机姿态估计中加速度补偿方法及系统 Download PDF

Info

Publication number
CN111221347B
CN111221347B CN202010314559.7A CN202010314559A CN111221347B CN 111221347 B CN111221347 B CN 111221347B CN 202010314559 A CN202010314559 A CN 202010314559A CN 111221347 B CN111221347 B CN 111221347B
Authority
CN
China
Prior art keywords
aerial vehicle
unmanned aerial
acceleration
acceleration compensation
shaft
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
Application number
CN202010314559.7A
Other languages
English (en)
Other versions
CN111221347A (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.)
Guangdong Innolux Technology Co ltd
Original Assignee
Guangdong Innolux Technology Co ltd
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 Guangdong Innolux Technology Co ltd filed Critical Guangdong Innolux Technology Co ltd
Priority to CN202010314559.7A priority Critical patent/CN111221347B/zh
Publication of CN111221347A publication Critical patent/CN111221347A/zh
Application granted granted Critical
Publication of CN111221347B publication Critical patent/CN111221347B/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
    • G05D1/0816Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20Instruments for performing navigational calculations

Landscapes

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

Abstract

本发明公开一种垂直起降固定翼无人机姿态估计中加速度补偿方法,该方法基于机体坐标系下的速度信息以及陀螺仪的角速度信息,计算得到三个方向上的线性加速度信息以及科氏加速度,选取两种加速度中影响较大的分量进行组合,再利用罚函数来最终确定无人机的加速度补偿量。本发明提供的加速度补偿方法同时兼顾了垂直起降固定翼无人机在旋翼飞行阶段和固定翼飞行阶段的姿态估计,可实现垂直起降固定翼全程飞行过程中加速度的平滑补偿,从而更准确地反映因姿态变化引起的重力加速度测量值变化,提高无人机姿态估计的精度;此外,可有效避免因传感器数据噪声和运算带来的影响。

Description

垂直起降固定翼无人机姿态估计中加速度补偿方法及系统
技术领域
本发明涉及飞行器导航技术领域,尤其是一种垂直起降固定翼无人机姿态估计中加速度补偿方法及系统。
背景技术
在多旋翼无人机的姿态估计中,常采用互补滤波、卡尔曼滤波等数据融合方法,对陀螺仪和加速度计的数据进行融合,从而得到无人机的姿态信息。
在利用加速度计数据进行姿态解算时,其误差来源于两个方面。其一是加速度计自身的噪声,另一方面是无人机运动造成的加速度叠加在加速度计的测量值上,使得加速度测量值不能准确反映因姿态改变引起的重力加速度计测量值(由于重力加速度是垂直向下的,所以在无人机上固定安装了一个三轴的加速度计,重力加速度在三个轴上的加速度分量可以反映无人机的姿态,从而可以将其用于无人机的姿态解算。如果无人机运动产生了加速度,加速度计测量值就是重力加速度与这个加速度的叠加值,但是姿态解算只希望测量得到的值为重力加速度的值,此时就会带来误差)。对于旋翼机而言,加速度计自身噪声可以通过低通滤波解决,由于飞行速度不快,加速也比较平稳,所以常忽略无人机自身运动引起的加速度,其误差可以通过数据融合算法来消除。
而上述处理方法在垂直起降固定翼无人机的姿态解算时遇到了问题。垂直起降固定翼无人机除了旋翼起飞和旋翼降落阶段,还有固定翼飞行阶段。在固定翼飞行阶段,由于无人机总是以空速作为控制量之一,当在有风环境下飞行时,需要不停地加减地速以保持空速,从而引起较大的加速度叠加在重力加速度的测量值上,无法忽略。另一方面,由于固定翼具有较大的前飞速度,在飞机转弯过程中容易产生科氏加速度,也会对加速度测量值产生干扰。因此,垂直起降固定翼无人机的姿态估计必须考虑以上因素,对加速度计测量数据进行补偿,以满足旋翼飞行和固定翼飞行的姿态估计。
发明内容
本发明提供一种垂直起降固定翼无人机姿态估计中加速度补偿方法及系统,用于克服现有技术中姿态估计方法无法同时满足垂直起降固定翼无人机在旋翼飞行阶段和固定翼飞行阶段的姿态估计要求等缺陷。
为实现上述目的,本发明提出一种垂直起降固定翼无人机姿态估计中加速度补偿方法,包括:
利用GPS获得当前时刻绝对坐标下无人机的速度信息;
获取上一时刻无人机的姿态转换矩阵,根据所述姿态转换矩阵和所述绝对坐标下无人机的速度信息,得到当前时刻机体坐标系下无人机的速度信息;所述机体坐标系包括x轴、y轴和z轴三个坐标轴,x轴为机头方向,z轴为垂直向下,y轴与x轴、z轴构成右手坐标系;
对所述机体坐标系下无人机的速度信息进行求导,得到无人机在x轴、y轴和z轴三个方向上的加速度信息;
根据垂直起降固定翼无人机在固定翼飞行时的状态和所述无人机在x轴、y轴和z轴三个方向上的加速度信息,获得第一加速度补偿量;
利用陀螺仪测量获得当前时刻无人机的角速度,对所述角速度进行低通滤波,得到平滑角速度;
根据科氏加速度计算法则,利用所述机体坐标系下无人机的速度信息和所述平滑角速度,计算获得科氏加速度;
根据科氏加速度对垂直起降固定翼无人机姿态的影响,获得第二加速度补偿量;
根据所述第一加速度补偿量和所述第二加速度补偿量,并利用罚函数,获得最终无人机的加速度补偿量。
为实现上述目的,本发明还提出一种垂直起降固定翼无人机姿态估计中加速度补偿系统,包括:
速度信息获取模块,用于利用GPS获得当前时刻绝对坐标下无人机的速度信息;获取上一时刻无人机的姿态转换矩阵,根据所述姿态转换矩阵和所述绝对坐标下无人机的速度信息,得到当前时刻机体坐标系下无人机的速度信息;所述机体坐标系包括x轴、y轴和z轴三个坐标轴,x轴为机头方向,z轴为垂直向下,y轴与x轴、z轴构成右手坐标系;
速度信息处理模块,用于对所述机体坐标系下无人机的速度信息进行求导,得到无人机在x轴、y轴和z轴三个方向上的加速度信息;根据垂直起降固定翼无人机在固定翼飞行时的状态和所述无人机在x轴、y轴和z轴三个方向上的加速度信息,获得第一加速度补偿量;
角速度信息获取模块,用于利用陀螺仪测量获得当前时刻无人机的角速度,对所述角速度进行低通滤波,得到平滑角速度;
角速度信息处理模块,用于根据科氏加速度计算法则,利用所述机体坐标系下无人机的速度信息和所述平滑角速度,计算获得科氏加速度;根据科氏加速度对垂直起降固定翼无人机姿态的影响,获得第二加速度补偿量;
加速度补偿模块,用于根据所述第一加速度补偿量和所述第二加速度补偿量,并利用罚函数,获得最终无人机的加速度补偿量。
为实现上述目的,本发明还提出一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现上述所述方法的步骤。
与现有技术相比,本发明的有益效果有:
本发明提供的垂直起降固定翼无人机姿态估计中加速度补偿方法,基于机体坐标系下的速度信息以及陀螺仪的角速度信息,计算得到三个方向上的线性加速度信息以及科氏加速度,选取两种加速度中影响较大的分量进行组合,再利用罚函数来最终确定无人机的加速度补偿量。本发明提供的加速度补偿方法同时兼顾了垂直起降固定翼无人机在旋翼飞行阶段和固定翼飞行阶段的姿态估计,可实现垂直起降固定翼全程飞行过程中加速度的平滑补偿,从而更准确地反映因姿态变化引起的重力加速度测量值变化,提高无人机姿态估计的精度;此外,可有效避免因传感器数据噪声和运算带来的影响。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图示出的结构获得其他的附图。
图1为本发明提供的垂直起降固定翼无人机姿态估计中加速度补偿方法流程图。
本发明目的的实现、功能特点及优点将结合实施例,参照附图做进一步说明。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
另外,本发明各个实施例之间的技术方案可以相互结合,但是必须是以本领域普通技术人员能够实现为基础,当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在,也不在本发明要求的保护范围之内。
本发明提出一种垂直起降固定翼无人机姿态估计中加速度补偿方法,如图1所示,包括:
101:利用GPS获得当前时刻绝对坐标下无人机的速度信息;
GPS,全球定位系统,是一种以空中卫星为基础的高精度无线电导航的定位系统,它在全球任何地方以及近地空间都能够提供准确的地理位置、车行速度及精确的时间信息。
绝对坐标,是一个固定的坐标位置,使用它输入的点坐标不会因参照物的不同而不同。
102:获取上一时刻无人机的姿态转换矩阵,根据所述姿态转换矩阵和所述绝对坐标下无人机的速度信息,得到当前时刻机体坐标系下无人机的速度信息;所述机体坐标系包括x轴、y轴和z轴三个坐标轴,x轴为机头方向,z轴为垂直向下,y轴与x轴、z轴构成右手坐标系;
机体坐标系,是指固定在飞行器或者飞机上的遵循右手法则的三维正交直角坐标系,其原点位于飞行器的质心,OX轴位于飞行器参考平面内平行于机身轴线并指向飞行器前方,OY轴垂直于飞行器参考面并指向飞行器右方,OZ轴在参考面内垂直于XOY平面,指向航空器下方。
103:对所述机体坐标系下无人机的速度信息进行求导,得到无人机在x轴、y轴和z轴三个方向上的加速度信息;
104:根据垂直起降固定翼无人机在固定翼飞行时的状态和所述无人机在x轴、y轴和z轴三个方向上的加速度信息,获得第一加速度补偿量;
垂直起降固定翼无人机在固定翼飞行时仅有前飞螺旋桨提供动力进行加速。
105:利用陀螺仪测量获得当前时刻无人机的角速度,对所述角速度进行低通滤波,得到平滑角速度;
陀螺仪,是用高速回转体的动量矩敏感壳体相对惯性空间绕正交于自转轴的一个或二个轴的角运动检测装置。
低通滤波,是一种过滤方式,规则为低频信号能正常通过,而超过设定临界值的高频信号则被阻隔、减弱。
106:根据科氏加速度计算法则,利用所述机体坐标系下无人机的速度信息和所述平滑角速度,计算获得科氏加速度;
在转动参考系中,物体在做牵连运动的同时沿旋转半径做相对运动,由牵连运动和相对运动交互耦合而形成的加速度为科氏加速度。
107:根据科氏加速度对垂直起降固定翼无人机姿态的影响,获得第二加速度补偿量;
科氏加速度主要影响无人机y轴向的加速度。
108:根据所述第一加速度补偿量和所述第二加速度补偿量,并利用罚函数,获得最终无人机的加速度补偿量。
罚函数用于判断最终无人机的加速度补偿量的大小。
步骤101~104,是为了获得第一加速度补偿量;
步骤105~107,是为了获得第二加速度补偿量;
本发明提供的垂直起降固定翼无人机姿态估计中加速度补偿方法中,第一加速度补偿量和第二加速度补偿量的获得是同时进行的,以提高方法效率。
在其中一个实施例中,对于步骤102,获取上一时刻无人机的姿态转换矩阵,根据所述姿态转换矩阵和所述绝对坐标下无人机的速度信息,得到当前时刻机体坐标系下无人机的速度信息,包括:
201:根据历史数据,获取上一时刻无人机的姿态转换矩阵;
202:将所述绝对坐标下无人机的速度信息乘以所述姿态转换矩阵,得到当前时刻机体坐标系下无人机的速度信息,为
Figure 496994DEST_PATH_IMAGE001
(4)
式中,
Figure 846067DEST_PATH_IMAGE002
为机体坐标系下无人机的速度信息;
Figure 549580DEST_PATH_IMAGE003
为绝对坐标下无人机的速度信息;
Figure 501356DEST_PATH_IMAGE004
为姿态变换矩阵;
Figure 628055DEST_PATH_IMAGE005
为机体坐标系;
Figure 7083DEST_PATH_IMAGE006
为惯性系;xyz为机体坐标系中三个坐标轴上的分量;
Figure 197893DEST_PATH_IMAGE007
为速度;
Figure 563147DEST_PATH_IMAGE008
为转置符号。
在另一个实施例中,对于步骤103,对所述机体坐标系下无人机的速度信息进行求导,得到无人机在x轴、y轴和z轴三个方向上的加速度信息,包括:
通过微分环节的传递函数对所述机体坐标系下无人机的速度信息进行求导,得到无人机在x轴、y轴和z轴三个方向上的加速度信息
Figure 421381DEST_PATH_IMAGE009
Figure 971311DEST_PATH_IMAGE010
(5)
式中,
Figure 790363DEST_PATH_IMAGE011
为速度求导;
Figure 818362DEST_PATH_IMAGE012
为机体坐标系;xyz为机体坐标系中三个坐标轴上的分量;
Figure 406469DEST_PATH_IMAGE013
为转置符号。
在某个实施例中,所述微分环节的传递函数为
Figure 392879DEST_PATH_IMAGE014
(1)
式中,
Figure 558281DEST_PATH_IMAGE015
为决定跟踪速度的速度因子;
Figure 29452DEST_PATH_IMAGE016
为自变量,是时间经过拉氏变换后的量。
在下一个实施例中,对于步骤104,考虑到垂直起降固定翼无人机在固定翼飞行时仅有前飞螺旋桨提供动力进行加速,其他两个方向上的加速度多为噪声,因此取加速度信息
Figure 331120DEST_PATH_IMAGE017
的第一个分量作为第一加速度补偿量
Figure 222853DEST_PATH_IMAGE018
,即
Figure 16496DEST_PATH_IMAGE019
(6)
式中,
Figure 651877DEST_PATH_IMAGE020
为速度求导;
Figure 948997DEST_PATH_IMAGE021
为机体坐标系;x为机体坐标系中x坐标轴上的分量;
Figure 277211DEST_PATH_IMAGE022
为转置符号。
在另一个实施例中,对于步骤105,利用陀螺仪测量获得当前时刻无人机的角速度,对所述角速度进行低通滤波,得到平滑角速度,包括:
501:利用陀螺仪测量获得当前时刻无人机的角速度;
502:通过低通滤波的传递函数对所述角速度进行低通滤波处理,得到平滑角速度
Figure 292571DEST_PATH_IMAGE023
Figure 731643DEST_PATH_IMAGE024
(7)
式中,
Figure 273483DEST_PATH_IMAGE025
Figure 883849DEST_PATH_IMAGE026
Figure 511139DEST_PATH_IMAGE027
分别为平滑角速度在机体坐标系中三个坐标轴上的分量;
Figure 222743DEST_PATH_IMAGE028
为转置符号。
在某个实施例中,所述低通滤波的传递函数为
Figure 494456DEST_PATH_IMAGE029
(2)
式中,
Figure 164472DEST_PATH_IMAGE030
为时间常数;k为放大系数;
Figure 13479DEST_PATH_IMAGE031
为自变量。
在下一个实施例中,对于步骤106,根据科氏加速度计算法则,利用所述机体坐标系下无人机的速度信息和所述平滑角速度,计算获得科氏加速度
Figure 404140DEST_PATH_IMAGE032
Figure 654993DEST_PATH_IMAGE033
(8)
式中,
Figure 495910DEST_PATH_IMAGE034
为平滑角速度;
Figure 973159DEST_PATH_IMAGE035
为机体坐标系下无人机的速度信息;
Figure 292145DEST_PATH_IMAGE036
为加速度;
Figure 36984DEST_PATH_IMAGE037
为机体坐标系;xyz为机体坐标系中三个坐标轴上的分量;
Figure 783224DEST_PATH_IMAGE038
为转置符号。
在下一个实施例中,对于步骤107,根据科氏加速度对垂直起降固定翼无人机姿态的影响,获得第二加速度补偿量,
考虑到科氏加速度
Figure 606823DEST_PATH_IMAGE039
主要影响无人机y轴方向的加速度,因此取科氏加速度
Figure 604866DEST_PATH_IMAGE039
的第二个分量作为第二加速度补偿量
Figure 830311DEST_PATH_IMAGE040
,即
Figure 747451DEST_PATH_IMAGE041
(9)
式中,
Figure 933713DEST_PATH_IMAGE042
为加速度;
Figure 860081DEST_PATH_IMAGE043
为机体坐标系;y为机体坐标系中y坐标轴上的分量;
Figure 815399DEST_PATH_IMAGE044
为转置符号。
在下一个实施例中,对于步骤108,根据所述第一加速度补偿量
Figure 434599DEST_PATH_IMAGE045
和所述第二加速度补偿量
Figure 967211DEST_PATH_IMAGE046
,并利用罚函数,获得最终无人机的加速度补偿量
Figure 808522DEST_PATH_IMAGE047
,包括:
801:将所述第一加速度补偿量
Figure 742980DEST_PATH_IMAGE048
和所述第二加速度补偿量
Figure 408447DEST_PATH_IMAGE049
相加,作为加速度补偿量
Figure 162777DEST_PATH_IMAGE050
由于传感器噪声影响,计算得到的加速度补偿量
Figure 430947DEST_PATH_IMAGE051
不可避免地存在误差,对于旋翼阶段而言,如果直接以所述加速度补偿量
Figure 219912DEST_PATH_IMAGE051
进行补偿容易引起额外的姿态估计误差,因此还需以无人机空速
Figure 525122DEST_PATH_IMAGE052
作为判断标准进行补偿,以进一步提高补偿方法的精确度。
802:根据无人机的固定翼飞行巡航速度设置第一空速阈值
Figure 32327DEST_PATH_IMAGE053
和第二空速阈值
Figure 104188DEST_PATH_IMAGE054
,通过空速计测量得到当前时刻无人机的空速
Figure 623025DEST_PATH_IMAGE055
Figure 754929DEST_PATH_IMAGE056
Figure 483851DEST_PATH_IMAGE057
为无人机固定翼飞行的巡航速度。
当空速
Figure 733304DEST_PATH_IMAGE058
大于阈值
Figure 496861DEST_PATH_IMAGE059
时,无人机处于固定翼飞行状态,加速度计测量值需要补偿;
当空速
Figure 409453DEST_PATH_IMAGE060
小于
Figure 891250DEST_PATH_IMAGE061
时,无人机处于旋翼飞行状态,加速度计测量值不需要补偿。
803:根据所述第一空速阈值
Figure 304914DEST_PATH_IMAGE061
、第二空速阈值
Figure 798343DEST_PATH_IMAGE062
和当前时刻无人机的空速
Figure 740892DEST_PATH_IMAGE063
,建立加速度补偿罚函数
Figure 444405DEST_PATH_IMAGE064
加速度补偿罚函数
Figure 396181DEST_PATH_IMAGE065
是为了使得当空速
Figure 9696DEST_PATH_IMAGE066
处于
Figure 654304DEST_PATH_IMAGE067
Figure 221945DEST_PATH_IMAGE068
之间时加速度补偿值能够平滑地过渡。
804:根据所述加速度补偿罚函数加速度补偿和所述加速度补偿量,计算获得最终无人机的加速度补偿量
Figure 977411DEST_PATH_IMAGE069
Figure 711012DEST_PATH_IMAGE070
(10)
式中,
Figure 260942DEST_PATH_IMAGE071
为加速度补偿罚函数;
Figure 673469DEST_PATH_IMAGE072
为加速度补偿量。
将加速度计测量得到的加速度减去
Figure 107992DEST_PATH_IMAGE073
,得到重力加速度在机体坐标系下的修正值,进而结合陀螺仪的数据对无人机姿态进行估计,可以得到比较准确的姿态数据。
在某个实施例中,所述加速度补偿罚函数为
Figure 555154DEST_PATH_IMAGE074
(3)
式中,
Figure 275986DEST_PATH_IMAGE075
为当前时刻无人机的空速;
Figure 582333DEST_PATH_IMAGE076
为第一空速阈值;
Figure 414023DEST_PATH_IMAGE077
为第二空速阈值。
本发明还提供一种垂直起降固定翼无人机姿态估计中加速度补偿系统,包括:
速度信息获取模块,用于利用GPS获得当前时刻绝对坐标下无人机的速度信息;获取上一时刻无人机的姿态转换矩阵,根据所述姿态转换矩阵和所述绝对坐标下无人机的速度信息,得到当前时刻机体坐标系下无人机的速度信息;所述机体坐标系包括x轴、y轴和z轴三个坐标轴,x轴为机头方向,z轴为垂直向下,y轴与x轴、z轴构成右手坐标系;
速度信息处理模块,用于对所述机体坐标系下无人机的速度信息进行求导,得到无人机在x轴、y轴和z轴三个方向上的加速度信息;根据垂直起降固定翼无人机在固定翼飞行时的状态和所述无人机在x轴、y轴和z轴三个方向上的加速度信息,获得第一加速度补偿量;
角速度信息获取模块,用于利用陀螺仪测量获得当前时刻无人机的角速度,对所述角速度进行低通滤波,得到平滑角速度;
角速度信息处理模块,用于根据科氏加速度计算法则,利用所述机体坐标系下无人机的速度信息和所述平滑角速度,计算获得科氏加速度;根据科氏加速度对垂直起降固定翼无人机姿态的影响,获得第二加速度补偿量;
加速度补偿模块,用于根据所述第一加速度补偿量和所述第二加速度补偿量,并利用罚函数,获得最终无人机的加速度补偿量。
本发明还提供一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现上述所述方法的步骤。
以上所述仅为本发明的优选实施例,并非因此限制本发明的专利范围,凡是在本发明的发明构思下,利用本发明说明书及附图内容所作的等效结构变换,或直接/间接运用在其他相关的技术领域均包括在本发明的专利保护范围内。

Claims (7)

1.一种垂直起降固定翼无人机姿态估计中加速度补偿方法,其特征在于,包括:
利用GPS获得当前时刻绝对坐标下无人机的速度信息;
获取上一时刻无人机的姿态转换矩阵,根据所述姿态转换矩阵和所述绝对坐标下无人机的速度信息,得到当前时刻机体坐标系下无人机的速度信息;所述机体坐标系包括x轴、y轴和z轴三个坐标轴,x轴为机头方向,z轴为垂直向下,y轴与x轴、z轴构成右手坐标系;
对所述机体坐标系下无人机的速度信息进行求导,得到无人机在x轴、y轴和z轴三个方向上的加速度信息;
根据垂直起降固定翼无人机在固定翼飞行时的状态和所述无人机在x轴、y轴和z轴三个方向上的加速度信息,获得第一加速度补偿量;
利用陀螺仪测量获得当前时刻无人机的角速度,对所述角速度进行低通滤波,得到平滑角速度;
根据科氏加速度计算法则,利用所述机体坐标系下无人机的速度信息和所述平滑角速度,计算获得科氏加速度;
根据科氏加速度对垂直起降固定翼无人机姿态的影响,获得第二加速度补偿量;
根据所述第一加速度补偿量和所述第二加速度补偿量,并利用罚函数,获得最终无人机的加速度补偿量,包括:
将所述第一加速度补偿量和所述第二加速度补偿量相加,作为加速度补偿量;
根据无人机固定翼巡航飞行速度设置第一空速阈值和第二空速阈值,通过空速计测量得到当前时刻无人机的空速;
根据所述第一空速阈值、第二空速阈值和当前时刻无人机的空速,建立加速度补偿罚函数;
根据所述加速度补偿罚函数和所述加速度补偿量,计算获得最终无人机的加速度补偿量。
2.如权利要求1所述的垂直起降固定翼无人机姿态估计中加速度补偿方法,其特征在于,获取上一时刻无人机的姿态转换矩阵,根据所述姿态转换矩阵和所述绝对坐标下无人机的速度信息,得到当前时刻机体坐标系下无人机的速度信息,包括:
根据历史数据,获取上一时刻无人机的姿态转换矩阵;
将所述绝对坐标下无人机的速度信息乘以所述姿态转换矩阵,得到当前时刻机体坐标系下无人机的速度信息。
3.如权利要求1所述的垂直起降固定翼无人机姿态估计中加速度补偿方法,其特征在于,对所述机体坐标系下无人机的速度信息进行求导,得到无人机在x轴、y轴和z轴三个方向上的加速度信息,包括:
通过微分环节的传递函数对所述机体坐标系下无人机的速度信息进行求导,得到无人机在x轴、y轴和z轴三个方向上的加速度信息。
4.如权利要求1所述的垂直起降固定翼无人机姿态估计中加速度补偿方法,其特征在于,利用陀螺仪测量获得当前时刻无人机的角速度,对所述角速度进行低通滤波,得到平滑角速度,包括:
利用陀螺仪测量获得当前时刻无人机的角速度;
通过低通滤波的传递函数对所述角速度进行低通滤波处理,得到平滑角速度。
5.如权利要求4所述的垂直起降固定翼无人机姿态估计中加速度补偿方法,其特征在于,所述低通滤波的传递函数为
Figure 835461DEST_PATH_IMAGE001
(2)
式中,
Figure 790779DEST_PATH_IMAGE002
为时间常数;k为放大系数;
Figure 144400DEST_PATH_IMAGE003
为自变量。
6. 如权利要求1所述的垂直起降固定翼无人机姿态估计中加速度补偿方法,其特征在于,所述加速度补偿罚函数为
Figure 817957DEST_PATH_IMAGE004
(3)
式中,
Figure 548016DEST_PATH_IMAGE005
为当前时刻无人机的空速;
Figure 918692DEST_PATH_IMAGE006
为第一空速阈值;
Figure 318581DEST_PATH_IMAGE007
为第二空速阈值。
7.一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1~6中任一项所述方法的步骤。
CN202010314559.7A 2020-04-21 2020-04-21 垂直起降固定翼无人机姿态估计中加速度补偿方法及系统 Active CN111221347B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010314559.7A CN111221347B (zh) 2020-04-21 2020-04-21 垂直起降固定翼无人机姿态估计中加速度补偿方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010314559.7A CN111221347B (zh) 2020-04-21 2020-04-21 垂直起降固定翼无人机姿态估计中加速度补偿方法及系统

Publications (2)

Publication Number Publication Date
CN111221347A CN111221347A (zh) 2020-06-02
CN111221347B true CN111221347B (zh) 2020-07-21

Family

ID=70813316

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010314559.7A Active CN111221347B (zh) 2020-04-21 2020-04-21 垂直起降固定翼无人机姿态估计中加速度补偿方法及系统

Country Status (1)

Country Link
CN (1) CN111221347B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2022072338A (ja) * 2020-10-29 2022-05-17 セイコーエプソン株式会社 慣性計測装置
CN112414365B (zh) * 2020-12-14 2022-08-16 广州昂宝电子有限公司 位移补偿方法和设备及速度补偿方法和设备
CN112762918B (zh) * 2020-12-28 2023-07-14 中国科学院长春光学精密机械与物理研究所 集成式数字角速度陀螺仪、预测滤波方法和电子设备
CN113359793B (zh) * 2021-06-01 2022-08-23 北京电子工程总体研究所 一种低速飞行器提高空速控制品质的补偿方法与装置
CN113485393B (zh) * 2021-06-22 2022-06-14 北京三快在线科技有限公司 飞行设备的控制方法、装置、存储介质及飞行设备

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102608596B (zh) * 2012-02-29 2013-06-05 北京航空航天大学 一种用于机载惯性/多普勒雷达组合导航系统的信息融合方法
CN102830414B (zh) * 2012-07-13 2014-12-24 北京理工大学 一种基于sins/gps的组合导航方法
CN107305393A (zh) * 2016-04-20 2017-10-31 比亚迪股份有限公司 无人机及其控制方法
CN106774378B (zh) * 2017-03-01 2019-08-09 杭州谷航科技有限公司 一种无人机飞行控制及定位方法
US10101749B1 (en) * 2017-03-21 2018-10-16 Bell Helicopter Textron Inc. Combined airspeed and inertial data for rotorcraft longitudinal control

Also Published As

Publication number Publication date
CN111221347A (zh) 2020-06-02

Similar Documents

Publication Publication Date Title
CN111221347B (zh) 垂直起降固定翼无人机姿态估计中加速度补偿方法及系统
RU2236697C2 (ru) Резервная система для индикации курса и пространственного положения на самолете
CN103363992B (zh) 基于梯度下降的四旋翼无人机姿态航向参考系统解算方法
CN103837151B (zh) 一种四旋翼飞行器的气动模型辅助导航方法
CN111551175B (zh) 一种航姿参考系统的互补滤波姿态解算方法
CN106444804A (zh) 基于互补滤波算法和串级pid的四旋翼飞行器设计方法
CN104503467A (zh) 基于双核架构的无人机自主起降飞行控制系统
CN111207745B (zh) 一种适用于大机动无人机垂直陀螺仪的惯性测量方法
CN202939489U (zh) 一种多旋翼自动平衡飞行控制器
CN110426032B (zh) 一种解析式冗余的飞行器容错导航估计方法
CN108592911B (zh) 一种四旋翼飞行器动力学模型/机载传感器组合导航方法
GB2444814A (en) Estimating inertial acceleration bias errors
CN107247459A (zh) 抗干扰飞行控制方法及装置
CN108759814B (zh) 一种四旋翼飞行器横滚轴角速度和俯仰轴角速度估计方法
CN111189442A (zh) 基于cepf的无人机多源导航信息状态预测方法
CN111976974A (zh) 一种飞行控制方法、无人机及存储介质
CN114018254B (zh) 一种激光雷达与旋转惯导一体化构架与信息融合的slam方法
CN108693372B (zh) 一种四旋翼飞行器的航向轴角速度估计方法
CN107063248A (zh) 基于旋翼转速的动力学模型辅助惯导的导航方法
CN107807375A (zh) 一种基于多gps接收机的无人机姿态追踪方法及系统
Emran et al. A cascaded approach for quadrotor's attitude estimation
EP4269044A1 (en) Slope location correction method and apparatus, robot and readable storage medium
Li et al. Unmanned aerial vehicle position estimation augmentation using optical flow sensor
Kwon et al. Real-time observability-aware inertia parameter estimation for quadrotors
CN207147501U (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