CN107509162B - 新型蜂窝下行通信能效优化方法 - Google Patents

Abstract

本发明公开了一种新型蜂窝下行通信能效优化方法，包括基站分配中继信道；基站辅助链路接收端节点调节发射端功率；确定数据接收模式。本发明方法能够将基站运行的集中式方式改进为基站辅助的分布式方式，有效减轻了基站负担，降低对用户设备配置的要求，避免发生一跳通信范围内的同频干扰，优化链路能效，降低寻找合适中继与构建传输路径的开销，保证了所建立的传输路径具有更好的持续服务能力，并且改善了通信质量。

Description

新型蜂窝下行通信能效优化方法

技术领域

本发明具体涉及一种新型蜂窝下行通信能效优化方法。

背景技术

随着新媒体、增强现实、虚拟现实、移动云等移动互联网业务的普及，万物互联的内在需求与物联网市场的高速发展，下一代移动通信网必须要注重能量效率和频谱效率的提升。D2D通信模式被广泛认同为提升蜂窝网络能量与频谱效率的一种可行技术，其主要作用表现如下：1)地理位置上邻近的通信对之间能够通过D2D模式直接传输数据，减轻了基站的邻近区域数据中转量，从而提升了为远距离通信对提供流量中继服务的能力；2)远离基站的蜂窝边缘节点或非视距范围内节点可以通过D2D中继节点改善与基站通信的链路质量，从而提升蜂窝上下行吞吐量。本发明内容属于后者范畴，且仅考虑单一蜂窝内一定资源约束下的下行通信能效优化问题。

D2D通信可分为带外D2D和带内D2D。带外D2D通信模式使用非授权的免费频谱，如WiFi或蓝牙技术等所采用的频段。因此，存在如下问题：1)非授权的免费频谱不受蜂窝基站控制，可能会导致很差的服务质量，2)至少具有两类无线接口(例如，LTE和WiFi)的设备才可以使用带外D2D通信，并涉及不同类型无线接口之间的电源管理、不同无线通信协议数据包的解码和编码、更复杂的D2D用户之间的发现和连接设置等问题，因此，对用户设备配置要求较高。带内D2D通信使用蜂窝运营商的授权频谱，又可以分为专用模式(dedicatedmode，在一些文献中也被称为overlay mode或者orthogonal mode)和共享模式(sharedmode，在一些文献中也被称为underlay mode或者non-orthogonal mode)。前者将授权频谱划分为两部分，一部分专供D2D用户使用，一部分给蜂窝用户使用，从而避免了相互干扰。后者允许D2D通信与蜂窝通信共享频谱资源，因此，存在同频干扰问题，需要有效的方法来控制干扰的影响。本专利采用的是共享模式。

当给定了一个蜂窝内的信道资源总量，并限定了基站的最大发射功率，能够为随机发生的大量用户请求提供的下行吞吐量的上限在哪里？通常这个问题属于NP难范畴，存在不同的近似求解方法，有文献提出的一种穷尽搜索的贪婪算法，其主要思路是由蜂窝基站为蜂窝内获得蜂窝信道的每个下行流量的接收节点确定是直接从基站接收？还是选择一个D2D中继节点帮它从基站转发？尽管该类方法实现简单，但存在的一些问题有待解决：1)由基站运行全网穷尽搜索的贪婪算法，导致基站负担过重；2)采用带外D2D通信对用户设备配置具有更高要求；3)规定每个下行数据接收节点至多选择一个D2D中继节点可能无法满足远离基站的蜂窝边缘节点或非视距范围内节点的下行吞吐量需求；4)未对发射节点的功率进行调节，蜂窝通信能效仍有提升空间；5)未约束D2D中继节点的能量储备水平，难以保障所建立的传输路径的持续服务能力。

发明内容

本发明的目的在于提供一种能够将实际执行的计算与通信开销在基站和数据的接收端节点之间进行合理分担，从而实现高效经济通信目标的新型蜂窝下行通行能效优化方法。

本发明提供的这种新型蜂窝下行通行能效优化方法，包括如下步骤：

基站分配中继信道的步骤：具体包括蜂窝中继节点信道分配的步骤和D2D中继链路信道分配的步骤；

基站辅助链路接收端节点调节发射端功率的步骤：具体包括调节基站对接收节点的发射功率的步骤、调节基站对D2D中继节点的发射功率的步骤、调节D2D中继链路发射端的功率的步骤和基站辅助UE进行功率调节的步骤，所述的UE为用户设备；

确定数据接收模式的步骤：具体包括接收节点确定数据接收模式的步骤和D2D中继节点确定数据接收模式的步骤。

所述的蜂窝中继节点信道分配的步骤，具体为所分配的信道用于从基站直接接收数据，具体包括如下步骤：

步骤1.0：将L^r _M的每个元素初始化为0；创建集合U_c和U_r，并初始化为空；进入步骤1.1；所述L^r _M为基站与用户之间的中继信道分配关系数组，其任一元素l^r _j的值为k且k∈{1,…,N}，则表示从基站到用户j之间的链路被分配了中继信道k，否则表示没有分配中继信道；j∈{1,…,M}；

步骤1.1：针对L^c _M,N中的每个元素l^c _j,k，若存在l^c _j,k为“真”，则将j节点加入到U_c集合中，否则加入到U_r集合中；进入步骤1.2；L^c _M,N表示用户直接与基站通信的信道分配关系矩阵，其任一元素l^c _j,k的值若为“真”，则表示用户j被分配了信道k用于与基站直接通信；若为“假”，则表示没有分配信道；

步骤1.2：针对U_r中的每个距离小于R_th且大于r_th的元素，用j表示，创建集合D并初始化为空；进入步骤1.3；

步骤1.3：针对U_r中的每个距离小于R_th且大于r_th的非j元素，用k表示，计算k与j之间关于基站的夹角

若

小于基站天线发射角

并且l^r _k大于0，则将l^r _k加入到D中；进入步骤1.4；

步骤1.4：针对U_c集合中的每个元素，用i表示，计算i与j之间关于基站的夹角

若

小于基站天线发射角

并且存在满足l^c _i,k为“真”的信道k，则将该信道k加入到D中；进入步骤1.5；

步骤1.5：创建集合E，其中的元素由集合C中所有不在集合D中的元素组成，从集合E中随机取出一个元素k，并将k赋值给l^r _j；重复步骤1.2，直至遍历完U_r中的每个元素。

所述的D2D中继链路信道分配的步骤，为分配的D2D中继链路信道用于接收D2D中继链路发射端发射的数据；具体包括如下步骤：

步骤2.0：将L^r _M,M中的每个元素初始化为0；L^r _M,M表示用户之间的中继信道分配关系矩阵，其任一元素l^r _i,j的值若为k，k∈{1,…,N}，表示从用户i到用户j之间的链路被分配了中继信道k，否则表示没有分配中继信道；

步骤2.1：针对U_r中的每个距离大于r_th的元素，创建集合D并初始化为空；进入步骤2.2。

步骤2.2：若l^r _j大于0，则将l^r _j加入到D中；进入步骤2.3。

步骤2.3：针对U_c中的每个元素，以i表示，若i是j的邻居集合V_j中元素，则在L^c _M,N中找到i与基站直接通信时使用的信道k，并将k加入到D中；进入步骤2.4；

步骤2.4：针对U_r中的每个距离小于R_th且大于r_th的元素，以i表示，若i是j的邻居集合V_j中元素且l^r _i大于0，则将l^r _i加入到D中；进入步骤2.5；

步骤2.5：针对j的邻居集合V_j内的每个元素，以i表示，若l^r _i,j大于0，则将l^r _i,j加入到D中；进入步骤2.6；

步骤2.6：针对j的邻居集合V_j内的每个元素，以k表示，针对k的邻居集合V_k内的每个元素，以l表示，若l^r _l,k大于0，则将l^r _l,k加入到D中；进入步骤2.7；

步骤2.7：针对j的邻居集合V_j内的每个元素，以i表示，若i到基站的距离大于j到基站的距离且l^r _j,i为0，则创建集合D`并初始化为空；进入步骤2.8；

步骤2.8：针对U_c中的每个元素，以m表示，计算m与i之间关于基站的夹角

若

小于基站天线发射角

则在L^c _M,N找到m与基站直接通信时候使用的信道k，并将k加入到D`中；进入步骤2.9；

步骤2.9：针对U_r中的每个距离小于R_th且大于r_th的元素，以m表示，计算m与i之间关于基站的夹角

若

小于基站天线发射角

并且l^r _m大于0，则将l^r _m加入到D`中；进入步骤2.10；

步骤2.10：针对i的邻居V_i内的每个元素，以k表示，针对k的邻居V_k内的每个元素，以m表示，若l^r _k,m大于0，则将l^r _k,m加入到D`中；进入步骤2.11；

