基于能量收集的UAV-D2D网络资源分配算法

王茜竹; 胡洪瑞; 徐勇军; 张海波; 周继华; 陈莉

doi:10.11999/JEIT210491

基于能量收集的UAV-D2D网络资源分配算法

doi: 10.11999/JEIT210491

王茜竹¹,
胡洪瑞¹,
徐勇军^{1, 2, ,},
张海波¹,
周继华³,
陈莉¹

1.
重庆邮电大学通信与信息工程学院重庆 400065
2.
西安邮电大学陕西省信息通信网络及安全重点实验室西安 710121
3.
航天新通科技有限公司重庆 401332

基金项目: 国家自然科学基金(61601071)，重庆市自然科学基金(cstc2019jcyj-xfkxX0002)，陕西省信息通信网络及安全重点实验室(ICNS201904)，重庆市留学人员回国创业创新支持计划项目(CX2020095)

详细信息

作者简介:
王茜竹：女，1975年生，正高级工程师，硕士生导师，研究方向为下一代移动通信、物联网技术等

胡洪瑞：男，1998年生，硕士生，研究方向为异构无线网络资源分配

徐勇军：男，1986年生，副教授，硕士生导师，研究方向为D2D通信、能量收集、异构无线网络资源分配等

张海波：男，1979年生，副教授，硕士生导师，研究方向为无线网络资源分配、车联网等

周继华：男，1979年生，研究员，博士生导师，研究方向为移动网络、无线通信、5G、6G等

陈莉：女，1979年生，工程师，研究方向为无线传输技术

通讯作者:
徐勇军　xuyj@cqupt.edu.cn

中图分类号: TN929.538
计量
- 文章访问数: 1118
- HTML全文浏览量: 577
- PDF下载量: 145
- 被引次数: 18
出版历程
- 收稿日期: 2021-06-01
- 修回日期: 2022-02-13
- 录用日期: 2022-02-16
- 网络出版日期: 2022-02-23
- 刊出日期: 2022-03-28

Resource Rllocation for UAV-assisted D2D Communications with Energy Harvesting

1.
School of Communication and Information Engineering, Chongqing University of Posts and Telecommunications, Chongqing 400065, China
2.
Shaanxi Key Laboratory of Information Communication Network and Security, Xi’an University of Posts & Telecommunications, Xi’ an 710121, China
3.
Aerospace New Generation Communications Co., Ltd, Chongqing 401332, China

Funds: The National Natural Science Foundation of China (61601071), The Natural Science Foundation of Chongqing (cstc2019jcyj-xfkxX0002), The Shaanxi Key Laboratory of Information Communication Network and Security (ICNS201904), The Chongqing Entrepreneurship and Innovation Program for the Returned Overseas Chinese Scholars (CX2020095)

摘要

摘要: 为更好地利用周围环境中的射频信号能量，提升终端直连(D2D)通信的运行时间和无人机(UAV)通信的频谱利用率，该文提出一种基于能量收集的UAV-D2D网络资源分配算法。考虑UAV最大发射功率和移动性约束，蜂窝用户和D2D用户的最小速率约束，建立了系统和速率最大化的多变量耦合资源分配问题。利用连续凸近似和变量替换方法将混合整数非线性规划问题转化为凸优化问题，并利用拉格朗日对偶方法获得闭式解。仿真结果表明，所提算法具有良好的收敛性能，并能够有效提升系统容量。
- 无人机 /
- 终端直通 /
- 能量收集 /
- 资源分配
Abstract: In order to utilize better the surrounding radio-frequency energy and improve the operation lifetime of Device-to-Device (D2D) communications as well as the spectrum efficiency of Unmanned Aerial Vehicle (UAV) communication, a resource allocation algorithm is proposed for UAV-D2D networks with energy harvesting. Considering the constraints of the maximum transmit power and the mobility of the UAV, the minimum rate requirements of both cellular users and D2D users, a multivariable coupling resource allocation problem is formulated to maximize the sum rates of both cellular users and D2D users. The mixed-integer nonlinear programming problem is transformed into a convex optimization problem by using the successive convex approximation and variable substitution methods, where the closed-form solutions are obtained by the using Lagrange dual method. Simulation results demonstrate that the proposed algorithm has good convergence performance and higher system capacity.
- Unmanned Aerial Vehicle (UAV) /
- Device to Device (D2D) /
- Energy harvesting /
- Resource allocation

HTML全文

1. 引言

逆合成孔径雷达(Inverse Synthetic Aperture Radar, ISAR)作为一种主动的微波探测手段，可实现针对非合作运动目标的全天时、全天候、远距离、高分辨成像探测。ISAR通过发射并接收大带宽信号获得距离向高分辨，通过在相干积累时间(Coherent Processing Interval, CPI)内合成一定目标视角变化的多脉冲回波获得方位向高分辨^[1,2]。典型的ISAR系统常放置于地面，用于监测空中威胁目标。由于空中目标相对空域稀疏，ISAR接收回波信号具有稀疏特性，可利用压缩感知(Compressed Sensing, CS)理论，以远低于奈奎斯特采样率的回波数据恢复ISAR成像结果^[3]。压缩感知理论由Donoho^[4]提出，该理论基于目标稀疏性假设，可实现图像或信号的稀疏重构，从而建立了由低维观测数据精确恢复高维目标信号的理论体系^[5]，在雷达信号处理领域应用广泛。

压缩感知理论框架下，在目标满足稀疏性的前提时，可通过有限的低维观测数据恢复高分辨的高维信号^[6,7]。其中一类经典的信号恢复算法为贪婪算法，该类算法主要利用最小均方误差的思想，通过人为设置阈值，在每个迭代步骤中选择局部最优解，并在算法收敛后获得全局最优稀疏解。文献[8]将贪婪算法应用于稀疏信号恢复，从理论方面验证了此类方法的适用性。尽管该类算法复杂程度较低，运算速度快，但人为设置阈值的方式并不总能保证算法收敛后的解为最优稀疏解，有时只是局部最优，计算结果稳定性欠佳。

另一类稀疏恢复算法为凸优化类算法，该类算法是通过利用优化策略对稀疏信号进行有效恢复，主要解决关于目标的 $\ell _0$ 范数最优化问题。已有数学推导证明 $\ell _0$ 范数优化问题为非确定性问题(Non-deterministic Problem, NP)^[9]，可用最小化 $\ell _1$ 范数方法进行近似计算， $\ell _1$ 范数可认为是最紧致的近似稀疏求解算法。典型的优化类算法包括凸优化(ConVeX, CVX)^[10]和迭代软阈值算法^[11]等，其中CVX算法基于不断的优化设计最终得到最优解，是一种标准的凸优化算法，在工程实验中应用广泛。但在高维数据情况下，典型的CVX算法面临着运行性能低，耗时过长的缺点。而迭代软阈值算法通过对目标函数进行近似替代，并将软阈值与优化策略相结合，提高了算法计算效率和收敛性能^[11]，但由于迭代过程中需要手动调节迭代步长，算法自主学习程度不高。文献[12]将正则化范数与凸优化算法结合，用于SAR超分辨成像，证明了该类方法在效率和性能方面的优势，但繁琐的手动调整正则项系数步骤严重限制了该类算法在超分辨成像的精度及实用性。

