1.能量获取D2D异构网络中高吞吐量资源分配方法，其特征在于，包括以下步骤：
步骤1：构建能量获取D2D异构网络中系统的数学模型；
步骤2：化简能量获取D2D异构网络的数学模型；
步骤3：基于变量松弛技术转换数学模型并证明其凹凸性；
步骤4：基于凸优化理论以及贪心策略求解数学模型；
步骤5：求满足约束条件的内点；
步骤6：基于可行方向探测与可变步长进行寻优。
2.根据权利要求1所述的能量获取D2D异构网络中高吞吐量资源分配方法，其特征在于，步骤1中构建能量获取D2D异构网络中系统的数学模型的步骤如下：
通信场景是以基站(Base Station，BS)为圆心位置，蜂窝用户(Cellular User 
Equipments，CUE)与终端直通通信(Device-to-Device，D2D)设备对随机分布于以BS为圆心r为半径的圆域内。系统中存在Nc个上行CUE用户，上行CUE链路的集合为C＝{1,2,...,Nc}，Nc个上行CUE链路对应有Nc个上行信道，信道集合为CH＝{1,2,...,Nc}，且按照编号顺序一一对应分配给CUE用户(各信道之间相互正交)。各CUE及使用的信道用用序号i∈C或i∈CH表示。系统中存在Nd对D2D设备，可以复用系统中上行信道进行D2D模式通信或者蜂窝模式通信。D2D链路的集合为D＝{1,2,...Nd}，D2D链路分别用序号j∈D表示。D2D设备通过从环境中收集能量为自己充能，能量获取速率服从泊松分布。在所考虑的通信系统模型中，系统传输总时槽数目为T，每时槽的时长为τt，各时槽用序号t表示，即t∈T。
D2D链路可以复用不同信道进行通信，信道分配变量为 当 时，表示t时槽D2D
链路j复用信道i进行通信， 为其它情况，公式如下：
D2D链路可以选择复用蜂窝信道进行D2D模式通信或者蜂窝模式通信， 表示模式选择， 时表示在时槽t内D2D链路j采用蜂窝模式进行通信， 时表示在时槽t内D2D链路j采用D2D模式进行通信， 的公式如下：
以基站为圆心建立笛卡尔直角坐标系，计算hi，hj，hi,j和dj，其中hi表示第i个上行CUE设备距离基站的距离，i∈C；hj表示第j对D2D设备对的发送设备距离基站的距离，j∈D；hi,j表示第j对D2D设备对的接收设备与CUE之间的距离；dj表示D2D设备对之间的距离。信号传输增益为距离倒数的3次方，设备受到的干扰为干扰源距离信号接收设备距离倒数的4次方。
在时槽t，蜂窝用户i的发送功率为 第j对D2D设备发送设备复用信道i的发送功率为第i个蜂窝用户使用i信道在t时槽内进行通信时的传输速率为 第i个蜂窝链路受到的干扰包括高斯白噪声N0B以及来自复用i信道进行通信的D2D发送设备带来的同频干扰根据香农公式可得其传输速率 的公式表示如下：
第j个D2D链路复用信道i在时槽t内进行蜂窝通信时的传输速率为 D2D链路j受到
的干扰包括高斯白噪声N0B以及蜂窝用户i所带来的同频干扰 其传输速率 的公
式如下：
第j个D2D链路复用信道i在时槽t内进行D2D通信时的传输速率为 D2D链路j受到的
干扰包括高斯白噪声N0B以及蜂窝用户i所带来的同频干扰 其传输速率 的公式
如下：
数学模型目标函数定义，我们考虑最大化T个时槽内全部D2D设备的总吞吐量，同一个D2D设备在同一个时槽内只能选择一种通信模式进行工作，目标是在满足CUE用户的QoS的同时，最大化T个时槽所有D2D设备的总的吞吐量，目标函数如下：
其中 表示t时槽所有D2D设备的总的吞吐量。
数学模型约束包括蜂窝用户的QoS约束；信道在被D2D设备复用时，同一时槽内信道i只能被一对D2D设备所复用；同一时槽内一对D2D设备也只复用一个信道i；D2D设备从环境中收集能量进行充能，每个时槽中D2D设备传输所消耗的能量始终不能大于D2D设备在本时槽起始时所剩余的能量；传输时间 不能大于时槽长度τt，数学模型的约束如下：
公式(7)表示为保障CUE用户的通信服务质量，蜂窝用户的最低传输速率为
公式(8)表示每个信道最多只能同时被一个D2D链路复用。
公式(9)表示每个D2D链路在同一时槽内只能复用一个蜂窝信道进行通信。
公式(10)表示能量获取约束，各个时槽内所消耗的总能量不得超过设备在本时槽已经获取的总能量，即在t时槽以及之前的时槽内所消耗的总能量不得超过在t时槽之前所总共获取的能量。D2D设备都能够从环境中获取能量为进行充能，用 表示D2D设备j在时间槽t内获取能量的速度，其值服从泊松分布， 代表D2D设备j在时间槽t内所获取到的总能量，则能量获取的公式如下；
表示D2D设备j在第1个时槽开始时所具有的初始能量。
公式(11)表示能量获取时间的约束，D2D设备因复用信道i进行工作，能量获取与数据传输不能同时进行，且在同一时槽t内需要满足数据传输时间与能量获取时间之和等于时槽长度τt。
公式(12)表示模型中优化变量的取值范围， 和 表示模式选择变量与
信道分配变量为二值变量； 和 表示CUE设备与D2D设备的传输功率为非负数；
表示数据传输时间为不大于时槽长度的非负数。
3.根据权利要求1所述的能量获取D2D异构网络中高吞吐量资源分配方法，其特征在于，步骤2中化简能量获取D2D异构网络中系统的数学模型的步骤为：
1)对于蜂窝用户的QoS约束如式(7)所示，即 其中 如式(3)所
示，令 则变量 需满足 通过观察目标函数可以
看出，目标函数是关于 的递减函数，若要使得目标函数取最大值，则变量 应该取最小值，即有
2)将变量 带入目标函式(6)，目标函数化简之后得到如下所示：
其中αi＝N0B(1+|hi|-1zi)，si,j＝|hi|3zi|hi,j|-4|hj|-4，ei,j＝|hi|-1|hj|-4zi+|hj|-3，ui,j＝|hi|-1zi|hj|-4，ki,j＝N0B(1+|hi|3zi|hi,j|-4)，mi,j＝|R2|-3+|hi|3zi|hi,j|-4|hj|-4，不难看出ei,j＞ui,j，mi,j＞si,j。
3)化简后的系统模型为：
4.根据权利要求1所述的能量获取D2D异构网络中高吞吐量资源分配方法，其特征在于，步骤3中基于变量松弛技术转换数学模型并证明其凹凸性的步骤为：
1)因为构建的数学模型是一个非线性混合整数约束优化问题，现利用变量松弛技术将二值变量 与 进行放缩为连续变量，即
2)根据复合函数凹凸性准则，对数函数凹凸性以及透视函数的性质，可以证明该转化的目标函数问题为凸函数，根据权利要求3中简化后的数学模型，简化后的目标函数如式(14)，设R(p,ρ,x)如下所示：
函数(16)是关于变量(p,ρ,x)的多元函数。先证明函数(16)是关于变量(p,x)的凹函数。若函数(16)在数学模型可行域内任意取不同两点：(p1,x1)、(p2,x2)，恒有式(17)成立，则函数(16)为凹函数，则目标函数(15)是关于变量 的凹函数。具体步骤为：
①将点(p1,x1)与(p2,x2)带入式(17)左边部分如下：
②将点(p1,x1)与(p2,x2)带入式(17)右边部分有
③将式(18)与式(19)相减并整理，则有如下式：
④根据式(21)证得函数(16)是关于变量(p,x)的凹函数，根据凸函数的可加性，可知目标函数(15)是关于变量 的凹函数。
⑤根据透视函数的性质，若f(x):Rn→R是凸函数，则
是凸函数。证明目标函数(15)
是关于变量 是凸函数的过程为：因目标函数(15)是关于变量 的凹函
数，又目标函数同时是关于变量 的函数，令 可知目标函数关
于变量 是凸的，根据透视函数的保凸性可知目标函数关于变量
是联合凸的。
5.根据权利要求1所述的能量获取D2D异构网络中高吞吐量资源分配方法，其特征在于，步骤4中基于凸优化理论以及贪心策略求解数学模型的步骤为：
1) 令 定义拉 格朗日 变量 矢量 其中
分别表示数学模型中的约束条件，如式(8)、(9)、(10)、(11)的拉格朗日乘
子。数学模型中的约束条件(12)在求解过程中会得到满足，则构造的拉格朗日函数为：
2)数学模型的对偶函数定义为 对偶问题为
3)对拉格朗日函数式(22)求偏导，如下所示:
4)根据KKT条件，联合约束条件(8)、(9)、(10)，求得 如下：
5)求解信道分配变量 与模式选择变量 因 与 是二值变量，经过变量松弛利用
KKT条件将问题初步求解后，采用贪心思想将松弛后的变量 与 进行求解使其二值化的值。利用式(24)所求得的传输功率 采用贪心思想，使信道i优先分配给能达到最大传输速率的D2D链路，记所有可供D2D设备复用的信道集合为SET_CH，所有未分配信道资源的D2D设备对的序号集合为SETρ_D，则信道分配变量 的二值化表达式如下：
6)求解模式选择变量 利用式(24)所求得的传输功率 采用贪心思想，令所有未选
择通信模式的D2D链路设为集合SETx_D，则变量 的二值化表达式如下：
式(27)中的 表示当前链路j复用了信道i进行工作。
7)求解传输时间变量 已求得传输功率变量 信道分配变量 与模式选择变量
的闭合表达式，将拉格朗日函数L(X,λ)对变量 求偏导，如下所示:
由KKT条件，可得出传输时间式变量 如下所示：
其中
6.根据权利要求1所述的能量获取D2D异构网络中高吞吐量资源分配方法，其特征在于，步骤5求满足约束条件的内点的具体步骤为：
1)根据步骤4中利用KKT条件解得的各变量拉格朗日对偶解，根据给定的 根
据式(24)、(25)、(27)和(30)依次计算变量
2)由于拉格朗日对偶问题关于对偶变量 是凸的，而约束条件(8)、(9)、(12)在
计算过程中将得到满足，我们以约束条件(10)、(11)来作为边界来确定 是否为内点，在各个时间槽中，当前确定的各个D2D链路的信道分配与模式选择时，由于凸函数的可加性，原问题可以按照D2D链路分解成Nd个子凸问题，令各D2D链路传输能量约束条件(10)用 表示，如下所示：
为 的拉格朗日乘子，各D2D链路传输时间约束条件(11)用 表示，
为 的拉格朗日乘子。βj为第j个链路寻找内点的步长系数，取值为小于1的正数，ε1为寻找内点步长阈值。内点必须满足约束条件(31)与(32)，即 与 必须非负。
3)确定内点迭代方向。迭代方向一定得使 与 增大，对所有D2D链路，如果D2D链路j有分配信道，则根据式(31)和式(32)计算 与 如果 或者 则说明
不满足约束条件，需要计算新的拉格朗日乘子向量Λm,k+1,m∈{1,2,3,4}，如下所示：
其中Λm＝1,k+1表示拉格朗日乘子 和 都增大，Λm＝2,k+1表示 增大， 减少，
Λm＝3,k+1表示 减少， 增大，Λm＝4,k+1表示 和 都减少。由计算得出的Λm,k+1,m∈{1,
2,3,4}，根据式(31)和式(32)再计算 其中 与 为四
个不同方向的值。
4)根据计算得出的 按照m的取值顺序依次判断式(34)是
否成立。
若第m个成立，则令 若m＝1,2,3,4时(34)均不成立，则令 并
且更新第j个链路寻找内点的步长系数，如下所示：
βj,k+1＝βj,k×βj,k   (35)
5)根据式(24)、(30)和(13)重新计算变量 如果j＜Nd，跳转到3)开始计算一
下个D2D链路j+1，否则跳转到6)。
6)当j＝Nd，判断是否所有的D2D链路 与 值都满足 若
所有链路j都满足约束条件，则说明已找到内点，跳转至步骤6开始进行寻找最优值；若仍然有不满足约束条件(31)与(32)的D2D链路，则计算不满足条件的链路的最小的步长系数如果最小步长系数 则认为当前信道分配方式在对应链路j上
不存在有效内点，且令对应链路j的功率 则根据式(25)计算信道分配变量，重新分配信道，根据式(26)计算模式选择变量，重新选择模式，根据式(27)计算传输时间，重置步长系数βj，并跳转到3)。
7.根据权利要求1所述的能量获取D2D异构网络中高吞吐量资源分配方法，其特征在于，步骤6基于可行方向探测与可变步长进行寻优的具体步骤为：
1)根据权利要求6中已求得的内点开始计算问题的最优值，具体包括迭代方向采用可行方向探测法，及可变步长方式进行计算，以使得算法能够快速达到收敛。因此时已经找到内点，目标函数为凹函数，最大值吞吐量，则使目标函数增大的方向即为可行方向。因为在寻优过程中，由于信道分配方式或模式选择的改变，可能导致目标值在迭代过程中有短暂缩小的情况，设寻优迭代的收敛条件为式 其中0＜ε2＜1为目标函数收敛阈
值。下面依次对不同D2D链路j迭代寻找最优解，j＝1。
2)对链路j，令 各链路j的迭代步长为stepj,n，步长更新系数设为α，α
为小于1的正数。计算新的拉格朗日乘子向量Λn+1,m,m∈{1,2,3,4}，如下所示：
并根据式(15)和Λn的值计算t时槽第n次迭代的目标值 并根据式(15)和Λn+1,m,m∈{1,2,3,4}的值探测计算t时槽第n+1次迭代的目标值
3)比较 与 中的最大值，若最大值为 则当前链路的迭代方向探
测失败，缩小当前链路的步长stepj,n+1＝(1-α)×stepj,n，令 若最大
值为 则当前链路迭代方向探测成功，增加当前链路的步长stepj,n+1＝(1+α)×
stepj,n，更新拉格朗日乘子 根据式(24)、(30)和(13)重新计算变量
如果j小于Nd，则j＝j+1，并跳转到2)，否则跳转4)。
4)根据式(24)、(25)、(27)、(30)和(13)依次计算变量 和
5)若本轮迭代中所有链路在当前步长下均探测不成功，则将各链路步长同时乘以α，j＝1，并跳转到2)继续迭代，重新开始探测各链路迭代方向，否则，说明存在有效探测，则根据式 (31) 和式(32) 计算 与 并判断是否所有D2D链路都满足
如果有不满足约束条件(31)和(32)的链路，则说明探测的点不在
内点，跳转到步骤5重新计算内点，否则判断是否达到收敛误差，即判断式(37)是否成立，若成立，说明已找到最优值并继续迭代下一时槽(t＝t+1)，若判断式(37)不成立，说明未达到收敛误差要求，则跳转到2)。