步骤2.11：创建集合E，其中的元素由集合C中所有不在集合D和D`中的元素组成，从集合E中随机取出一个元素k，或从集合E中取出第一个元素k，并将k赋值给l^r _j,i；重复步骤2.7，直至遍历完V_j内的每个元素；进入步骤2.12；

步骤2.12：重复步骤2.1，直至遍历完U_r中的每个元素。

所述的调节基站对接收节点的发射功率的步骤，具体包括如下步骤：

步骤0.0：事先获得全网所有节点的最大发射功率参数；

步骤1.0：针对从基站直接接收自己数据的节点，以j表示，初始化p_0j为0，若从基站接收到T_s ^max,P_s ^max,T_s ^half,P_s ^half，则判断从基站到自身的距离d_0,j是否小于R_th，如果小于则进入步骤1.1，否则不进行功率调节；所述的P_s ^max和P_s ^half表示包括潜在D2D链路和直接通信链路在内的整个网络潜在通信链路的发射端的总功耗最大值和对半值，而T_s ^max和T_s ^half表示整个网络中潜在通信链路的吞吐能力的最大值和对半值；

步骤1.1：计算基站对自身采用发射功率p₀ ^half下的效益μ_j ^half和信噪比，以及最大功率下p₀ ^max的效益μ_j ^max,若μ_j ^half大于等于μ_j ^max并且在一半功率情况下的信噪比大于等于预设的信噪比阈值γ_th，将p_0j赋值为p₀ ^half，否则将p_0j赋值为p₀ ^max，并将p_0j发送给基站，进入步骤1.2；p₀ ^half为最大发射功率p₀ ^max的一半；

步骤1.2：若从基站接收到T_s和P_s的更新信息，则进入步骤1.3；若从基站接收到终止信息，则结束运行；若未接收到任何信息，则继续等待；

步骤1.3：若j与基站的距离d_0,j小于R_th，则计算基站采用p_0j减去步长ε后的值作为发射功率时的效益μ`<sub>j</sub>和信噪比，以及采用p<sub>0j</sub>作为发射功率时的效益μ<sub>j</sub>，若μ`_j大于μ_j且信噪比大于预设的信噪比阈值γ_th，则将p_0j的值更新为减去了步长ε后的值，并将更新后的p_0j发送给基站，并重复步骤1.2。

所述的调节基站对D2D中继节点的发射功率的步骤，具体包括如下步骤：

步骤0.0：事先获得全网所有节点的最大发射功率参数；

步骤2.0：针对从基站接收其他节点数据的中继节点，以j表示。初始化p_0j为0，若从基站接收到T_s ^max,P_s ^max,T_s ^half,P_s ^half，则判断从基站到自身的距离d_0,j是否小于R_th且大于r_th，若满足条件则进入步骤2.1，否则不进行功率调节；

步骤2.1：计算基站对自身采用发射功率p₀ ^half下的效益μ_j ^half和信噪比，以及最大功率下p₀ ^max的效益μ_j ^max,若μ_j ^half大于等于μ_j ^max并且在一半功率情况下的信噪比大于等于预设的信噪比阈值γ_th，将p_0j赋值为p₀ ^half，否则将p_0j赋值为p₀ ^max，并将p_0j发送给基站，进入步骤2.2；

步骤2.2：若从基站接收到T_s和P_s的更新信息，则进入步骤2.3；若从基站接收到终止信息，则结束运行；若未接收到任何信息，则继续等待；

步骤2.3：若j与基站的距离d_0,j小于R_th且大于r_th，则计算基站采用p_0j减去步长ε后的值作为发射功率时的效益μ`<sub>j</sub>和信噪比，以及采用p<sub>0j</sub>作为发射功率时的效益μ<sub>j</sub>，若μ`_j大于μ_j且信噪比大于预设的信噪比阈值γ_th，则将p_0j的值更新为减去了步长ε后的值，并将更新后的p_0j发送给基站，并重复步骤2.2。

所述的调节D2D中继链路发射端的功率的步骤，具体包括如下步骤：

步骤0.0：事先获得全网所有节点的最大发射功率参数；

步骤3.0：针对任一D2D中继节点，以j表示，若从基站接收到T_s ^max,P_s ^max,T_s ^half,P_s ^half，则针对自己邻居集合V_j中所有已分配中继信道的邻居，以i表示，初始化p_ij为0；进入步骤3.1；

步骤3.1：计算i对j采用发射功率p_i ^half时的效益μ_j ^half和信噪比，以及最大功率下p_i ^max的效益μ_j ^max，若μ_j ^half大于等于μ_j ^max并且在一半功率情况下的信噪比大于等于预设的信噪比阈值γ_th，则将p_ij赋值为p_i ^half，否则将p_ij赋值为p_i ^max；进入步骤3.2；

步骤3.2：若邻居集合V_j中的元素未遍历完，则返回步骤3.1；否则将发射功率集合{p_ij|i∈V_j}中的所有功率值发送给基站，然后进入步骤3.3；

步骤3.3：若从基站接收到T_s和P_s的更新信息，则进入步骤3.4；若从基站接收到终止信息，则结束运行；若未接收到任何信息，则继续等待；

步骤3.4：针对自己邻居集合V_j中所有已分配中继信道的邻居，以i表示，计算i采用p_ij减去步长ε后的值作为发射功率时的效益μ`<sub>j</sub>和信噪比，以及采用p<sub>ij</sub>作为发射功率时的效益μ<sub>j</sub>，若μ`_j大于μ_j且信噪比大于预设的信噪比阈值γ_th，则将p_ij的值更新为减去了步长ε后的值；在遍历完邻居集合V_j中的所有元素后，将发射功率集合{p_ij|i∈V_j}中的所有功率值发送给基站，并重复步骤3.3。

所述的基站辅助UE进行功率调节的步骤，具体包括如下步骤：

步骤0.0：基站事先获得全网所有节点的最大发射功率参数；

步骤4.0：计算所有发射端都采用最大功率时的吞吐量之和T_s ^max和功率之和P_s ^max，以及所有发射端都采用最大功率的一半时的吞吐量之和T_s ^half和功率之和P_s ^half，然后向所有节点广播T_s ^max,P_s ^max,T_s ^half,P_s ^half；进入步骤4.1；

步骤4.1：初始化功率改变标志flag为“假”，将定时器t_Δ设置为Δ，在定时器t_Δ没有超时的情况下，若基站收到来自任一UE的功率报告，则检验功率有无发生变化，若发生了变化，则将功率改变标志flag置为“真”；进入步骤4.2；

步骤4.2：若功率改变标志flag为“真”，则更新T_s和P_s值，并将更新的T_s和P_s广播给所有的UE，并重复步骤4.1；否则，向所有节点广播功率调节终止包，结束算法的运行。

所述的接收节点确定数据接收模式的步骤，具体包括如下步骤：

步骤1.0：针对任一获得蜂窝信道的从基站接收自己数据的节点，以UE i表示，若收到基站的功率调节终止信息，则将α_0j、α_ji、T^ji _in、T^ji _out初始化为0,将集合D_i初始化为“空”，将标志位flag初始化为“假”；进入步骤1.1；

步骤1.1：针对i的邻居集合V_i中的每个邻居，以j表示；若j到基站的距离d_0,j小于R_th，则计算T^ji _in，否则计算T^ji _out；对所有的T^ji _in和T^ji _out分别按从大到小的顺序排序；进入步骤1.2；

步骤1.2：若标志位flag为“假”，则找到当前最大的T^ji _in并赋值给T^ji _max，然后将当前最大的T^ji _in值置为0，否则，找到当前最大的T^ji _out并赋值给T^ji _max，然后将当前最大的T^ji _out值置为0；进入步骤1.3；

步骤1.3：若T^ji _max大于0，则说明选择中继比直接通信效果更好，因此，继续判断路径延时与中继能量约束是否满足要求；若d_ji小于延迟阈值d_th并且j的剩余能量e_j大于能量阈值e_th，进入步骤1.4，否则，返回步骤1.2；若T^ji _max不大于0且标志位flag为“假”，则将标志位flag置为“真”，返回步骤1.2，否则，说明当前节点找不到合适的中继，结束运行；

步骤1.4：给邻居j发送中继请求包；将定时器t_Δ设置为Δ，在定时器t_Δ没有超时的情况下，若收到了j的响应中继请求的响应包，则将α_ji和α_0j分别置为1，将j加入集合给D_i，退出等待接收响应包的操作；若定时器t_Δ超时且未选中中继(即α_ji为0)，则返回步骤1.2；若选中了中继(即α_ji为1)，则结束运行。

所述的D2D中继节点确定数据接收模式的步骤，具体包括如下步骤：

步骤2.0：针对任一D2D中继节点，以j表示，若接收到任一邻居i的中继请求包，查看自身是否已充当了中继，若没有为任何邻居充当中继，则向i发送中继请求的响应包，否则，不发送任何响应；并进入步骤2.1；

步骤2.1：将α_0k、α_kj、T^kji _in、T^kji _out分别初始化为0，将集合D_j初始化为“空”，将标志位flag初始化为“假”；进入步骤2.2；

步骤2.2：针对j的邻居集合V_j中的每个邻居，以k表示；若k到基站的距离d_0,k小于R_th且大于r_th，则计算T^kji _in，否则，若k到基站的距离d_0,k不小于R_th，则计算T^kji _out；对所有的T^kji _in和T^kji _out分别按从大到小的顺序排序；进入步骤2.3；

步骤2.3：若标志位flag为“假”，则找到当前最大的T^kji _in并赋值给T^kji _max，然后将当前最大的T^kji _in值置为0，否则，找到当前最大的T^kji _out并赋值给T^kji _max，然后将当前最大的T^kji _out值置为0；进入步骤2.4；

步骤2.4：若T^kji _max大于0，则说明选择中继比直接通信效果更好，因此，继续判断路径延时与中继能量约束是否满足要求；若d_kji小于延迟阈值d_th并且k的剩余能量e_k大于能量阈值e_th，进入步骤2.5，否则，返回步骤2.3；若T^kji _max不大于0且标志位flag为“假”，则将标志位flag置为“真”，返回步骤2.3，否则，说明当前节点找不到合适的中继，结束运行；

步骤2.5：给邻居k发送中继请求包；将定时器t_Δ设置为Δ，在定时器t_Δ没有超时的情况下，若收到了k的中继响应包，则将α_kj和α_0k分别置为1，将k加入集合给D_j，不再接收其它响应包；若定时器t_Δ超时且未选中中继，则返回步骤2.3；若选中了中继，则结束运行。

本发明提供的这种新型蜂窝下行通行能效优化方法，在保留了现有方法实现的简单有效特征的同时，将基站运行的集中式方式改进为基站辅助的分布式方式，有效减轻了基站负担；采用带内共享模式的D2D通信以提升频谱利用率、避免带外D2D通信模式中干扰的不可控性，同时也可以降低对用户设备配置(尤其是无线接口类型)的要求；由基站为潜在的D2D中继节点集中规划蜂窝信道的复用规则，尽可能确保被分配者复用其一跳通信范围以外的信道资源，从而尽可能避免发生一跳通信范围内的同频干扰；通过引入博弈机制，依据提出的新效用函数，潜在的D2D中继链路发射端独立调节其发射功率，以达到本链路能效的优化目标；采用以接收节点为中心的一跳通信范围内的局部搜索，以确定接收节点是采用蜂窝链路从基站直接接收数据还是将蜂窝链路与D2D链路串联成数据传输路径来接收蜂窝下行数据，相比现有技术的全局部搜索，降低寻找合适中继与构建传输路径的开销；在选用D2D中继时，通过约束D2D中继节点的能量储备水平以保障所建立的传输路径具有更好的持续服务能力；最后，本发明规定D2D中继节点必须执行接收节点同样的数据接收模式决策过程，从而有可能为接收节点扩展接收数据的路径跳数、进一步改善远离基站的蜂窝边缘节点或非视距范围内节点与基站之间的通信质量。

附图说明

图1为本发明方法的方法流程图。

图2为蜂窝下行平均传输路径能量效率随固定区域中节点数量的变化趋势示意图。

图3为蜂窝下行平均传输路径持续服务能力随固定区域中节点数量的变化趋势示意图。

具体实施方式

如图1所示为本发明方法的方法流程图：

本发明中涉及的计算公式与相关参数阐述如下：

对任一用户设备(User Equipment,UE)(如i)，若它采用从基站直接接收数据的通信方式，则可获得的吞吐量可使用香农公式估算，具体表示如下：

T_0i＝b_i·log₂(1+γ_i)(1)

在公式(1)中，T_0i和γ_i分别是从基站到UE i的通信链路上的吞吐量和信号干扰噪声比(Signal to Interference Noise Ratio,SINR)；b_i是该链路上的信道带宽。γ_i可由公式(2)估算。

在公式(2)中，N_i表示UE i感知到的环境噪声功率；F_i表示当基站通过蜂窝信道向UE i发送数据时，UE i感知到的来自其他链路的干扰，p_i是基站向UE i发射数据时采用的发射功率(通常取10瓦)；g_i表示从基站到UE i的链路上的信道衰减系数；p_i与g_i的乘积表示UE i感知到的接收功率强度；g_i的值涉及路径损耗、多径衰落、阴影衰落等因素，通常较难获得。既然p_i与g_i的乘积表示UE i感知到的接收功率强度，我们可采取如下近似公式直接得到接收功率。

在公式(3)和(4)中，G_t和G_r分别是发射天线和接收天线增益；h_t和h_r分别是发射天线高度和接收天线高度；λ和L分别是载波信号波长和***损耗因子；d_u,v表示发射节点u与接收节点v之间的距离；当发射节点u向接收节点v发射数据时，p^t _u,v和p^r _u,v分别表示发射节点u的发射功率和接收节点v的接收功率。依据任一对通信节点间的距离与节点交叉距离d_crossover之间的关系，确定使用公式(3)和(4)中的哪一个。对任一链路u→v，若其链路长度小于d_crossover，则使用公式(3)，否则使用公式(4)。交叉距离d_crossover由公式(5)得到。

任一UE(如i)通过中继UE(如j)的数据转发来提升下行吞吐量，UE j从基站获得的吞吐量可由公式(6)估算。

T_0j＝b_j·log₂(1+γ_j) (6)

在公式(6)中，T_0j和γ_j分别表示从基站到UE j的直接链路上的吞吐量和信号干扰噪声比。

当中继UE j将数据转发给UE i时，其D2D链路的吞吐量可由公式(7)估算。

T_ji＝b_ji·log₂(1+γ_ji) (7)

在公式(7)中，T_ji和γ_ji分别是从UE j到UE i的D2D链路上的吞吐量和信号干扰噪声比；b_ji是该D2D链路上的信道带宽。γ_ji由公式(8)估算。

在公式(8)中，p_ji是UE j向UE i发射数据时采用的发射功率；g_ji表示从UE j到UEi的链路上的信道衰减系数；既然p_ji与g_ji的乘积表示UE i感知到的接收功率强度，我们可采用公式(3)或(4)近似得到；F_i表示当UE i使用某个蜂窝信道接收数据时，UE i感知到的来自其他同频链路的干扰，包括来自基站的干扰，以及来自邻居范围内采用同频信道的UE的干扰。由公式(9)估算。

在公式(9)中，I_i表示UE i所有的干扰源集合，而UE k表示该集合中的一个元素；p_k是发射端k采用的发射功率；g_ki表示从发射端k到UE i的链路上的信道衰减系数；p_k与g_ki的乘积表示UE i感知到的干扰功率强度，同样我们可采用公式(3)或(4)近似得到之。

UE i通过信道状态更优的中继UE j进行数据转发，达到优化该通信路径下行吞吐量的目的，综合公式(6)和(7)，其吞吐量值可由公式(10)估算。

为进一步优化下行链路吞吐量，UE i通过信道状态更优的中继UE k和UE j依次进行数据转发，同理，该通信路径下行吞吐量值可由公式(11)估算。

V_i,V_j分别表示UE i与UE j的邻居节点集合，对UE i来说，V_i集合内的任一节点以给定的发射功率(如p_ue，通常取值为0.1瓦)向UE i发送数据都应保证接收端UE i的误码率不高于设定的BER值(可取10^-10)。V_i集合内节点的分布范围(即接收节点UE i的覆盖半径D_i)可根据公式(12)估算得到。

在公式(12)中，γ_th是与设定的BER值对应的信号干扰噪声比值，其值可以通过公式(13)得到。

γ_th＝-2lnBER (13)

由于任一接收节点(如UE i)无法获得实际的干扰节点数目、分布、干扰节点采用的发射功率等估算干扰大小的相关信息，因此，F_i是一个不确定的值。在本发明中，通过合理的信道分配，尽量避免同频信道干扰的存在，因此，在计算覆盖半径D_i时，我们将F_i置为0。同理，基站以最大功率(可取10瓦)发射且满足最低误码率情况下的覆盖半径R_th的计算公式如下。

在(14)中，p₀ ^max表示基站的最大发射功率，而N_avg表示蜂窝内的平均噪声功率。该R_th值说明，在不考虑干扰的情况下，若某节点到基站的距离小于R_th，在基站以最大功率对其传输数据，必然满足最低误码率要求，因此，也可能存在调节功率的余地。若该范围内(即以基站为中心的半径R_th范围内)的节点未获得蜂窝信道，则无法作为接收节点，但可以充当其他接收节点(通常位于以基站为中心的半径R_th范围外)的中继，为其转发来自基站的数据，因此，需要为这样的中继节点分配中继信道。为了避免这样的中继离基站太近，本发明限制了其距离基站的距离大于r_th(可取值为120米)才给其分配中继信道。

本发明中假定基站采用方向天线发射信息，因此，处于基站方向天线角度范围内且与基站使用同频信道的链路的接收端才可能受到基站发射信号的干扰。在基站、信道的复用者(如UE_j)、被复用者(如UE_k)所形成的三角形区域中，以基站为中心的夹角(如

)必须大于基站向信道复用者发送数据的发射角

(可取60度)，否则，为了避免同频干扰，该信道不能被该复用者复用。若已知任意两个节点(如UE_j、UE_k)以及基站的平面坐标(x_j,y_j)、(x_k,y_k)和(x₀,y₀)，则以基站为中心的夹角

可以通过下列一组公式估算。

在本发明关注的问题场景(即只考虑单一蜂窝内的干扰，或假定相邻蜂窝间的同频干扰能得到很好管控)下，对一个蜂窝信道来说，能够方便地确定保障接收端接收误码率不低于接收比特误码率阈值(该阈值间接反映了用户的网络应用体验)的发射功率，只要知道用户设备的接收比特误码率阈值(如10^-10)，接收端的环境噪声功率，信道带宽、信道衰减系数(基于接收端测量的CSI信息进行量化得到)即可。若接收比特误码率阈值表示为BER，结合公式(13)，相应的吞吐量和发射功率可以通过公式(16)，(17)估算得到。

若从基站到任一接收端i之间的蜂窝信道的发射端功率不小于按式(17)所得的pⁱ _th，则可确保接收端接收误码率不低于接收比特误码率阈值BER。若发射端的最大发射功率达不到pⁱ _th，则最好采用最大发射功率以尽量保障接收质量。

对任一采用蜂窝信道的D2D通信链路j→i来说，为确保接收端i的比特误码率不高于BER，发射端j的发射功率不应低于

其值可以通过下式估算。

实际上，接收端i的潜在干扰源的发射功率可能会随时被调整，因此，F_i的值只能基于当前所已知的信息做估算，而获得的

也可能再被进一步调整。由于发射端j也可能是其它同频链路接收端的潜在干扰源，因此，它的功率调整也可能触发其它同频链路发射端作出新的功率调整。当所有同频链路的发射端都不再调整其发射功率时，整个***才能稳定。因此，本发明引入潜博弈(ordinal potential game(OPG))理论来建模解决该功率调整问题，包括效用函数(utility function)与博弈决策算法(game decision processalgorithm)的设计。在该模型中，规定博弈的参与者(players)为所有的潜在通信链路(即被分配了信道的通信对)，每个参与者拥有一个评估自身收益的效用函数和一组行为规则集。对任一D2D通信链路j→i或直接通信链路0→i的接收端i,其效用函数表示如下。

在(19)，右边第一项表示整个网络通信潜能的改善给个体(例如，D2D通信链路j→i或直接通信链路0→i，这里0表示基站)带来的效用，而右边第二项表示个体的效用，在这里，效用的具体含义是能量效率，即单位能耗所能传输的数据量。

U＝{1,2,…M}；μ_i(P)表示链路j→i(或链路0→i)的发射端j(或基站)使用功率p_ji(或p_0i)发射时获得的效用；w_g和w_l是权值系数，两者取值范围都为(0,1)，且两者之和为1；P_s表示整个网络潜在通信链路(包括潜在D2D链路和直接通信链路)的发射端的总功耗，而T_s为表示整个网络中潜在通信链路的吞吐能力，其函数表达式如下。

参与者(即某条通信链路)的行为规则集可具体化为它对自身发射端功率的调节级别(即预先定义的一个步长)，每次决策时，通过使用当前发射功率减去一个步长得到调整后的功率，再依据效用函数判断，若更新后的功率能带来更高的收益，则确定使用更新后的功率，否则，仍使用更新前的功率。

根据已存在的一种潜博弈决策过程算法(即better response algorithm^[2])，能使得各自参与者的发射功率值分布更为合理，从而获得更高网络整体效用。但是，该算法从高到底顺序调整发射功率的决策过程导致***收敛慢(即到达每个参与者都不再更新发射功率的均衡点所需时间更长)，尤其在步长太小的情况下，网络应用将无法容忍极慢的收敛速度。若步长太大，则达到均衡点时的整体效用会不够理想。因此，基于该算法思想，本发明提出一种更快的博弈决策过程算法，其基本思路是，先通过折半查找缩短顺序调整发射功率的行动集空间，然后从高到底顺序调整发射功率。

在本发明中，若UE i选择UE j作为其从基站接收数据的中继，则优先选择与基站的距离小于R_th的中继。因为这类中继与基站之间信道条件更好，且已经被分配了从基站到自身的中继信道，能够直接与基站通信；而与基站的距离大于R_th的中继没有分配与基站直接通信的中继信道，因此，它要么需要寻找离基站更近的邻居作为中继，要么复用基站分配给UE i的信道以达到与基站直接通信的目的。

当中继UE j位于以基站为中心的半径R_th范围内时，从基站通过UE j中转到达UE i的最大吞吐量，相比于从基站到UE i的直接链路的最大吞吐量，其提升幅度值表示为T^ji _in；若UE j位于以基站为中心的半径R_th范围外时，其提升幅度值表示为T^ji _out。相应的度量由公式(21)～(22)估算。

类似地，若中继UE j继续选择任一节点(如UE k)作为中继，则相应的吞吐量提升度量公式由(23)～(24)估算。

本发明算法描述中用到的一些数据结构描述如下：为方便算法描述，定义如下数据结构：1)L^c _M,N表示用户直接与基站通信的信道分配关系矩阵，其任一元素

的值若为“真”，则表示某用户j(即j∈{1,…,M})被分配了某信道k(即k∈{1,…,N})用于与基站直接通信；若为“假”，则表示没有分配信道；该数据结构保存的信息作为本发明算法执行的输入参数。2)L^r _M,M表示用户之间的中继信道分配关系矩阵，其任一元素l^r _i,j的值若为k(即k∈{1,…,N})，则表示从某用户i(即i∈{1,…,M})到某用户j(即j∈{1,…,M})之间的链路被分配了某中继信道k(即k∈{1,…,N})，否则表示没有分配中继信道。3)L^r _M表示基站与用户之间的中继信道分配关系数组，其任一元素l^r _j的值若为k(即k∈{1,…,N}，则表示从基站到某用户j(即j∈{1,…,M})之间的链路被分配了某中继信道k(即k∈{1,…,N})，否则表示没有分配中继信道。4)V_i表示第i个用户以最大发射功率所覆盖的邻域节点集合，

5)集合C＝{c_j|j∈{1,…,N}}表示所有的信道集合；该集合保存的信息作为本发明算法执行的输入参数。

另外，若UE j作为中继从基站接收别的UE的数据，则将α_0j置为1来记录这种情况，同样，若UE i选择中继UE j帮它转发数据，则将α_ji置为1来记录这种情况。d_ji表示数据经过UE j的中转从基站传输到接收端UE i的传输延时；d_kji表示数据依次经过UE k、UE j的中转从基站传输到接收端UE i的传输延时；d_th表示延时阈值(可取110纳秒/比特)。

通信链路j→i的帧成功率f(γ_i)可由公式(25)计算得到。

在公式(25)中，l和γ_i分别是数据包的帧格式长度和接收节点i的SINR。基于计算得到的链路帧成功率和本实施方案给定的延时定义，通信链路j→i的数据传输延时值主要取决于接收端i的转发能力，若接收端i为传输路径上最终的接收者，则其不再转发，故不再计算其转发延时。通过下列公式(26)可得到接收端i的转发延时，其中，t^b _i是UE i的数据包转发能力。

从基站经过中继UE j到达接收端UE i的传输延时由公式(27)近似得到。

从基站依次经过中继UE k和UE j到达接收端UE i的传输延时由公式(28)近似得到。

使用蜂窝下行平均传输路径数据量来衡量蜂窝下行平均传输路径持续服务能力，可由公式(29)计算得到。

AVG_D＝Sum_D/Amount_D2D (29)

其中AVG_D表示平均传输路径数据量，Sum_D表示所有完整的D2D路径的传输路径数据量之和，Amount_D2D表示***中存在的所有D2D路径的数目，不包括与基站直接通信的节点的路径数目。某条路径的传输路径数据量，其值(如下行传输路径k→j→i)可由公式(30)计算得到。

在(30)中，D_kji表示接收节点i通过下行传输路径k→j→i从基站接收的数据包的数量；t_r是接收节点i的持续服务时长，其值取决于中继节点k和j能量储备和功耗；RBL_k和RBL_j分别是中继节点k，j的剩余能量储备；p_k和p_j则分别是中继节点k，j作为链路发射端时的发射功耗；中继节点k，j处理每比特数据的电子元器件耗能，相比于其作为发射端每比特的数据传输耗能，其值可以忽略不计。

本发明提供的这种新型蜂窝下行通行能效优化方法，总共分为三个阶段，共8个部分，具体步骤如下：

第1阶段：基站分配中继信道

算法1-1：蜂窝中继节点的信道(用于从基站直接接收数据)分配方法

步骤1.0：将L^r _M的每个元素，以l^r _j为例，初始化为0；创建集合U_c和U_r，并初始化为空。进入步骤1.1。

步骤1.1：针对L^c _M,N中的每个元素，以

为例，若存在

为“真”，则将j节点加入到U_c集合中，否则加入到U_r集合中。进入步骤1.2。

步骤1.2：针对U_r中的每个距离小于R_th且大于r_th的元素，以j为例，创建集合D并初始化为空。进入步骤1.3。

步骤1.3：针对U_r中的每个距离小于R_th且大于r_th的非j元素，以k为例，根据公式(15)计算k与j之间关于基站的夹角

若

小于基站天线发射角

并且l^r _k大于0，则将l^r _k加入到D中。进入步骤1.4。

步骤1.4：针对U_c集合中的每个元素，以i为例，根据公式(15)计算i与j之间关于基站的夹角

若

小于基站天线发射角

并且存在满足l^c _i,k为“真”的某个信道k，则将该信道k加入到D中。进入步骤1.5。

步骤1.5：创建集合E，其中的元素由集合C中所有不在集合D中的元素组成，从集合E中随机取出一个元素k，并将k赋值给l^r _j。重复步骤1.2，直至遍历完U_r中的每个元素。

算法1-2：D2D中继链路信道(用于接收D2D中继链路发射端发射的数据)分配方法

步骤2.0：将L^r _M,M中的每个元素，以l^r _i,j为例，初始化为0。

步骤2.1：针对U_r中的每个距离大于r_th的元素，以j为例，创建集合D并初始化为空。进入步骤2.2。

步骤2.2：若l^r _j大于0，则将l^r _j加入到D中。进入步骤2.3。

步骤2.3：针对U_c中的每个元素，以i为例，若i是j的邻居集合V_j中元素，则在L^c _M,N中找到i与基站直接通信时使用的信道k，并将k加入到D中。进入步骤2.4。

步骤2.4：针对U_r中的每个距离小于R_th且大于r_th的元素，以i为例，若i是j的邻居集合V_j中元素且l^r _i大于0，则将l^r _i加入到D中。进入步骤2.5。

步骤2.5：针对j的邻居集合V_j内的每个元素，以i为例，若l^r _i,j大于0，则将l^r _i,j加入到D中。进入步骤2.6。

步骤2.6：针对j的邻居集合V_j内的每个元素，以k为例，针对k的邻居集合V_k内的每个元素，以l为例，若l^r _l,k大于0，则将l^r _l,k加入到D中。进入步骤2.7。

步骤2.7：针对j的邻居集合V_j内的每个元素，以i为例，若i到基站的距离大于j到基站的距离且l^r _j,i为0，则创建集合D`并初始化为空。进入步骤2.8。

