Processing math: 12%
  • ISSN 0258-2724
  • CN 51-1277/U
  • EI Compendex
  • Scopus 收录
  • 全国中文核心期刊
  • 中国科技论文统计源期刊
  • 中国科学引文数据库来源期刊

基于混合式磁浮平台的解耦及控制分析

邹圣楠 刘畅 邓舒同 刘英 陈鹏荣

张先玉, 陈勇, 张余, 杨华. 无人机通信中的资源分配及部署位置联合优化[J]. 西南交通大学学报, 2024, 59(4): 917-924. doi: 10.3969/j.issn.0258-2724.20230400
引用本文: 邹圣楠, 刘畅, 邓舒同, 刘英, 陈鹏荣. 基于混合式磁浮平台的解耦及控制分析[J]. 西南交通大学学报, 2022, 57(3): 540-548. doi: 10.3969/j.issn.0258-2724.20210750
ZHANG Xianyu, CHEN Yong, ZHANG Yu, YANG Hua. Joint Optimization of Resource Allocation and Deployment Location in Unmanned Aerial Vehicle-Assisted Communication[J]. Journal of Southwest Jiaotong University, 2024, 59(4): 917-924. doi: 10.3969/j.issn.0258-2724.20230400
Citation: ZOU Shengnan, LIU Chang, DENG Shutong, LIU Ying, CHEN Pengrong. Decoupling and Control Stability Analysis Based on Hybrid Repulsion Maglev Platform[J]. Journal of Southwest Jiaotong University, 2022, 57(3): 540-548. doi: 10.3969/j.issn.0258-2724.20210750

基于混合式磁浮平台的解耦及控制分析

doi: 10.3969/j.issn.0258-2724.20210750
基金项目: 国家自然科学基金(51807192);福建省中科院科技服务网络计划(2021T3034);泉州市高层次人才创新创业项目(2021C024R)
详细信息
    作者简介:

    邹圣楠(1988—),女,博士,研究方向为电磁场仿真与超导电工学,E-mail:shengnan.zou@fjirsm.ac.cn

    通讯作者:

    邓舒同(1988—),男,博士,研究方向为电磁场计算分析与电磁检测,E-mail:shutong.deng@fjirsm.ac.cn

  • 中图分类号: TH134