还有一类稀疏信号恢复算法是贝叶斯类算法^[13,14]，该类算法主要结合贝叶斯机器学习的思想对信号进行重建。典型的贝叶斯类算法包括变分贝叶斯(Variational Bayes, VB)推断算法和期望最大化(Expectation-Maximization, EM)算法等。文献[15]将贝叶斯机器学习结合应用于SAR地面动目标成像，选取拉普拉斯分布进行先验建模，利用VB-EM算法进行贝叶斯推理计算，并通过实测数据获得了高精度的恢复结果。贝叶斯类算法的优点在于建模过程应用灵活，而且能提供参数的区间估计。但传统基于贝叶斯机器学习的成像算法在贝叶斯推理过程中面临复杂的参数调整步骤和高维的矩阵逆运算的问题，计算效率较低。

基于以上分析，本文提出一种基于贝叶斯Lasso的稀疏ISAR超分辨成像算法，研究建立凸优化Lasso模型与贝叶斯概率模型之间的等价联系，解决凸优化稀疏算法中面临的繁琐手动调节正则项系数的问题。本文通过在分层贝叶斯框架下实现对稀疏Laplace先验的概率建模，同时考虑对优化过程中相应 ${\ell _1}$ 范数正则项稀疏调节系数建立条件概率依赖模型，并利用马尔科夫链蒙特卡洛(Markov Chain Monte Carlo, MCMC)采样算法从数据中自动学习产生正则项系数，避免繁琐手动调节步骤，最终实现对目标超分辨ISAR成像的精确恢复。针对高维情况下超参数难以求解的问题，本文采用吉布斯(Gibbs)采样的方法对建模的参数进行统计采样估计。Gibbs采样方法是一种经典的MCMC算法，在统计机器学习方面应用广泛^[16]。本文将贝叶斯Lasso模型与Gibbs采样方法结合，避免了复杂的参数调整步骤，实现了算法自动化程度提升，增强了稀疏恢复算法的稳健性，提高了算法实用性和恢复精度。本文实验部分应用多种ISAR仿真和实测数据，通过与稀疏贝叶斯学习(Sparse Bayesian Learning, SBL)算法和CVX算法进行对比实验，从抑噪和超分辨性能等不同角度定性的分析了所提算法在性能上的优势，验证所提算法在ISAR图像恢复方面的高精度性和实用性。

2. 回波信号模型

ISAR成像过程中，雷达保持静止状态，非合作目标通过相对运动形成虚拟孔径。根据ISAR成像原理可将目标相对雷达视线方向的运动分为平动和转动，其中平动会造成成像时目标散焦，需要进行平动补偿，转动提供了方位向分辨率，是成像过程中需要利用的关键部分。经过包络对齐和相位自校准后，ISAR成像可等效为转台几何模型如图1所示，其中目标飞机从A位置运动到B位置。在目标运动过程中，目标观测角发生变化，用 $\theta \left( t \right)$ 表示，其中 $t$ 为慢时间。坐标系的原点O与目标中心重合，Y轴与雷达观测目标视线一致，X轴与Y轴垂直，雷达与目标中心点之间的距离以 $R(t)$ 表示，初始距离设为 ${R_0}$ ，则目标散射点与雷达之间的距离可表示为

图 1 ISAR成像示意图

下载: 全尺寸图片幻灯片

$R\left( t \right){\rm{ = }}{R_0} + {X_l}\sin \theta \left( t \right) + {Y_l}\cos \theta \left( t \right)$

(1)

其中 $\left( {{X_l},{Y_l}} \right)$ 为目标散射点的位置坐标。

当CPI较大时，ISAR成像分辨率高，但由于高阶运动分量的影响，目标运动复杂多变，处于非匀速运动状态，使得在长CPI内对ISAR的成像分析面临很大的困难。当CPI较小时，ISAR成像积累时间短，目标可认为处于匀速运动状态，目标运动状态简单，因此成像过程易于分析。但由于CPI小时，雷达接收回波数量过少，会造成成像分辨率低，难以满足ISAR超分辨的需求。针对此问题，本文利用压缩感知理论，在有限CPI内进行稀疏成像，实现ISAR超分辨，保证最终成像分辨率^[17,18]。在CPI较小时雷达观测角可表示为 $\theta \left( t \right) = \omega t$ ，其中 $\omega$ 为雷达转动角速度。由于有限CPI内雷达回波积累时间相对较短，此时式(1)可进一步表示为

$R(t) \approx {R_0} + {X_l}\omega t + {Y_l}$

(2)

通常ISAR发射信号为宽带信号，在对接收回波进行距离向去斜处理，包络对齐和相位自校准补偿^[19,20]后，散射点回波信号为

$\begin{split} S\left( {{t_r},t} \right)=\,&\sum\limits_{l = 1}^L {A_l}{\rm{sinc}} \left( {B\left( {{t_r} - \frac{{2\left( {{X_l}\omega t + {Y_l}} \right)}}{{\rm{c}}}} \right)} \right)\\ & \cdot \exp \left( { - {\rm j}{\rm\pi} \frac{{4{f_{\rm{c}}}}}{{\rm{c}}}\left( {{X_l}\omega t + {Y_l}} \right)} \right) \\[-15pt] \end{split}$

(3)

其中 $L$ 表示散射点数目， ${A_l}$ 为散射点后向散射系数， $B$ 为发射信号带宽， ${f_{\rm{c}}}$ 为载波频率， ${\rm{c}}$ 为光速， ${\rm{sinc}}$ 函数表示距离压缩后的响应函数。

如式(3)所示的回波信号可以简化表示成矩阵形式，如式(4)所示

${{Y}}{\rm{ = }}{{AX}} + {{N}}$

(4)

其中 ${{Y}} \in {{{C}}^{K \times N}}$ 表示雷达接收回波复数据，对应式(3)中的 $S\left( {{t_r},t} \right)$ , ${{N}} \in {{{C}}^{K \times N}}$ 表示加性噪声或杂波， ${{X}} \in {{{C}}^{M \times N}}$ 为需要进行重构的成像目标， ${{A}}$ = ${\left\{ {{a_{km}}} \right\}_{K \times M}}$ 为方位向傅里叶字典，可构建为

$\left. \begin{aligned} & {{A}} = \left\{ {{{a}}\left( {{f_{\rm{d}}}\left( 1 \right)} \right)\; ··· \;{{a}}\left( {{f_{\rm{d}}}\left( M \right)} \right)} \right\} \\ & {{a}}\left( {{f_{\rm{d}}}\left( m \right)} \right) = {\left\{ {{{\rm{e}}^{ - {\rm{j}}2{\rm\pi} {f_{\rm{d}}}\left( {\rm{m}} \right){t_1}}}\; ··· \;{{\rm{e}}^{ - {\rm{j}}2{\rm\pi} {f_{\rm{d}}}\left( m \right){t_K}}}} \right\}^{\rm{T}}} \\ \end{aligned} \right\}$

(5)

其中 ${f_{\rm{d}}}$ 表示多普勒频率， $t$ 表示慢时间， $1 \le m \le M$ 。