步骤2.8：针对U_c中的每个元素，以m为例，根据公式(15)计算m与i之间关于基站的夹角

若

小于基站天线发射角

则在L^c _M,N找到m与基站直接通信时候使用的信道k，并将k加入到D`中。进入步骤2.9。

步骤2.9：针对U_r中的每个距离小于R_th且大于r_th的元素，以m为例，根据公式(15)计算m与i之间关于基站的夹角

若

小于基站天线发射角

并且l^r _m大于0，则将l^r _m加入到D`中。进入步骤2.10。

步骤2.10：针对i的邻居V_i内的每个元素，以k为例，针对k的邻居V_k内的每个元素，以m为例，若l^r _k,m大于0，则将l^r _k,m加入到D`中。进入步骤2.11。

步骤2.11：创建集合E，其中的元素由集合C中所有不在集合D和D`中的元素组成，从集合E中随机取出一个元素k，(或从集合E中取出第一个元素k)，并将k赋值给l^r _j,i。重复步骤2.7，直至遍历完V_j内的每个元素。进入步骤2.12。

步骤2.12：重复步骤2.1，直至遍历完U_r中的每个元素。

第2阶段：链路接收端节点在基站的辅助下调节其发射端功率

算法2-1：调节基站对接收节点的发射功率

该算法由下行流量的接收节点(即在L^c _M,N中记录了已获得蜂窝信道的节点)执行，该类节点需要事先获得全网所有节点的最大发射功率参数。