Decoupling and Control Stability Analysis Based on Hybrid Repulsion Maglev Platform

  • 摘要:

    为了研究磁悬浮平台系统中存在的多自由度耦合问题,提出一种利用永磁体之间的被动受力来减少竖直方向上主动控制的设计思路,给出一种混合斥力式磁浮平台的结构设计,该磁悬浮结构的定子由永磁体和电磁线圈共同组成,由永磁体提供主要的悬浮力,电磁线圈提供水平方向的驱动力,以此减少负责主动悬浮的线圈数量,减小线圈功耗及产热. 基于磁荷模型推导出磁标量势满足的拉普拉斯方程,利用分离变量法求出磁标量势的解析表达式,并对浮子在整个磁场中的受力进行精确的计算;充分研究探讨了定子与动子永磁体之间被动悬浮力的稳定区域,简化忽略了竖直方向上力的解耦,建立被控对象的数学模型,并研制了以微控制单元为中心的数字集成控制器,通过试验研究了平台的悬浮性能. 研究结果表明:本文所提出的混合斥力式磁浮平台在悬浮高度23 mm水平范围 ± 4 mm内,能够实现稳定的水平运动,并且浮子在垂直方向的位移变化不超过0.2 mm.

     

  • 近几十年来,借助于信息技术的快速迭代和制造工艺的飞速发展,无人机在军事、运输、农业等领域广泛应用[1-3]. 无人机具有价格低廉、部署方便、移动灵活、控制简单等优势,将其融入移动通信系统已成为提升通信服务质量的重要举措,得到了学界和工业部门的广泛关注[4-6]. 无人机既可作为移动用户,利用现有的地面基础通信网络对其提供指令、传输数据等业务,提高无人机通信的可靠性和有效性;同时,无人机也可充当空中基站或中继节点,为特定地区或时段提供不间断、宽覆盖的通信服务,并利用良好的空域传输条件满足高速率、低延迟的通信需求. 无人机通信作为现有通信系统的重要补充,克服了地面基础设施的诸多限制,在提高网络覆盖率、改善网络服务质量、增进系统可靠性等方面具有重要意义. 然而,无人机通信也面临诸多挑战. 一方面,受尺寸、载荷、功率等限制,无人机通信的频谱和发射功率有限,需要有效分配通信资源;另一方面,无人机的灵活性和移动性对节点部署提出了高要求,需要合理设计轨迹或位置,以提升网络性能[7-9].

    无人机通信系统的资源分配和轨迹/位置优化问题备受关注. 文献[10]研究了多无人机作为移动基站的通信系统中用户分配和无人机轨迹联合优化方法;文献[11]针对不同的用户需求,以提升用户公平性为目标,研究无人机通信系统中的带宽、功率和轨迹联合优化问题;文献[12]利用博弈论研究任务驱动下的无人机通信系统频谱分配;文献[13-14]探讨了无人机作为中继节点的网络资源分配和轨迹优化等问题. 正交频分多址接入(orthogonal frequency division multiple access, OFDMA)作为高效的多载波通信技术,已成为多个商用系统的标准多用户接入方式[15-16]. 文献[17]以无人机作为空中基站工作在OFDMA模式下,研究无人机轨迹和网络资源联合优化方法,提升网络公平性; 文献[18]研究无人机辅助蜂窝通信系统中的资源分配和轨迹联合优化问题,提出迭代优化算法实现网络容量的最大化;文献[19]基于智能超表面(intelligent reflecting surface, IRS)辅助的OFDMA无人机通信系统,研究无人机轨迹、网络资源、IRS配置的联合优化问题. 上述文献均是基于理论容量进行优化分析,实际应用中,主流的无人机系统和无线通信系统均为数字通信方式,数据速率和可靠性是系统设计的重要依据. 文献[20]将一定数据速率和误码率作为OFDMA无人机网络的限制条件,以最小化传输功率为目标,研究子载波分配、调制模式选择和轨迹联合优化问题,但尚未考虑功率分配等因素,也未进一步考虑提升网络速率、公平性等指标.

    基于上述问题,本文以最大化用户最小传输速率为目标,在一定的误码率要求下,考虑网络资源分配和无人机位置联合优化问题. 所建问题为混合整数优化问题,直接求解难度较大. 为解决变量间的耦合性,将联合优化问题分解为2个子优化问题,并通过适当变换转化为0-1线性优化和凸优化问题,利用迭代优化方法进行求解. 最后,通过实验仿真和算法对比验证所提算法的有效性和公平性.

    考虑将无人机作为空中基站的场景,系统模型如图1所示. 在通信服务区域内,存在一个空中悬停无人机作为空中基站,同时为M个地面用户提供通信服务. 为描述方便,设定用户索引集合为M={1,2,,m,,M}. 特别地,为提升频谱效率,设定无人机基站工作在OFDMA工作模式,系统总带宽为B,均分为N个子载波,单个子载波带宽Δf=B/N,子信道索引集合为N={1,2,,n,N}. 基站会根据用户的需求及信道质量分配一定数量的子信道,且同一子信道仅能分配到至多1个用户. 因此,不同用户之间不存在信道间干扰. 定义子信道分配变量αmn,若子信道n分配给用户mαmn=1,反之,αmn=0. 另外,为提升系统性能,基站可根据一定的算法自适应调整各子信道的功率分配和调制方式. 不失一般性,本文考虑正交幅度调制(quadrature amplitude modulation,QAM),选定的调制维度集为cn{1,2,4,6},即系统可自适应选择BPSK、QPSK、16QAM和64QAM 4种调制方式. 在实际应用中,可根据具体情况选定合适的维度集或其他调制方式.

    图  1  所考虑的无人机通信系统模型示意
    Figure  1.  System model of considered UAV-assisted communication

    设定无人机悬停高度为H,其平面坐标q=(xU,yU). 类似地,假定地面用户的垂直高度为0,相应地,用户m的平面坐标wm=(xm,ym).

    相比于地面通信系统,无人机通信具有良好的视距到达特性,无人机—地面信道极少会受到障碍物及非视距分量的影响,信道传播质量主要取决于直达链路. 类似于文献[21-22],本文选择直达链路为信道传播模型. 基于自由空间信道传播模型,UAV与用户m之间的信号功率衰减系数为[21-24]

    hm=β0H2+ (1)

    式中: \;{\beta _0} 为单位距离(1 m)处的信号功率衰减量.

    假定各子信道独立选择QAM调制方式,用户 m 的接收信号误码率为 P _m^{\mathrm{\;e}} ,若 {\alpha _{mn}} = 1 ,根据文献[25],用户 m 利用子信道 n 实现可靠接收信息符号所需的信噪比为

    {f_m}\left( {{c_n}} \right) = \left\{ {\begin{array}{*{20}{l}} {\dfrac{1}{2}{{\left[ {{Q^{ - 1}}\left( {P _m^{\mathrm{\;e}}} \right)} \right]}^2},\quad{c_n} = 1}, \\ {\dfrac{1}{3}{{\left[ {{Q^{ - 1}}\left( {\dfrac{{P _m^{\mathrm{\;e}}}}{4}} \right)} \right]}^2}\left( {{2^{{c_n}}} - 1} \right),\quad{c_n} = 2,4,6} , \end{array}} \right. (2)

    式中: {Q^{ - 1}}\left( x \right) 为式(3)所示Q函数的逆函数.

    Q\left( x \right) = \frac{1}{{\sqrt {2{\text{π}} } }}\int_x^\infty {{{\mathrm{e}}^{ - \tfrac{{{\;\;t^2}}}{2}}}{\mathrm{d}}t}. (3)

    进一步地,UAV利用信道 n 服务于用户 m 所需的功率下界可表示为

    p_{mn}^{\lim } = {\alpha _{mn}} {\frac{{{f_m}\left( {{c_n}} \right)\delta _0^2\Delta f}}{{{h_m}}}} , (4)

    式中: \delta _0^2 为用户 m 处的噪声功率谱密度.

    网络公平性是衡量通信系统性能的重要指标,能够确保稳定可靠的通信服务. 本文将提升网络的公平性作为系统设计的主要指标. 为得到优化的系统设计方案,首先需要对问题进行建模:将最大化用户最小信息速率作为优化目标,以有限的发射功率、子信道、调制方式作为限制条件,实现对子信道分配、功率分配、调制方式选择及无人机位置的联合优化. 为表示方便,优化变量可表示为: {{A}} = \{ {\alpha _{mn}}, \forall m \in \mathcal{M},\forall n \in \mathcal{N} \} {{C}} = \left\{ {{c_n},\forall n \in \mathcal{N}} \right\} {{P}} = \left\{ {{p_n},\forall n \in \mathcal{N}} \right\} {q} = \left( {{x_{\mathrm{U}}},{y_{\mathrm{U}}}} \right) pn为第n个子信道功率. 相应地,联合优化问题可建模为

    \begin{split} &\text{P0}:\underset{\left\{A,C,P,q\right\}}{\text{max}}\text{min}\;{R}_{m},\\ & \text{ s}\text{.t}\text{. }\quad {\displaystyle \sum _{n=1}^{N}{\alpha }_{mn}{c}_{n}}={R}_{m},\quad\forall m\in \mathcal{M},\\ &\qquad\quad {\displaystyle \sum _{m=1}^{M}{\alpha }_{mn}}\leqslant 1,\quad\forall n\in \mathcal{N},\\ &\qquad\quad {P} _{m}^{{\mathrm{\;e}}}\leqslant {P} _{\mathrm{max}}^{{\mathrm{\;e}}},\quad\forall m\in \mathcal{M},\\ &\qquad\quad {\displaystyle \sum _{n=1}^{N}{p}_{n}}\leqslant {P}_{{\mathrm{T}}},\\ &\qquad\quad {\alpha }_{mn}\in \left\{0,1\right\},\quad\forall m\in \mathcal{M},\;\forall n\in \mathcal{N},\\ &\qquad\quad {c}_{n}\in \left\{1,2,4,6\right\},\quad\forall n\in \mathcal{N}, \end{split}

    式中: {R_m} 为用户 m 获取的信息速率, P _{\max }^{\mathrm{\;e}} 为误码率上限, {P_{\mathrm{T}}} 为UAV总发射功率限制.

    明显地,所建立的优化问题P0为一个混合整数非线性优化问题,并且 {\alpha _{mn}} {c_n} {c_n} {q} 之间互相耦合. 因此,很难直接求解该问题.

    为解决优化问题P0中的变量耦合问题,采用迭代优化的方法进行求解. 利用梯度下降法将优化问题分解为2个子问题:子信道分配及调制模式识别优化问题、功率分配和无人机位置优化问题. 然后,利用迭代优化求解,得到优化的系统设计方案.

    子信道分配系数 {\alpha _{mn}} 和调制模式选择系数 {c_n} 之间互耦,直接联合优化两类参数难度很大. 由于 {c_n} 取值为有限的4个离散值,因此,可将表达式 {\alpha _{mn}}{f_m}\left( {{c_n}} \right) 变换为式(5).

    {q_m}\left( {b_{mn}} \right) = {\rho _{m1}}b_{mn} + 3{\rho _{m2}}b_{mn}^2 + 15{\rho _{m2}}b_{mn}^3 + 63{\rho _{m2}}b_{mn}^4, (5)
    {\rho _{m1}} = \frac{1}{2}{\left[ {{Q^{ - 1}}\left( {P_m^{\mathrm{e}}} \right)} \right]^2}, (6)
    {\rho _{m2}} = \frac{1}{3}{\left[ {{Q^{ - 1}}\left( {\frac{{P_m^{\mathrm{e}}}}{4}} \right)} \right]^2}, (7)
    0 \leqslant \sum\limits_{s = 1}^4 {b_{mn}^s} \leqslant 1, (8)

    式中:\rho _{m1}\rho _{m2}为系数; b_{mn} \in \left\{ {0,1} \right\} ,为二进制变量,且应满足式(8)所示约束条件.

    通过上述变换,可将优化变量 {\alpha _{mn}} {c_n} b_{mn}^s 替代表示,因此,子信道分配及调制模式选择优化问题可表示为

    \begin{split} &\text{P1}:\;\underset{\left\{A,C\right\}}{\text{max }}\text{min }{R}_{m},\\ & \text{s}\text{.t}\text{.}\quad {\displaystyle \sum _{n=1}^{N}{r}_{mn}}={R}_{m},\quad\forall m\in \mathcal{M},\\ &\qquad {r}_{mn} = {b}_{mn} + 2{b}_{mn}^{2} + 4{b}_{mn}^{3} + 6{b}_{mn}^{4},\;\;\\ &\qquad\quad \forall m\in \mathcal{M},\;\forall n\in \mathcal{N},\\ &\qquad 0\leqslant {\displaystyle \sum _{s=1}^{4}{b}_{mn}^{s}}\leqslant 1,\quad\forall m\in \mathcal{M},\; \forall n\in \mathcal{N},\\ &\qquad 0\leqslant {\displaystyle \sum _{m=1}^{M}{\displaystyle \sum _{s=1}^{4}{b}_{mn}^{s}}}\leqslant 1,\quad\forall n\in \mathcal{N},\\ &\qquad{p}_{mn}\geqslant \dfrac{{q}_{m}\left({b}_{mn}^{s}\right){\delta }_{0}^{2}\Delta f}{{h}_{m}},\quad\forall m\in \mathcal{M},\;\forall n\in \mathcal{N}. \end{split}

    优化问题P1为一个0-1整数线性规划问题,可通过成熟的算法(CPLEX或者MOSEK等)予以求解.

    从信息速率的形式上看,在满足功率、误码率等约束条件下,用户的信息速率由子信道选择系数和调制模式选择系数决定. 因此,子问题1决定了用户的信息速率,在分配方案不变、满足相应约束条件情况下,用户信息速率不会随着功率及位置发生变化. 另一方面,由信道传输模型可知,无人机和用户节点之间的距离直接影响了用户的信噪比,进而通过调整调制模式改变接收的信息速率. 因此,功率分配及无人机位置联合优化问题可将最小化最大传输距离作为优化目标,通过调整功率和无人机位置,尽量提升最低信道衰减值,同时又满足子问题1分配方案所需的功率及误码率约束. 基于上述考虑,功率分配及无人机位置联合优化子问题可建模为

    \text{P2}:\;\underset{\left\{P,q\right\}}{\text{min }}\text{max }{D}_{m},
    \begin{array}{l} \text{s}\text{.t}\text{.}\quad {H}^{2} + {\Vert q-{w}_{m}\Vert }^{2}={D}_{m},\quad\forall m\in \mathcal{M},\\ \qquad {p}_{n}\geqslant \dfrac{{q}_{m}\left({b}_{mn}\right){\delta }_{0}^{2}\Delta f\left({H}^{2} + {\Vert q-{w}_{m}\Vert }^{2}\right)}{{\beta }_{0}},\quad \\ \qquad\quad \forall m\in \mathcal{M},\;\forall n\in \mathcal{N},\\ \qquad {\displaystyle \sum _{n=1}^{N}{p}_{n}}\leqslant {P}_{{\mathrm{T}}},\\ \qquad {\alpha }_{mn}\in \left\{0,1\right\},\quad\forall m\in \mathcal{M},\;\forall n\in \mathcal{N},\\ \qquad {c}_{n}\in \left\{1,2,4,6\right\},\quad\forall n\in \mathcal{N}.\end{array}

    进一步,引入松弛变量 \eta = {\text{max }} {D_m} ,将问题变换为

    {\text{P}}3:\;\mathop {{\text{min }}}\limits_{\left\{ {{{P,q}}} \right\}} \eta,
    \begin{gathered} {\text{s}}{\text{.t}}{\text{.}}\quad {H^2} + {\left\| {{q} - {w_m}} \right\|^2} \leqslant \eta ,\quad\forall m \in M , \\ \qquad {p_n} \geqslant \frac{{{q_m}\left( {b_{mn}} \right)\delta _0^2\Delta f\left( {{H^2} + {{\left\| {{q} - {w_m}} \right\|}^2}} \right)}}{{{\beta _0}}}, \quad\\ \qquad\quad \forall m \in \mathcal{M},\; \forall n \in \mathcal{N} , \\ \qquad \sum\limits_{n = 1}^N {{p_n}} \leqslant {P_{\mathrm{T}}} , \\ \qquad {\alpha _{mn}} \in \left\{ {0,1} \right\},\quad\forall m \in \mathcal{M},\;\forall n \in \mathcal{N} , \\ \qquad {c_n} \in \left\{ {1,2,4,6} \right\},\quad\forall n \in \mathcal{N} . \\ \end{gathered}

    明显可见,问题P3为凸优化问题,因此,可通过标准算法(CVX软件包等)进行求解.

    总体而言,直接求解通信资源分配及位置联合优化问题难度较大,可通过梯度下降法将其转变为2个子优化问题,即3.1节和3.2节所述的问题,进而可通过迭代优化的方式进行求解[26],直至算法收敛. 相关算法流程为:

    1) 输入:无人机通信网络设定的误码率要求P_m^{\mathrm{e}},无人机最大传输功率限制 {P_{\mathrm{T}}} ,调制模式选择类型 {c_n} ;最大迭代次数 N_{{\mathrm{umi}}} 和收敛门限 \varepsilon ;无人机初始位置 {{{q}}_0} = \left( {{x_{\mathrm{U}}},{y_{\mathrm{U}}}} \right) 和初始功率分配方案 {{{P}}_0} = \left\{ {{p_n},\forall n \in \mathcal{N}} \right\} ,子信道分配方案 {{{A}}_0} = \left\{ {{\alpha _{mn}},\forall m \in \mathcal{M},\forall n \in \mathcal{N}} \right\} ,调制模式选择方案{{{C}}_0} = \{ {c_n}, \forall n \in \mathcal{N} \} .

    2) 基于参数 {{{A}}_0} {{{C}}_0} {{{P}}_0} {{{q}}_0} 计算得到用户速率初值 R_{m,0} 及可获取的最小用户速率{\phi _0} = {\text{min}}\; R_{m,0}.

    3) for i = 1:N_{\mathrm{umi}}

    4) 基于第i−1次迭代功率分配方案 {{{P}}_{i - 1}} 、无人机位置 {{{q}}_{i - 1}} ,建立子信道分配及调制模式选择优化模型,利用二进制线性规划优化算法求解得到第i次迭代的子信道分配方案 {{{A}}_i} 和调制模式选择方案 {{{C}}_b} .

    5) 基于 {{{A}}_i} {{{C}}_i} ,建立功率分配及无人机位置联合优化模型,利用凸优化算法求解得到优化变量 {{{P}}_i} {{{q}}_i} .

    6) 基于参数 {{{A}}_i} {{{C}}_i} {{{P}}_i} {{{q}}_i} 计算得到第i次迭代的用户速率 R_{m,i} 及可获取的最小用户速率 {\phi_i} = {\text{min}}\; R_{m,i} .

    7) 计算判断:若 i = N_{{\mathrm{umi}}} 或者 {{\left( {{\phi _i} - {\phi _{i - 1}}} \right)} / {{\phi _{i - 1}}}} \leqslant \varepsilon 迭代停止退出;否则, i = i + 1 ,迭代继续.

    8) End For

    9) 输出:优化的资源分配及无人机位置方案( {{{A}}_i} {{{C}}_i} {{{P}}_i} {{{q}}_i} ).

    为验证所提算法的有效性,本节通过设计不同条件下的仿真,利用结果比对验证所提算法的性能. 不失一般性,设定实验仿真条件为:服务区域为半径500 m的圆形区域,无人机以一定的高度悬停在空中,作为空中基站为服务区域内的 K 个用户提供通信服务,且用户均匀分布在服务区域内;另外,无人机的可用子载波个数、总带宽及最大发射功率均受限. 关键参数设置如表1所示. 仿真时,选定具有代表性的均匀分配算法作为对比对象:无人机的水平位置设置在所有服务用户的中心点,即 {x_{\mathrm{U}}} = {{\displaystyle\sum {{x_m}} } / M} , {y_{\mathrm{U}}} = {{\displaystyle\sum{{y_m}} } / M};子信道平均分配至 M 个用户,即每个用户占用的子信道数目为 {N \mathord{\left/ {\vphantom {N M}} \right. } M} ;各子信道功率一致,均为 {p_n} = {{{P_{\mathrm{T}}}} \mathord{/ {\vphantom {{{P_{\mathrm{T}}}} N}} } N};各子信道均选定能够满足误码率要求的最高阶调制方式. 所提算法的迭代参数可设定为 N_{{\mathrm{umi}}}= 20 次, \varepsilon = 0.01 .

    表  1  部分关键系统仿真参数
    Table  1.  Some key system simulation parameters
    参数 数值
    M /个 20
    H /m 100
    N /个 100
    B /MHz 10
    {\beta _0} /dB −50
    \delta _0^2 /(dBm·Hz−1) −169
    P_{\max }^{\mathrm{e}} {10^{ - 4}}
    {P_{\mathrm{T}}} /dBm 22
    下载: 导出CSV 
    | 显示表格

    网络各节点的平面部署情况示意如图2所示,从图中可见,20个地面用户均匀分布在服务区域内,假设无人机初始位置均选择部署在网络中心区域内. 在UAV部署位置优化中,本文设定的优化目标是:尽量减小最大传输距离,降低信号的最大衰减量. 图2中用户9更靠近服务区域的边界,UAV会朝着更靠近用户9的方向部署,同时受到其他用户的限制,UAV的部署不会偏离中心位置太远. 因此,所提算法能够考虑到网络公平性,能够平衡各用户的通信服务需求.

    图  2  网络节点部署位置示意
    Figure  2.  Deployment location of network nodes

    相应地,图3为所提算法得到的子载波分配方案. 图中,蓝色矩形表示“未分配”,黄色矩形代表“分配”,余图同. 明显可见,相比均匀分配算法,本文所提算法对子信道的分配方案实际调整较少,仅用户2和用户12增加了1个子信道,各分配6个子信道;用户14减少了2个子信道,仅分配3个子信道,其余用户仍为5个子信道. 这是由于在本系统设定中子信道数量相对于用户数量较少,可调整的自由度不足,为满足用户的公平性,需要子信道数目保持一定的稳定性. 用户14与UAV空中基站的距离最近,具有最优的信道质量,因此,利用较少数量的子信道既可达到较高的信息传输速率.

    图  3  无人机各子信道的优化分配情况
    Figure  3.  Optimal sub-channel allocation of UAV

    图4表示均匀分配算法及本文所提优化算法所采用的调制方式选择方案. 从图中可以发现:2种算法均没有选择BPSK调制,仅有用户14选择64QAM调制. 因此,所选仿真系统设定中,大量的用户适用于QPSK、16QAM调制方式. 基站可通过调整各子信道功率,使得较多的子信道适用于高阶调制方式. 2种算法的功率分配方案如图5所示. 由于本文设定功率为连续可调变量,功率域的自由度较大,系统可通过合理分配各信道功率使得较多的子信道适用于更高阶的数字调制方式,进而提升用户速率. 为提升所提算法的适用性,也可通过适当的调整,将功率设定为固定可调档位选择,得到离散的功率调整方案. 另外,从图5可以发现,对于同一用户的不同子信道,系统分配的功率不同. 这是因为本文所提算法中功率分配仅是一类变量,当算法获得优化结果后即停止迭代. 事实上,若将同一用户的子信道功率一致增加为约束条件,即可获得均匀的功率分配方案.

    图  4  调制方式选择
    Figure  4.  Modulation mode selection

    利用得到的通信资源分配及无人机位置优化方案,通过计算得到各用户信息速率. 2种算法下用户获取速率的对比如图6所示. 由于本文所提的优化算法综合考虑了频域、功率域和空间域,变量的自由度较大,能够综合利用各空间的参数改进系统的总体性能. 因此,本文所提算法能够有效提高用户的最小获取速率,改善网络的公平性.

    图  5  功率分配方案
    Figure  5.  Power allocation schemes
    图  6  用户获取速率
    Figure  6.  Achievable user rate

    图7为不同误码率要求及不同发射功率限制下的系统最小获取用户速率性能对比,横坐标表示误码率上限及最大功率限制,即 ( {P_{\max }^{\mathrm{e}},{P_{\mathrm{T}}}}) ,其中,PT的单位为dBW. 由图可见,系统最小获取信息速率会随着设定误码率的增大呈增长趋势,也会随着最大发射功率的增大呈增长趋势;所提优化算法能极大提升系统可获取的最小信息速率. 另一方面,最小获取用户速率呈阶梯状增加,这是由于在数字调制方式中,信息速率仅与调制方式及信噪比有关,而信噪比在一定范围内,所选的调制方式可保持一定的稳定性,因此,速率变化为离散变化,而非连续变化.

    图  7  不同误码率及发射功率下最小获取用户速率
    Figure  7.  Minimum achievable user rates with different BERs and transmitted powers

    无人机的高度也是系统设计中需要考虑的重要因素,实际情况中无人机高度会直接影响系统的覆盖范围和可获取的信息速率. 本文中仅考虑通信直达链路,信号仅与传播距离有关,因此,未考虑非直达链路造成的性能损失. 图8显示了系统获取的最小信息速率与无人机高度的关系曲线. 由图8可见,均匀分配算法下,利用设定的系统参数获取的最小速率未发生变化,这是因为均匀分配会使得最小速率用户获得足够的带宽和功率,确保信道在一定变换范围内保持稳定的信息获取速率. 相应地,本文优化算法中,系统最小获取速率会随着无人机高度的提升呈下降趋势,这与预期结果一致:随着无人机高度的增加,信号传输衰减会随之增大,从而造成信息速率的降低. 同时,最小获取信息速率阶梯状下降,这与图7中的变化趋势相似:信息速率在一定信噪比变换范围内保持稳定.

    图  8  不同无人机高度下的用户最小获取信息速率
    Figure  8.  Minimum achievable user rates under different UAV heights

    事实上,以上仿真实验仅是针对某种用户分布情况下做出的针对性分析对比. 下一步,可针对不同的用户节点分布展开分析比对,或者利用传输信道的瞬时或统计特性进行研究,使本文所提算法更具针对性和普适性.

    本文解决了基于OFDMA模式的无人机通信系统中的通信资源及无人机位置联合优化问题. 首先,以提升系统用户的公平性为目标,将联合优化问题建模为子信道、功率、误码率限制下的最大-最小问题;之后,利用迭代优化的方法将联合优化分解为2个子优化问题,通过适当的变换将子优化问题转变为标准的0-1整数线性规划和凸优化问题分别求解. 仿真结果证实所提算法能有效提升用户最小获取速率,改善系统的服务质量.

  • 图 1  混合斥力式磁浮平台简易结构图

    Figure 1.  Simple structure diagram of hybrid repulsion maglev platform

    图 2  动态平衡下的悬浮

    Figure 2.  Suspension under dynamic equilibrium

    图 3  圆环与圆柱永磁体之间的等价关系

    Figure 3.  Equivalence relation between ring and cylindrical permanent magnet

    图 4  沿轴向磁化的圆柱永磁体

    Figure 4.  Axially magnetized cylindrical permanent magnet

    图 5  不同高度处圆环永磁体的磁感应强度分布

    Figure 5.  Magnetic flux density distribution of ring permanent magnets at different heights

    图 6  悬浮力随磁体之间轴向距离的变化

    Figure 6.  Variation of levitation force with axial distance between magnets

    图 7  被动悬浮力在不同高度处随水平位移的变化

    Figure 7.  Variation of passive suspension force with horizontal displacement at different heights

    图 8  浮子受力分析

    Figure 8.  Force analysis of float

    图 9  浮子垂直方向受力与水平位移的关系

    Figure 9.  Relationship between vertical force and horizontal displacement of float

    图 10  数字控制器结构

    Figure 10.  Digital controller structure diagram

    图 11  主程序流程

    Figure 11.  Main program flow diagram

    图 12  实验平台结构

    Figure 12.  Experimental platform structure diagram

    图 13  位移为1.0 mm的阶跃轨迹跟踪

    Figure 13.  Step trajectory tracking with 1.0 mm final value

    图 14  步长0.5 mm、位移为2.0 mm的轨迹跟踪

    Figure 14.  Trajectory tracking with 0.5 mm step length and 2.0 mm final value

  • [1] 翟婉明,赵春发. 现代轨道交通工程科技前沿与挑战[J]. 西南交通大学学报,2016,51(2): 209-226. doi: 10.3969/j.issn.0258-2724.2016.02.001

    ZHAI Wanming, ZHAO Chunfa. Frontiers and challenges of sciences and technologies in modern railway engineering[J]. Journal of Southwest Jiaotong University, 2016, 51(2): 209-226. doi: 10.3969/j.issn.0258-2724.2016.02.001
    [2] 徐飞,罗世辉,邓自刚. 磁悬浮轨道交通关键技术及全速度域应用研究[J]. 铁道学报,2019,41(3): 40-49. doi: 10.3969/j.issn.1001-8360.2019.03.006

    XU Fei, LUO Shihui, DENG Zigang. Study on key technologies and whole speed range application of maglev rail transport[J]. Journal of the China Railway Society, 2019, 41(3): 40-49. doi: 10.3969/j.issn.1001-8360.2019.03.006
    [3] 邓自刚,张勇,王博,等. 真空管道运输系统发展现状及展望[J]. 西南交通大学学报,2019,54(5): 1063-1072. doi: 10.3969/j.issn.0258-2724.20180204

    DENG Zigang, ZHANG Yong, WANG Bo, et al. Present situation and prospect of evacuated tube transportation system[J]. Journal of Southwest Jiaotong University, 2019, 54(5): 1063-1072. doi: 10.3969/j.issn.0258-2724.20180204
    [4] 潘毅,周盟,郭瑞,等. 基于磁悬浮技术的建筑隔震研究进展[J]. 西南交通大学学报,2019,54(3): 475-482. doi: 10.3969/j.issn.0258-2724.20170413

    PAN Yi, ZHOU Meng, GUO Rui, et al. Research progress on building isolation based on electromagnetic levitation techniques[J]. Journal of Southwest Jiaotong University, 2019, 54(3): 475-482. doi: 10.3969/j.issn.0258-2724.20170413
    [5] 熊嘉阳,邓自刚. 高速磁悬浮轨道交通研究进展[J]. 交通运输工程学报,2021,21(1): 177-198.

    XIONG Jiayang, DENG Zigang. Research progress of high-speed maglev rail transit[J]. Journal of Traffic and Transportation Engineering, 2021, 21(1): 177-198.
    [6] 秦新燕,雷金. 磁悬浮定位平台的研究综述[J]. 机床与液压,2012,40(21): 160-166. doi: 10.3969/j.issn.1001-3881.2012.21.043

    QIN Xinyan, LEI Jin. Research review on magnetic levitation positioning stage[J]. Machine Tool & Hydraulics, 2012, 40(21): 160-166. doi: 10.3969/j.issn.1001-3881.2012.21.043
    [7] KIM W J, TRUMPER D L. High-precision magnetic levitation stage for photolithography[J]. Precision Engineering, 1998, 22(2): 66-77. doi: 10.1016/S0141-6359(98)00009-9
    [8] CHEN M Y, LIN T B, HUNG S K, et al. Design and experiment of a macro—micro planar maglev positioning system[J]. IEEE Transactions on Industrial Electronics, 2012, 59(11): 4128-4139. doi: 10.1109/TIE.2011.2174531
    [9] LAHDO M, STROHLA T, KOVALEV S. Magnetically levitated planar positioning systems based on Lorentz forces[C]//2017 11th International Symposium on Linear Drives for Industry Applications (LDIA). Osaka: IEEE, 2017: 1-6.
    [10] LAHDO M, STRÖHLA T, KOVALEV S. Design and implementation of an new 6-DoF magnetic levitation positioning system[J]. IEEE Transactions on Magnetics, 2019, 55(12): 1-7.
    [11] ZHANG H, KOU B Q, ZHOU Y H. Analysis and design of a novel magnetic levitation gravity compensator with low passive force variation in a large vertical displacement[J]. IEEE Transactions on Industrial Electronics, 2020, 67(6): 4797-4805. doi: 10.1109/TIE.2019.2924858
    [12] SHAKIR H, KIM W J. Nanoscale path planning and motion control with maglev positioners[J]. IEEE/ASME Transactions on Mechatronics, 2006, 11(5): 625-633. doi: 10.1109/TMECH.2006.882995
    [13] 蒋启龙,梁达,阎枫. 数字单周期电流控制在电磁悬浮系统中的应用[J]. 西南交通大学学报,2019,54(1): 1-8,22.

    JIANG Qilong, LIANG Da, YAN Feng. Application of digital one-cycle control for current in electromagnetic suspension system[J]. Journal of Southwest Jiaotong University, 2019, 54(1): 1-8,22.
    [14] 周振雄,杨建东,曲永印,等. 基于自抗扰控制器的磁浮平台水平推力控制[J]. 机械工程学报,2008,44(9): 193-199,204. doi: 10.3321/j.issn:0577-6686.2008.09.032

    ZHOU Zhenxiong, YANG Jiandong, QU Yongyin, et al. Horizontal thrust control of magnetic suspension platform based on active disturbance rejection controller[J]. Chinese Journal of Mechanical Engineering, 2008, 44(9): 193-199,204. doi: 10.3321/j.issn:0577-6686.2008.09.032
    [15] 李黎川,丁玉成,卢秉恒. 超精密磁悬浮工作台及其解耦控制[J]. 机械工程学报,2004,40(9): 84-88,94. doi: 10.3321/j.issn:0577-6686.2004.09.018

    LI Lichuan, DING Yucheng, LU Bingheng. High-precision magnetically suspended table and its decoupling control[J]. Chinese Journal of Mechanical Engineering, 2004, 40(9): 84-88,94. doi: 10.3321/j.issn:0577-6686.2004.09.018
    [16] 张生果,朱煜,尹文生,等. 动圈式磁悬浮工件台线圈阵列实时电流分配法则[J]. 机械工程学报,2011,47(6): 180-185. doi: 10.3901/JME.2011.06.180

    ZHANG Shengguo, ZHU Yu, YIN Wensheng, et al. Coil array real-time commutation law for magnetically levitated stage with moving-coils[J]. Journal of Mechanical Engineering, 2011, 47(6): 180-185. doi: 10.3901/JME.2011.06.180
    [17] 孙立军,张涛,赵兵. 永磁磁轴承数学模型的研究[J]. 机械工程学报,2005,41(4): 69-74. doi: 10.3321/j.issn:0577-6686.2005.04.014

    SUN Lijun, ZHANG Tao, ZHAO Bing. Study of mathematical model of permanent magnet bearings[J]. Chinese Journal of Mechanical Engineering, 2005, 41(4): 69-74. doi: 10.3321/j.issn:0577-6686.2005.04.014
    [18] RAVAUD R, LEMARQUAND G, LEMARQUAND V, et al. Analytical calculation of the magnetic field created by permanent-magnet rings[J]. IEEE Transactions on Magnetics, 2008, 44(8): 1982-1989. doi: 10.1109/TMAG.2008.923096
    [19] 苟晓凡,杨勇,郑晓静. 矩形永磁体磁场分布的解析表达式[J]. 应用数学和力学,2004,25(3): 271-278. doi: 10.3321/j.issn:1000-0887.2004.03.008

    GOU Xiaofan, YANG Yong, ZHENG Xiaojing. Analytic expression of magnetic field distribution of rectangular permanent magnets[J]. Applied Mathematics and Mechanics, 2004, 25(3): 271-278. doi: 10.3321/j.issn:1000-0887.2004.03.008
    [20] TANG W B, XIAO L Y, XIA D, et al. 2-D and 3-D analytical calculation of the magnetic field and levitation force between two halbach permanent magnet arrays[J]. IEEE Transactions on Magnetics, 2021, 57(4): 1-8.
    [21] 王瑞凯,左洪福,吕萌. 环形磁铁空间磁场的解析计算与仿真[J]. 航空计算技术,2011,41(5): 19-23. doi: 10.3969/j.issn.1671-654X.2011.05.005

    WANG Ruikai, ZUO Hongfu, LYU Meng. Analytical calculation and simulation for magnetic field distribution of ring magnet[J]. Aeronautical Computing Technique, 2011, 41(5): 19-23. doi: 10.3969/j.issn.1671-654X.2011.05.005
    [22] 吴崇试. 均匀带电圆盘的静电势问题[J]. 大学物理,2000,19(11): 1-4. doi: 10.3969/j.issn.1000-0712.2000.11.001

    WU Chongshi. Electrostatic potential of uniformly charged disc[J]. College Physics, 2000, 19(11): 1-4. doi: 10.3969/j.issn.1000-0712.2000.11.001
    [23] ROBERTSON W. Modelling and design of magnetic levitation systems for vibration isolation[D]. Australia: The University of Adelaide , 2013.
    [24] VERMA S, KIM W J, GU J. Six-axis nanopositioning device with precision magnetic levitation technology[J]. IEEE/ASME Transactions on Mechatronics, 2004, 9(2): 384-391. doi: 10.1109/TMECH.2004.828648
    [25] ZHANG C, LU Y H, LIU G C, et al. Research on one-dimensional motion control system and method of a magnetic levitation ball[J]. The Review of Scientific Instruments, 2019, 90(11): 115005.1-115005.9.
  • 期刊类型引用(5)

    1. 李铮,郑泽熙,李亮,王萌,邓晶雪. 匈牙利铁路平交道口封锁方案设计与实现. 铁道通信信号. 2024(05): 92-97 . 百度学术
    2. 郑泽熙,李亮,邓晓鹏,邓晶雪,段鹏宇. 匈牙利铁路区间道口与公路路口联动方案研究. 铁道通信信号. 2024(06): 85-90 . 百度学术
    3. 赵大地,王恪铭. 基于状态图转形式化B模型的安全苛求系统开发方法. 计算机工程. 2024(11): 173-186 . 百度学术
    4. 王霞,王恪铭,徐扬,唐伟健. 基于形式化方法的平交道口控制系统安全设计. 西南交通大学学报. 2023(01): 109-116 . 本站查看
    5. 王恪铭,王霞,程鹏,刘宁,张传东. 车站联锁系统行为验证与数据确认的形式化方法. 西南交通大学学报. 2021(03): 587-593+613 . 本站查看

    其他类型引用(11)

  • 加载中
图(14)
计量
  • 文章访问数:  616
  • HTML全文浏览量:  148
  • PDF下载量:  36
  • 被引次数: 16
出版历程
  • 收稿日期:  2021-09-28
  • 修回日期:  2022-03-02
  • 刊出日期:  2022-03-11

目录

/

返回文章
返回