摘要
无人机通信面临路径损耗和组间干扰等问题.为了满足离散用户通信需求,实现无人机组网的静态部署并最大化能效,本文针对多比凹凸函数分式规划问题进行研究,提出一种凸优化协同群体智能优化策略,该策略将原问题解耦为功率控制与高度优化问题并迭代求解.首先,引入视距概率路径损耗模型,通过俯仰角研究部署高度与水平距离之间关系,将部署问题三维化.其次,利用二次变换解耦原问题,旨在提升视距概率链路下的系统能效.最后,提出快速反馈粒子群算法对高度进行部署,以解决复杂多目标协同优化问题.仿真结果表明,在本文模型下,该策略能够实现算法复杂度与准确性之间平衡,对无人机基站进行高效准确部署.
无线通信技术的迅速发展引领了移动通信、互联网普及、物联网应用等革命性变革,连接了全球人类的交流与互动,塑造了智能城市、医疗创新等领域,为未来科技与生活持续演进开辟了广阔空
在城市或障碍物密集的区域,UAV难以保证始终处于LoS状态,视距概率描述了在无线通信中信号传输路径与接收器之间没有被障碍物阻挡的可能性,尤其在UAV通信或毫米波通信等场景下,它是影响信号质量的关键因素.路径损耗在LoS状态下通常远低于非视距链路(non line of sight, NLoS)状态下,因此引入视距概率模型有助于更精确地计算UAV通信中的路径损耗,进而优化能效.为了准确预测在不同部署高度和水平距离条件下,UAV与地面用户之间存在视距概率,文献[
针对UAV部署高度,文献[
针对系统能效提升,文献[
然而,上述文献仅根据覆盖边缘信干噪比(signal to interference plus noise ratio, SINR)来确定UAV最佳部署高度,未考虑离散用户情况下UAV部署高度问题,且缺少UAV部署高度对能效优化影响的研究,限制了模型优化过程中达到较为理想能效.此外,在提升能效过程中多采用凸优化策略去寻求变量最优解,然而凸函数约束条件限制了问题建模能力,复杂多目标问题难以直接表示为凸优化问题.因此,本文在上述研究基础上,针对提高UAV在城市环境下应急通信的网络能效,做出以下贡献:
1) 改进视距概率路径损耗模型,将UAV部署高度引入能效优化模型中,研究部署高度与能效之间的关系,以找到最佳部署策略.
2) 提出一种凸优化协同群体智能算法,对包含高度部署与功率控制的非凸问题进行求解,提升系统能效.
3) 权衡算法效率与准确性,减少计算资源需求和执行时间,提出适用于特定问题的优化策略,确保算法在实际应用中具有可行性和实用性.
1 系统模型
1.1 随机用户生成
考虑在城市环境下,固定生成在个UAV服务区域,UAV水平位置固定在服务区域的中心,令最大服务范围为,每个UAV服务个随机分布在服务范围内的用户,用户总数为.此服务范围内UAV对用户发射功率平均分配,在固定水平位置下调整UAV高度以及总发射功率.将UAV合集与用户合集定义为与,服务用户合集表示为且.定义用户的服务关联向量,表示该用户在与所有UAV通信节点服务关联,存在以下服务对应关系:
(1) |
上式表示用户仅由单个UAV节点提供通信服务,接收信号为1,其他UAV信号为干扰信号0,其中为向量1-范数.
如
(2) |