步骤1.0：针对从基站直接接收自己数据的节点，以j为例，初始化p_0j为0，若从基站接收到T_s ^max,P_s ^max,T_s ^half,P_s ^half，则判断从基站到自身的距离d_0,j是否小于R_th，如果小于则进入步骤1.1，否则不进行功率调节。

步骤1.1：根据公式(19)计算基站对自身采用发射功率p₀ ^half下(即最大发射功率p₀ ^max的一半)的效益μ_j ^half和信噪比，以及最大功率下p₀ ^max的效益μ_j ^max,若μ_j ^half大于等于μ_j ^max并且在一半功率情况下的信噪比大于等于预设的信噪比阈值γ_th，将p_0j赋值为p₀ ^half，否则将p_0j赋值为p₀ ^max，并将p_0j发送给基站，进入步骤1.2。

步骤1.2：若从基站接收到更新的信息(即更新的T_s和P_s)，则进入步骤1.3；若从基站接收到终止信息，则结束运行；若未接收到任何信息，则继续等待。

步骤1.3：若j与基站的距离d_0,j小于R_th，则根据公式(19)计算基站采用p_0j减去步长ε后的值作为发射功率时的效益μ`<sub>j</sub>和信噪比，以及采用p<sub>0j</sub>作为发射功率时的效益μ<sub>j</sub>，若μ`_j大于μ_j且信噪比大于预设的信噪比阈值γ_th，则将p_0j的值更新为减去了步长ε后的值，并将更新后的p_0j发送给基站。重复步骤1.2。

