用于辅助通信系统高能效的无人机的三维轨迹优化方法技术方案

本发明专利技术公开了一种用于辅助通信系统高能效的无人机的三维轨迹优化方法，用于规划单个无人机完成一个通信节点的数据收集或传播的三维轨迹，使得可以显著提高无人机辅助通信的效率。本发明专利技术在进行无人机三维轨迹优化时，考虑了无人机的平动以及转动的动力学方程，使得优化出的轨迹更符合无人机的动力学特性；能量消耗模型在建模时，没有对无人机运动做过多限制，因此能耗模型适合无人机二维平面运动、三维空间运动、变加速度运动等运动状态的能量消耗估计；在优化建模时，优化的能量结果直接体现无人机的电能消耗，相比现有仅考虑有效能量的方法，本发明专利技术可以考虑了电能转化为有效能量的转化效率，因此优化结果更能体现最终能量消耗。耗。耗。

【技术实现步骤摘要】
用于辅助通信系统高能效的无人机的三维轨迹优化方法


[0001]本专利技术涉及无人机动力消耗
，具体地说，是涉及用于辅助通信系统高能效的无人机的三维轨迹优化方法。

技术介绍

[0002]在“后5G”(Beyond 5G)或6G时代，建立并实现“空
‑
天
‑
陆
‑
海”全维度网络架构，是通信领域的研究热点。无人机作为空中载体，有望作为空基网络的载体之一。由于无人机受有限机载能量的束缚，所以提高无人机能量利用效率，是提高无人机辅助通信实际应用可能性的关键技术之一。在此背景下，如何准确评估无人机能量消耗、如何在无人机完成通信任务中最大限度节省能量，是评价此类方法优劣性的重要指标。
[0003]目前，在无人机作为动态空中基站的应用场景中，关于无人机辅助通信的能量模型以及三维调度研究尚不充分。
[0004]文献“Energy
‑
Efficient UAV Communication With Trajectory Optimization”和文献“Energy Minimization for Wireless Communication With Rotary
‑
Wing UAV”中Y.Zeng等人基于P＝FV公式分别推导了一种固定翼无人机和旋翼无人机有效功率计算方法，但是其并不能直接用来估计无人机的原始能量(如燃油、电能)消耗；其二，他们的能量消耗模型是基于无人机做平面运动的假设推导的，不便于推广到三维空间运动。r/>[0005]文献“Optimal 3D
‑
Trajectory Design and Resource Allocation for Solar
‑
Powered UAV Communication Systems”Y.Sun等人将旋翼无人机能量消耗分解成水平面运动做功和垂直方向做功之和，用于估计无人机在三维空间的运动耗能，然而功率并不满足叠加原理，因此此方法的理论可行性不足。

技术实现思路

