3D Parameters Estimation of Helicopter with Constant Speed Using Single Hydrophone
-
摘要: 针对空中匀速飞行运动目标所激发的水声信号特征，该文将传统的2维平面内估计目标飞行高度、速度等参数的问题扩展到3维空间，可以求解飞行时偏航距离，更符合实际情况，解决了空中快速飞行目标状态3维参数估计问题。该文首先以直升机离散线谱为特征声源，建立其在空气-水两层介质中声学多普勒的3维传播模型，考虑了目标的飞行速度、高度和偏离水听器的偏航距离。然后根据多普勒频移曲线及其1阶、2阶导数的不对称性，推导出水下探空应用中飞行器的3维参数估计方法。最后，通过分析单水听器接收的实测信号，验证了文章构建3维空间多普勒频移飞行参数估计模型的合理性及APP-LMS算法相较于短时傅里叶瞬时频率估计算法能够更准确反演直升机的航行参数。Abstract: The three-dimensional parameter estimation algorithm of the helicopter with constant speed flight from the underwater acoustic data with single hydrophone, which extended the traditional two-dimension flight parameters estimation is proposed. Firstly, the helicopter line spectrum is used as the exciting sound source, and its three-dimensional Doppler propagation model in two-layer air-water medium, including altitude, speed and deviation distance of the helicopter, is established. The asymmetry of the Doppler frequency curve and its first- and second-order derivatives is related with the three-dimensional motion parameter of the helicopter, which can be estimated from the received data. Finally, with the measured data, the rationality of the three-dimensional Doppler shift flight model is verified and the result is compared with short-time Fourier instantaneous frequency estimation algorithm, APP-LMS algorithm can more accurately retrieve the flight parameters such as natural frequency, velocity, altitude and yaw distance of the helicopter.
-
1. 引言
未来无线通信网络将面临高密度用户、高速率及低时延业务等网络需求的挑战[1]。传统的通信系统常采用正交多址接入(Orthogonal Multiple Access, OMA)技术，如时分多址接入(Time-Division Multiple Access, TDMA)、频分多址接入(Frequency-Division Multiple Access, FDMA)、码分多址接入(Code-Division Multiple Access, CDMA)、正交频分复用(Orthogonal Frequency-Division Multiple Access, OFDMA)等技术。正交多址接入技术将网络资源在时域、频域或码域划分成正交的资源块，从而避免了用户之间的干扰。但是，正交的资源块也限制了单个网络的接入用户数。为了应对高密度用户的通信网络需求，在有限的频谱资源下，亟需更高效的多址接入技术来提升接入用户容量和频谱利用率。因此，非正交多址接入(Non-Orthogonal Multiple Access, NOMA)已成为5G通信系统中的研究热点之一。
NOMA技术在功率域为不同的用户提供不同的增益，并采用叠加编码(Superposition Coding, SC)实现多路信号在相同时域、频域或码域的资源块上传输[2]。同时，接收端通过串行干扰消除(Successive Interference Cancellation, SIC)技术来区别多路用户信号，并提取本接收端所需的信息[3]。因此，NOMA技术能够成倍地提升通信网络系统的接入用户容量，并显著地提升频谱资源利用率。在前人研究工作中，文献[4]考虑了多用户功率分配、信令负载、SIC错误率和高移动性等网络因素，通过系统级仿真评估了NOMA技术的性能，实验结果表明NOMA技术的优越性。文献[5]分析了两用户非正交多址接入场景下网络容量与功率控制参数的变化关系，并验证了采用NOMA技术的网络总能够为用户提供不差于采用OMA技术时的性能。此外，由于NOMA技术需要在功率域区分不同用户的信号并通过SIC解调用户信号，因此NOMA技术的应用需要配合高效的功率控制算法。在不完美信道状态信息的场景下，文献[6]考虑了信道状态信息误差、用户服务质量需求等因素，以最小化发送功率为目标，研究了联合功率控制、速率控制、用户调度和SIC技术的多维度资源分配算法；文献[7]考虑以最大化网络能效为优化目标，提出了一种次优的联合用户调度与功率控制算法；文献[8]基于统计的信道状态信息提出了联合速率与功率分配的资源管理方法，在满足各用户的吞吐量需求的条件下最小化网络的总传输功率。基于完美信道状态信息的情况下，文献[9]提出了以最大化能量效率为目标的最优功率控制策略；文献[10]通过松弛方法研究了最小化功率消耗的方法；文献[11]提出了以最大化网络吞吐量为目标的联合用户配对与功率分配方法；文献[12]考虑传统电网与可再生能源的融合电能，研究了联合用户接入与功率分配的资源管理方法；文献[13]研究了异构无线网络下以最大化微小区用户的总吞吐量为目标的资源分配方法。但是，上述资源管理方法均需要通过循环迭代的方式来获得功率控制的次优解或最优解，通常具有较高的时间复杂度。因此，本文致力于研究简单可行的资源管理方法，通过推导功率分配的最优解的闭式表达式，使得基站能够在不需要循环迭代的情况下将各个时隙的瞬时信道状态信息与瞬时队列状态信息直接映射为该时隙下的功率控制结果，从而极大地降低了基站的计算开销。
此外，随着视频分享网络、在线游戏、社交网络等时延敏感性业务的飞速发展，数据流的时延性能逐渐地成为了用户体验的重要指标，单纯最优化吞吐量无法提供良好的时延性能。例如，网络中某个节点具有较差的信道条件和很长的数据队列，仅基于信道条件的资源管理方法倾向于给该节点分配少量的网络资源以提升全网的资源利用率，进一步恶化了该节点的时延性能。因此，高效的资源管理方法应联合信道状态信息和队列状态信息，从而均衡节点排队时延并优化全网吞吐量。基于马尔科夫决策过程和贝尔曼方程，文献[14]设计了时延感知的联合信道分配和功率控制算法来最小化D2D终端的平均时延与平均丢弃率之和。基于李雅普诺夫优化方法(Lyapunov optimization)，文献[15]提出了一种在线的功率控制算法来处理能量效率与排队时延之间的均衡关系。文献[16]研究了联合时延均衡和串行干扰消除的多能源融合网络资源管理方法。文献[17]研究了大规模多小区中联合时延均衡和功率控制的分布式资源管理方法。
本文致力于研究在采用NOMA技术的单小区网络中联合时延均衡和功率控制的资源管理算法。由于功率控制在多条数据流之间的耦合性以及时延均衡和功率控制在长时间维度上的相互影响，联合时延均衡和功率控制的资源管理算法通常具有较高的复杂度。本文的主要贡献如下：
(1) 考虑信道衰落的随机性、业务到达的动态性以及控制变量的耦合性，本文建模了联合时延均衡和功率控制的随机优化问题，在保障各用户的队列稳定性与基站的功率约束下，最大化网络的总吞吐量。其次，针对目标函数的时间平均特性和约束条件的参数在不同时隙的随机特性，本文采用李雅普诺夫优化方法将随机优化问题等效地分解为单时隙的时延均衡问题与功率控制问题。
(2) 本文设计了一种低复杂度的最优时延均衡和功率控制的资源管理方法(DEPC-NOMA)。虽然单时隙的功率控制问题仍为非凸的优化问题，但是本文利用该优化问题的KKT(Karush-Kuhn-Tucker)条件的数学特征，通过理论推导获得了最优解在每个时隙的闭式表达式，极大地降低了功率控制算法的复杂度。
(3) 本文通过仿真对比了DEPC-NOMA方法与其它方法的性能。仿真结果表明了本文提出的DEPC-NOMA方法能够有效地均衡各节点的排队时延并且显著地提升网络总吞吐量。
2. 系统模型
本文研究单小区下行链路，一个基站采用NOMA技术在一个信道上同时服务两个用户，用户集合定义为
U={1,2} 。网络系统在时域上采用时隙结构，其中时隙t∈{1,2,···,∞} 。定义时隙t 的信道状态为G(t)={gu(t)}u∈U ，其中，gu(t) 的取值受影响于路径损耗、多径快衰落和阴影慢衰落。在信道的频谱带宽小于信道的相干带宽的情况下，本文考虑信道状态G(t) 在不同时隙上是独立同分布的，且G(t) 在单个时隙内保持恒定而在时隙边界上可能变化。在每个时隙t 对g1(t) 和g2(t) 进行排序，并定义gl(t) 和gs(t) 分别为集合G(t)={g1(t),g2(t)} 中较大值的元素和较小值的元素，其下标l 和s 对应地分别等于g1(t) 和g2(t) 中较大值的下标和较小值的下标。通过在基站侧采用叠加编码(Superposition Coding, SC)，在任意的功率控制分配结果下，基站利用自适应编码调制保证其给用户s 传输的数据包能够被用户s 成功地解调，则具有较优信道条件的用户l 必定能够成功解调基站给用户s 传输的数据包，并在用户l 采用串行干扰消除技术剥离用户s 的数据包带来的干扰[3]。因此，基站在时隙t 给用户l 和用户s 的传输数据量可分别表示为Rl(t)=wτlog2(1+Pl(t)gl(t)σ2) (1) Rs(t)=wτlog2(1+Ps(t)gs(t)Pl(t)gs(t)+σ2) (2) 其中，
w 表示信道带宽，τ 表示时隙长度，Pu(t) 为基站在时隙t 给用户u 发送数据的发射功率，σ2 为白噪声功率。由于设备能力和功率节约的约束，基站的发射功率存在瞬时功率约束与时间平均的功率约率。基站的瞬时发射功率需大于等于零并且小于等于其允许的最大发射功率Pmax ，而基站的时间平均发射功率需小于等于门限值Pmean ，即Pu(t)≥0,∀u∈U,∀t (3) ∑u∈UPu(t)≤Pmax,∀t (4) limsupT→∞1T T−1∑t=0E{∑u∈UPu(t)}≤Pmean (5) 为了均衡数据排队时延，本文采用接纳控制方法来动态调整数据的接纳过程。考虑动态业务到达模型，将时隙
t 内到达基站的用户u 的到达数据量表示为Au(t) ，而实际填入数据队列的接入数据量表示为ru(t) 。由于接入数据量必定小于到达数据量，可得ru(t)∈[0,Au(t)],∀u∈U,∀t (6) 基站为每个用户
u 维护一个缓存接入数据的队列。令Qu(t) 表示发送节点u 在时隙t 的数据队列长度，则节点u 的数据队列更新方程可定义为Qu(t+1)=max{Qu(t)−Ru(t),0}+ru(t),∀u∈U (7) 其中，
ru(t) 和Ru(t) 的取值受约束于业务产生速率和调制编码方案，因此两者分别存在上界值rmaxu 和Rmaxu 。此外，数据队列的稳定性需要保证时间平均上的数据队列长度为有限值，从而保障接入的业务数据能够在有限时间内离开数据队列，即limsupT→∞1T T−1∑t=0E{Qu(t)}<∞,∀u∈U (8) 用户
u 的吞吐量是数据队列Qu 在单位时间内输出的数据量。在有限排队时延的条件下，数据队列将趋于稳定状态，从而在时间平均上数据队列的输入数据量等于输出数据量。因此，由数据队列更新式(7)可得，在队列稳定状态下，网络u 的吞吐量定义为limT→∞1T T−1∑t=0E{min{Ru(t),Qu(t)}}=limT→∞1T T−1∑t=0E{ru(t)}≜ˉru (9) 综上，在保障基站功率限制和数据队列稳定性的约束下，以最大化网络总吞吐量为目标的联合时延均衡和功率控制的随机优化问题可建模为
\begin{align} {{\rm P}1:} & \ {\mathop {\max }\limits_{_{\{ r(t),{{P}}(t)\} }} } \quad{\sum\limits_{u \in {\cal U}} {{{\bar r}_u}} }\\ {} & \ {{\rm{s}}{\rm{.t}}{\rm{.}}} \quad{{{式}}(3) {—} {{式}}(6),{{式}}(8)} \end{align} (10) 其中，
{{r}}(t) = {[{r_u}(t)]_{u \in \cal{U}}} 和{{P}}(t) = {[{P_u}(t)]_{u \in \cal{U}}} 为随机优化问题的控制变量。首先，任意节点的在任意时隙的传输数据量不但取决于自身的信道条件，而且受其他节点的功率控制结果的影响。其次，由于网络变量的随机性以及时间平均上的约束条件式(5)和式(8)，控制变量\{ {{r}}(t),{{P}}(t)\} 在不同时隙间的耦合性给求解优化问题P1带来了挑战。 虽然本文采用李雅普诺夫优化方法[18]将问题P1转换为单时隙的时延均衡与功率控制问题，然而功率控制优化问题的非凸目标函数使得求解该优化问题具有较高的复杂度。3. 单时隙的时延均衡与功率控制
本节通过李雅普诺夫优化方法分析随机稳定性的特点，依据瞬时的网络状态进行每个时隙的控制决策，设计了低复杂度的最优联合时延均衡与功率控制方法，并推导了最优解的闭式表达式。
为了处理时间平均上的功率约束式(5)，本节定义虚拟队列
Z(t) 的队列长度更新方程定义为{Z(t + 1)} { = \max \left\{ {Z(t) + \sum\limits_{u \in \cal{U}} {{P_u}} (t) - {P_{\rm{mean}}},0} \right\}} (11) 依据队列稳定理论(Rate stability theorem)[18]，通过保证虚拟队列
Z(t) 在时间平均上的队列长度是有限的，即可确保虚拟队列Z(t) 在时间平均上流入速率小于等于流出速率，从而满足时间平均上的功率约束式(5)。令
{{S}}(t) = [{Q_u}(t),Z(t)] 表示时隙t 的队列状态信息，依据Lyapunov Drift定理[18]，定义李雅普诺夫函数(Lyapunov function)和单时隙的条件李雅普诺夫偏移(One-slot conditional Lyapunov drift)分别为L({{S}}(t)) = \left[\displaystyle\sum\nolimits_{u \in \cal{U}} {Q_u^2(t) + {Z^2}(t)} \right]/2 和\Delta ({{S}}(t)) = {\rm E}\{ L({{S}}(t \!+\! 1)) \!-\! L({{S}}(t))|{{S}}(t)\} 。随后，依据Lyapunov优化定理[18]，优化问题\rm{P}1 的李雅普诺夫优化函数(Lyapunov drift-plus-penalty function)构建为\Delta ({{S}}(t)) - V\; {\rm E}\left\{ {\sum\limits_{u \in \cal{U}} {{r_u}} (t)|{{S}}(t)} \right\}, (12) 其上界可通过定理1给出。
定理1 假设
{{G}}(t) 在不同时隙上是独立同分布的，则对于任意的联合时延均衡和功率控制的资源管理方法，所有可能的{{S}}(t) 和任意的大于零的控制参数V\; ，优化问题P1的李雅普诺夫优化函数式(12)的上界为\begin{align} {F_{\rm upper}} =& {\rm E}\left\{ \sum\limits_{u \in \cal{U}} {({Q_u}(} t) - V\,){r_u}(t) - \sum\limits_{u \in \cal{U}} {{Q_u}} (t){R_u}(t)\right. \\ {\rm{}}& \left.+ \sum\limits_{u \in \cal{U}} Z (t){P_u}(t)|{{S}}(t) \right\} \!+ B + Z(t){P_{\rm{mean}}}\quad\quad \end{align} (13) 其中，
{B = } {\sum\limits_{u \in \cal{U}} {\frac{{{{(R_u^{\max })}^2} + {{(r_u^{\max })}^2}}}{2}} + \frac{{P_{\max }^2 - P_{\rm{mean}}^2}}{2}} (14) 证明 首先，对数据队列更新式(7)与虚拟队列更新式(11)分别在等式两边进行取平方、移位、累加等操作可得
\begin{align} {\rm{}}& \sum\limits_{u \in \cal{U}} {\frac{{Q_u^2(t + 1) - Q_u^2(t)}}{2}} \le \sum\limits_{u \in \cal{U}} {\frac{{R_u^2(t) + r_u^2(t)}}{2}} \\ {\rm{}}&\quad\quad - \sum\limits_{u \in \cal{U}} {{Q_u}} (t)({R_u}(t) - {r_u}(t)) \end{align} (15) \begin{align} {\rm{}}& \frac{{{Z^2}(t + 1) - {Z^2}(t)}}{2} \le \frac{{{{\left( {\displaystyle\sum\limits_{u \in \cal{U}} {{P_u}} (t)} \right)}^2} - P_{\rm{mean}}^2}}{2} \\ {\rm{}}& \quad\quad\ + Z(t)\left( {\sum\limits_{u \in \cal{U}} {{P_u}} (t) - {P_{\rm{mean}}}} \right) \end{align} (16) 其中，不等式符号源自于不等式
{(\max [Q - R,0] + r)^2} \le{Q^2} + {R^2} + {r^2} + 2Q(R - r) 。进一步，将不等式(15)和不等式(16)相加可得\begin{align} {\rm{}}& L({{S}}(t + 1)) - L({{S}}(t)) \\ {\rm{}}& \le \sum\limits_{u \in {\cal U}} {\frac{{R_u^2(t) + r_u^2(t)}}{2}} + \frac{{{{\left( {\displaystyle\sum\limits_{u \in {\cal U}} {{P_u}} (t)} \right)}^2} - P_{{\rm{mean}}}^2}}{2}\\ {\rm{}}&\quad+ \sum\limits_{u \in {\cal U}} {{Q_u}} (t){r_u}(t) - \left[ \sum\limits_{u \in {\cal U}} {{Q_u}} (t){R_u}(t)\right. \\ {\rm{}} & \quad\left.+ Z(t)\left( {\sum\limits_{u \in {\cal U}} {{P_u}} (t) - {P_{{\rm{mean}}}}} \right) \right] \end{align} (17) 给定
{{S}}(t) ，对不等式(17)两边同时取条件期望 并减去V\; {\rm E}\left\{ {\sum\limits_{u \in \cal{U}} {{r_u}} (t)|{{S}}(t)} \right\} 可得\begin{align} \Delta ({{S}}(t)) &- V\; {\rm E}\left\{ {\sum\limits_{u \in \cal{U}} {{r_u}} (t)|{{S}}(t)} \right\} \\ {\rm{}}& \le B + {\rm E}\left\{ {\sum\limits_{u \in \cal{U}} {({Q_u}(} t) - V){r_u}(t)|{{S}}(t)} \right\} \\ {\rm{}}& - {\rm E}\left\{ \left[ \sum\limits_{u \in \cal{U}} {{Q_u}} (t){R_u}(t)- Z(t)\right.\right. \\ {\rm{}}& \left.\left.\left( {\sum\limits_{u \in \cal{U}} {{P_u}} (t) - {P_{\rm{mean}}}} \right) \right]\Biggr|{{S}}(t) \right\} \end{align} (18) 其中，
{B}= {\displaystyle\sum\limits_{u \in \cal{U}} {\dfrac{{{{(R_u^{\max })}^2} + {{(r_u^{\max })}^2}}}{2}} + \dfrac{{P_{\max }^2 - P_{\rm{mean}}^2}}{2}} 。证毕
依据Min Drift-Plus-Penalty定理[18]，对于给定的
{{S}}(t) ，优化问题P1可以转换为最小化上界函数式(13)的优化问题，并受限于瞬时约束式(3)、式(4)和式(6)。研究上界函数式(13)的结构可得：(1) 给定
{{S}}(t) ，该上界函数的最后2项在时隙t 内是常数；(2) 给定
{{S}}(t) ，最小化该上界函数的期望等同于根据当前的信道状态信息和队列状态信息的情况来最小化该期望内部的函数；(3) 给定
{{S}}(t) ，该上界函数中期望内部的函数可以分解为2个相互独立的函数并独立地进行最小化操作。具体地，最小化上界函数式(13)中期望内部的函数可独立地分解为：(a)最小化该期望中的第1项对应的时延均衡问题，该时延均衡问题依据各用户的数据队列长度对当前到达的数据进行接纳控制，从而有效地均衡数据传输的时延；(b)最小化该期望中的第2项对应的功率控制问题，该功率控制问题依据各用户的数据队列长度与信道条件建模功率控制对吞吐量的增益并维持关于功率消耗的虚拟队列Z(t) 的有限队列长度。3.1 时延均衡
最小化上界函数式(12)的期望中的第1项可建模为如下优化问题P2
\left. \begin{aligned} {{{\rm P}}2:} & {\mathop {\min }\limits_{_{r(t)}} } \quad{\:\sum\limits_{u \in {\cal U}} {({Q_u}(} t) - V){r_u}(t)}\\ {} & {{\rm{s}}{\rm{.t}}{\rm{.}}} \quad {{{式}}(6)} \end{aligned}\right\} (19) 优化问题P2为线性规划问题[19]，其在时隙
t 的最优解{{{r}}^ * }(t) 为{{r}_u^ * (t)} { = {A_u}(t) \cdot {{{1}}_{\{ {Q_u}(t) < V\,\} }},\forall u \in \cal{U}} (20) 最优解
{{{r}}^ * }(t) 表明，调整控制参数V\; 可以限制数据队列的最大队列长度，从而达到均衡排队时延的目的。3.2 功率控制
在每个时隙
t 下，功率控制的决策是从信道状态信息{{G}}(t) 和队列状态信息{{S}}(t) 到发射功率向量{{P}}(t) 的映射。在两用户的场景下，最小化上界函数式(12)的期望中的第2项可建模为如下优化问题P3\begin{align} {\rm{}}& {{{\rm P}}3:} \\ {\rm{}}& {\mathop {\min }\limits_{_{{\hat{{P}}}(t)}} } \ Z(t)({P_l}(t) \!+\! {P_s}(t)) \\ {\rm{}}& \quad\quad- {Q_l}(t)w\tau {{\log }_2}\left( {1 \!+\! \frac{{{P_l}(t){g_l}(t)}}{{{\sigma ^2}}}} \right)\\ {\rm{}}{\rm{}}& \quad\quad { - {Q_s}(t)w\tau {{\log }_2}\left( {1 + \frac{{{P_s}(t){g_s}(t)}}{{{P_l}(t){g_s}(t) + {\sigma ^2}}}} \right)}\\ {\rm{}}& {{\rm{s}}{\rm{.t}}{\rm{.}}} \quad\quad\quad\quad {{{式}}(3),{{式}}(4)} \end{align} (21) 其中，
{\hat{{P}}}(t) = \{ {P_l}(t),{P_s}(t)\} 。需要注意的是，{R_u}(t) 的信干噪比中发射功率的相互耦合关系使得优化问题P3的目标函数为非凸函数。如定理2所示，本节分析了优化问题P3的KKT条件[19]，并利用KKT条件的函数特性推导了优化问题P3的最优解的闭式表达式，极大地减少了功率控制方法的复杂度。定理2：如果优化问题P3取得了最优解
{\hat {{P}}^*}(t) ，则最优解{\hat {{P}}^*}(t) 必定为最优解集合\tilde {{P}}(t) 中的元素\begin{align} \tilde {{P}}(t) =& \left\{ \left(P_l^{(1)}(t),P_s^{(1)}(t)\right),\left(P_l^{(2)}(t),P_s^{(2)}(t)\right),\right.\\ {\rm{}}& \left.\left(P_l^{(3)}(t),P_s^{(3)}(t)\right)\right\} \end{align} (22) 其中，
P_l^{(1)}(t) = \frac{{({Q_s}(t){g_s}(t) - {Q_l}(t){g_l}(t)){\sigma ^2}}}{{({Q_l}(t) - {Q_s}(t)){g_l}(t){g_s}(t)}} \hspace{90pt} (23) \begin{align} P_s^{(1)}(t) =& \min \biggr\{ \frac{{w\tau {Q_s}(t)}}{{Z(t)}} - \frac{{{\sigma ^2}}}{{{g_s}(t)}} - P_l^{(1)}(t), \\ {\rm{}}& {P_{\max }} - P_l^{(1)}(t) \biggr\}\hspace{160pt} \end{align} (24) \begin{align} {\rm{}}&P_l^{(2)}(t) = \max \left\{ {0,\min \left\{ {\frac{{w\tau {Q_l}(t)}}{{Z(t)}} - \frac{{{\sigma ^2}}}{{{g_l}(t)}},{P_{\max }}} \right\}} \right\}\quad\\ {\rm{}}& \hspace{212pt}(25) \end{align} P_s^{(2)}(t) = 0 \hspace{183pt} (26) P_l^{(3)}(t) = 0 \hspace{183pt} (27) \begin{align} {\rm{}}& P_s^{(3)}(t) = \max \left\{ {0,\min \left\{ {\frac{{w\tau {Q_s}(t)}}{{Z(t)}} - \frac{{{\sigma ^2}}}{{{g_s}(t)}},{P_{\max }}} \right\}} \right\} \ \ \\ {\rm{}}& \hspace{211pt}(28) \end{align} 证明 依据KKT条件定义[19]，优化问题P3的最优解满足KKT条件式(29)—式(33)
\begin{align} {\rm{}}& \frac{{{Q_l}(t)w\tau {g_l}(t)}}{{{P_l}(t){g_l}(t) + {\sigma ^2}}} - \frac{{{Q_s}(t)w\tau {g_s}(t)}}{{{P_l}(t){g_s}(t) + {\sigma ^2}}} \\ {\rm{}}& \quad+ \frac{{{Q_s}(t)w\tau {g_s}(t)}}{{({P_l}(t) + {P_s}(t)){g_s}(t) + {\sigma ^2}}}\! -\! Z(t) \!+\! {\lambda _l} \!-\! \mu = 0\quad\quad\quad\; \end{align} (29) \frac{{{Q_s}(t)w\tau {g_s}(t)}}{{({P_l}(t) + {P_s}(t)){g_s}(t) + {\sigma ^2}}} - Z(t) + {\lambda _s} - \mu = 0 \quad\quad\ (30) {\lambda _l}{P_l}(t) = 0 \hspace{180pt} (31) {\lambda _s}{P_s}(t) = 0 \hspace{180pt} (32) \mu ({P_l}(t) + {P_s}(t) - {P_{\max }}) = 0 \hspace{90pt} (33) 其中，
{\lambda _l} 与{\lambda _s} 的取值可分为以下3种情况：(1) 当
{\lambda _l} = {\lambda _s} = 0 时，优化问题P3的最优解满足条件式(34)—式(36)\begin{align} {\rm{}}& \frac{{{Q_l}(t)w\tau {g_l}(t)}}{{{P_l}(t){g_l}(t) + {\sigma ^2}}} - \frac{{{Q_s}(t)w\tau {g_s}(t)}}{{{P_l}(t){g_s}(t) + {\sigma ^2}}} \\ {\rm{}}& + \frac{{{Q_s}(t)w\tau {g_s}(t)}}{{({P_l}(t) + {P_s}(t)){g_s}(t) + {\sigma ^2}}} - Z(t) - \mu = 0 \end{align} (34) \frac{{{Q_s}(t)w\tau {g_s}(t)}}{{({P_l}(t) + {P_s}(t)){g_s}(t) + {\sigma ^2}}} - Z(t) - \mu = 0 \quad\ (35) \mu ({P_l}(t) + {P_s}(t) - {P_{\max }}) = 0 \hspace{65pt} (36) 求解式(34)—式(36)可得
P_l^{(1)}(t) = \frac{{({Q_s}(t){g_s}(t) - {Q_l}(t){g_l}(t)){\sigma ^2}}}{{({Q_l}(t) - {Q_s}(t)){g_l}(t){g_s}(t)}} \hspace{25pt} (37) \begin{align} P_s^{(1)}(t) =& \min \left\{ \frac{{w\tau {Q_s}(t)}}{{Z(t)}} - \frac{{{\sigma ^2}}}{{{g_s}(t)}} - P_l^{(1)}(t),\right.\\ {\rm{}}& {P_{\max }} - P_l^{(1)}(t) \biggr\} \end{align} (38) (2) 当
{\lambda _s} \ne 0 时，则优化问题P3的最优解满足P_s^{(2)}(t) = 0 且\frac{{{Q_l}(t)w\tau {g_l}(t)}}{{{P_l}(t){g_l}(t) + {\sigma ^2}}} - Z(t) + {\lambda _l} - \mu = 0 (39) {\lambda _l}{P_l}(t) = 0 \hspace{105pt} (40) \mu ({P_l}(t) - {P_{\max }}) = 0 \hspace{67pt} (41) 求解式(39)—式(41)可得
\begin{align} {\rm{}}&P_l^{(2)}(t) = \max \left\{ {0,\min \left\{ {\frac{{w\tau {Q_l}(t)}}{{Z(t)}} - \frac{{{\sigma ^2}}}{{{g_l}(t)}},{P_{\max }}} \right\}} \right\}\\ {\rm{}}& \hspace{212pt}(42) \end{align} (3) 当
{\lambda _l} \ne 0 时，则优化问题P3的最优解满足P_l^{(3)}(t) = 0 且\frac{{{Q_s}(t)w\tau {g_s}(t)}}{{{P_s}(t){g_s}(t) + {\sigma ^2}}} - Z(t) + {\lambda _s} - \mu = 0 (43) {\lambda _s}{P_s}(t) = 0 \hspace{107pt} (44) \mu ({P_s}(t) - {P_{\max }}) = 0 \hspace{72pt} (45) 求解式(43)—式(45)可得
\begin{align} {\rm{}}&P_s^{(3)}(t) = \max \left\{ {0,\min \left\{ {\frac{{w\tau {Q_s}(t)}}{{Z(t)}} - \frac{{{\sigma ^2}}}{{{g_s}(t)}},{P_{\max }}} \right\}} \right\}\\ {\rm{}}& \hspace{212pt}(46) \end{align} 综上，优化问题
{\rm P}3 的最优解必定属于最优解集合\tilde {{P}}(t) 的元素。 证毕4. 仿真结果与数值分析
依据Min Drift-Plus-Penalty定理[18]，数据队列式(7)和虚拟队列式(11)是稳定的，从而满足了时间平均上的有限队列长度约束式(8)和功率约束式(5)。本节通过MATLAB仿真验证本文提出的最优时延均衡和功率控制的非正交多址接入方法(DEPC-NOMA)在网络吞吐量和数据队列长度上的性能。本文采用如下2种方法做性能对比：
(1) DEEP-NOMA：该方法的时延均衡方法与DEPC-NOMA方法相同。其次，该方法的功率控制采用均分的方式(Equal Power, EP)，基站在每个时隙
t 给两个用户分配相同的发射功率，并保证满足功率约束。(2) DEPC-TDMA：该方法的时延均衡方法与DEPC-NOMA方法相同。其次，两用户通过时分多址接入(TDMA)的方式共享频率与功率资源，则该方法的功率控制的最优解
{\hat {{P}}^*}(t) 属于集合式(47)\begin{align} {\rm{}}&{\tilde {{P}}_{\rm TDMA}}(t) \!=\! \left\{ \left(P_l^{(2)}(t),P_s^{(2)}(t)),(P_l^{(3)}(t),P_s^{(3)}(t)\right)\right\} \\ {\rm{}}& \hspace{212pt}(47) \end{align} 4.1 仿真设置
考虑一个半径为250 m的圆形小区，基站位于中心位置，用户1和用户2分别随机分布在距离基站125 m和200 m的圆环上。信道带宽
w 为10 MHz。每个时隙长度\tau 为1 ms。基站的瞬时功率约束{P_{\max }} 与时间平均的功率约束{P_{\rm{mean}}} 分别为20 dBm与18 dBm。噪声功率谱密度{\sigma ^2} 为–174 dBm/Hz。每个时隙内的平均数据到达量为50 kbit/slot。路径损耗模型采用PL = 128.1 + 37.6{\log _2}d ，其中，路径距离d 单位为km。多径快衰落服从相互独立的均值为1的指数分布。阴影慢衰落模型采用标准方差为8 dB的独立对数正态分布。在仿真中，数值结果为50次拓扑生成且每个拓扑运行2000个时隙所获得的平均值结果。4.2 仿真结果
图1描述了参数V的取值范围为[1, 400]时平均数据队列长度与网络总吞吐量的权衡关系。在平均数据队列长度较小时，随着平均数据队列长度逐渐增加，用户有充足的缓存数据来保证网络传输机会得到高效的利用，从而网络总吞吐量迅速地上升。而在平均数据队列长度较大时，网络传输机会的利用趋于饱和状态，从而网络总吞吐量随着平均数据队列长度的增加而缓慢上升。此外，对比DEEP-NOMA方法和DEPC-TDMA方法，本文提出的DEPC-NOMA方法能够显著地提升网络吞吐量并降低排队时延。
图2和图3分别说明了各用户的吞吐量和平均数据队列长度与控制参数
V\; 的关系。各用户的吞吐量随着控制参数V\; 的增大而增大，但是各用户的平均数据队列长度及排队时延也随着控制参数V\; 的增大而增大。相对于DEPC-TDMA方法，DEPC-NOMA方法能够创造更多的传输机会，从而有效地提升用户的吞吐量性能，同时保持良好的排队时延性能。此外，在DEEP-NOMA方法中，由于其功率控制未考虑时延均衡的需求，虽然信道条件较好的用户1能够获得很高的吞吐量与很低的排队时延，但是信道条件较差的用户2却忍受着极低的吞吐量与极高的排队时延，用户间的公平性极差。在本文提出的DEPC-NOMA方法中，通过联合时延均衡与功率控制，虽然用户1的吞吐量与排队时延性能稍微下降，但是用户2的吞吐量与排队时延能够获得极大地性能提升。5. 结束语
针对采用非正交多址接入技术的两用户单小区网络，考虑信道衰落的随机性、业务到达的动态性以及控制变量的耦合性，本文提出了一种低复杂度的最优时延均衡和功率控制的资源管理方法(DEPC-NOMA)，以期最大化网络总吞吐量，并满足数据队列稳定性与基站功率约束。仿真结果表明，本文提出的DEPC-NOMA方法能够显著地提升网络总吞吐量并均衡用户的排队时延。
-
表 1 各参数在穿过CPA点前后的变化规律
位置 从负无穷至CPA CPA 从CPA至正无穷 t (s) - \infty \to 0 0 0 \to + \infty \alpha \left( t \right) (s) 0 \to {{\pi} \mathord{\left/ {\vphantom {{\pi} 2}} \right. } 2} {{\pi} \mathord{\left/ {\vphantom {{\pi} 2}} \right. } 2} {{\pi} \mathord{\left/ {\vphantom {{\pi} 2}} \right. } 2} \to {\pi} {\theta _{\text{I}}}\left( t \right) (s) \arcsin {\text{n}} \to \dfrac{{\pi}}{{\text{2}}} - \arctan \dfrac{h}{{{w_{{\text{amin}}}}}} \dfrac{{\pi}}{2} - \arctan \dfrac{h}{{{w_{{\text{amin}}}}}} \dfrac{ {\pi} }{2} - \arctan \dfrac{h}{ { {w_{ {\text{amin} } } } } } \to \arcsin {{n} } {\theta _{\text{T}}}\left( t \right) (s) \dfrac{{\pi}}{2} \to \dfrac{{\pi}}{2} - \arctan \dfrac{d}{{{w_{{\text{wmin}}}}}} \dfrac{{\pi}}{2} - \arctan \dfrac{d}{{{w_{{\text{wmin}}}}}} \dfrac{{\pi}}{2} - \arctan \dfrac{d}{{{w_{{\text{wmin}}}}}} \to \dfrac{{\pi}}{2} 
下载: 导出CSV
-
[1] 修建娟, 张敬艳, 董凯. 基于动力学模型约束的空间目标精确跟踪算法研究[J]. 电子学报, 2021, 49(4): 781–787. doi: 10.12263/DZXB.20200336XIU Jianjuan, ZHANG Jingyan, and DONG Kai. Precise tracking algorithm of space target based on dynamic model[J]. Acta Electronica Sinica, 2021, 49(4): 781–787. doi: 10.12263/DZXB.20200336 [2] 穆森, 李京华, 张恒, 等. 基于谐波集检测的飞行目标水下声探测算法研究[J]. 兵工学报, 2019, 40(5): 1050–1057. doi: 10.3969/j.issn.1000-1093.2019.05.018MU Sen, LI Jinghua, ZHANG Heng, et al. Research on underwater acoustic detection algorithm of aerial targets based on harmonic set detection[J]. Acta Armamentarii, 2019, 40(5): 1050–1057. doi: 10.3969/j.issn.1000-1093.2019.05.018 [3] PENHALE M and BARNARD A. Direction of arrival estimation in practical scenarios using moving standard deviation processing for localization and tracking with acoustic vector sensors[J]. Applied Acoustics, 2020, 168: 107421. doi: 10.1016/j.apacoust.2020.107421 [4] BUCKINGHAM M J, GIDDENS E M, POMPA J B, et al. Sound from a light aircraft for underwater acoustics experiments?[J]. Acta Acustica United with Acustica, 2002, 88(5): 752–755. [5] BUCKINGHAM M J, GIDDENS E M, SIMONET F, et al. Propeller noise from a light aircraft for low-frequency measurements of the speed of sound in a marine sediment[J]. Journal of Computational Acoustics, 2002, 10(4): 445–464. doi: 10.1142/S0218396X02001760 [6] FERGUSON B G. A ground-based narrow-band passive acoustic technique for estimating the altitude and speed of a propeller-driven aircraft[J]. The Journal of the Acoustical Society of America, 1992, 92(3): 1403–1407. doi: 10.1121/1.403934 [7] 陈韶华, 陈川, 郑伟. 单矢量水听器线谱多目标分辨研究[J]. 电子与信息学报, 2010, 32(5): 1253–1256. doi: 10.3724/SP.J.1146.2009.00460CHEN Shaohua, CHEN Chuan, and ZHENG Wei. Multi-sources distinguishing by exploring line spectra with single acoustic vector sensor[J]. Journal of Electronics &Information Technology, 2010, 32(5): 1253–1256. doi: 10.3724/SP.J.1146.2009.00460 [8] 王彪, 陈宇, 徐千驰, 等. 非理想条件下基于矢量水听器阵列的一种快速方位估计算法[J]. 电子与信息学报, 2021, 43(3): 745–751. doi: 10.11999/JEIT200541WANG Biao, CHEN Yu, XU Qianchi, et al. A fast direction estimation algorithm based on vector hydrophone array under non-ideal conditions[J]. Journal of Electronics &Information Technology, 2021, 43(3): 745–751. doi: 10.11999/JEIT200541 [9] LO K W and FERGUSON B G. Broadband passive acoustic technique for target motion parameter estimation[J]. IEEE Transactions on Aerospace and Electronic Systems, 2000, 36(1): 163–175. doi: 10.1109/7.826319 [10] LO K W. Flight parameter estimation using instantaneous frequency and direction of arrival measurements from a single acoustic sensor node[J]. The Journal of the Acoustical Society of America, 2017, 141(3): 1332–1348. doi: 10.1121/1.4976091 [11] 刘凯悦, 彭朝晖, 张灵珊, 等. 水下水平阵对空中运动声源的线谱探测[J]. 声学学报, 2019, 44(4): 566–575. doi: 10.15949/j.cnki.0371-0025.2019.04.017LIU Kaiyue, PENG Zhaohui, ZHANG Lingshan, et al. Line spectrum detection of airborne moving source by underwater horizontal array[J]. Acta Acustica, 2019, 44(4): 566–575. doi: 10.15949/j.cnki.0371-0025.2019.04.017 [12] 安春莲, 杨古月, 杨延菊. 基于中值滤波预处理的强冲击噪声背景测向方法[J]. 电子学报, 2021, 49(6): 1159–1166. doi: 10.12263/DZXB.20200392AN Chunlian, YANG Guyue, and YANG Yanju. DOA estimation under strong impulsive noise based on median value filtering[J]. Acta Electronica Sinica, 2021, 49(6): 1159–1166. doi: 10.12263/DZXB.20200392 [13] URICK R J. Noise signature of an aircraft in level flight over a hydrophone in the sea[J]. The Journal of the Acoustical Society of America, 1972, 52(1A): 172. doi: 10.1121/1.1982074 [14] FERGUSON B G, CULVER R L, and GEMBA K L. International student challenge problem in acoustic signal processing 2019[J]. Acoustics Today, 2019, 15(1): 71–73. doi: 10.1121/AT.2019.15.1.73 [15] DIPASSIO III J. Tre DiPassio’s solution to the 2019 international student challenge problem in acoustic signal processing[EB/OL]. https://static1.squarespace.com/static/5d362ec264480900010504be/t/5d765d6527a38b3a5d5be508/1568038249817/TreDiPassio_ChallengeProblem2019Solution.pdf, 2019. [16] SUN Weitao, WANG Huigang, GU Qingyue, et al. Exact frequency estimation in the i. i. d. noise via KL divergence of accumulated power[J]. IEEE Communications Letters, 2021, 25(8): 2574–2578. doi: 10.1109/LCOMM.2021.3077315 [17] SUN Weitao, WANG Huigang, GU Qingyue, et al. Exact and robust time-frequency estimation via accumulation of phase-difference power on multiple log-sum[EB/OL]. https://doi.org/10.13140/RG.2.2.21036.59523, 2021. [18] SADEGHI M, BEHNIA F, and AMIRI R. Window selection of the Savitzky-Golay filters for signal recovery from noisy measurements[J]. IEEE Transactions on Instrumentation and Measurement, 2020, 69(8): 5418–5427. doi: 10.1109/TIM.2020.2966310 期刊类型引用(18)
1. 寇子明，李腾宇，吴娟. 矿井提升机天轮智能润滑系统设计及控制策略. 煤炭学报. 2024(02): 1230-1239 . 
百度学术2. 崔胜胜，李汐，牟颖莹，李振. 间接接入式直流电能表异常计量数据识别算法. 微型电脑应用. 2024(07): 130-133 . 百度学术3. 毛冬，张辰，陈又咏，刘永清，焦艳斌. 基于国产CPU环境的国产数据库历史数据迁移技术. 太赫兹科学与电子信息学报. 2024(10): 1154-1160 . 百度学术4. 潘海霞，曹宁. 基于SDN的数据中心端口故障检测方法设计. 自动化与仪器仪表. 2023(02): 56-60 . 百度学术5. 刘浪，时宏伟. 基于注意力机制的CNN-LSTM的ADS-B异常数据检测. 计算机系统应用. 2023(04): 94-103 . 百度学术6. 卢宁. 基于加权广义匹配滤波联合SVDD的FOD检测算法. 现代雷达. 2023(05): 112-119 . 百度学术7. 何佳佳，安大庆，叶美丽. 面向学生体质测量的数据异常可视化检测方法. 信息技术. 2023(06): 83-87+93 . 百度学术8. 王鹏，董方正，张伟，汪克念，赵长啸. 基于TriLSTM-SVDD的ADS-B欺骗数据检测方法. 电光与控制. 2023(10): 77-81+119 . 百度学术9. 李妮锶，林琳，姚成. ADS-B系统安全问题解决策略对比分析. 信息技术. 2022(03): 30-35+41 . 百度学术10. 杨章勇，李欢. 工业以太网通信协议关键方法的优化与仿真. 计算机仿真. 2022(03): 377-380+395 . 百度学术11. 符士侃，夏元轶，杜钰，石廷川. 基于概率主题模型的大数据平台隐私泄露自动检测方法. 自动化与仪器仪表. 2022(04): 115-118+123 . 百度学术12. 曾雄飞. 基于粒子群算法优化BP神经网络的PID控制算法. 电子设计工程. 2022(11): 69-73+78 . 百度学术13. 汤双霞. 一种毫米波雷达机场跑道异物检测算法. 激光与红外. 2022(06): 820-826 . 百度学术14. 暴佳伟，田小平，刘宇娜，邹长宽，张雨晴. ADS-B安全问题研究综述. 现代防御技术. 2022(05): 28-35 . 百度学术15. 李娟，李海瑛，边玲，刘秋红. 基于跨模态深度度量学习的时序数据异常检测. 计算机仿真. 2022(10): 533-537 . 百度学术16. 刘明，张弘. AI交互终端异常数据入侵识别与仿真. 计算机仿真. 2022(11): 467-471 . 百度学术17. 张维，杜兰. 一种集成式Beta过程最大间隔一类分类方法. 电子与信息学报. 2021(05): 1219-1227 . 
本站查看18. 郭峰. 远距离多信道光纤通信网络异常数据检测研究. 激光杂志. 2021(09): 98-103 . 百度学术其他类型引用(13)
-
下载:
百度学术
计量
- 文章访问数: 411
- HTML全文浏览量: 253
- PDF下载量: 73
- 被引次数: 31
下载:
