CN113050453B - 一种固定翼无人机通信过程能量优化方法 - Google Patents
一种固定翼无人机通信过程能量优化方法 Download PDFInfo
- Publication number
- CN113050453B CN113050453B CN202110321632.8A CN202110321632A CN113050453B CN 113050453 B CN113050453 B CN 113050453B CN 202110321632 A CN202110321632 A CN 202110321632A CN 113050453 B CN113050453 B CN 113050453B
- Authority
- CN
- China
- Prior art keywords
- unmanned aerial
- aerial vehicle
- equation
- max
- drone
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B17/00—Systems involving the use of models or simulators of said systems
- G05B17/02—Systems involving the use of models or simulators of said systems electric
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种固定翼无人机通信过程能量优化方法,首先通过无人机坐标系的转换推导出无人机在高度固定与高度可变两种情况下的模型方程;然后建立一个无人机通信过程中的信道方程;建立无人机发动机能耗方程;所需要的方程建立完成后,将研究的问题描述为具有连续状态不等式约束和等式约束的非线性最优控制问题;最后将最优控制问题方程进行转换,将无人机无线通信问题转化为非线性最优参数选择的典型数学问题,之后通过梯度的方法进行求解,来最小化无人机的能量损耗。与时间优化相比,优化后的能量轨迹更加稳定,能耗降低了很多。
Description
技术领域
本发明涉及无人机技术领域,尤其涉及一种固定翼无人机通信过程能量优化方法。
背景技术
近来,无人飞行器(UAV)无线通信已被广泛使用,引起了越来越多的关注。使用UAV进行通信的优点是,与使用固定基础结构进行通信相比,它具有显着的性能,例如灵活的部署,强大的视线(LoS)链接,并且在大多数情况下都可以建立链接,从而使其成为可能通过直接交流。它的应用也相当广泛,并且在民用和商业领域出现了许多新应用,包括天气监视,森林火灾检测,货运,紧急搜寻和救援,通讯中继等。因此,无人机无线通信是当今最热门的通信手段之一。如何使无人机在通信过程中的能量消耗最小化成为当前比较关注的一个问题,在以上的背景下,需要进行无人机的轨迹优化,以减少无人机在飞行过程中所损耗的能量。
在“An energy efficient design for UAV communication with mobile edgecomputing”这篇文章中提出了一种带有移动边缘计算(MEC)的无人机通信系统。同时优化无人机的轨迹和任务分配以及CPU的计算速度,以最小化整个系统的能耗。通过应用离散线性状态空间近似技术获得近似优化问题,然后通过线性化将非凸约束转化为凸约束。最后,提出了一种基于凹凸过程的算法来逼近优化问题。在“Energy-Efficient UAVCommunication With Trajectory Optimization”整篇文章中,综合考虑了无人机的通信吞吐量和能耗,得出了固定翼无人机推进能耗的理论模型。在此基础上,定义了无人机通信的能效,并进行了有效的设计,使用离散线性状态空间逼近和顺序凸优化来最大化无人机的能效。
在“基于移动边缘计算的无人机能耗优化研究”这篇文章中,通过优化固定翼无人机的飞行轨迹、速度、加速度和数据卸载率以及CPU的周期频率以使得整个网络的能耗最小化。导出了数据计算的能耗模型,然后结合了数据传输能耗模型和固定翼无人机的飞行消耗模型,最后对整个问题进行了建模。由于该问题存在无穷多个优化变量且部分约束是非凸的,难以直接求解,所以首先对优化问题进行离散化,其次导出目标函数的上界并对耦合的优化变量进行解耦合,接着线性化非凸约束,最后提出了一种基于凹凸过程的算法。
以上文章虽然都通过轨迹优化来最小化无人机的能耗,大都是采用凸优化的方法来进行求解,但是可以发现,在他们的文章中都未曾考虑到无人机系统的动力学和运动学方程,也就是说,他们缺少对无人机自身的姿态进行研究。因此,本发明提出了一种控制参数化的方法,将问题转化为一个具有连续状态不等式约束和终端等式约束的非线性最优控制问题,通过梯度的方法进行求解,通过仿真,可以发现本发明的方法不仅可以有效的减少无人机在飞行过程中所产生的能耗,还可以使时间尽可能地缩小。
发明内容
本发明的目的是要提供一种固定翼无人机通信过程能量优化方法,首先通过无人机的运动学和动力学方程推导了固定高度和可变高度条件下的固定翼无人机模型。然后将该问题描述为具有连续状态不等式约束和等式约束的非线性最优控制问题,并通过控制参数化,约束转录方法和局部平滑技术等一系列方法加以解决,来最小化无人机通信所损耗的能量。最后提出一个只优化时间的情况进行对比,通过仿真证明本发明的可行性。
为达到上述目的,本发明是按照以下技术方案实施的:
本发明包括以下步骤:
S1:通过无人机坐标系的转换推导出无人机在高度固定与高度可变两种情况下的模型方程;
S2:建立一个无人机通信过程中的信道方程;
S3:建立无人机发动机能耗方程;
S4:所需要的方程建立完成后,将研究的问题描述为具有连续状态不等式约束和等式约束的非线性最优控制问题;
S5:将最优控制问题方程进行转换,将无人机无线通信问题转化为非线性最优参数选择的典型数学问题,之后通过梯度的方法进行求解,来最小化无人机的能量损耗。
本发明的有益效果是:
本发明是一种固定翼无人机通信过程能量优化方法,与现有技术相比,本发明考虑到了无人机的自身姿态问题,将无人机的运动学方程与动力学方程融入到问题之中,将问题转化成了一个具有连续状态不等式约束和等式约束的非线性最优控制问题,并提供了一个非线性最优控制问题的解决方法。通过时间变换的方法,在优化能量的同时,尽可能减少了无人机在飞行过程中所消耗的时间。通过仿真仿真结果表明,随着最小通信吞吐量的增加,无人机沿着一条近似“8”形的路径盘旋。并且考虑了一个时间优化的对比,与时间优化相比,优化后的能量轨迹更加稳定,能耗降低了很多。
同时,本发明所使用的方法同样适用于旋翼无人机上,只需要将系统模型,信道模型,能量模型改变为相对应的旋翼无人机的模型即可。
附图说明
图1是本发明的时间变化方法示意图;
图2是本发明的控制参数化离散化示意图;
图3是本发明高度不变情况下的约束转换方法示意图;
图4是本发明高度不变情况下的无人机轨迹图;
图5是本发明高度不变情况下的花费时间对比图;
图6是本发明高度不变情况下的无人机能量损耗对比图;
图7是本发明高度可变情况下的无人机轨迹图;
图8是本发明高度可变情况下的花费时间对比图;
图9是本发明高度可变情况下的无人机能量损耗对比图。
具体实施方式
下面结合附图以及具体实施例对本发明作进一步描述,在此发明的示意性实施例以及说明用来解释本发明,但并不作为对本发明的限定。
首先,通过在坐标系之间进行转换获得了无人机的动力学和运动学方程,得到了无人机水平飞行与高度可变两种情况的模型。建立了无人机的通信通道模型与无人机发动机能耗模型。本发明将解决的问题描述为固定翼无人机从起点飞往终点,在确保与基站之间的最小吞吐量的同时优化无人飞行器的能耗。为此,本发明添加了一些等式约束和不等式约束,并将该问题写为具有连续状态不等式约束和终端等式约束的非线性最优控制问题。为了解决该问题,针对一些难点,本发明提出了以下解决方案:(1)由于本发明的终端时间是自由的,因此本发明使用时间缩放方法,将时间作为系统参数,并从[0,T]开始转换为[0,1];(2)问题中的控制u是连续的,这对本发明来说不方便解决,因此本发明提出了一种控制参数化方法来离散控制u;(3)为了处理连续的状态不等式约束,本发明使用了约束转录方法和局部平滑技术;最后,为了更清楚地观察结果,本发明还提出了在满足无人机最小吞吐量的同时优化时间的问题。通过仿真模拟,将固定高度和可变高度条件下的能量优化与时间优化进行了比较,得出了无人机完成任务的轨迹,证明了该方法可以有效地优化飞行过程中无人机的能耗。具体如下:
设有一个任务,本发明设置其起点和终点,并在途中建立基站,让固定翼无人机满足与基站之间通信的最小吞吐量的同时,最小化无人机飞行过程中产生的能耗。
首先通过无人机坐标系的转换推导出无人机在高度固定与高度可变两种情况下的动力学方程与运动学方程。公式1为无人机在高度可变下的模型,其中γ为爬升角,χ为偏航角,μ为滚转角。T为推力,D为阻力,n代表的是无人机的负载系数。公式中的m代表包括无人机所有载荷的总质量,而g是以m/s^2为单位表示的重力加速度。公式2为推导出的无人机在高度固定下的模型。
之后,本发明建立出一个无人机通信过程中的信道模型,如公式3
其中,设基站处于坐标系原点处,B是信道带宽,ρ0为距离基站上空1m处参考接收到的信噪比。S是平均信号功率,N是平均噪声功率,S/N是信噪比(SNR),此处的信噪比通常用分贝表示,分贝数Fb=10×log10(S/N)。
由于本发明所有研究的是固定翼无人机,对于一架固定翼无人机来说,在飞行过程中所损耗的能量主要是发动机所产生的能耗,通信产生的能耗和其他能耗相对较小,因此本发明可以忽略其他能耗,而仅考虑发动机能耗,所以本发明可以将无人机的能耗模型写为
公式4中,F为无人机发动机所产生的力,用于使无人机向前移动并克服前进产生的阻力,V为无人机飞行速度。
本发明所需要的模型建立完成后,便可以将本发明所研究的问题描述为具有连续状态不等式约束和等式约束的非线性最优控制问题,如下所示
x(0)=x0
x(T)=xt
y(0)=y0
y(T)=yt
z(0)=z0
z(T)=zt
V(0)=v0
χ(0)=χ0
γ(0)=γ0
Vmin≤V≤Vmax
z≥zmin
Fmin≤F≤Fmax
μmin≤μ≤μmax
nmin≤n≤nmax
x(0)=x0
x(T)=xt
y(0)=y0
y(T)=yt
V(0)=v0
χ(0)=χ0
Vmin≤V≤Vmax
Fmin≤F≤Fmax
μmin≤μ≤μmax
P1为高度可变下的问题描述,P2为高度不变下的问题描述,其中x0,y0和z0是无人机的初始位置,xt,yt和zt是无人机的端点,v0,χ0和γ0分别是无人机的初始速度,初始航向角和初始俯仰角。σ是无人机在飞行任务中必须满足的最小吞吐量。考虑到无人机本身并为防止固定翼无人机因低速而失速,本发明设置了最小无人机速度Vmin和最大无人机速度Vmax。并且限制了无人机的最小飞行高度zmin。并且由于发动机的硬件,本发明的推力F也受到限制,Fmin,Fmax分别为最小推力和最大推力。为了防止UAV在固定翼UAV的盘旋或升空飞行中失速,也必须限制UAV的最小滚转角μmin,最大滚转角μmax和最小载荷系数nmin,最大载荷系数nmax。
对于这样一个最优控制问题,可以发现还有许多难点,并不好进行处理,为此本发明要将此优化问题方程进行进一步的转换。
首先,由于问题中断时间并不固定,为此本发明采用了时间变换的方法,将时间由0-T,转换为0-1。
如图1所示,可以得到
这样,便可以将原来的方程写为
其次,发现控制u是连续的,这并不方便进行求解,为此,本发明通过控制参数化的方法,将其进行离散化,如图2所示。
本发明在时间段[0,1]中将u分为P个部分,每个部分都是一个常数,并且u的值从一个常数跳到另一个常数,本发明可以得到分段常量的近似函数
公式7中sk=k/p,k=0,1,……,p,并且
本发明便可以将控制u写为
最后,对于本发明的不等式约束,将约束写为
如图3所示,发现并不平滑,这意味着无法找到公式10的梯度来解决问题,因此,通过约束转换和平滑技术来将本发明的不等式约束进行处理。
本发明需要将约束近似为一个光滑函数
公式11中,ε>0,本发明可以将公式10写为
如图3所示,可以看出公式12现在是一个平滑函数,但是由于它是一个受约束的近似值,因此更精确地表示为ε→0,因此这种近似值可能导致无法满足所有约束条件。为了克服这些困难,对于每个i,本发明引入了一个新的约束
这样,本发明通过一系列的变换,将无人机无线通信问题转化为非线性最优参数选择的典型数学问题,之后便可以通过梯度的方法进行求解,来最小化无人机的能量损耗。
图4为本发明在高度不变情况下最小吞吐量500Mbits时的无人机轨迹;
图5为本发明在高度不变情况下时间优化与能量优化在不同最小吞吐量花费时间对比数据;
图6为本发明在高度不变情况下时间优化与能量优化在不同最小吞吐量能量损耗对比数据;
图7为为本发明在高度可变情况下最小吞吐量500Mbits时的无人机轨迹;
图8为本发明在高度可变情况下时间优化与能量优化在不同最小吞吐量花费时间对比数据;
图9为本发明在高度可变情况下时间优化与能量优化在不同最小吞吐量能量损耗对比数据。
本发明的技术方案不限于上述具体实施例的限制,凡是根据本发明的技术方案做出的技术变形,均落入本发明的保护范围之内。
Claims (4)
1.一种固定翼无人机通信过程能量优化方法,其特征在于,包括以下步骤:
S1:通过无人机坐标系的转换推导出无人机在高度固定与高度可变两种情况下的模型方程;公式(1)为无人机在高度可变下的模型方程,公式(2)为推导出的无人机在高度固定下的模型方程;
其中γ为爬升角,χ为偏航角,μ为滚转角,F为推力,D为阻力,n代表的是无人机的负载系数,m代表包括无人机所有载荷的总质量,而g是以m/s^2为单位表示的重力加速度;
S2:建立一个无人机通信过程中的信道方程;
S3:建立无人机发动机能耗方程;
S4:所需要的方程建立完成后,将研究的问题描述为具有连续状态不等式约束和等式约束的非线性最优控制问题;
所述步骤S4中将研究的问题描述为具有连续状态不等式约束和等式约束的非线性最优控制问题具体为:
x(0)=x0
x(T)=xt
y(0)=y0
y(T)=yt
z(0)=z0
z(T)=zt
V(0)=v0
χ(0)=χ0
γ(0)=γ0
Vmin≤V≤Vmax
z≥zmin
Fmin≤F≤Fmax
μmin≤μ≤μmax
nmin≤n≤nmax
x(0)=x0
x(T)=xt
y(0)=y0
y(T)=yt
V(0)=v0
χ(0)=χ0
Vmin≤V≤Vmax
Fmin≤F≤Fmax
μmin≤μ≤μmax
P1为高度可变下的问题描述,P2为高度不变下的问题描述,其中x0,y0和z0是无人机的初始位置,xt,yt和zt是无人机的端点,v0,χ0和γ0分别是无人机的初始速度,初始航向角和初始俯仰角;σ是无人机在飞行任务中必须满足的最小吞吐量;Vmin是最小无人机速度,Vmax是最大无人机速度,zmin是无人机的最小飞行高度,F是推力,Fmin,Fmax分别为最小推力和最大推力,μmin是最小滚转角,μmax是最大滚转角,nmin是最小载荷系数,nmax是最大载荷系数;B是信道带宽,ρ0为距离基站上空1m处参考接收到的信噪比;V为无人机飞行速度;
S5:将最优控制问题方程进行转换,将无人机无线通信问题转化为非线性最优参数选择的典型数学问题,之后通过梯度的方法进行求解,来最小化无人机的能量损耗。
4.根据权利要求1所述的固定翼无人机通信过程能量优化方法,其特征在于:所述步骤S5的具体步骤为:
首先采用时间变换的方法将时间由0-T,转换为0-1;得到
将原来的方程写为
通过控制参数化的方法,将控制u进行离散化,
在时间段[0,1]中将u分为P个部分,每个部分都是一个常数,并且u的值从一个常数跳到另一个常数,得到分段常量的近似函数
公式(7)中sk=k/p,k=0,1,……,p,并且
将控制u写为:
最后,对于不等式约束,将约束写为
通过约束转换和平滑技术来将不等式约束进行处理;
将约束近似为一个光滑函数
公式11中,ε>0,将公式10写为
公式12中表示为ε→0,对于每个i,引入一个新的约束:
这样,通过一系列的变换,将无人机无线通信问题转化为非线性最优参数选择的典型数学问题,之后通过梯度的方法进行求解,来最小化无人机的能量损耗。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110321632.8A CN113050453B (zh) | 2021-03-25 | 2021-03-25 | 一种固定翼无人机通信过程能量优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110321632.8A CN113050453B (zh) | 2021-03-25 | 2021-03-25 | 一种固定翼无人机通信过程能量优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113050453A CN113050453A (zh) | 2021-06-29 |
CN113050453B true CN113050453B (zh) | 2022-02-01 |
Family
ID=76515690
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110321632.8A Active CN113050453B (zh) | 2021-03-25 | 2021-03-25 | 一种固定翼无人机通信过程能量优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113050453B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2018016203A (ja) * | 2016-07-28 | 2018-02-01 | 株式会社テクノアクセルネットワークス | 通信システム及び通信パスの設定方法 |
CN109067490A (zh) * | 2018-09-29 | 2018-12-21 | 郑州航空工业管理学院 | 蜂窝网联下多无人机协同移动边缘计算系统资源分配方法 |
CN109286913A (zh) * | 2018-09-29 | 2019-01-29 | 郑州航空工业管理学院 | 基于蜂窝网联的无人机移动边缘计算系统能耗优化方法 |
CN109682380A (zh) * | 2019-01-16 | 2019-04-26 | 北京邮电大学 | 一种通信无人机路径优化方法及设备 |
CN109839955A (zh) * | 2019-02-25 | 2019-06-04 | 南京航空航天大学 | 一种无人机与多个地面终端进行无线通信的轨迹优化方法 |
CN109857143A (zh) * | 2019-04-09 | 2019-06-07 | 大连大学 | 一种吞吐量最大化的无人机轨迹规划方法 |
CN110364031A (zh) * | 2019-07-11 | 2019-10-22 | 北京交通大学 | 地面传感器网络中无人机集群的路径规划和无线通信方法 |
CN110856191A (zh) * | 2019-10-24 | 2020-02-28 | 广东工业大学 | 一种基于无线通信的无人机轨迹优化方法 |
-
2021
- 2021-03-25 CN CN202110321632.8A patent/CN113050453B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2018016203A (ja) * | 2016-07-28 | 2018-02-01 | 株式会社テクノアクセルネットワークス | 通信システム及び通信パスの設定方法 |
CN109067490A (zh) * | 2018-09-29 | 2018-12-21 | 郑州航空工业管理学院 | 蜂窝网联下多无人机协同移动边缘计算系统资源分配方法 |
CN109286913A (zh) * | 2018-09-29 | 2019-01-29 | 郑州航空工业管理学院 | 基于蜂窝网联的无人机移动边缘计算系统能耗优化方法 |
CN109682380A (zh) * | 2019-01-16 | 2019-04-26 | 北京邮电大学 | 一种通信无人机路径优化方法及设备 |
CN109839955A (zh) * | 2019-02-25 | 2019-06-04 | 南京航空航天大学 | 一种无人机与多个地面终端进行无线通信的轨迹优化方法 |
CN109857143A (zh) * | 2019-04-09 | 2019-06-07 | 大连大学 | 一种吞吐量最大化的无人机轨迹规划方法 |
CN110364031A (zh) * | 2019-07-11 | 2019-10-22 | 北京交通大学 | 地面传感器网络中无人机集群的路径规划和无线通信方法 |
CN110856191A (zh) * | 2019-10-24 | 2020-02-28 | 广东工业大学 | 一种基于无线通信的无人机轨迹优化方法 |
Non-Patent Citations (1)
Title |
---|
基于hp自适应RPM的无人机群编队重组最优控制;刘高旗等;《2020中国自动化大会(CAC2020)论文集》;20201106;第537页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113050453A (zh) | 2021-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Muthusamy et al. | Real-time adaptive intelligent control system for quadcopter unmanned aerial vehicles with payload uncertainties | |
CN109286913B (zh) | 基于蜂窝网联的无人机移动边缘计算系统能耗优化方法 | |
WO2018232196A1 (en) | Autonomuos aircraft health systems and methods | |
Maqsood et al. | Optimization of transition maneuvers through aerodynamic vectoring | |
CN113254188B (zh) | 调度优化方法和装置、电子设备及存储介质 | |
CN112965523B (zh) | 一种绳系连接多飞行器离线轨迹生成方法 | |
CN111722533A (zh) | 针对含有传感器故障的多时滞系统滑模预测容错控制方法 | |
AU2022283385A1 (en) | Roll-biased skid-to-turn terminal guidance with rudder integrator feedback | |
WO2022242468A1 (zh) | 任务卸载方法、调度优化方法和装置、电子设备及存储介质 | |
CN116027809A (zh) | 一种DoS攻击下的多四旋翼无人机编队控制方法 | |
CN113050453B (zh) | 一种固定翼无人机通信过程能量优化方法 | |
Tysse et al. | Lyapunov-based damping controller with nonlinear MPC control of payload position for a knuckle boom crane | |
CN113900449B (zh) | 多无人机轨迹规划方法、装置、无人机及存储介质 | |
CN111953397A (zh) | 一种面向自由信道的全双工无人机中继通信方法 | |
CN114564044B (zh) | 一种输入限幅事件触发的无人机有限时间编队控制方法 | |
CN117042003A (zh) | 混合式ris辅助uav移动边缘计算系统的能量效率优化方法 | |
CN115968010A (zh) | 基于稳定性分析模型的多智能体系统的网络拓扑优化方法 | |
CN114003052B (zh) | 一种基于动态补偿系统的固定翼无人机纵向运动鲁棒自适应控制方法 | |
CN114943168A (zh) | 一种水上浮桥组合方法及系统 | |
Yeh et al. | Fuzzy sliding-mode control for a Mini-UAV | |
Jia et al. | Control of an airship using particle swarm optimization and neural network | |
CN113378451A (zh) | 飞行器着陆点预测方法、设备及存储器 | |
CN112491458A (zh) | 一种变加速无人机移动中继的通信方法 | |
Wang et al. | Design of Adaptive Time-Varying Sliding Mode Controller for Underactuated Overhead Crane Optimized via Improved Honey Badger Algorithm | |
Dong et al. | Quantitative Study of Load Stability of Quadrotor Based on Lyapunov Exponents |
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 |