算法2-2：调节基站对D2D中继节点的发射功率

该算法由蜂窝中继节点(即在L^c _M,N中记录了未获得蜂窝信道，但在L^r _M中记录了已获得中继信道的节点)执行，该类节点需要事先获得全网所有节点的最大发射功率参数。

步骤2.0：针对从基站接收其他节点数据的中继节点，以j为例。初始化p_0j为0，若从基站接收到T_s ^max,P_s ^max,T_s ^half,P_s ^half，则判断从基站到自身的距离d_0,j是否小于R_th且大于r_th，若满足条件则进入步骤2.1，否则不进行功率调节。

步骤2.1：根据公式(19)计算基站对自身采用发射功率p₀ ^half下(即最大发射功率p₀ ^max的一半)的效益μ_j ^half和信噪比，以及最大功率下p₀ ^max的效益μ_j ^max,若μ_j ^half大于等于μ_j ^max并且在一半功率情况下的信噪比大于等于预设的信噪比阈值γ_th，将p_0j赋值为p₀ ^half，否则将p_0j赋值为p₀ ^max，并将p_0j发送给基站，进入步骤2.2。

步骤2.2：若从基站接收到更新的信息(即更新的T_s和P_s)，则进入步骤2.3；若从基站接收到终止信息，则结束运行；若未接收到任何信息，则继续等待。

步骤2.3：若j与基站的距离d_0,j小于R_th且大于r_th，则根据公式(19)计算基站采用p_0j减去步长ε后的值作为发射功率时的效益μ`<sub>j</sub>和信噪比，以及采用p<sub>0j</sub>作为发射功率时的效益μ<sub>j</sub>，若μ`_j大于μ_j且信噪比大于预设的信噪比阈值γ_th，则将p_0j的值更新为减去了步长ε后的值，并将更新后的p_0j发送给基站。重复步骤2.2。

算法2-3：调节D2D中继链路发射端的功率