图1 系统模型图
Fig.1 System model diagram
1.2 视距概率路径损耗模型
基于理想介质中电磁波扩散的自由传播模型,载波波长 为定值, 为载波频率,为真空中光的速率,考虑UAV发射天线与用户接收天线均为全向天线,将其增益置为1,自由传播损耗表示为
(3) |
由于密集障碍物所产生的阴影和散射会减弱信号强度,不存在完全的菲涅耳区域,并且UAV在与用户通信的过程当中无法保证所有用户均处于LoS通信状态.忽略地面路径带来的散射、衍射等影响因素,主要考虑阴影衰落引起的传播损耗,对过度路径损耗建模.阴影衰落受城市中建筑高度、密度等因素影响,因素随机性使得阴影衰落效应也呈现随机性,常被建模为正态分布,即与, 与分别为LoS与NLoS下阴影衰落均值,方差取决于时间、位置或其他因素环境,考虑阴影衰落的强度相对稳定,不会产生显著变化.有关与间仰角的视距概率函数表示
(4) |
式中:与为环境参数;主要与UAV和用户间俯仰角有关.与之间通信过度路径损耗表示为
(5) |
平均路径损耗在通信中受到自由传播损耗和过度路径损耗影响,表示为
(6) |
尽管提升UAV的部署高度可以增加视距概率,但会增大自由路径损耗,从而导致QoS恶化;而在降低UAV的部署高度以增加QoS时,会降低视距概率,进一步导致QoS下降.将式(3)~
(7) |
在UAV静态部署中,电磁波传播路径较为稳定,信道增益可以近似认为不变,在模型下,仅受到平均路径损耗影响,表示为
(8) |
1.3 组间干扰模型
UAV服务多用户时,将频带分为若干正交子载波分配给用户,UAV服务范围内各用户间不存在干扰,用户在通信过程中容易受到其他UAV正常通信带来的组间干扰.因此,服务范围内接收用户的SINR表示为
(9) |
式中:为接收用户处加性高斯白噪声功率.令UAV覆盖范围中所有用户都使用相同频带进行通信,为通信连接状态下向平均分配的发射功率,即,其中表示总发射功率, 为同频UAV干扰功率.
每架UAV服务范围内用户平均分配传输带宽,表示总传输带宽.区域范围内下行链路总传输速率表示为
(10) |
当UAV与地面站或终端设备之间距离很远时,传输质量较差,适当提升传输功率才能保证信号强度和可靠性,超过一定发射功率限制可能会干扰其他合法无线电通信,并可能对UAV的控制和安全造成威胁.设定最大发射功率阈值, 发射功率满足.系统能效表示为
(11) |
为分析UAV高度对QoS的影响,本文根据门限阈值定义最低SINR,即 .考虑UAV悬停功率恒定不变,在此模型下,续航时间仅与变量发射功率有关.对UAV高度与发射功率进行联合优化,以最大化系统能效的问题表示为
(12) |
2 算法设计
由于上述问题属于分式规划(fractional planning, FP)问题,是一类涉及比率项的优化问题,难以利用传统优化方法求解,本文利用二次变换(quadratic transformation, QT)协同快速反馈粒子群优化(rapid feedback particle swarm optimization, RFPSO)算法,解决UAV部署中高度优化和功率控制问题.QT-RFPSO是一种交替迭代优化算法,结合了传统凸优化算法、反馈机制与智能群体算
QT-RFPSO算法中包括2个循环迭代部分,后文分别用与分别表示算法1与算法4,即主迭代与内部算法的迭代次数,且最大迭代次数分别为与,QT-RFPSO算法伪代码如
算法1:QT-RFPSO算法 |
---|
输入:UAV水平位置矩阵、用户位置矩阵; 输出:功率分配向量、UAV位置矩阵; 初始化UAV位置矩阵、迭代次数以及功率分配向量; Repeat ; 执行算法2,得到功率分配向量; 执行算法4,得到UAV位置矩阵; 通过式(7)~ Until 能效收敛或. |
2.1 高度与发射功率初始化
本文考虑在区域内部署架UAV,UAV坐标矩阵,每个UAV为位终端用户提供服务,通过
在服务范围内产生用户坐标矩阵.由
(13) |
为简化计算,取区间内若干离散值,以上方法可得实际值,通过最小二乘法最小化预测值和之间的误差平方和来寻找区间最佳拟合函数,即最小化误差函数
(14) |
可得最小二乘法下最佳拟合函数.将代入拟合函数,得到满足平均路径损耗最小化的个体最优高度矩阵,在后续算法迭代中保持不变.
由于UAV在服务范围内与多用户进行通信,无法保证UAV服务高度均为期望值,引入高度权重矩阵随迭代次数或变化,每次迭代更新后通过权重系数矩阵对UAV部署高度向量进行适应性调整,通过重复扩充得到迭代过程中群体最优高度矩阵,表示为
(15) |
2.2 功率控制
固定UAV水平位置以及高度,并保证服务对应关系不变,对UAV发射功率进行优化,给出功率控制问题模型,将P1变换表示为
(16) |
此时P2为最大化能效的功率控制问题.令,,处理分子部分为凹函数,分母部分为凸函数的凹凸函数分式规划问题(concave-convex FP,CCFP),传统方法是将问题的分子和分母解耦进行联合优化,以满足用户通信速率需求.为解决CCFP问题,常用经典技术如Dinkelbach变换,用一个新的辅助变量进行迭代更新,表示
(17) |
式中:为的迭代值.上述变换应用到多比率FP问题时,不能保证最优解性质.相比之下,QT可以将单一比率问题推广到涉及多个比率情况,单链路CCFP问题的情况下,P3的QT表示
(18) |
式中:为迭代时的更新值,对于固定的,最优存在封闭解,表示为
(19) |
在固定时,通过凸优化可以有效地得到最优.通过
(20) |
式中:用户接收速率向量,由
(21) |
(22) |
式中:为全1矩阵, 为信道增益矩阵,由
(23) |
在交替迭代中,当其他变量被固定时,辅助变量存在最优值表示为
(24) |
当更新完,辅助变量存在最优值表示为
(25) |
同时固定和,P6对于是标准凸二次规划问题.CVX是用于构造和求解DCP(disciplined convex program)的建模系统,本文利用CVX工具,通过标准凸优化技术求解该问
算法2: QT算法 |
---|
输入:功率向量、UAV位置矩阵、用户位置矩阵; 输出:功率分配向量; 根据 根据 利用CVX求解P6解得. |
2.3 高度优化
在UAV确定当前最优发射功率情况下,调整UAV高度实现能效最大化,由于
(26) |
RFPSO是一种群体智能算法,能够有效解决多目标问题,在搜索空间中通过反馈机制引导个体变化,从而实现目标搜索.在算法迭代过程中,个体会根据适应度反馈更新突变量,同时参考pbest与gbest来指导更新
算法3:适应度算法 |
---|
输入:个体最优高度矩阵pbest、群体最优高度矩阵gbest以及高度权重矩阵; 输出:适应度向量; 通过 更新UAV位置矩阵; 通过 计算适应度向量. |
在RFPSO迭代过程中产生趋于个体最佳高度与群体最佳高度差的随机突变向量,为零向量,初始化,更新随机突变向量,表示为
(27) |
式中:为惯性权重,为认知权重,为随机阶矩阵.引入惯性反馈权重向量,其中为反馈逻辑值,对中每个元素进行逻辑判断,用于评价迭代中突变量优劣,进而产生对应元素为0与1的向量.对粒子速度和位置进一步调整可提高算法的有效性,不同UAV根据在随机过程中产生的进行反馈调节,更新权重系数矩阵,表示为
(28) |
令,根据
算法4:RFPSO算法 |
---|
输入:UAV位置矩阵、个体最优高度矩阵与群体最优高度 矩阵; 输出: UAV位置矩阵、群体最优高度矩阵; 初始化迭代次数、、随机突变向量、群适应度向量以及惯性反馈权重向量 ; Repeat ; 执行算法3计算下适应度; 根据 通过 执行算法3计算下适应度; 根据 令; Until适应度收敛或; 通过 更新UAV位置矩阵. |
3 仿真分析
在本节中,对UAV网络视距概率路径损耗模型与QT-FEPSO算法进行仿真.利用MATLAB 2021a进行仿真,实验CPU为i5-12600k,内存为16 GB,本实验采用的算法收敛阈值为1
定义 | 符号/单位 | 参数值 |
---|---|---|
最小部署高度 | 100 | |
最大部署高度 | 200 | |
最大发射功率 | 30 | |
视距概率参数 | 12.08 | |
视距概率参数 | 0.29 | |
SINR阈值 | -10 | |
LoS阴影衰落均值 | 0.2 | |
NLoS阴影衰落均值 | 20 | |
加性高斯白噪声功率 | -105 |
3.1 应用场景与可行性分析
仿真环境设定为城市或障碍物密集区域,部署若干UAV,为其对应服务范围内的用户提供应急通信服务,保证该区域下正常通信覆盖.假定UAV在下行链路多播场景下,采用2.4 GHz频段进行通信.UAV同时向多个地面用户发送相同的数据,采用正交多址方式为地面用户提供通信服务.国内小型通信UAV飞行高度通常限制在200 m左右,例如中国电信研发的系留无人机最大飞行高度为200 m,移动DG-M20无人机应急通信高空基站常部署在高度200 m的空中.城市环境下较低的飞行高度会导致视距概率降低、电磁信号干扰增加、覆盖范围受限、安全风险增加等系列问题,故设定最低飞行高度为100 m.
由