3. 贝叶斯Lasso

观察如式(4)所示的ISAR信号模型，目标成像问题可以解释为关于散射场 ${{X}}$ 的线性方程组求解问题。由于ISAR接收回波有时是不完整的，在强噪声和干扰情况下关于目标 ${{X}}$ 的线性方程组求解的数学模型一般是非齐次欠定的，进而造成了解的非唯一性和计算复杂性。针对这个问题，可利用压缩感知技术提供解决途径，在稀疏先验假设条件下，将非齐次欠定线性方程组的求解转变为对Lasso问题进行求解，数学模型如式(6)所示

${\hat{ X}} = \mathop {\arg \min }\limits_{{X}} \frac{1}{2}\left\| {{{Y}} - {{AX}}} \right\|_{\rm{F}}^2 + \lambda {\left\| {{X}} \right\|_1}$

(6)

其中 ${\left\| \cdot \right\|_1}$ 和 ${\left\| \cdot \right\|_{\rm{F}}}$ 分别表示 $\ell _1$ 范数和Frobenius范数，等号右侧第1项为保真项，表征恢复结果精度，第2项为 $\ell _1$ 范数正则项，代表引入的稀疏先验条件， $\lambda$ 为正则项约束调节系数，通常 $\lambda > 0$ ，调节正则化系数 $\lambda$ 可以控制在模型中施加正则项的程度，即稀疏度。但在不同降采样率和噪声条件下正则项系数变化很大，繁琐的手动调整步骤严重限制了算法的实用性和成像的精度。为实现参数的自动调整，本文深入挖掘凸优化与统计学习理论间的等价性，通过将Lasso问题与贝叶斯概率模型结合，建立贝叶斯Lasso模型，并将正则项系数引入分层框架内，求取正则项系数依赖的条件概率分布，最后利用Gibbs采样进行贝叶斯推理，有效提高了算法的稳定性和自适应性。

通常，距离压缩后不同距离单元之间独立同分布，且服从圆对称复高斯分布。此时，由式(4)可得回波数据似然函数服从式(7)分布

$\begin{split} p({{Y}}\left| {{{X}},\beta } \right.)\,&= \prod\limits_{n = 1}^N {{\cal C}{\cal N}\left( {{{{Y}}_{:n}}\left| {{{A}}{{{X}}_{:n}},{\beta ^{{\rm{ - }}1}}{{{I}}_M}} \right.} \right)}\\ & = {\left( {\frac{\beta }{{2{\rm\pi} }}} \right)^{KN/2}}\exp \left( { - \frac{\beta }{2}\left\| {{{Y}} - {{AX}}} \right\|_{\rm{F}}^2} \right) \end{split}$

(7)

其中 $\beta$ 表示噪声精度，等于噪声方差的倒数，且 $\beta > 0$ 。 ${\cal C}{\cal N}\left( \cdot \right)$ 代表复高斯函数分布， ${{{I}}_M}$ 表示 $M \times M$ 的单位矩阵。

为得到散射场 ${{X}}$ 的后验概率分布，需要选取符合目标特征的先验分布。由于ISAR目标具有稀疏性的特点，本文利用拉普拉斯先验分布对目标场景进行先验建模。拉普拉斯分布具有良好的重尾特性，可以有效地模拟目标的稀疏特征^[7]。但由于拉普拉斯先验分布与式(7)的高斯似然函数不共轭，因此无法直接计算后验分布的闭合解析解。针对这个问题，本文采用层级(hierarchical)贝叶斯模型对先验进行概率建模，通过引入超参数构建层级之间变量的紧密联系，进而实现ISAR成像目标 ${{X}}$ 的有效估计计算。

由于ISAR成像回波数据为复数形式，本文在分层模型的第1步将 ${{X}}$ 建模成依赖于超参数 ${{\alpha}}$ 的复高斯分布，建模过程如式(8)所示

$\begin{split} p({{X}}\left| {{\alpha}} \right.)=\,&\prod\limits_{m = 1}^M {\prod\limits_{n = 1}^N {{\cal C}{\cal N}\left( {{{{X}}_{mn}}\left| {0,{{\alpha}} _{mn}^{{\rm{ - }}1}} \right.} \right)} } {\rm{ = }}{\left( {\frac{{{{{\alpha}} _{mn}}}}{{2{\rm\pi} }}} \right)^{MN}}\\ & \cdot\prod\limits_{m = 1}^M {\prod\limits_{n = 1}^N {\exp \left( { - \frac{{{{{\alpha}} _{mn}}}}{2}{{\left| {{{{X}}_{mn}}} \right|}^2}} \right)} } \\[-20pt] \end{split}$

(8)

其中 ${{\alpha}}$ 表示方差倒数矩阵，且 ${{\alpha}}$ 的第 $mn$ 个元素 ${{{\alpha}} _{mn}} > 0$ 。由于高斯函数固有的平滑特性，式(8)不属于稀疏分布范畴。为促进先验的稀疏特性并简化后验分布的推导，本文选取与高斯分布成对共轭的伽马(Gamma)分布作为超参数 ${{\alpha}}$ 的先验，并在分层模型的第2步将超参数 ${{\alpha}}$ 建模成如式(9)的形式

$\begin{split} p\left( {{{\alpha}} \left| \lambda \right.} \right) =\,& \prod\limits_{m = 1}^M {\prod\limits_{n = 1}^N {{\rm{Gamma}}\left( {{{{\alpha}} _{mn}}\left| {\eta ,\lambda } \right.} \right)} } \\ = \,&\prod\limits_{m = 1}^M \prod\limits_{n = 1}^N \frac{{{\lambda ^\eta }}}{{\Gamma (\eta )}}{{\left( {{{{\alpha}} _{mn}}} \right)}^{\eta - 1}}\\ & \cdot \exp \left( { - \lambda {{{\alpha}} _{mn}}} \right){{\cal I}_{{\mathbb R} + }}\left( {{{{\alpha}} _{mn}}} \right) \end{split}$

(9)

其中 ${{\cal I}_{{\mathbb R} + }}$ 代表指示函数，当 ${{{\alpha}} _{mn}} > 0$ 时， ${{\cal I}_{{\mathbb R} + }}\left( {{{{\alpha}} _{mn}}} \right) = 1$ ，当 ${{{\alpha}} _{mn}} \le 0$ 时， ${{\cal I}_{{\mathbb R} + }}\left( {{{{\alpha}} _{mn}}} \right)$ =0。根据式(8)和式(9)可计算边缘分布，如 $p({{X}}\left| \lambda \right.)$ = $\displaystyle\int p ({{X}}\left| {{\alpha}} \right.)p\left( {{{\alpha}} \left| \lambda \right.} \right){\rm{d}}{{\alpha}}$ 。经过验证，当 $\eta$ 的值为10^–3时，该边缘分布可服从预期的Laplace分布，从而可以保证稀疏性的先验假设。此外，式(9)中的超参数 $\lambda$ 与式(6)的Lasso模型中正则项系数的作用是等效的。通过对 $\lambda$ 赋予相应的无信息先验分布，可在贝叶斯Lasso框架下将正则项系数引入到分层建模框架内，通过推导建立相应的后验概率分布可实现参数自学习的目的。