该算法由D2D中继节点(即在L^c _M,N中记录了未获得蜂窝信道，但在L^r _M,M中记录了已获得中继信道的节点)执行，该类节点需要事先获得全网所有节点的最大发射功率参数。

步骤3.0：针对任一D2D中继节点，以j为例，若从基站接收到T_s ^max,P_s ^max,T_s ^half,P_s ^half，则针对自己邻居集合V_j中所有已分配中继信道的邻居，以i为例，初始化p_ij为0。进入步骤3.1。

步骤3.1：根据公式(19)计算i对j采用发射功率p_i ^half(即最大发射功率p_i ^max的一半)时的效益μ_j ^half和信噪比，以及最大功率下p_i ^max的效益μ_j ^max，若μ_j ^half大于等于μ_j ^max并且在一半功率情况下的信噪比大于等于预设的信噪比阈值γ_th，则将p_ij赋值为p_i ^half，否则将p_ij赋值为p_i ^max。进入步骤3.2。

步骤3.2：若邻居集合V_j中的元素未遍历完，则返回步骤3.1；否则将发射功率集合{p_ij|i∈V_j}中的所有功率值发送给基站，然后进入步骤3.3。

步骤3.3：若从基站接收到更新的信息(即更新的T_s和P_s)，则进入步骤3.4；若从基站接收到终止信息，则结束运行；若未接收到任何信息，则继续等待。

步骤3.4：针对自己邻居集合V_j中所有已分配中继信道的邻居，以i为例，根据公式(19)计算i采用p_ij减去步长ε后的值作为发射功率时的效益μ`<sub>j</sub>和信噪比，以及采用p<sub>ij</sub>作为发射功率时的效益μ<sub>j</sub>，若μ`_j大于μ_j且信噪比大于预设的信噪比阈值γ_th，则将p_ij的值更新为减去了步长ε后的值；在遍历完邻居集合V_j中的所有元素后，将发射功率集合{p_ij|i∈V_j}中的所有功率值发送给基站。重复步骤3.3。

算法2-4：基站辅助UE进行功率调节

该算法基站执行，且基站需要事先获得全网所有节点的最大发射功率参数。

步骤4.0：先根据公式(20)计算所有发射端都采用最大功率时的吞吐量之和T_s ^max和功率之和P_s ^max，以及所有发射端都采用最大功率的一半时的吞吐量之和T_s ^half和功率之和P_s ^half，然后向所有节点广播T_s ^max,P_s ^max,T_s ^half,P_s ^half。进入步骤4.1。

步骤4.1：初始化功率改变标志flag为“假”，将定时器t_Δ设置为Δ，在定时器t_Δ没有超时的情况下，若基站收到来自任一UE的功率报告，则检验功率有无发生变化，若发生了变化，则将功率改变标志flag置为“真”。进入步骤4.2。

步骤4.2：若功率改变标志flag为“真”，则根据公式(20)更新T_s和P_s值，并将更新的T_s和P_s广播给所有的UE，并重复步骤4.1。否则，向所有节点广播功率调节终止包，结束算法的运行。

第3阶段：数据接收模式的确定

算法3-1：接收节点确定数据接收模式

该算法由下行流量的接收节点(即在L^r _M,M中记录了已获得蜂窝信道的节点)执行。

步骤1.0：针对任一获得蜂窝信道的从基站接收自己数据的节点，以UE i为例，若收到基站的功率调节终止信息，则将α_0j、α_ji、T^ji _in、T^ji _out初始化为0,将集合D_i初始化为“空”，将标志位flag初始化为“假”。进入步骤1.1。

步骤1.1：针对i的邻居集合V_i中的每个邻居，以j为例。若j到基站的距离d_0,j小于R_th，则根据公式(21)计算T^ji _in，否则根据公式(22)计算T^ji _out。对所有的T^ji _in和T^ji _out分别按从大到小的顺序排序。进入步骤1.2。

步骤1.2：若标志位flag为“假”，则找到当前最大的T^ji _in并赋值给T^ji _max，然后将当前最大的T^ji _in值置为0，否则，找到当前最大的T^ji _out并赋值给T^ji _max，然后将当前最大的T^ji _out值置为0。进入步骤1.3。

步骤1.3：若T^ji _max大于0，则说明选择中继比直接通信效果更好，因此，继续判断路径延时与中继能量约束是否满足要求；若d_ji小于延迟阈值d_th并且j的剩余能量e_j大于能量阈值e_th，进入步骤1.4，否则，返回步骤1.2；若T^ji _max不大于0且标志位flag为“假”，则将标志位flag置为“真”，返回步骤1.2，否则，说明当前节点找不到合适的中继，结束运行。

步骤1.4：给邻居j发送中继请求包；将定时器t_Δ设置为Δ，在定时器t_Δ没有超时的情况下，若收到了j的响应中继请求的响应包，则将α_ji和α_0j分别置为1，将j加入集合给D_i，退出等待接收响应包的操作；若定时器t_Δ超时且未选中中继(即α_ji为0)，则返回步骤1.2；若选中了中继(即α_ji为1)，则结束运行。

算法3-2：D2D中继节点确定数据接收模式

该算法由D2D中继节点(即在L^c _M,N中记录了已获得中继信道的节点)执行。

步骤2.0：针对任一D2D中继节点，以j为例，若接收到任一邻居i的中继请求包，查看自身是否已充当了中继，若没有为任何邻居充当中继，则向i发送中继请求的响应包，否则，不发送任何响应。并进入步骤2.1。

步骤2.1：将α_0k、α_kj、T^kji _in、T^kji _out分别初始化为0，将集合D_j初始化为“空”，将标志位flag初始化为“假”。进入步骤2.2。

步骤2.2：针对j的邻居集合V_j中的每个邻居，以k为例。若k到基站的距离d_0,k小于R_th且大于r_th，则根据公式(23)计算T^kji _in，否则，若k到基站的距离d_0,k不小于R_th，则根据公式(24)计算T^kji _out。对所有的T^kji _in和T^kji _out分别按从大到小的顺序排序。进入步骤2.3。

步骤2.3：若标志位flag为“假”，则找到当前最大的T^kji _in并赋值给T^kji _max，然后将当前最大的T^kji _in值置为0，否则，找到当前最大的T^kji _out并赋值给T^kji _max，然后将当前最大的T^kji _out值置为0。进入步骤2.4。

步骤2.4：若T^kji _max大于0，则说明选择中继比直接通信效果更好，因此，继续判断路径延时与中继能量约束是否满足要求；若d_kji小于延迟阈值d_th并且k的剩余能量e_k大于能量阈值e_th，进入步骤2.5，否则，返回步骤2.3；若T^kji _max不大于0且标志位flag为“假”，则将标志位flag置为“真”，返回步骤2.3，否则，说明当前节点找不到合适的中继，结束运行。

步骤2.5：给邻居k发送中继请求包；将定时器t_Δ设置为Δ，在定时器t_Δ没有超时的情况下，若收到了k的中继响应包，则将α_kj和α_0k分别置为1，将k加入集合给D_j，不再接收其它响应包；若定时器t_Δ超时且未选中中继(即α_kj为0)，则返回步骤2.3；若选中了中继(即α_kj为1)，则结束运行。

以下以一个仿真实例说明本发明方法与现有技术的差别：

本发明(为方便表述，称为DEEO-GPA)和接收节点直接与基站通信(为方便表述，称为ORGIN)在平均传输路径能量效率上进行了比较。平均传输路径能量效率特指接收节点所在下行传输路径上的数据接收能力与付出的功率代价(包括基站和中继节点的发射功率之和)之比的平均值；本发明也与未对中继进行能量约束的方案(为方便表述，称为DEEO-NO-ETH)，在平均传输路径持续服务能力方面进行了比较。可用平均传输路径数据量来衡量平均传输路径持续服务能力，其特指接收节点所在下行传输路径断连之前，接收节点平均从基站接收的数据包的个数。

基本仿真参数设置如下：仿真网络是半径为500m的圆形区域，一个基站位于该区域中央，拥有的可分配蜂窝信道数量N为500，该区域内UE数量M变化范围从1000到2800，延时约束阈值d_th、能量约束阈值e_th分别设置为110纳秒/比特,0.1焦耳，信道带宽设置为2兆赫兹。其它参数取值见表1。

表1仿真参数设置表