[0006]本专利技术的目的在于提供一种用于辅助通信系统高能效的无人机的三维轨迹优化方法，用于规划单个无人机完成一个通信节点的数据收集或传播的三维轨迹，使得可以显著提高无人机辅助通信的效率。
[0007]为实现上述目的，本专利技术采用的技术方案如下：
[0008]用于辅助通信系统高能效的无人机的三维轨迹优化方法，包括以下步骤：
[0009](1)建立地面坐标系和无人机的机体坐标系，根据已知的十字形旋翼动力学模型，建立无人机的平移运动方程和旋转运动方程；
[0010](2)建立自由空间的LoS信道模型，根据已知的环境噪声、发射功率、以及载波波长，确定每个时刻接收端信道容量；
[0011](3)根据无人机所用电动机性能参数，以及无人机的通信系统额定功率，得到无人机的无刷电机功率消耗及无人机的总体消耗功率；
[0012](4)根据无人机的平移运动方程和旋转运动方程，信道容量方程和功率消耗方程，构建无人机的状态空间方程；
[0013](5)根据状态空间方程，分别以时间最优性和能量最优性为目标，添加无人机安全飞行状态约束、通信任务约束，建立优化模型，获得无人机的三维优化轨迹。
[0014]进一步地，在步骤(1)中，所述地面坐标系为以地面O
e
为原点东北天为坐标轴正方向的O
e
‑
XYZ；所述机体坐标系为以无人机质心O
b
为原点的坐标系O
b
‑
xyz。
[0015]进一步地，所述无人机的平移运动方程和旋转运动方程：
[0016][0017]其中m是无人机机体总质量，g是当地重力加速度，L是四旋翼的悬臂长度；J
x
,J
y
,J
z
分别是绕机体x，y，z轴的转动惯量，J
m
是电机转动惯量；x,y,z是无人机在地面坐标中的位置坐标，表示x,y,z的一阶导数即无人机各坐标方向的速度；表示x,y,z的二阶导数即无人机各坐标方向的加速度；F
i
,i＝1,2,3,4表示各螺旋桨产生的拉力，且F
i
(t)＝C
t
ω
i2
(t)，ω
i
为相应螺旋桨的旋转角速度；φ,θ,ψ是欧拉角，分别表示无人机的滚转角、俯仰角和偏航角；表示φ,θ,ψ的一阶导数，即各欧拉角的角速度；表示φ,θ,ψ的二阶导数，即各欧拉角的角加速度；C
t
是电机拉力系数，C
m
是电机扭矩系数；C
dx
,C
dy
,C
dz
是机体x，y，z轴方向的阻力系数；C
dxm
,C
dym
,C
dzm
是机体x，y，z轴方向的阻尼力矩系数。
[0018]进一步地，所述LoS信道模型每个时刻接收端信道容量：
[0019][0020]其中，W是通信带宽，P是通信发射端发射功率，σ2表示信道噪声功率，β(d(t))表示路径损耗；γ0＝Pβ0/σ2表示参考距离为1米处的信噪比，β0表示参考距离为1米时路径损耗；p
x
,p
y
,p
z
表示信息接收端位置坐标。
[0021]进一步地，所述无人机的无刷电机功率消耗：
[0022]P
m
(ω)＝c4ω4+c3ω3+c2ω2+c1ω+c0[0023]其中
[0024][0025][0026]C
m
表示单桨综合力矩系数，K
T
是转矩常数，I
m0
表示电机的空载电流，表示反电动势常数，K
T
与K
E
之间有关系式：K
T
＝9.55K
E
，U
m0
表示电机的空载电压，R
m0
表示空载等效电阻，K
V
表示电机标称KV值；
[0027]从而得到无人机总体消耗功率：
[0028][0029]其中，P
c
表示通信设备发射端功率。
[0030]进一步地，所述状态空间方程的计算中，
[0031]无人机的状态向量为：
[0032][0033]无人机的控制信号向量为：
[0034]u＝[u1,u2,u3,u4]T
；
[0035]得到状态空间方程为：
[0036][0037]其中
[0038][0039][0040][0041][0042]状态空间方程简记为
[0043][0044]进一步地，所述时间最优性目标约束模型为：
[0045]P1:
[0046]s.t.C0:
[0047]C1:0≤u1本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.用于辅助通信系统高能效的无人机的三维轨迹优化方法，其特征在于，包括以下步骤：(1)建立地面坐标系和无人机的机体坐标系，根据已知的十字形旋翼动力学模型，建立无人机的平移运动方程和旋转运动方程；(2)建立自由空间的LoS信道模型，根据已知的环境噪声、发射功率、以及载波波长，确定每个时刻接收端信道容量；(3)根据无人机所用电动机性能参数，以及无人机的通信系统额定功率，得到无人机的无刷电机功率消耗及无人机的总体消耗功率；(4)根据无人机的平移运动方程和旋转运动方程，信道容量方程和功率消耗方程，构建无人机的状态空间方程；(5)根据状态空间方程，分别以时间最优性和能量最优性为目标，添加无人机安全飞行状态约束、通信任务约束，建立优化模型，获得无人机的三维优化轨迹。2.根据权利要求1所述的用于辅助通信系统高能效的无人机的三维轨迹优化方法，其特征在于，在步骤(1)中，所述地面坐标系为以地面O
e
为原点东北天为坐标轴正方向的O
e
‑
XYZ；所述机体坐标系为以无人机质心O
b
为原点的坐标系O
b
‑
xyz。3.根据权利要求2所述的用于辅助通信系统高能效的无人机的三维轨迹优化方法，其特征在于，所述无人机的平移运动方程和旋转运动方程：其中m是无人机机体总质量，g是当地重力加速度，L是四旋翼的悬臂长度；J
x
,J
y
,J
z
分别是绕机体x，y，z轴的转动惯量，J
m
是电机转动惯量；x,y,z是无人机在地面坐标中的位置坐标，表示x，y，z的一阶导数即无人机各坐标方向的速度；表示x,y,z的二阶导数即无人机各坐标方向的加速度；F
i
,i＝1,2,3,4表示各螺旋桨产生的拉力，且F
i
(t)＝C
t
ω
i2
(t)，ω
i
为相应螺旋桨的旋转角速度；φ,θ,ψ是欧拉角，分别表示无人机的滚转角、俯仰角和偏航角；表示φ,θ,ψ的一阶导数，即各欧拉角的角速度；表示φ,θ,ψ的二阶导数，即各欧拉角的角加速度；C
t
是电机拉力系数，C
m
是电机扭矩系数；C
dx
,C
dy
,C
dz
是机体x，y，z轴方向的阻力系数；C
dxm
,C
dym
,C
dzm
是机体x，y，z轴方向的阻尼力矩系数。4.根据权利要求3所述的用于辅助通信系统高能效的无人机的三维轨迹优化方法，其特征在于，所述LoS信道模型每个时刻接收端信道容量：
其中，W是通信带宽，P是通信发射端发射功率，σ2表示信道噪声功率，β(d(t))表示路径损耗；γ0＝Pβ0/σ2表示参考距离为1米处的信噪比，β0表示参考距离为1米时路径损耗；p
x
,p
y
,p
z
表示信息接收端位置坐标。5.根据权利要求4所述的用于辅助通信系统高能效的无人机的三维轨迹优化方法，其特征在于，所述无人机的无刷电机功率消耗：P
m
(ω)＝c4ω4+c3ω3+c2ω2+c1ω+c0其中其中C
m
表示单桨综合力矩系数，K
T
是转矩常数，I
m0
表示电机的空载电流，表示反电动势常数，K
T
与K

【专利技术属性】
技术研发人员：李彬，李清亮，宁召柯，张凯，
申请(专利权)人：四川大学，
类型：发明
国别省市：