为提高ISAR成像精度和分辨率，并实现参数的自动调整，本文通过对正则化系数 $\lambda$ 设定一个先验分布并进行分层贝叶斯建模，将惩罚项系数 $\lambda$ 作为超参数引入分层贝叶斯中。同时求取 $\lambda$ 相应的后验概率密度分布，进而将惩罚系数与后验概率密度分布相联系，实现参数自主学习。为减少外界人为因素的影响，本文将 $\lambda$ 的先验建模为Gamma分布，概率模型如式(10)所示

$\begin{split} p(\lambda \left| {a,b} \right.) =\,& {\rm{Gamma}}\left( {\lambda \left| {a,b} \right.} \right) = \frac{{{b^a}}}{{\varGamma (a)}}{\left( \lambda \right)^{a - 1}}\\ & \cdot \exp \left( { - b\lambda } \right){{\cal I}_{{\mathbb R} + }}\left( \lambda \right) \end{split}$

(10)

其中 $a$ 和 $b$ 的值均固定为10^-3， ${{\cal I}_{{\mathbb R} + }}$ 代表指示函数。考虑到噪声功率的随机性，高斯噪声精度 $\beta$ 的先验分布也采用Gamma分布，建模如式(11)所示

$\begin{split} p(\beta \left| {c,d} \right.) = \,&{\rm{Gamma}}\left( {\beta \left| {c,d} \right.} \right) = \frac{{{d^c}}}{{\varGamma (c)}}{\left( \beta \right)^{c - 1}}\\ & \cdot \exp \left( { - d\beta } \right){{\cal I}_{{\mathbb R} + }}\left( \beta \right) \end{split}$

(11)

其中 $c$ 和 $d$ 分别为Gamma分布的形状参数和尺度参数，一般将 $c$ 和 $d$ 设置为很小的值以保证先验的无信息性，通过交叉实验验证，本文将 $c$ 和 $d$ 的值均固定为10^–5。

综上所述，可将 ${{X}}$ 和各参数的联合后验概率密度分布表示为

$\begin{split} & p({{X}},{{\alpha}} ,\lambda ,\beta \left| {{Y}} \right.) \\ & \quad = \frac{{p({{Y}}\left| {{{X}},\beta } \right.)p({{X}}\left| {{\alpha}} \right.)p({{\alpha}} \left| \lambda \right.)p(\lambda )p(\beta )}}{{p({{Y}})}} \end{split}$

(12)

由于 $p\left( {{Y}} \right)$ 与 ${{X}}$ 无关，因此

$\begin{split} & p({{X}},{{\alpha}} ,\lambda ,\beta \left| {{Y}} \right.) \propto p({{Y}}\left| {{{X}},\beta } \right.)p({{X}}\left| {{\alpha}} \right.)p({{\alpha}} \left| \lambda \right.)\\ & \quad \cdot p(\lambda )p(\beta ) \end{split}$

(13)

符号“ $\propto$ ”代表“正比例于”函数关系。将式(7)—式(11)代入式(13)可得

$\begin{split} & p({{X}},{{\alpha}} ,\lambda ,\beta \left| {{Y}} \right.) \propto {\left( {\frac{\beta }{{2{\rm\pi} }}} \right)^{KN/2}}\\ & \quad \cdot \exp \left( { - \frac{\beta }{2}\left\| {{{Y}} - {{AX}}} \right\|_{\rm{F}}^2} \right) \\ & \quad \cdot \prod\limits_{m = 1}^M {\prod\limits_{n = 1}^N {{\cal C}{\cal N}\left( {{{{X}}_{mn}}\left| {0,{{\alpha}} _{mn}^{{\rm{ - }}1}} \right.} \right)} } \\ & \quad \cdot \prod\limits_{m = 1}^M {\prod\limits_{n = 1}^N {{\rm{Gamma}}\left( {{{{\alpha}} _{mn}}\left| {\eta ,\lambda } \right.} \right)} } \\ & \quad \cdot {\rm{Gamma}}\left( {\lambda \left| {a,b} \right.} \right) \cdot {\rm{Gamma}}\left( {\beta \left| {c,d} \right.} \right) \end{split}$

(14)

上述贝叶斯分层中参数与超参数的关系可用有向无环图(DAG)进行表示，对应贝叶斯模型如图2所示。其中 ${{Y}}$ 表示ISAR接收回波数据， ${{X}}$ 为待恢复的成像目标， ${{\alpha}}$ 和 $\lambda$ 为引入的超参数。 ${{X}}$ 和 ${{\alpha}}$ 的紧密结合能够促进先验的稀疏性，而对 $\lambda$ 的概率建模则将正则项系数引入层级模型，为参数自学习提供了条件，虚线方框中的 $a$ , $b$ , $c$ , $d$ 和 $\eta$ 均为固定参数。

图 2 贝叶斯层级模型(有向无环图，DAG)

下载: 全尺寸图片幻灯片

常见文献中对式(12)后验分布的求解方法可分为两种，一种是利用VB的方法进行近似求解，另外一种是利用基于统计采样的MCMC采样方法进行估计计算。在高维情况下，采用VB的方法进行后验求解面临繁杂的参数调整和矩阵逆运算，算法自学习能力不强，需要人工手动介入，计算效率较低。而MCMC统计采样方法依赖马尔科夫链进行递归采样，在进行贝叶斯推理的过程依赖各参数的条件后验概率分布，是一个数据驱动的过程，无需过多的参数调整步骤，在算法自学习方面具有很大的优势。为提高算法性能和成像精度，减少复杂的人为介入参数调整，本文在贝叶斯层级模型下建立各参数依存的概率分布关系，并利用经典的Gibbs采样算法进行后验概率估计，最终达到对目标散射场 ${{X}}$ 的高精度恢复。

4. Gibbs采样算法

Gibbs采样算法是MCMC方法中一种经典的采样方式，其以变量的后验概率分布为基础，构造具有马尔可夫统计特性的转移矩阵，并以转移矩阵为采样核构建采样链，进而对变量及其统计特征进行采样估计。Gibbs采样算法通过将高维度计算问题进行分解，转换为低维度问题进行运算，并且能够对随机采样样本充分保留，具有高接受率。因此本文利用Gibbs采样算法对所提出的贝叶斯Lasso层级模型中的变量和超参数进行统计采样计算，并对参数依赖的后验分布进行估计，实现对待恢复目标的高精度求解。

在应用Gibbs采样算法进行采样计算之前，需要对各参数服从的后验概率密度分布进行推导。首先根据式(11)中噪声精度 $\beta$ 的先验分布与式(7)中的似然函数可得精度 $\beta$ 的后验分布如式(15)所示

$\begin{split} & p\left( {\beta \left| {{{X}},{{Y}},{{\alpha}} ,\lambda } \right.} \right) \propto p\left( {{{Y}}\left| {{{X}},\beta } \right.} \right)p\left( \beta \right)\\ & \quad ={\rm{Gamma}}\left( {\beta \left| {c + {{KN} / 2},d + \frac{1}{2}\left\| {{{Y}} - {{AX}}} \right\|_{\rm{F}}^2} \right.} \right) \end{split}$