采用OMNeT++4.6网络仿真器得到图2至图3所示结果。从图2我们可以发现，无论蜂窝内节点数量M如何变化，DEEO-GPA的平均传输路径能量效率始终优于ORGIN。其主要原因是，对于远离基站的蜂窝边缘节点来说，DEEO-GPA通过中继链路辅助其接收基站数据，高比特误码率得到了有效地控制，而直接通信中，基站即使以最大功率发射，也无法满足蜂窝边缘处节点的比特误码率。而根据公式(6)～(13)可知，低比特误码率提升了下行数据传输路径的吞吐量，从而也提升了单位耗能所获得的吞吐量增益。另外，DEEO-GPA通过改进的博弈决策过程算法将中继节点的发射功率控制在合理的范围内，在下行传输路径吞吐量没有降低的基础上，降低了中继节点的发射功率，从而进一步地提升了蜂窝下行通信的能量效率。

从图3可观察到类似于图2的变化趋势，无论蜂窝内节点数量M如何变化，DEEO-GPA的平均传输路径持续服务能力始终优于未对中继的剩余能量进行约束的DEEO-NO-ETH方案。虽然DEEO-NO-ETH算法可能会有更多的接收节点选择中继，但是实质上吞吐量的增长率比不上中继数目的增长，并且，当UE数目达到一定数目后，由于UE的邻居也随之增多，可选择的中继也增多了，最终导致两个方案选择的中继几乎一样多，DEEO-NO-ETH方案未对中继的剩余能量进行约束，导致其可能会选择更多的较低能量储备的中继。因此由公式(29)和(30)可知，在平均传输路径持续服务能力方面，DEEO-GPA会高于DEEO-NO-ETH的方案。

Claims (8)

1.一种新型蜂窝下行通信能效优化方法，包括如下步骤：

基站分配中继信道的步骤：具体包括蜂窝中继节点信道分配的步骤和D2D中继链路信道分配的步骤；所述的蜂窝中继节点信道分配的步骤，具体为所分配的信道用于从基站直接接收数据，具体包括如下步骤：

步骤1.0：将L^r _M的每个元素初始化为0；创建集合U_c和U_r，并初始化为空；进入步骤1.1；所述L^r _M为基站与用户之间的中继信道分配关系数组，其任一元素

的值为k且k∈{1,…,N}，则表示从基站到用户j之间的链路被分配了中继信道k，否则表示没有分配中继信道；j∈{1,…,M}；

步骤1.1：针对L^c _M,N中的每个元素l^c _j,k，若存在l^c _j,k为“真”，则将j节点加入到U_c集合中，否则加入到U_r集合中；进入步骤1.2；L^c _M,N表示用户直接与基站通信的信道分配关系矩阵，其任一元素l^c _j,k的值若为“真”，则表示用户j被分配了信道k用于与基站直接通信；若为“假”，则表示没有分配信道；

步骤1.2：针对U_r中的每个距离小于R_th且大于r_th的元素，用j表示，创建集合D并初始化为空；进入步骤1.3；

步骤1.3：针对U_r中的每个距离小于R_th且大于r_th的非j元素，用k表示，计算k与j之间关于基站的夹角

若

小于基站天线发射角

并且l^r _k大于0，则将l^r _k加入到D中；进入步骤1.4；

步骤1.4：针对U_c集合中的每个元素，用i表示，计算i与j之间关于基站的夹角

若

小于基站天线发射角

并且存在满足l^c _i,k为“真”的信道k，则将该信道k加入到D中；进入步骤1.5；

步骤1.5：创建集合E，其中的元素由集合C中所有不在集合D中的元素组成，从集合E中随机取出一个元素k，并将k赋值给

重复步骤1.2，直至遍历完U_r中的每个元素；

基站辅助链路接收端节点调节发射端功率的步骤：具体包括调节基站对接收节点的发射功率的步骤、调节基站对D2D中继节点的发射功率的步骤、调节D2D中继链路发射端的功率的步骤和基站辅助UE进行功率调节的步骤，所述的UE为用户设备；

确定数据接收模式的步骤：具体包括接收节点确定数据接收模式的步骤和D2D中继节点确定数据接收模式的步骤。

2.根据权利要求1所述的新型蜂窝下行通行能效优化方法，其特征在于所述的D2D中继链路信道分配的步骤，为分配的D2D中继链路信道用于接收D2D中继链路发射端发射的数据；具体包括如下步骤：

步骤2.0：将L^r _M,M中的每个元素初始化为0；L^r _M,M表示用户之间的中继信道分配关系矩阵，其任一元素l^r _i,j的值若为k，k∈{1,…,N}，表示从用户i到用户j之间的链路被分配了中继信道k，否则表示没有分配中继信道；

步骤2.1：针对U_r中的每个距离大于r_th的元素，创建集合D并初始化为空；进入步骤2.2；

步骤2.2：若

大于0，则将

加入到D中；进入步骤2.3；

步骤2.3：针对U_c中的每个元素，以i表示，若i是j的邻居集合V_j中元素，则在L^c _M,N中找到i与基站直接通信时使用的信道k，并将k加入到D中；进入步骤2.4；

步骤2.4：针对U_r中的每个距离小于R_th且大于r_th的元素，以i表示，若i是j的邻居集合V_j中元素且l^r _i大于0，则将l^r _i加入到D中；进入步骤2.5；

步骤2.5：针对j的邻居集合V_j内的每个元素，以i表示，若l^r _i,j大于0，则将l^r _i,j加入到D中；进入步骤2.6；

步骤2.6：针对j的邻居集合V_j内的每个元素，以k表示，针对k的邻居集合V_k内的每个元素，以l表示，若l^r _l,k大于0，则将l^r _l,k加入到D中；进入步骤2.7；

步骤2.7：针对j的邻居集合V_j内的每个元素，以i表示，若i到基站的距离大于j到基站的距离且

为0，则创建集合D`并初始化为空；进入步骤2.8；

步骤2.8：针对U_c中的每个元素，以m表示，计算m与i之间关于基站的夹角

若

小于基站天线发射角

则在L^c _M,N找到m与基站直接通信时候使用的信道k，并将k加入到D`中；进入步骤2.9；

步骤2.9：针对U_r中的每个距离小于R_th且大于r_th的元素，以m表示，计算m与i之间关于基站的夹角

若

小于基站天线发射角

并且l^r _m大于0，则将l^r _m加入到D`中；进入步骤2.10；

步骤2.10：针对i的邻居V_i内的每个元素，以k表示，针对k的邻居V_k内的每个元素，以m表示，若l^r _k,m大于0，则将l^r _k,m加入到D`中；进入步骤2.11；

