CN106774400A - 一种基于逆动力学的无人机三维轨迹制导方法 - Google Patents
一种基于逆动力学的无人机三维轨迹制导方法 Download PDFInfo
- Publication number
- CN106774400A CN106774400A CN201611233604.6A CN201611233604A CN106774400A CN 106774400 A CN106774400 A CN 106774400A CN 201611233604 A CN201611233604 A CN 201611233604A CN 106774400 A CN106774400 A CN 106774400A
- Authority
- CN
- China
- Prior art keywords
- unmanned plane
- angle
- flight path
- instruction
- error
- 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 46
- 238000004458 analytical method Methods 0.000 claims abstract description 11
- 238000006073 displacement reaction Methods 0.000 claims description 34
- 230000007704 transition Effects 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- NCGICGYLBXGBGN-UHFFFAOYSA-N 3-morpholin-4-yl-1-oxa-3-azonia-2-azanidacyclopent-3-en-5-imine;hydrochloride Chemical compound Cl.[N-]1OC(=N)C=[N+]1N1CCOCC1 NCGICGYLBXGBGN-UHFFFAOYSA-N 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000004088 simulation Methods 0.000 description 4
- 238000006467 substitution reaction Methods 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 206010021703 Indifference Diseases 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
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)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种基于逆动力学的无人机三维轨迹制导方法,属于无人机导航制导与控制技术领域;首先,针对某无人机的某个轨迹段,通过将解析法和数值迭代法相结合,依次对无人机的位置状态方程进行一级逆动力学解算,以及对无人机的地速和航迹角状态方程进行二级逆动力学解算,得到无人机的指令推力、指令迎角和指令航迹滚转角,再将这几个指令量输入已设计好的姿态控制回路,即可控制无人机的实际地速和位置,同时采用PID控制,使无人机的实际位置收敛于航路点之间的参考轨迹上;本发明可同时精确的控制无人机的地速和三维位置,在线调整无人机的期望飞行轨迹,且控制器计算成本低。
Description
技术领域
本发明属于无人机导航制导与控制技术领域,具体涉及一种基于逆动力学的无人机三维轨迹制导方法。
背景技术
无人机又称无人驾驶飞行器,广泛应用于军用和民用领域;无人机制导是指通过指令程序,使无人机沿给定的轨迹飞行。随着无人机的执行任务日趋多样,人们对无人机机动性能的要求也日益提高,传统的二维轨迹制导律已无法满足需求,实现三维轨迹的高精度跟踪控制,能让无人机完成地形回避,编队飞行以及自主空中加油等特殊任务,具有重要的意义。
在无人机的轨迹跟踪控制中,实现无人机直接沿给定航路点飞行具有重要的实用价值,目前很多无人机的参考轨迹设计方法是将各航路点以空间三维曲线进行拟合,或将轨迹拆分成若干段分别进行拟合,这类方法虽然可得到较为光滑的参考轨迹,但算法计算量较大,特别是面对复杂参考轨迹时,其拟合难度将直线上升,这对于机载计算机设备性能较差的低成本无人机来说是非常不利的。
另外,从无人机的飞行动力学角度来看,制导回路所属的各控制量之间存在高度耦合的关系。其中,无人机的位置方程为
式中,(x,y,h)为无人机的三维坐标;
地速及航迹角状态方程为
(Vg,χ,γ)分别为无人机的地速、无人机的航迹偏角和无人机的航迹倾角;(α,β,μ)分别为无人机的迎角、无人机的侧滑角和无人机的航迹滚转角;(T,FD,FY,FL)分别为无人机的发动机推力、无人机的阻力、无人机的侧力和无人机的升力;m和g分别为无人机质量和重力加速度。
由式(1)(2)可见,(Vg,χ,γ)和(T,α,β,μ)均为相互耦合的控制量,若采用传统PID控制,控制效果必然会因耦合作用而受到严重影响。
基于上述分析,所设计的三维轨迹制导方法应满足如下要求:
(1)能够在飞行包线约束条件下,同时控制无人机的地速和位置,精确跟踪三维轨迹,并实现无人机直接沿给定航路点飞行任务;
(2)期望飞行轨迹便于线生成和调整,且所需的时间和计算成本较小;
(3)制导方法物理意义明确,形式简洁,参数整定方便,易于工程实现。
发明内容
本发明的目的是为了解决上述问题,并满足上述要求,提出一种基于逆动力学的无人机三维轨迹制导方法,具体包括如下步骤:
步骤一、针对某个无人机,根据期望轨迹将无人机飞行经过的航路点设定坐标,并划分轨迹段;
期望轨迹由若干航路点联结成的航路段构成;航路段共有n(n≥1)段,对应有n+1个航路点。
步骤二、针对每个航路点,根据无人机的机动性能和所设定好的航路情况,设定该航路点对应的转换半径为d;
步骤三、针对第k段轨迹段,利用无人机的期望轨迹分别计算该轨迹段的实际速度误差和实际位移误差,采用PID控制律加和构成该轨迹段期望的横侧向位移状态的微分项和纵向位移状态的微分项
初始k=1;
首先,无人机的实际速度误差包括横侧向速度误差和纵向速度误差计算如下:
是无人机在第k段轨迹段上设定的指令跟踪横侧向速度;是无人机当前在第k段轨迹段上实际的横侧向速度;是无人机在第k段轨迹段上设定的指令跟踪纵向速度;是无人机当前在第k段轨迹段上实际的纵向速度;Vgk是第k段轨迹段上为无人机设定的指令地速;γk是第k个轨迹段所指向方向的航迹倾角;χk是第k个轨迹段所指向方向的航迹偏角;Vg是无人机当前在第k段轨迹段上实际的地速;γ是无人机当前在第k段轨迹段上实际的航迹倾角;χ是无人机当前在第k段轨迹段上实际的航迹偏角。
然后,无人机的实际位移误差包括横侧向位移误差(yline-y)和纵向位移误差(hline-h);计算如下:
从无人机当前实际位置处作此时对应直线轨迹的空间垂线,得到垂足坐标(xline,yline,hline),该垂线段即为无人机的位置误差总和,分解为横侧向位移误差(yline-y)和纵向位移误差(hline-h)。
最后,PID控制律加和具体为:对无人机的位移误差(yline-y)和(hline-h)采用PI控制,将速度误差和作为微分项,并将PI控制与微分项求和,PID的参数即为速度误差和位移误差在飞行制导律中分别所占的权重,得到期望的横侧向位移状态的微分项和纵向位移状态的微分项
步骤四、将期望的微分项和微分项分别带入横侧向位移状态和纵向位移状态的微分方程中,采用解析法和数值迭代法相结合的方法对微分方程进行第一级逆动力学解算,输出指令航迹偏角χc和航迹倾角γc;
首先,将期望的微分项带入纵向位移状态的微分方程中,采用解析法解算出指令航迹倾角γc,计算如下:
然后,将期望的微分项和指令航迹倾角γc带入到横侧向位移状态的微分方程中,利用数值迭代法解算出指令航迹偏角χc,计算如下:
步骤五、将指令航迹偏角χc和指令航迹倾角γc以及无人机在该航路段上的指令地速Vgk作为下一级逆动力学解算的输入指令值,采用解析法和数值迭代法相结合的方法对无人机地速、航迹倾角和航迹偏角的状态方程进行第二级逆动力学解算,输出指令推力Tc、指令迎角αc和指令航迹滚转角μc;
具体步骤如下:
步骤501、将指令航迹偏角χc,指令航迹倾角γc和指令地速Vgk分别与当前无人机实际的航迹偏角χ,航迹倾角γ和地速Vg作差,并采用PD控制,分别得到期望的航迹偏角微分值期望的航迹倾角微分值和期望的地速微分值
步骤502、利用当前无人机的实际地速Vg和航迹倾角γ,以及期望的航迹偏角微分值和期望的航迹倾角微分值计算指令航迹滚转角μc的解析解;
步骤503、利用指令航迹滚转角μc的解析解,结合期望地速状态方程,代入期望航迹倾角状态方程中,得到关于指令迎角αc的非线性一元方程;
步骤504、用数值迭代法对关于指令迎角αc的非线性一元方程进行解算,得到指令迎角αc。
步骤505、将指令迎角αc代入期望地速状态方程中,解算出指令推力Tc;
期望地速状态方程如下:
步骤六、将输出的指令推力Tc、指令迎角αc、侧滑角0°和航迹滚转角μc作为无人机姿态控制回路的输入,即可实现对无人机三维轨迹的跟踪控制;
步骤七、当无人机飞行至第k段轨迹段对应的转换半径d的范围内时,则无人机飞到所对应的航路点位置,继续将跟踪下一段轨迹,返回至步骤三,直至追踪至最后一个航路点。
本发明的优点在于:
(1)一种基于逆动力学的无人机三维轨迹制导方法,可同时精确的控制无人机的地速和三维位置。
(2)一种基于逆动力学的无人机三维轨迹制导方法,可在线调整无人机的期望飞行轨迹,且控制器计算成本低。
(3)一种基于逆动力学的无人机三维轨迹制导方法,控制结构简单,各部分物理意义明确,便于参数整定。
附图说明
图1是本发明整个制导控制器的结构框图;
图2是本发明一种基于逆动力学的无人机三维轨迹制导方法流程图;
图3是本发明采用的三维制导方法中存在的几何关系示意图;
图4是本发明采用解析法和数值迭代法相结合进行第二级逆动力学解算的流程图;
图5a是无人机实际位置与对应直线航路段之间的距离误差图;
图5b是无人机地速误差图;
图5c是无人机XY平面轨迹跟踪效果图;
图5d是无人机XZ平面轨迹跟踪效果图;
图5e是无人机YZ平面轨迹跟踪效果图。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明一种基于逆动力学的无人机三维轨迹制导方法,通过将解析法和数值迭代法相结合,依次对无人机的位置状态方程进行一级逆动力学解算,以及对无人机的地速和航迹角状态方程进行二级逆动力学解算,得到无人机的指令推力、指令迎角和指令航迹滚转角,再将这几个指令量输入已设计好的姿态控制回路,即可控制无人机的实际地速和位置,同时采用PID控制,使无人机的实际位置收敛于航路点之间的参考轨迹上。
整个制导控制器的结构框图如图1,其整体工作流程为:
针对某个无人机,首先设定好各航路点的坐标,根据指令航路点坐标,根据无人机当前位置(x,y,h)获取计算出期望航迹角,并设定好无人机在各航路段上的指令地速;根据指令地速和期望航迹角,可计算出无人机在横侧向与纵向的指令跟踪地速分量和令其分别于无人机实际的地速分量和分别作差,获取无人机的实际速度误差和再乘以系数KD,作为PID控制器的微分项;
同时根据指令航路点坐标,生成对应的各直线轨迹段,从无人机的当前位置向当前对应的直线航路段作垂线,垂足坐标为(xline,yline,hline),令yline和hline分别与无人机的实际横侧向坐标y与纵向高度h作差,计算无人机在各直线轨迹段上的实际位移误差(yline-y)和(hline-h);并并采用PI控制律,分别乘以系数KP和KI;
接下来要采用PID控制,构成期望的横侧向和纵向位移状态的微分值和作为第一级逆动力学解算的输入,将PI控制项与微分项分别对应加和构成和作为第一级逆动力学解算的输入,采用解析法与数值迭代法相结合的逆动力学方法解算出指令航迹角,作为第二级逆动力学解算的输入,解算出指令推力Tc、指令迎角αc和指令航迹滚转角μc。将Tc输入至无人机动力控制模块,将αc和μc输入至无人机的姿态控制回路中,进一步生成无人机的相应舵角指令,从而同时控制无人机的姿态与地速。
如图2所示,具体包括如下步骤:
步骤一、针对某个无人机,根据期望轨迹将无人机飞行经过的航路点设定坐标,并划分轨迹段;
根据设定的无人机飞行所需要经过的航路点坐标,构成无人机的期望轨迹,并将轨迹划分为若干轨迹段;无人机的期望轨迹由若干航路点所联结成的航路段构成,设航路点按序号可联结为n(n≥1)段直线航路,对应有n+1个航路点。
步骤二、针对每个航路点,根据无人机的机动性能和所设定好的航路情况,设定该航路点对应的转换半径为d;
本实施例中将每个航路点的转换半径均设定为d。
步骤三、针对第k段轨迹段,利用无人机的期望轨迹分别计算该轨迹段的实际速度误差和实际位移误差,采用PID控制律加和构成该轨迹段期望的横侧向位移状态的微分项和纵向位移状态的微分项
将无人机基于航路点的期望轨迹跟踪问题划分为速度误差收敛问题和位置误差收敛问题,分别设计相应的PID控制策略,再将所设计的PID控制律加和构成期望横侧向位移状态和纵向位移状态的微分项;
具体为:
(1)无人机速度误差收敛问题
如图3所示,当无人机UAV收敛于第k(1≤k≤n)段航路上时,无人机UAV的三维坐标为(x,y,h),第k-1个航路点WPk-1(xk-1,yk-1,hk-1)与第k个航路点WPk(xk,yk,hk)之间的航迹偏角χk∈[0,2π)和航迹倾角γk∈(-π/2,π/2)即为此时无人机的期望航迹偏角和航迹倾角,其表达式为:
其中yk-1和yk起到了指令横侧向位移的作用,hk-1和hk起到了指令纵向位移的作用,χk和γk即为所要跟踪的虚拟航迹角;χk是第k个轨迹段所指向方向的航迹偏角;γk是第k个轨迹段所指向方向的航迹倾角;
注意:当k=1时,第k-1个航路点为(xk-1,yk-1,hk-1)=(x0,y0,h0),该点被定义为飞行的起点。
将第k个轨迹段所要跟踪的虚拟航迹角χk和γk以及无人机此时设定的指令地速Vgk代入式(1),得到此时的指令跟踪横侧向速度和纵向速度其与实际速度的误差分别为和表达式为
(2)无人机位置误差收敛问题
如图3所示,从无人机当前实际位置处作此时所对应第k段航路的垂线,可得垂足坐标FP(xline,yline,hline),该垂线段即为无人机的位置误差总和,可分解为横侧向位移误差(yline-y)和纵向位移误差(hline-h)。其中,yline起着指令横侧向位移的作用,hline起着指令纵向位移的作用。
也就是说,yk-1,hk-1,yk,hk,yline和hline均起到了指令位置的作用,其中,yk和hk是驱使无人机收敛到指令横侧向和纵向速度的位置指令,yline和hline是驱使无人机收敛到指令轨迹的位置指令,制导方法中存在的几何关系如图1所示,此时需要对位置指令进行加权求和处理,本发明引入了PID控制的思路来解决这一问题。
对于指令轨迹的收敛问题,其位移误差(yline-y)和(hline-h)采用PI控制,从而使无人机能够无差收敛到指令轨迹上;对于指令速度的收敛问题,通过式(3)和(4)已经求得了横侧向和纵向速度误差和该项可直接作为微分项。将PI控制的结构与微分项求和,其中PID的参数即为速度误差收敛问题和位置误差收敛问题分别在飞行制导律中所占的权重,可得期望的横侧向位移状态的微分项和纵向位移状态的微分项
步骤四:将步骤三中所得的期望状态的微分项和分别带入横侧向位移状态的微分方程和纵向位移状态的微分方程中,采用解析法和数值迭代法相结合的方法对无人机的横侧向位移和纵向位移状态方程进行第一级逆动力学解算,输出指令航迹偏角χc和航迹倾角γc。
用期望横侧向位移状态的微分项和纵向位移状态的微分项替代式(1)中的对应的微分项和
首先,针对期望纵向位移的状态方程,通过解析法解算出指令航迹倾角γc的解析解:
然后,将解算出的γc带入到期望横侧向位移的微分状态方程中,可得仅包含未知数χc的一元方程:
利用数值方法,例如牛顿迭代法解算出指令航迹偏角χc。
步骤五:将步骤四解算出的航迹偏角χc和航迹倾角γc以及设定的无人机指令地速Vgk作为下一级逆动力学解算的输入,采用解析法和数值迭代法相结合的方法对无人机地速、航迹倾角和航迹偏角的状态方程进行第二级逆动力学解算,输出指令推力Tc、指令迎角αc和指令航迹滚转角μc。
如图4所示,具体步骤如下:
步骤501、将指令航迹偏角χc,指令航迹倾角γc和指令地速Vgk分别与当前无人机实际的航迹偏角χ,航迹倾角γ和地速Vg作差,并采用PD控制,分别得到期望的航迹偏角微分值期望的航迹倾角微分值和期望的地速微分值
步骤502、利用当前无人机的实际地速Vg和航迹倾角γ,以及期望的航迹偏角微分值和期望的航迹倾角微分值计算指令航迹滚转角μc的解析解;
由于无人机采用BTT方式转弯,侧滑角β约等于0,此时侧力FY也约等于0,因此式(2)的三个方程可化简为
上式中的微分值和被其期望值 和所替换,联立式(8)和(9),可得指令航迹滚转角的解析解:
步骤503、利用指令航迹滚转角μc的解析解,结合期望地速状态方程,代入期望航迹倾角状态方程中,得到关于指令迎角αc的非线性一元方程;
将解算出的指令航迹滚转角μc和式(7)代入式(9)中,得只含未知数指令迎角αc的非线性一元方程:
注:阻力FD和升力FL在此时均可视为变量αc的函数;
步骤504、用数值迭代法对关于指令迎角αc的非线性一元方程进行解算,得到指令迎角αc。
步骤505、将指令迎角αc代入期望地速状态方程中,解算出指令推力Tc;
期望地速状态方程如下:
将αc代入式(7),即可解算出指令推力Tc。
步骤六、将输出的指令推力Tc、指令迎角αc、侧滑角0°和航迹滚转角μc作为无人机姿态控制回路的输入,即可实现对无人机三维轨迹的跟踪控制;
步骤七:当无人机飞行至距离目前所跟踪轨迹段的对应航路点d的位置时,即可认为无人机已飞到所对应的航路点位置,此时无人机将跟踪下一段轨迹,整个制导算法循环至步骤三,直至追踪至最后一个航路点。
本实施例中,设定无人机以200m/s的地速飞行,初始高度为7010m;采用本发明方法的风扰条件下,获取的无人机三维轨迹跟踪效果图;
如图5a所示,是无人机实际位置与对应直线航路段之间的距离误差图,可以看出无人机在进行恒速爬升及转弯动作时,位置不可避免的存在一定程度的超调,但很快位置误差就能收敛到接近0的程度;
如图5b所示,是无人机地速误差图,可以看出无人机在进行恒速爬升及转弯动作时,会不可避免的存在一定的速度损失,但很快速度误差就能收敛到接近0的程度,且速度损失不大;
如图5c,5d和5e所示,分别是XY平面轨迹跟踪效果图,XZ平面轨迹跟踪效果图和YZ平面轨迹跟踪效果图,这三个图直观的体现了无人机的实际飞行过程,可以看出:在爬升-拉平阶段,纵向高度上一开始存在约100m左右的超调,但是很快无人机就恢复到指令轨迹上了;在45°大角度水平转弯时,高度变化仅在10m-20m左右,横侧向位置超调达到约400m后,能较快收敛到指令位置,相对于无人机200m/s的高速,该位置超调是较为可观的。
Claims (4)
1.一种基于逆动力学的无人机三维轨迹制导方法,其特征在于,包括如下步骤:
步骤一、针对某个无人机,根据期望轨迹将无人机飞行经过的航路点设定坐标,并划分轨迹段;
期望轨迹由若干航路点联结成的航路段构成;航路段共有n(n≥1)段,对应有n+1个航路点;
步骤二、针对每个航路点,根据无人机的机动性能和所设定好的航路情况,设定该航路点对应的转换半径为d;
步骤三、针对第k段轨迹段,利用无人机的期望轨迹分别计算该轨迹段的实际速度误差和实际位移误差,采用PID控制律加和构成该轨迹段期望的横侧向位移状态的微分项和纵向位移状态的微分项
初始k=1,1≤k≤n;
步骤四、将期望的微分项和微分项分别带入横侧向位移状态和纵向位移状态的微分方程中,采用解析法和数值迭代法相结合的方法对微分方程进行第一级逆动力学解算,输出指令航迹偏角χc和航迹倾角γc;
步骤五、将指令航迹偏角χc和指令航迹倾角γc以及无人机在该航路段上的指令地速Vgk作为下一级逆动力学解算的输入指令值,采用解析法和数值迭代法相结合的方法对无人机地速、航迹倾角和航迹偏角的状态方程进行第二级逆动力学解算,输出指令推力Tc、指令迎角αc和指令航迹滚转角μc;
步骤六、将输出的指令推力Tc、指令迎角αc、侧滑角0°和航迹滚转角μc作为无人机姿态控制回路的输入,实现对无人机三维轨迹的跟踪控制;
步骤七、当无人机飞行至第k段轨迹段对应的转换半径d的范围内时,则无人机飞到所对应的航路点位置,继续将跟踪下一段轨迹,返回至步骤三,直至追踪至最后一个航路点。
2.如权利要求1所述的一种基于逆动力学的无人机三维轨迹制导方法,其特征在于,所述的步骤三具体为:
首先,无人机的实际速度误差包括横侧向速度误差和纵向速度误差计算如下:
是无人机在第k段轨迹段上设定的指令跟踪横侧向速度;是无人机当前在第k段轨迹段上实际的横侧向速度;是无人机在第k段轨迹段上设定的指令跟踪纵向速度;是无人机当前在第k段轨迹段上实际的纵向速度;Vgk是第k段轨迹段上为无人机设定的指令地速;γk是第k个轨迹段所指向方向的航迹倾角;χk是第k个轨迹段所指向方向的航迹偏角;Vg是无人机当前在第k段轨迹段上实际的地速;γ是无人机当前在第k段轨迹段上实际的航迹倾角;χ是无人机当前在第k段轨迹段上实际的航迹偏角;
然后,无人机的实际位移误差包括横侧向位移误差(yline-y)和纵向位移误差(hline-h);计算如下:
从无人机当前实际位置处作此时对应直线轨迹的空间垂线,得到垂足坐标(xline,yline,hline),该垂线段即为无人机的位置误差总和,分解为横侧向位移误差(yline-y)和纵向位移误差(hline-h);
最后,PID控制律加和具体为:对无人机的位移误差(yline-y)和(hline-h)采用PI控制,将速度误差和作为微分项,并将PI控制与微分项求和,PID的参数即为速度误差和位移误差在飞行制导律中分别所占的权重,得到期望的横侧向位移状态的微分项和纵向位移状态的微分项
3.如权利要求1所述的一种基于逆动力学的无人机三维轨迹制导方法,其特征在于,所述的步骤四具体为:
首先,将期望的微分项带入纵向位移状态的微分方程中,采用解析法解算出指令航迹倾角γc,计算如下:
然后,将期望的微分项和指令航迹倾角γc带入到横侧向位移状态的微分方程中,利用数值迭代法解算出指令航迹偏角χc,计算如下:
4.如权利要求1所述的一种基于逆动力学的无人机三维轨迹制导方法,其特征在于,所述的步骤五具体为:
具体步骤如下:
步骤501、将指令航迹偏角χc,指令航迹倾角γc和指令地速Vgk分别与当前无人机实际的航迹偏角χ,航迹倾角γ和地速Vg作差,并采用PD控制,分别得到期望的航迹偏角微分值期望的航迹倾角微分值和期望的地速微分值
步骤502、利用当前无人机的实际地速Vg和航迹倾角γ,以及期望的航迹偏角微分值和期望的航迹倾角微分值计算指令航迹滚转角μc的解析解;
步骤503、利用指令航迹滚转角μc的解析解,结合期望地速状态方程,代入期望航迹倾角状态方程中,得到关于指令迎角αc的非线性一元方程;
FL为无人机的升力,m和g分别为无人机的质量和重力加速度,FD为无人机的阻力;
步骤504、用数值迭代法对关于指令迎角αc的非线性一元方程进行解算,得到指令迎角αc;
步骤505、将指令迎角αc代入期望地速状态方程中,解算出指令推力Tc;
期望地速状态方程如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611233604.6A CN106774400B (zh) | 2016-12-28 | 2016-12-28 | 一种基于逆动力学的无人机三维轨迹制导方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611233604.6A CN106774400B (zh) | 2016-12-28 | 2016-12-28 | 一种基于逆动力学的无人机三维轨迹制导方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106774400A true CN106774400A (zh) | 2017-05-31 |
CN106774400B CN106774400B (zh) | 2020-02-07 |
Family
ID=58921346
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611233604.6A Active CN106774400B (zh) | 2016-12-28 | 2016-12-28 | 一种基于逆动力学的无人机三维轨迹制导方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106774400B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107992069A (zh) * | 2017-11-29 | 2018-05-04 | 上海无线电设备研究所 | 一种无人机路径跟踪控制的制导律设计方法 |
CN108873929A (zh) * | 2018-05-24 | 2018-11-23 | 北京航空航天大学 | 一种固定翼飞机自主着舰方法及*** |
CN109189079A (zh) * | 2018-11-05 | 2019-01-11 | 南京理工大学 | 基于gps定位的移动机器人导航控制方法 |
CN109828602A (zh) * | 2018-01-19 | 2019-05-31 | 北京航空航天大学 | 一种基于观测补偿技术的航迹回路非线性模型变换方法 |
CN110879611A (zh) * | 2019-11-01 | 2020-03-13 | 中国电子科技集团公司电子科学研究院 | 无人机集群三维曲线路径跟踪方法及装置 |
CN111487966A (zh) * | 2020-04-13 | 2020-08-04 | 哈尔滨工程大学 | 一种基于航路点的水面无人艇自适应路径跟踪控制方法 |
CN111580542A (zh) * | 2019-02-15 | 2020-08-25 | 北京京东尚科信息技术有限公司 | 动态无人机编队控制方法、装置及存储介质 |
WO2020238966A1 (zh) * | 2019-05-29 | 2020-12-03 | 深圳市道通智能航空技术有限公司 | 一种参考航点数据的传输方法、装置及无人机 |
CN112099528A (zh) * | 2020-09-18 | 2020-12-18 | 中国航空无线电电子研究所 | 基于化学反应优化算法的制导方法 |
CN112686921A (zh) * | 2021-01-08 | 2021-04-20 | 西安羚控电子科技有限公司 | 一种基于轨迹特征的多干扰无人机检测跟踪方法 |
CN113419557A (zh) * | 2021-06-17 | 2021-09-21 | 哈尔滨工业大学 | 运动无人机音频合成方法 |
CN114326793A (zh) * | 2021-12-09 | 2022-04-12 | 北京机电工程研究所 | 基于plc的无人飞行器飞行控制装置、***及无人飞行器 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103135546A (zh) * | 2013-01-28 | 2013-06-05 | 北京航空航天大学 | 基于航段距离和相对位置向量点积的无人机飞行航段切换方法 |
CN105159305A (zh) * | 2015-08-03 | 2015-12-16 | 南京航空航天大学 | 一种基于滑模变结构的四旋翼飞行控制方法 |
CN105259760A (zh) * | 2015-10-23 | 2016-01-20 | 天津工业大学 | 一种小型四轴飞行器控制器设计方法 |
-
2016
- 2016-12-28 CN CN201611233604.6A patent/CN106774400B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103135546A (zh) * | 2013-01-28 | 2013-06-05 | 北京航空航天大学 | 基于航段距离和相对位置向量点积的无人机飞行航段切换方法 |
CN105159305A (zh) * | 2015-08-03 | 2015-12-16 | 南京航空航天大学 | 一种基于滑模变结构的四旋翼飞行控制方法 |
CN105259760A (zh) * | 2015-10-23 | 2016-01-20 | 天津工业大学 | 一种小型四轴飞行器控制器设计方法 |
Non-Patent Citations (4)
Title |
---|
HAI CHEN等: "《Attitude control of UAV based on PI dynamic inversion》", 《2016 35TH CHINESE CONTROL CONFERENCE》 * |
QING LIN等: "《Adaptive Flight Control Design for Quadrotor UAV Based on Dynamic Inversion and Neural Networks》", 《2013 THIRD INTERNATIONAL CONFERENCE ON INSTRUMENTATION, MEASUREMENT, COMPUTER, COMMUNICATION AND CONTROL》 * |
卢晓东等: "《基于改进动态逆的四旋翼无人机三维路径跟踪控制》", 《仪表技术》 * |
王永寿等: "《采用逆动力学法的无人机跟踪飞行速度控制***设计》", 《控制与制导》 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107992069A (zh) * | 2017-11-29 | 2018-05-04 | 上海无线电设备研究所 | 一种无人机路径跟踪控制的制导律设计方法 |
CN107992069B (zh) * | 2017-11-29 | 2021-01-19 | 上海无线电设备研究所 | 一种无人机路径跟踪控制的制导律设计方法 |
CN109828602A (zh) * | 2018-01-19 | 2019-05-31 | 北京航空航天大学 | 一种基于观测补偿技术的航迹回路非线性模型变换方法 |
CN108873929A (zh) * | 2018-05-24 | 2018-11-23 | 北京航空航天大学 | 一种固定翼飞机自主着舰方法及*** |
CN109189079B (zh) * | 2018-11-05 | 2021-07-23 | 南京理工大学 | 基于gps定位的移动机器人导航控制方法 |
CN109189079A (zh) * | 2018-11-05 | 2019-01-11 | 南京理工大学 | 基于gps定位的移动机器人导航控制方法 |
CN111580542A (zh) * | 2019-02-15 | 2020-08-25 | 北京京东尚科信息技术有限公司 | 动态无人机编队控制方法、装置及存储介质 |
WO2020238966A1 (zh) * | 2019-05-29 | 2020-12-03 | 深圳市道通智能航空技术有限公司 | 一种参考航点数据的传输方法、装置及无人机 |
CN110879611A (zh) * | 2019-11-01 | 2020-03-13 | 中国电子科技集团公司电子科学研究院 | 无人机集群三维曲线路径跟踪方法及装置 |
CN111487966A (zh) * | 2020-04-13 | 2020-08-04 | 哈尔滨工程大学 | 一种基于航路点的水面无人艇自适应路径跟踪控制方法 |
CN111487966B (zh) * | 2020-04-13 | 2022-09-09 | 哈尔滨工程大学 | 一种基于航路点的水面无人艇自适应路径跟踪控制方法 |
CN112099528A (zh) * | 2020-09-18 | 2020-12-18 | 中国航空无线电电子研究所 | 基于化学反应优化算法的制导方法 |
CN112099528B (zh) * | 2020-09-18 | 2023-02-28 | 中国航空无线电电子研究所 | 基于化学反应优化算法的制导方法 |
CN112686921A (zh) * | 2021-01-08 | 2021-04-20 | 西安羚控电子科技有限公司 | 一种基于轨迹特征的多干扰无人机检测跟踪方法 |
CN112686921B (zh) * | 2021-01-08 | 2023-12-01 | 西安羚控电子科技有限公司 | 一种基于轨迹特征的多干扰无人机检测跟踪方法 |
CN113419557A (zh) * | 2021-06-17 | 2021-09-21 | 哈尔滨工业大学 | 运动无人机音频合成方法 |
CN113419557B (zh) * | 2021-06-17 | 2022-07-19 | 哈尔滨工业大学 | 运动无人机音频合成方法 |
CN114326793A (zh) * | 2021-12-09 | 2022-04-12 | 北京机电工程研究所 | 基于plc的无人飞行器飞行控制装置、***及无人飞行器 |
CN114326793B (zh) * | 2021-12-09 | 2023-09-15 | 北京机电工程研究所 | 基于plc的无人飞行器飞行控制装置、***及无人飞行器 |
Also Published As
Publication number | Publication date |
---|---|
CN106774400B (zh) | 2020-02-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106774400A (zh) | 一种基于逆动力学的无人机三维轨迹制导方法 | |
Sujit et al. | Unmanned aerial vehicle path following: A survey and analysis of algorithms for fixed-wing unmanned aerial vehicless | |
CN104615144B (zh) | 基于目标规划的高超声速飞行器再入轨迹在线优化方法 | |
CN103488814B (zh) | 一种适用于再入飞行器姿态控制的闭环仿真*** | |
Singh et al. | Automatic path planning and control design for autonomous landing of UAVs using dynamic inversion | |
CN108388269A (zh) | 基于四旋翼飞行器的无人机编队飞行控制方法 | |
Bulka et al. | Automatic control for aerobatic maneuvering of agile fixed-wing UAVs | |
CN111399531B (zh) | 高超声速飞行器滑翔段制导与姿态控制一体化设计方法 | |
CN102289207B (zh) | 一种可变飞行模态无人机广义指令生成器及其指令生成方法 | |
CN104199303B (zh) | 一种基于矢量场制导的平流层卫星平面路径跟踪控制方法 | |
CN110471456A (zh) | 高超声速飞行器俯冲段制导、姿控、变形一体化控制方法 | |
CN105629734A (zh) | 一种近空间飞行器的轨迹跟踪控制方法 | |
CN106371312A (zh) | 基于模糊控制器的升力式再入预测‑校正制导方法 | |
CN104019701A (zh) | 一种直接力气动力复合控制方法与前向拦截制导方法 | |
CN104536448B (zh) | 一种基于Backstepping法的无人机姿态***控制方法 | |
CN105425812A (zh) | 一种基于双模型下的无人机自动着舰轨迹控制方法 | |
CN109828602A (zh) | 一种基于观测补偿技术的航迹回路非线性模型变换方法 | |
Manjunath et al. | Application of virtual target based guidance laws to path following of a quadrotor UAV | |
CN105759609A (zh) | 一种基于显式非线性模型预测控制的舰载机自主着舰方法 | |
Brezoescu et al. | Straight-line path following in windy conditions | |
CN113093789B (zh) | 一种基于路径点优选的飞行器禁飞区规避轨迹规划方法 | |
Karimoddini et al. | Multi-layer flight control synthesis and analysis of a small-scale uav helicopter | |
CN108958278A (zh) | 一种空天飞行器巡航段快速抗干扰制导方法 | |
CN105094144A (zh) | 一种无人飞艇自适应抗风路径跟踪控制方法 | |
Xu et al. | Multi-eye guidance method for uavs path following |
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 |