(15)

式(15)为标准的伽马分布。由式(9)和式(10)可将正则化系数 $\lambda$ 的后验概率密度分布构建为

$\begin{split} & p\left( {\lambda \left| {{{X}},{{Y}},{{\alpha}} ,\beta } \right.} \right) \propto p\left( {{{Y}}\left| {{{X}},{{\alpha}} ,\beta } \right.} \right)p({{\alpha}} \left| \lambda \right.)p(\lambda ) \\ & \quad \propto p({{\alpha}} \left| \lambda \right.)p(\lambda ){\rm{ = }}{\rm{Gamma}}\left( {\lambda \left| {a+{{MN} / 2},b+{{\left\| {{\alpha}} \right\|_2^2} / 2}} \right.} \right) \end{split}$

(16)

由式(8)和式(9)可将超参数 ${{\alpha}}$ 的后验分布推导如式(17)所示

$\begin{split} p\left( {{{\alpha}} \left| {{{X}},{{Y}},\lambda ,\beta } \right.} \right) \,& \propto p\left( {{{Y}}\left| {{{X}},{{\alpha}} ,\beta ,\lambda } \right.} \right)p({{X}}\left| {{\alpha}} \right.)p\left( \lambda \right) \\ & \propto p({{X}}\left| {{\alpha}} \right.)p\left( {{{\alpha}} \left| \lambda \right.} \right)\\[-10pt] \end{split}$

(17)

假设 ${{\alpha}}$ 的元素相互独立，则 ${{{\alpha}} _{mn}}$ 的后验分布可如式(18)表示

$\begin{split} & p\left( {{{{\alpha}} _{mn}}\left| {{{X}},\lambda } \right.} \right) \propto {\rm{Inv}} - {\rm{Gauss}}\left( {{{{\alpha}} _{mn}}\left| {\rho ,{\lambda ^\prime }} \right.} \right)\; \\ & \quad =\sqrt {\dfrac{{{\lambda ^\prime }}}{{2{\rm\pi} }}} {\left( {{{{\alpha}} _{mn}}} \right)^{ - 3/2}}\exp \left\{ { - \frac{{{\lambda ^\prime }{{({{{\alpha}} _{mn}} - \rho )}^2}}}{{2{\rho ^2}{{{\alpha}} _{mn}}}}} \right\} \end{split}$

(18)

其中 $\rho = \sqrt {\dfrac{{{\lambda ^2}}}{{X_{mn}^2}}}$ , ${\lambda ^\prime } = {\lambda ^2}$ , ${\rm{Inv}} - {\rm{Gauss}}\left( \cdot \right)$ 表示逆高斯分布。由式(7)和式(8)可得 ${{X}}$ 的条件后验分布为

$\begin{split} & p\left( {{{X}}\left| {{{Y}},{{\alpha}} ,\lambda ,\beta } \right.} \right) \propto p\left( {{{Y}}\left| {{{X}},{{\alpha}} ,\beta ,\lambda } \right.} \right)p({{X}}\left| {{\alpha}} \right.) \\ & \quad = {\cal C}{\cal N}\left( {{{X}}\left| {{{\mu}} ,{{\varSigma}} } \right.} \right) \end{split}$

(19)

其中

$\left. \begin{aligned} & {{\mu}} = \beta {{\varSigma}} {{{A}}^{\rm{H}}}{{Y}},\\ & {{\varSigma}} = {\left( {\beta {{{A}}^{\rm{H}}}{{A}} + {\rm{diag}}\left( {\frac{1}{{{{{\alpha}}_1}}}\; ··· \;\frac{1}{{{{{\alpha}}_{MN}}}}} \right)} \right)^{ - 1}} \end{aligned}\right\}$

(20)

其中 ${{\mu}}$ 和 ${{\varSigma}}$ 分别代表式(19)的均值与协方差矩阵。综上所述，提出本文所采用的贝叶斯Lasso算法计算流程，如表1所示。

表 1 贝叶斯Lasso算法流程

　(1) 初始化

${{{X}}^{\left( 0 \right)}}$ ,

${\beta ^{\left( 0 \right)}}$ ,

${{{\alpha}} ^{\left( 0 \right)}}$ 和

${\lambda ^{\left( 0 \right)}}$ 并构造字典

${{A}}$ ；

　(2) 设置老化期门限

$T$ 和迭代次数

$k$ =

$1,2, ··· ,K$ ，其中

$K$ 为迭代总次数且

$T < K$ ，开始循环采样；

　(3) 根据式(15)噪声精度

$\beta$ 的条件后验概率密度分布

$p\left( {\beta \left| {{{X}},{{Y}},{{\alpha}} ,\lambda } \right.} \right)$ 对

$\beta$ 进行采样；

　(4) 根据式(16)正则化系数

$\lambda$ 所服从的条件后验概率密度分布

$p\left( {\lambda \left| {{{X}},{{Y}},{{\alpha}} ,\beta } \right.} \right)$ 对

$\lambda$ 进行采样；

　(5) 根据式(17)超参数

${{\alpha}}$ 所服从的条件后验概率密度分布

$p\left( {{{\alpha}} \left| {{{X}},{{Y}},\lambda ,\beta } \right.} \right)$ 对

${{\alpha}}$ 进行采样；

　(6) 根据式(19)目标

${{X}}$ 所服从的条件后验概率密度分布

$p\left( {{{X}}\left| {{{Y}},{{\alpha}} ,\lambda ,\beta } \right.} \right)$ 对

${{X}}$ 进行采样；

　(7) 若

$k > T$ ，收集

${{X}}$ 的采样样本

$\left[ {{{{X}}^{\left( T \right)}},{{{X}}^{\left( {T + 1} \right)}}, ··· ,{{{X}}^{\left( K \right)}}} \right]$ ，重复步骤(3)—步骤(6)直至收敛。

下载: 导出CSV

| 显示表格

由于采样算法的马尔可夫链需要一定的迭代次数才能达到收敛状态，收敛后最终得到的样本就是服从后验概率分布的样本估计值，而收敛之前的迭代过程称为老化期(burn-in period)。处于老化期的采样结果具有较大的误差，在进行估算过程中需对这部分数据进行舍弃。表1中步骤(2)设置了老化期门限 $T$ ，经过一定的迭代后由Gibbs采样器得到的样本逐渐收敛到目标分布，此后由采样得到的样本可用于参数和变量的估计计算。为降低误差的影响，对收敛后的采样样本进行均值求解可得出最终 ${{X}}$ 的采样估计值，即为待恢复的ISAR成像目标散射场 ${{X}}$ 。

5. 实验验证

5.1 仿真数据验证