图2 不同水平距离与高度下平均路径损耗
Fig.2 Average path loss at different horizontal distances and heights

图3 平均路径损耗等高线与最佳高度
Fig.3 Average path loss contour and optimal height
3.2 初始化函数拟合分析
在
由

(a) 线性拟合残差

(b) 4阶多项式残差

(c) 5阶多项式残差

(d) 6阶多项式残差
图4 拟合残差
Fig.4 Fitting residual
通过蒙特卡洛实验(Monte Carlo experiment),设定多次总实验次数,在上述范围内生成个均匀分布的随机数,通过
(29) |
如

(a) N=100

(b) N=1 000

(c) N=10 000
图5 蒙特卡洛实验
Fig.5 Monte Carlo experiment
为进一步验证拟合的可靠性,计算拟合曲线均方根误差(root mean square error, RMSE):
(30) |
以及决定系数(R-squared):
(31) |
计算在=500 m以及制约条件下IntProp,即
(32) |
不同拟合曲线拟合参数,结果如
拟合函数 | RMSE | R-squared | IntProp/% |
---|---|---|---|
线性 | 12.525 8 | 0.912 9 | 53.914 0 |
4阶多项式 | 3.923 7 | 0.991 5 | 54.149 7 |
5阶多项式 | 1.396 1 | 0.998 9 | 54.156 6 |
6阶多项式 | 1.064 5 | 0.999 4 | 54.208 0 |
3.3 优化效果分析
在服务范围内随机生成若干用户,并将服务范围的中心设定为UAV的水平位置.假设UAV数量为7,每个UAV服务10名用户,最大服务范围为500 m作为标准场景.
通过调整UAV数量、服务用户数量以及最大服务范围,得到了如