步骤2.11：创建集合E，其中的元素由集合C中所有不在集合D和D`中的元素组成，从集合E中随机取出一个元素k，或从集合E中取出第一个元素k，并将k赋值给

重复步骤2.7，直至遍历完V_j内的每个元素；进入步骤2.12；

步骤2.12：重复步骤2.1，直至遍历完U_r中的每个元素。

3.根据权利要求2所述的新型蜂窝下行通行能效优化方法，其特征在于所述的调节基站对接收节点的发射功率的步骤，具体包括如下步骤：

步骤0.0：事先获得全网所有节点的最大发射功率参数；

步骤1.0：针对从基站直接接收自己数据的节点，以j表示，初始化p_0j为0，若从基站接收到T_s ^max,P_s ^max,T_s ^half,P_s ^half，则判断从基站到自身的距离d_0,j是否小于R_th，如果小于则进入步骤1.1，否则不进行功率调节；所述的P_s ^max和P_s ^half表示包括潜在D2D链路和直接通信链路在内的整个网络潜在通信链路的发射端的总功耗最大值和对半值，而T_s ^max和T_s ^half表示整个网络中潜在通信链路的吞吐能力的最大值和对半值；

步骤1.1：计算基站对自身采用发射功率p₀ ^half下的效益μ_j ^half和信噪比，以及最大功率下p₀ ^max的效益μ_j ^max,若μ_j ^half大于等于μ_j ^max并且在一半功率情况下的信噪比大于等于预设的信噪比阈值γ_th，将p_0j赋值为p₀ ^half，否则将p_0j赋值为p₀ ^max，并将p_0j发送给基站，进入步骤1.2；p₀ ^half为最大发射功率p₀ ^max的一半；

步骤1.2：若从基站接收到T_s和P_s的更新信息，则进入步骤1.3；若从基站接收到终止信息，则结束运行；若未接收到任何信息，则继续等待；

步骤1.3：若j与基站的距离d_0,j小于R_th，则计算基站采用p_0j减去步长ε后的值作为发射功率时的效益μ`<sub>j</sub>和信噪比，以及采用p<sub>0j</sub>作为发射功率时的效益μ<sub>j</sub>，若μ`_j大于μ_j且信噪比大于预设的信噪比阈值γ_th，则将p_0j的值更新为减去了步长ε后的值，并将更新后的p_0j发送给基站，并重复步骤1.2。

4.根据权利要求3所述的新型蜂窝下行通行能效优化方法，其特征在于所述的调节基站对D2D中继节点的发射功率的步骤，具体包括如下步骤：

步骤0.0：事先获得全网所有节点的最大发射功率参数；

步骤2.0：针对从基站接收其他节点数据的中继节点，以j表示；初始化p_0j为0，若从基站接收到T_s ^max,P_s ^max,T_s ^half,P_s ^half，则判断从基站到自身的距离d_0,j是否小于R_th且大于r_th，若满足条件则进入步骤2.1，否则不进行功率调节；

步骤2.1：计算基站对自身采用发射功率p₀ ^half下的效益μ_j ^half和信噪比，以及最大功率下p₀ ^max的效益μ_j ^max,若μ_j ^half大于等于μ_j ^max并且在一半功率情况下的信噪比大于等于预设的信噪比阈值γ_th，将p_0j赋值为p₀ ^half，否则将p_0j赋值为p₀ ^max，并将p_0j发送给基站，进入步骤2.2；

步骤2.2：若从基站接收到T_s和P_s的更新信息，则进入步骤2.3；若从基站接收到终止信息，则结束运行；若未接收到任何信息，则继续等待；

步骤2.3：若j与基站的距离d_0,j小于R_th且大于r_th，则计算基站采用p_0j减去步长ε后的值作为发射功率时的效益μ`<sub>j</sub>和信噪比，以及采用p<sub>0j</sub>作为发射功率时的效益μ<sub>j</sub>，若μ`_j大于μ_j且信噪比大于预设的信噪比阈值γ_th，则将p_0j的值更新为减去了步长ε后的值，并将更新后的p_0j发送给基站，并重复步骤2.2。

5.根据权利要求4所述的新型蜂窝下行通行能效优化方法，其特征在于所述的调节D2D中继链路发射端的功率的步骤，具体包括如下步骤：

步骤0.0：事先获得全网所有节点的最大发射功率参数；

步骤3.0：针对任一D2D中继节点，以j表示，若从基站接收到T_s ^max,P_s ^max,T_s ^half,P_s ^half，则针对自己邻居集合V_j中所有已分配中继信道的邻居，以i表示，初始化p_ij为0；进入步骤3.1；

步骤3.1：计算i对j采用发射功率p_i ^half时的效益μ_j ^half和信噪比，以及最大功率下p_i ^max的效益μ_j ^max，若μ_j ^half大于等于μ_j ^max并且在一半功率情况下的信噪比大于等于预设的信噪比阈值γ_th，则将p_ij赋值为p_i ^half，否则将p_ij赋值为p_i ^max；进入步骤3.2；

步骤3.2：若邻居集合V_j中的元素未遍历完，则返回步骤3.1；否则将发射功率集合{p_ij|i∈V_j}中的所有功率值发送给基站，然后进入步骤3.3；

步骤3.3：若从基站接收到T_s和P_s的更新信息，则进入步骤3.4；若从基站接收到终止信息，则结束运行；若未接收到任何信息，则继续等待；

步骤3.4：针对自己邻居集合V_j中所有已分配中继信道的邻居，以i表示，计算i采用p_ij减去步长ε后的值作为发射功率时的效益μ`<sub>j</sub>和信噪比，以及采用p<sub>ij</sub>作为发射功率时的效益μ<sub>j</sub>，若μ`_j大于μ_j且信噪比大于预设的信噪比阈值γ_th，则将p_ij的值更新为减去了步长ε后的值；在遍历完邻居集合V_j中的所有元素后，将发射功率集合{p_ij|i∈V_j}中的所有功率值发送给基站，并重复步骤3.3。

6.根据权利要求5所述的新型蜂窝下行通行能效优化方法，其特征在于所述的基站辅助UE进行功率调节的步骤，具体包括如下步骤：

步骤0.0：基站事先获得全网所有节点的最大发射功率参数；

步骤4.0：计算所有发射端都采用最大功率时的吞吐量之和T_s ^max和功率之和P_s ^max，以及所有发射端都采用最大功率的一半时的吞吐量之和T_s ^half和功率之和P_s ^half，然后向所有节点广播T_s ^max,P_s ^max,T_s ^half,P_s ^half；进入步骤4.1；

步骤4.1：初始化功率改变标志flag为“假”，将定时器t_Δ设置为Δ，在定时器t_Δ没有超时的情况下，若基站收到来自任一UE的功率报告，则检验功率有无发生变化，若发生了变化，则将功率改变标志flag置为“真”；进入步骤4.2；

步骤4.2：若功率改变标志flag为“真”，则更新T_s和P_s值，并将更新的T_s和P_s广播给所有的UE，并重复步骤4.1；否则，向所有节点广播功率调节终止包，结束算法的运行。

7.根据权利要求6所述的新型蜂窝下行通行能效优化方法，其特征在于所述的接收节点确定数据接收模式的步骤，具体包括如下步骤：

步骤1.0：针对任一获得蜂窝信道的从基站接收自己数据的节点，以UE i表示，若收到基站的功率调节终止信息，则将α_0j、α_ji、T^ji _in、T^ji _out初始化为0,将集合D_i初始化为“空”，将标志位flag初始化为“假”；进入步骤1.1；

步骤1.1：针对i的邻居集合V_i中的每个邻居，以j表示；若j到基站的距离d_0,j小于R_th，则计算T^ji _in，否则计算T^ji _out；对所有的T^ji _in和T^ji _out分别按从大到小的顺序排序；进入步骤1.2；

步骤1.2：若标志位flag为“假”，则找到当前最大的T^ji _in并赋值给T^ji _max，然后将当前最大的T^ji _in值置为0，否则，找到当前最大的T^ji _out并赋值给T^ji _max，然后将当前最大的T^ji _out值置为0；进入步骤1.3；

步骤1.3：若T^ji _max大于0，则说明选择中继比直接通信效果更好，因此，继续判断路径延时与中继能量约束是否满足要求；若d_ji小于延迟阈值d_th并且j的剩余能量e_j大于能量阈值e_th，进入步骤1.4，否则，返回步骤1.2；若T^ji _max不大于0且标志位flag为“假”，则将标志位flag置为“真”，返回步骤1.2，否则，说明当前节点找不到合适的中继，结束运行；

步骤1.4：给邻居j发送中继请求包；将定时器t_Δ设置为Δ，在定时器t_Δ没有超时的情况下，若收到了j的响应中继请求的响应包，则将α_ji和α_0j分别置为1，将j加入集合给D_i，退出等待接收响应包的操作；若定时器t_Δ超时且未选中中继，则返回步骤1.2；若选中了中继，则结束运行。

8.根据权利要求7所述的新型蜂窝下行通行能效优化方法，其特征在于所述的D2D中继节点确定数据接收模式的步骤，具体包括如下步骤：

步骤2.0：针对任一D2D中继节点，以j表示，若接收到任一邻居i的中继请求包，查看自身是否已充当了中继，若没有为任何邻居充当中继，则向i发送中继请求的响应包，否则，不发送任何响应；并进入步骤2.1；

步骤2.1：将α_0k、α_kj、T^kji _in、T^kji _out分别初始化为0，将集合D_j初始化为“空”，将标志位flag初始化为“假”；进入步骤2.2；

步骤2.2：针对j的邻居集合V_j中的每个邻居，以k表示；若k到基站的距离d_0,k小于R_th且大于r_th，则计算T^kji _in，否则，若k到基站的距离d_0,k不小于R_th，则计算T^kji _out；对所有的T^kji _in和T^kji _out分别按从大到小的顺序排序；进入步骤2.3；

步骤2.3：若标志位flag为“假”，则找到当前最大的T^kji _in并赋值给T^kji _max，然后将当前最大的T^kji _in值置为0，否则，找到当前最大的T^kji _out并赋值给T^kji _max，然后将当前最大的T^kji _out值置为0；进入步骤2.4；

步骤2.4：若T^kji _max大于0，则说明选择中继比直接通信效果更好，因此，继续判断路径延时与中继能量约束是否满足要求；若d_kji小于延迟阈值d_th并且k的剩余能量e_k大于能量阈值e_th，进入步骤2.5，否则，返回步骤2.3；若T^kji _max不大于0且标志位flag为“假”，则将标志位flag置为“真”，返回步骤2.3，否则，说明当前节点找不到合适的中继，结束运行；

步骤2.5：给邻居k发送中继请求包；将定时器t_Δ设置为Δ，在定时器t_Δ没有超时的情况下，若收到了k的中继响应包，则将α_kj和α_0k分别置为1，将k加入集合给D_j，不再接收其它响应包；若定时器t_Δ超时且未选中中继，则返回步骤2.3；若选中了中继，则结束运行。

Images