为进一步验证贝叶斯Lasso算法的优越性和通用性，本文利用仿真数据Mig-25对所提模型及算法的成像性能进行验证。通过在不同降采样率和信噪比下对目标场景进行恢复重建，并与由距离多普勒(Range-Doppler, RD), CVX和SBL算法得到的成像结果进行对比，从而验证所提算法的优越性和稳健性。图3给出了信噪比为5 dB降采样率变化时，通过不同方法所得的实验结果。其中图3第1列和第2列分别表示降采样率为0.25和未进行降采样时分别利用RD, SBL, CVX和本文贝叶斯Lasso算法进行处理后的恢复结果。从图3中可以看出利用RD方法进行成像时，成像结果背景噪声大，目标不明显。而SBL, CVX和贝叶斯Lasso 3种算法则实现了对噪声的抑制，以及对成像效果的增强，其中又以贝叶斯Lasso算法的效果最好，从而证明了本文贝叶斯Lasso算法的实用性。

图 3 信噪比5 dB不同降采样率成像结果

下载: 全尺寸图片幻灯片

为验证本文贝叶斯Lasso算法在抑噪方面的优越性能，图4给出了降采样率为0.5条件下信噪比变化时通过不同算法恢复的结果。其中图4第1列和第2列分别表示信噪比为0 dB和信噪比为5 dB时分别利用RD, SBL ,CVX和该文贝叶斯Lasso算法恢复的成像结果。从图4中可以看出RD重建的结果背景噪声很大，目标不明显，利用SBL和CVX算法恢复的图像虽然有部分抑噪效果，但目标丢失了多数特征。与这3种算法相比，采用该文贝叶斯Lasso算法恢复的目标散射场飞机轮廓特征突出，背景噪声消除干净，性能最优。

图 4 降采样率0.5不同信噪比成像结果

下载: 全尺寸图片幻灯片

5.2 实测数据验证

实测部分采用Yak-42数据^[12]进行实验分析，图5给出了信噪比5 dB条件下降采样率分别为0.25和未降采样的恢复结果。从图5中可以看出随着采样率的降低，RD恢复的目标淹没在噪声下，SBL算法恢复的结果目标丢失部分特征，CVX算法对噪声的抑制较弱，恢复的目标背景噪声斑点较多。而本文贝叶斯Lasso算法在低降采样率下保留目标大部分特征的同时，对噪声也起到了很好的抑制作用，总体性能优于SBL和CVX算法。

图 5 信噪比5 dB不同降采样率成像结果

下载: 全尺寸图片幻灯片

为验证本文贝叶斯Lasso算法在抑噪方面的优势，图6给出了降采样为0.5信噪比分别为0 dB和5 dB时不同算法重建的结果。其中在信噪比为0 dB时，RD算法得到的结果目标基本识别不清，并且背景噪声很大，SBL恢复的目标丢失大部分特征，并且机头部分基本无法识别，而CVX算法抑噪能力不足。从图6中可以看出，本文采用的贝叶斯Lasso算法在抑制噪声的同时保留了目标的多数特征，对目标进行了有效恢复，整体性能明显优于其他3种算法。

图 6 降采样0.5不同信噪比成像结果

下载: 全尺寸图片幻灯片

5.3 地面动目标实验验证

本部分通过对地面动目标进行成像分析，验证贝叶斯Lasso算法在实测数据的性能优势。地面动目标成像实验采用Gotcha-GMTI数据集^[21]，该数据集由美国空军实验室发布，经由三通道机载雷达获取，雷达工作于X波段，作用距离10 km，发射信号带宽约640 MHz, PRF约2100 Hz。由于在对地面动目标进行成像过程中目标的运动特性会引起成像的偏移与散焦，本文利用杂波相消并采用吕氏变换方法(LV’s Distribution, LVD)^[22]对调频发射信号中的多普勒频率和线性调频率进行协同估计，以改善动目标成像的聚焦和抑噪性能，提高成像分辨率。

图7(a)给出了采用RD变换后的参考成像，由图可见背景噪声抑制效果不明显。图7(b)、图7(c)分别为利用CVX和SBL算法的成像结果，从图中可以看出两种算法在抑噪方面均起到一定作用，但由于CVX算法难以适应不同距离单元的不同稀疏度，在成像过程中失去了部分目标特征，而SBL算法求解精度有限，目标恢复的效果不佳。图7(d)为采用本文贝叶斯Lasso算法得到的成像结果，与其他3种算法相比，由于贝叶斯Lasso算法无需人工干预参数选择，能够自主地适应空变的稀疏度，因此得到的成像结果信噪比更高，目标散射特征更清晰，从而验证了本文贝叶斯Lasso算法在成像性能方面的优越性。

图 7 地面动目标成像

下载: 全尺寸图片幻灯片

5.4 相变图实验分析

为了定性验证本文贝叶斯Lasso算法性能的优越性，本文利用Donoho等人^[23]提出的相位变化图(Phase Transition Diagram, PTD)进行实验分析。PTD是对算法精确性和稳定性进行评估的一种经典方式，其通过多次蒙特卡洛独立重复实验，获得恢复图像与参考图像相关量的平均值，并利用均方误差比较确定算法恢复的准确率，从而对算法的恢复精度进行定性分析。PTD将相变区域分为不可重建和可重建区域，并利用不同区域对算法的恢复性能进行直观的描述。本文在对算法进行评估分析时，通过100次蒙特卡洛实验得到重建图像与目标图像的平均相关值，并利用归一化均方误差(Normalized Mean Square Error, NMSE)进行判别图像恢复的准确率，NMSE的表达式为

${\rm{NMSE}}\; = \;\left\| {{{\left( {{{{X}}^k} - {{X}}} \right)} / {{X}}}} \right\|_2^2$

(21)

其中 ${{{X}}^k}$ 代表第 $k$ 次蒙特卡洛实验所得的恢复结果， ${{X}}$ 表示参考目标图像，阈值设定为0.0001，当 ${\rm{NMSE < 0}}{\rm{.0001}}$ 时认为图像重建成功。

图8为利用相变图对SBL,CVX和贝叶斯Lasso算法进行相变分析所得的结果。图8中各子图的右上方区域代表可重建区域，左下方区域代表不可重建区域。从图8中可以看出利用贝叶斯Lasso算法得到的相变图可重建区域要明显大于利用SBL算法得到的结果，并且与CVX得到的相变图非常接近，从而可以得出本文贝叶斯Lasso算法对目标的成像恢复能力要明显优于SBL算法，并且能够达到与CVX标准相近的性能。尽管CVX算法损耗大，运算时间长，并且高维情况下算法实用性不高，但由于CVX是一种标准的凸优化恢复算法，其可作为衡量其他算法性能优劣的参考方法。从图8(b)和图8(c)中可以看出由贝叶斯Lasso算法和CVX算法得到可重建区域非常接近，表明两者在算法恢复性能方面具有很高的相似度，从而验证了贝叶斯Lasso算法的优越性。

图 8 相变图

下载: 全尺寸图片幻灯片

6. 结束语

本文针对凸优化ISAR成像算法繁琐的参数调整过程，研究建立凸优化Lasso模型与贝叶斯概率模型之间的联系。通过结合贝叶斯理论，建立贝叶斯Lasso模型，并采用分层方法对参数建立依赖的概率分布关系，在求取条件后验分布的同时利用Gibbs采样算法进行参数后验估计，实现了ISAR目标高分辨成像。实验利用相变图进行定性比较分析，并利用仿真和实测数据进行算法性能验证，证明了本文采用的贝叶斯Lasso算法的有效性和优越性。但是，在研究过程中，笔者发现标量正则项系数不能很好地模拟目标场景，有必要研究矩阵化的正则项，并实现矩阵正则项系数的自动调整，这将是后续工作研究的重点。