(a) 标准场景
(b) UAV数量调整为3

(c) 服务用户数量调整为5
(d) 最大服务范围调整为800 m
图6 不同场景初始位置
Fig.6 Initial position of different scenario
由

图7 UAV功率与高度优化对比
Fig.7 Comparison of UAV power and height optimization
通过调整UAV数量、服务用户数量以及最大服务范围,考虑到用户生成的随机性,每个场景统计若干次能效优化情况并求平均值,得到了如

图8 不同场景下能效优化对比
Fig.8 Comparison of energy efficiency optimization in different scenarios
3.4 性能与复杂度分析
对标准场景下QT-RFPSO算法进一步分析,由

图9 QT-RFPSO算法内部不同迭代次数能效对比
Fig.9 Comparison of different iteration energy efficiencies within the QT-RFPSO algorithm
由

图10 QT-RFPSO与协同传统算法能效对比
Fig.10 Energy efficiency comparison between QT-RFPSO and collaborative traditional algorithms
对QT-RFPSO算法进行时间与空间复杂度分析.该算法主要包括迭代高度优化和功率控制两个核心步骤, 算法核心是一个至多执行次的循环,为时间复杂度主要决定因素.在主循环内部,算法执行系列计算和优化步骤:
1) CVX工具箱利用内点法涉及矩阵运算,涉及大小为的功率变量,维度为的信道增益矩阵,时间复杂度包括多项式级别的矩阵运算,故时间复杂度表示为.内点法需存储与问题规模相对应的功率控制变量、位置信息与信道增益矩阵,其中信道增益矩阵占用空间最多,其空间复杂度为.
2) 在内部RFPSO算法中,时间复杂度为多项式级别,取决于迭代次数、粒子数与比较维度,在每次迭代中需要计算粒子的适应度,则时间复杂度可以表示为,需要存储位置信息以及信道增益矩阵,信道增益矩阵的空间复杂度为,高度优化算法复杂度对比如
内部算法 | 时间复杂度 | 空间复杂度 |
---|---|---|
RFPSO | ||
PSO | ||
GA | ||
SA |
综上,由于内部迭代次数与用户数量的不确定性,QT-RFPSO算法总体时间复杂度主要取决于同级别两种内部算法时间复杂度之和,表示为.由于主循环内部两部分空间复杂度相同,均取决于信道增益矩阵的存储空间大小,故总体空间复杂度为.通过
4 结 论
本文针对城市环境应急通信服务,改进了视距概率路径损耗模型,并提出了能效优化复杂分式规划问题,将QT算法结合智能群体算法,提出了QT-RFPSO算法对UAV高度与发射功率联合优化.实验证明,视距概率路径损耗模型对城市环境下的信道状态以及路径损耗描述具有一定的实用性.在上述模型下通过已知UAV与用户的水平距离,寻求该模型下使平均路径损耗最小的部署高度,文中多项式拟合函数在RMSE、R-squared与IntProp分别为1.396 1、0.998 9与54.156 6%,经不同数量蒙特卡洛实验计算,原始数IntProp估计值为54.40%,拟合结果与原始数据高度吻合,具有一定可靠性,在后续算法中一定程度简化了求解部署高度的复杂度.在标准场景下,QT-RFPSO算法对UAV发射功率和部署高度均做出适应性调整,实现了系统性能显著优化,能效提升率可达16.98%.系统能效会随服务用户数量、UAV数量增加而提高,QT-RFPSO算法在不同服务用户数量、UAV数量下能效提升率存在一定幅度波动,在UAV数量较小的情况下能效提升率较低. 在UAV服务范围在300 m附近时,该模型算法提升率能达到最佳;随着服务范围提升至1 000 m, 算法提升率降低至2.08%,算法受到一定限制.将算法内部迭代次数设定在20时,可以平衡优化效果和计算复杂度.QT-RFPSO算法虽然复杂度略高于传统算法,但在收敛效果以及收敛速度方面优于大部分传统算法.对比类似算法,QT-RFPSO在迭代开始阶段能对UAV进行更为精确的部署,提升幅度更为明显,实现复杂度与能效性能之间的权衡.
未来的工作考虑通过采用聚类分簇与有效算法对UAV水平位置进行更精准的调整,进一步提升部署的灵活性和能效优化的效果.采用更优的组内功率控制方式,提高算法在更大服务范围的能效提升效率,同时提高系统中用户的公平性.
参考文献
张海君, 陈安琪, 李亚博, 等. 6G移动网络关键技术[J]. 通信学报, 2022, 43(7): 189-202. [百度学术]
ZHANG H J, CHEN A Q, LI Y B, et al. Key technologies of 6G mobile network[J]. Journal on Communications, 2022, 43(7):189-202.(in Chinese) [百度学术]
LIU Z X,QI J X, SHEN Y Y, et al. Maximizing energy efficiency in UAV-assisted NOMA-MEC networks[J]. IEEE Internet of Things Journal,2023,10(24):22208-22222. [百度学术]
LIU X, FENG J H,LI F, et al. Downlink energy efficiency maximization for RSMA-UAV assisted communications[J]. IEEE Wireless Communications Letters, 2024, 13(1): 98-102. [百度学术]
SU Z J, FENG W M,TANG J,et al.Energy-efficiency optimization for D2D communications underlaying UAV-assisted industrial IoT networks with SWIPT[J].IEEE Internet of Things Journal, 2023, 10(3): 1990-2002 [百度学术]
FENG Q X,MCGEEHAN J,TAMEH E K,et al.Path loss models for air-to-ground radio channels in urban environments[C]//2006 IEEE 63rd Vehicular Technology Conference. Melbourne,VIC,Australia. IEEE, 2006: 2901-2905. [百度学术]
AL-HOURANI A, KANDEEPAN S, JAMALIPOUR A. Modeling air-to-ground path loss for low altitude platforms in urban environments[C]//2014 IEEE Global Communications Conference. Austin, TX, USA. IEEE, 2014: 2898-2904. [百度学术]
AL-HOURANI A, KANDEEPAN S, LARDNER S. Optimal LAP altitude for maximum coverage[J]. IEEE Wireless Communi cations Letters, 2014, 3(6): 569-572. [百度学术]
CHEN L,ZHANG W T,KISHK M A,et al.Correlation of line-of-sight probabilities in aerial-terrestrial communications:modeling,analysis,and application[J].IEEE Transactions on Vehicular Technology,2024,73(5):7049-7065. [百度学术]
KHUWAJA A A,ZHENG G,CHEN Y F,et al.Optimum deployment of multiple UAVs for coverage area maximization in the presence of co-channel interference[J].IEEE Access,2019,7:85203-85212. [百度学术]
XU X Z,TAO R T,LI S K,et al.Collaborative UAV deployment and task allocation for environment sensing in multi-UAV networks[C]//2022 IEEE International Conference on Unmanned Systems (ICUS). Guangzhou, China. IEEE, 2022: 738-743. [百度学术]
PLIATSIOS D,GOUDOS S K,LAGKAS T, et al.Drone-base-station for next-generation internet-of-things:a comparison of swarm intelligence approaches[J]. IEEE Open Journal of Antennas and Propagation, 2021, 3: 32-47. [百度学术]
韩宾, 邓冬梅, 江虹. 基于功率控制的认知无线网络能效研 究[J].湖南大学学报(自然科学版), 2019, 46(4): 115-120. [百度学术]
HAN B,DENG D M,JIANG H.Study on energy-efficient of cognitive wireless networks based on power control[J].Journal of Hunan University (Natural Sciences), 2019, 46(4):115-120.(in Chinese) [百度学术]
MASARACCHIA A,NGUYEN L D,DUONG T Q,et al.Energy-efficient and throughput fair resource allocation for TS-NOMA UAV-assisted communications[J]. IEEE Transactions on Communications, 2020, 68(11): 7156-7169. [百度学术]
LIU Z X,QI J X,SHEN Y Y, et al.Maximizing energy efficiency in UAV-assisted NOMA-MEC networks[J].IEEE Internet of Things Journal, 2023, 10(24): 22208-22222. [百度学术]
ARAFAT M Y,MOH S. Localization and clustering based on swarm intelligence in UAV networks for emergency communications[J]. IEEE Internet of Things Journal, 2019, 6(5): 8958-8976. [百度学术]
李国权, 林金朝, 徐勇军, 等.无人机辅助的NOMA网络用户分组与功率分配算法[J].通信学报, 2020, 41(9):21-28. [百度学术]
LI G Q,LIN J Z,XU Y J,et al.User grouping and power allocation algorithm for UAV-aided NOMA network[J].Journal on Communications, 2020, 41(9): 21-28.(in Chinese) [百度学术]
OMONIWA B,GALKIN B,DUSPARIC I. Optimizing energy efficiency in UAV-assisted networks using deep reinforcement learning[J].IEEE Wireless Communications Letters, 2022, 11(8): 1590-1594. [百度学术]
杨青青, 韩卓廷, 彭艺, 等.基于多密度流聚类的UAV-NOMA用户分簇与功率分配算法[J].湖南大学学报(自然科学版), 2024, 51(6): 86-97. [百度学术]
YANG Q Q, HAN Z T, PENG Y, et al.User clustering and power allocation algorithm for UAV-NOMA based on multi-density stream clustering[J]. Journal of Hunan University (Natural Sciences), 2024, 51(6): 86-97.(in Chinese) [百度学术]
LIU L L, WANG A M, SUN G,et al.Maximizing data gathering and energy efficiency in UAV-assisted IoT:a multi-objective optimization approach[J].Computer Networks,2023,235:109986. [百度学术]
ABDEL-RAZEQ S, SHAKHATREH H, ALENEZI A, et al.PSO-based UAV deployment and dynamic power allocation for UAV-enabled uplink NOMA network[J]. Wireless Communications and Mobile Computing, 2021, 2021(1): 2722887. [百度学术]
AL-HOURANI A, KANDEEPAN S, LARDNER S. Optimal LAP altitude for maximum coverage[J]. IEEE Wireless Communi cations Letters, 2014, 3(6): 569-572. [百度学术]
DINKELBACH W. On nonlinear fractional programming[J].Management Science, 1967, 13(7): 492-498. [百度学术]
SHEN K M,YU W.Fractional programming for communication systems:part Ⅰ:power control and beamforming[J]. IEEE Transactions on Signal Processing,2018,66(10):2616-2630. [百度学术]
NEMIROVSKI A S, TODD M J. Interior-point methods for optimization[J]. Acta Numerica, 2008, 17: 191-234. [百度学术]
PANG J W,LI X H,HAN S.PSO with mixed strategy for global optimization[J].Complexity,2023, 2023(1): 7111548. [百度学术]
MA W,FANG Y,FU W,et al.Cooperative localisation of UAV swarm based on adaptive SA-PSO algorithm[J]. The Aeronautical Journal, 2023, 127(1307): 57-75. [百度学术]