图 1 基于能量收集的UAV-D2D网络系统模型

下载: 全尺寸图片幻灯片

图 2 速率收敛性能

下载: 全尺寸图片幻灯片

图 3 系统和速率与UAV最大发射功率的关系

下载: 全尺寸图片幻灯片

图 4 系统和速率与D2D最低速率要求的关系

下载: 全尺寸图片幻灯片

图 5 系统和速率与UAV飞行高度的关系

下载: 全尺寸图片幻灯片

图 6 系统和速率与用户簇数量关系图

下载: 全尺寸图片幻灯片

图 7 UAV飞行轨迹图

下载: 全尺寸图片幻灯片

图 8 不同算法下系统和速率与UAV最大发射功率的关系

下载: 全尺寸图片幻灯片

图 9 不同算法下系统和速率与UAV飞行周期的关系

下载: 全尺寸图片幻灯片

表 1 基于交替迭代的资源分配算法

初始化系统参数：用户簇个数、用户簇内D2D对数、D2D发射　机功率和UAV轨迹 ${{\mathbf{q}}_n}$ ；设置迭代收敛精度为 $\varepsilon$ ，最大迭代次数　 ${L_{\rm max}}$ ；穷尽搜索 $N$ 个时隙UAV与用户簇连接情况。
(1) 循环；
(2) 迭代次数更新 $l = l + 1$ ；
(3) 通过问题式(18)得到 ${\tau _{k,n}}$ 和 $p_{k,n}^{\rm C}$ ；
(4) 通过问题式(35)得到 $p_{k,n,m}^{\rm D}$ ；
(5) 通过问题式(47)得到 ${{\mathbf{q}}_n}$ ；
(6) 直到算法满足收敛条件或达到最大迭代次数，即　　 $R\left( {l + 1} \right) - R\left( l \right) \le \varepsilon$ 或 $l = {L_{{\rm{max}}} }$ 。

下载: 导出CSV

参考文献(27)

[1]	MOZAFFARI M, SAAD W, BENNIS M, et al. A tutorial on UAVs for wireless networks: Applications, challenges, and open problems[J]. IEEE Communications Surveys & Tutorials, 2019, 21(3): 2334–2360. doi: 10.1109/COMST.2019.2902862
[2]	ZENG Yong, ZHANG Rui, and LIM T J. Wireless communications with unmanned aerial vehicles: Opportunities and challenges[J]. IEEE Communications Magazine, 2016, 54(5): 36–42. doi: 10.1109/MCOM.2016.7470933
[3]	ZHANG Haijun, ZHANG Jianmin, and LONG Keping. Energy efficiency optimization for NOMA UAV network with imperfect CSI[J]. IEEE Journal on Selected Areas in Communications, 2020, 38(12): 2798–2809. doi: 10.1109/JSAC.2020.3005489
[4]	XU Yongjun, GUI Guan, GACANIN H, et al. A survey on resource allocation for 5G heterogeneous networks: Current research, future trends, and challenges[J]. IEEE Communications Surveys & Tutorials, 2021, 23(2): 668–695. doi: 10.1109/COMST.2021.3059896
[5]	ZHAO Nan, LU Weidang, SHENG Min, et al. UAV-assisted emergency networks in disasters[J]. IEEE Wireless Communications, 2019, 26(1): 45–51. doi: 10.1109/MWC.2018.1800160
[6]	WANG Feiran, XU Chen, SONG Lingyang, et al. Energy-efficient resource allocation for device-to-device underlay communication[J]. IEEE Transactions on Wireless Communications, 2015, 14(4): 2082–2092. doi: 10.1109/TWC.2014.2379653
[7]	徐勇军, 刘子腱, 李国权, 等. 基于NOMA的无线携能D2D通信鲁棒能效优化算法[J]. 电子与信息学报, 2021, 43(5): 1289–1297. doi: 10.11999/JEIT200175 XU Yongjun, LIU Zijian, LI Guoquan, et al. Robust energy efficiency optimization algorithm for NOMA-based D2D communication with simultaneous wireless information and power transfer[J]. Journal of Electronics &Information Technology, 2021, 43(5): 1289–1297. doi: 10.11999/JEIT200175
[8]	XU Yongjun, XIE Hao, LIANG Chengchao, et al. Robust secure energy-efficiency optimization in SWIPT-aided heterogeneous networks with a nonlinear energy-harvesting model[J]. IEEE Internet of Things Journal, 2021, 8(19): 14908–14919. doi: 10.1109/JIOT.2021.3072965
[9]	ZHANG Haijun, FENG Mengting, LONG Keping, et al. Energy efficient resource management in SWIPT enabled heterogeneous networks with NOMA[J]. IEEE Transactions on Wireless Communications, 2020, 19(2): 835–845. doi: 10.1109/TWC.2019.2948874
[10]	HUANG Jun, XING Congcong, and GUIZANI M. Power allocation for D2D communications with SWIPT[J]. IEEE Transactions on Wireless Communications, 2020, 19(4): 2308–2320. doi: 10.1109/TWC.2019.2963833
[11]	YANG H H, LEE J, and QUEK T Q S. Heterogeneous cellular network with energy harvesting-based D2D communication[J]. IEEE Transactions on Wireless Communications, 2016, 15(2): 1406–1419. doi: 10.1109/TWC.2015.2489651
[12]	ZHONG Xijian, GUO Yan, LI Ning, et al. Joint optimization of relay deployment, channel allocation, and relay assignment for UAVs-aided D2D networks[J]. IEEE/ACM Transactions on Networking, 2020, 28(2): 804–817. doi: 10.1109/TNET.2020.2970744
[13]	WANG Haichao, WANG Jinlong, DING Guoru, et al. Resource allocation for energy harvesting-powered D2D communication underlaying UAV-assisted networks[J]. IEEE Transactions on Green Communications and Networking, 2018, 2(1): 14–24. doi: 10.1109/TGCN.2017.2767203
[14]	JI Jiequ, ZHU Kun, NIYATO D, et al. Joint trajectory design and resource allocation for secure transmission in cache-enabled UAV-relaying networks with D2D communications[J]. IEEE Internet of Things Journal, 2021, 8(3): 1557–1571. doi: 10.1109/JIOT.2020.3013647
[15]	JI Jiequ, ZHU Kun, NIYATO D, et al. Probabilistic cache placement in UAV-assisted networks with D2D connections: Performance analysis and trajectory optimization[J]. IEEE Transactions on Communications, 2020, 68(10): 6331–6345. doi: 10.1109/TCOMM.2020.3006908
[16]	HUANG Wenhuan, YANG Zhaohui, PAN Cunhua, et al. Joint power, altitude, location and bandwidth optimization for UAV with underlaid D2D communications[J]. IEEE Wireless Communications Letters, 2019, 8(2): 524–527. doi: 10.1109/LWC.2018.2878706
[17]	WANG Baoji, ZHANG Rongqing, CHEN Chen, et al. Graph-based file dispatching protocol with D2D-enhanced UAV-NOMA communications in large-scale networks[J]. IEEE Internet of Things Journal, 2020, 7(9): 8615–8630. doi: 10.1109/JIOT.2020.2994549
[18]	HUQ K M S, MUMTAZ S, ZHOU Zhenyu, et al. Energy-efficiency maximization for D2D-enabled UAV-aided 5G networks[C]. 2020 IEEE International Conference on Communications, Dublin, Ireland, 2020: 1–6.
[19]	LIU Xiaonan, LI Zan, ZHAO Nan, et al. Transceiver design and multihop D2D for UAV IoT coverage in disasters[J]. IEEE Internet of Things Journal, 2019, 6(2): 1803–1815. doi: 10.1109/JIOT.2018.2877504
[20]	NGUYEN H T, TUAN H D, DUONG T Q, et al. Joint D2D assignment, bandwidth and power allocation in cognitive UAV-enabled networks[J]. IEEE Transactions on Cognitive Communications and Networking, 2020, 6(3): 1084–1095. doi: 10.1109/TCCN.2020.2969623
[21]	XIE Lifeng, XU Jie, and ZHANG Rui. Throughput maximization for UAV-enabled wireless powered communication networks[J]. IEEE Internet of Things Journal, 2019, 6(2): 1690–1703. doi: 10.1109/JIOT.2018.2875446
[22]	NGUYEN M N, NGUYEN L D, DUONG T Q, et al. Real-time optimal resource allocation for embedded UAV communication systems[J]. IEEE Wireless Communications Letters, 2019, 8(1): 225–228. doi: 10.1109/LWC.2018.2867775
[23]	ZHANG Jing, ZHANG Yanxia, XIANG Lin, et al. Robust energy-efficient transmission for wireless-powered D2D communication networks[J]. IEEE Transactions on Vehicular Technology, 2021, 70(8): 7951–7965. doi: 10.1109/TVT.2021.3095626
[24]	BOSHKOVSKA E, NG D W K, ZLATANOV N, et al. Robust resource allocation for MIMO wireless powered communication networks based on a non-linear EH model[J]. IEEE Transactions on Communications, 2017, 65(5): 1984–1999. doi: 10.1109/TCOMM.2017.2664860
[25]	徐勇军, 谷博文, 陈前斌, 等. 基于能效最大的无线供电反向散射网络资源分配算法[J]. 通信学报, 2020, 41(10): 202–210. doi: 10.11959/j.issn.1000-436x.2020132 XU Yongjun, GU Bowen, CHEN Qianbin, et al. Energy efficiency maximization resource allocation algorithm in wireless-powered backscatter communication network[J]. Journal on Communications, 2020, 41(10): 202–210. doi: 10.11959/j.issn.1000-436x.2020132
[26]	BOYD S and VANDENBERGHE L. Convex Optimization[M]. Cambridge, UK: Cambridge University Press, 2004.
[27]	XU Yongjun, HU Yuan, LI Guoquan, et al. Robust resource allocation for heterogeneous wireless network: A worst-case optimisation[J]. IET Communications, 2018, 12(9): 1064–1071. doi: 10.1049/iet-com.2017.0626

施引文献

期刊类型引用(8)

1.	杨子聪，李小勇，白雪茹. 基于即插即用2D-FISTA的高分辨ISAR成像方法. 空军工程大学学报. 2023(01): 38-44 . 百度学术
2.	杨磊，周弘昊，黄博，廖仙华，夏亚波. 机载合成孔径雷达高度计高程参数贝叶斯估计. 电子与信息学报. 2023(04): 1254-1264 . 本站查看
3.	舒婷，罗幼喜，李翰芳. 面板数据贝叶斯双惩罚分位回归方法研究. 广西师范大学学报(自然科学版). 2022(01): 150-165 . 百度学术
4.	潘小义，艾夏，刘佳琪，陈吉源，肖顺平. 自卫干扰位置可控式ISAR假目标合成方法. 太赫兹科学与电子信息学报. 2022(02): 140-146+153 . 百度学术
5.	王樾，白雪茹，周峰. 基于稀疏步进调频信号的低信噪比逆合成孔径雷达成像. 电子与信息学报. 2022(03): 1034-1043 . 本站查看
6.	杜文静，刘庆华，欧阳缮. 一种复块稀疏贝叶斯探地雷达成像算法. 现代雷达. 2022(05): 33-39 . 百度学术
7.	杨敏佳，白雪茹，刘士豪，曾磊，周峰. 基于高斯原型网络的小样本逆合成孔径雷达目标识别. 电子与信息学报. 2022(10): 3566-3573 . 本站查看
8.	侯颖妮，谢洁. 稀疏采样数据大转角高分辨ISAR成像技术研究. 雷达科学与技术. 2021(05): 604-609+624 . 百度学术

其他类型引用(10)

资源附件(0)

访问统计

图(9) / 表(1)

计量

文章访问数: 1118
HTML全文浏览量: 577
PDF下载量: 145
被引次数: 18

1. 引言
2. 回波信号模型
3. 贝叶斯Lasso
4. Gibbs采样算法
5. 实验验证
5.1 仿真数据验证
5.2 实测数据验证
5.3 地面动目标实验验证
5.4 相变图实验分析
6. 结束语

姓名
邮箱
手机号码
标题
留言内容
验证码

留言板

基于能量收集的UAV-D2D网络资源分配算法

doi: 10.11999/JEIT210491

通讯作者:
徐勇军　xuyj@cqupt.edu.cn

计量

Resource Rllocation for UAV-assisted D2D Communications with Energy Harvesting

1. 引言

2. 回波信号模型

3. 贝叶斯Lasso

4. Gibbs采样算法

5. 实验验证

5.1 仿真数据验证

5.2 实测数据验证

5.3 地面动目标实验验证

5.4 相变图实验分析

6. 结束语

期刊类型引用(8)

其他类型引用(10)

计量

目录

1. 引言

2. 回波信号模型

3. 贝叶斯Lasso

4. Gibbs采样算法

5. 实验验证

5.1 仿真数据验证

5.2 实测数据验证

5.3 地面动目标实验验证

5.4 相变图实验分析

6. 结束语

留言板

基于能量收集的UAV-D2D网络资源分配算法

doi: 10.11999/JEIT210491

通讯作者: 徐勇军 xuyj@cqupt.edu.cn

计量

出版历程

Resource Rllocation for UAV-assisted D2D Communications with Energy Harvesting

1. 引言

2. 回波信号模型

3. 贝叶斯Lasso

4. Gibbs采样算法

5. 实验验证

5.1 仿真数据验证

5.2 实测数据验证

5.3 地面动目标实验验证

5.4 相变图实验分析

6. 结束语

期刊类型引用(8)

其他类型引用(10)

计量

出版历程

目录

1. 引言

2. 回波信号模型

3. 贝叶斯Lasso

4. Gibbs采样算法

5. 实验验证

5.1 仿真数据验证

5.2 实测数据验证

5.3 地面动目标实验验证

5.4 相变图实验分析

6. 结束语

通讯作者:
徐勇军　xuyj@cqupt.edu.cn