一种基于NOMA的车联网广播通信资源分配方法与流程

一种基于noma的车联网广播通信资源分配方法
技术领域
1.本发明属于无线通信技术领域，尤其涉及一种基于noma的车联网广播通信资源分配方法。


背景技术：

2.车联网资源分配是提升车联网性能的重要因素，对车联网的优化具有重要的作用。基站bs集中式资源分配建立在基站可获知覆盖范围内的所有用户的信道状态信息(csi)，而车联网场景中车辆信道是快速时变信道，集中式资源分配方法不再适用。此外，由于无线频谱资源与功率资源有限，传统oma下，车联网中的网络接口规模受限。因此车联网资源分配受到多方面挑战。
3.车联网资源分配对频谱资源、功率资源存在优化需求，且在车辆用户逐日增多的情况下，车联网用户入网接入需求逐渐增多，车联网安全性受到挑战。因此有必要开发一种基于noma的车联网广播通信资源分配方法。现有文献大多数研究的是基站bs集中式资源分配，常见的有基于d2d的车联网资源分配方法、基于图论的v2v资源分配方法等。基站bs集中式资源分配在车联网通信场景中，因其难以获取快速时变信道csi而难以实施。基于d2d的车联网资源分配方法除了上述原因，还会在基站bs端占用大量的计算开销，不适用于v2v通信场景，现有的基于图论的v2v资源分配方法文献多数对图的规模有限制而难以适配车联网大连接需求场景。此外，现有车联网资源分配方案多集中于车联网系统吞吐量的提升，对车联网用户公平性的研究较为欠缺。


技术实现要素：

4.本发明目的在于提供一种基于noma的车联网广播通信资源分配方法,以解决上述的技术问题。
5.为解决上述技术问题，本发明的具体技术方案如下：
6.一种基于noma的车联网广播通信资源分配方法，包括以下步骤：
7.步骤1、在半静态调度(sps)的开始阶段，基站bs根据覆盖范围内的车辆与车辆通信用户类型中的发送方(v2v tx)用户的地理位置，对v2v tx归并成簇，同簇内的用户数最多为θ，分簇条件为：各个簇内的v2v tx用户彼此的距离不低于η且分簇过程使得各簇内v2v tx的最小间距最大化；
8.同簇内的v2v tx用户可共享同一子信道，且每个蜂窝车辆用户v2i用户只与一个v2v tx簇共享子信道：
9.步骤2、根据步骤1中产生的簇数与基站覆盖范围内的v2i用户数，对kuhn
‑
munkres算法做出调整，将各v2v tx簇与v2i用户的距离矩阵算出，记为d
v2i_cluster
，将调整后的kuhn
‑
munkres算法应用于距离矩阵d
v2i_cluster
的负矩阵，即
‑
d
v2i_cluster
，以获取v2i用户与各v2i tx簇子信道共享结果；
10.步骤3、在sps的数据传输阶段，v2v tx用户发送广播信息给d≤r(其中，r为v2v tx
用户的最大发送距离)范围内的v2v rx用户，v2v rx用户根据sic解码条件及noma下的接收各v2v tx的信干噪比情况，以及v2v tx对于v2i用户干扰受限的条件下，将pso算法应用于分布式场景中，以pf函数作为适应度函数，设定迭代次数及精确度获取v2v tx用户的发射功率最优解。
11.进一步的，步骤1具体包括以下步骤：
12.步骤1.1、定义迭代变量t，令t＝1，对基站bs覆盖范围内的所有v2v tx建立标识，记为tx1、tx2、
…
、tx
p
，基站bs覆盖范围内存在p个v2v tx用户，采用生成随机数的方法，将p个v2v tx用户的标识自动填入一个空簇中，记单个簇的最大v2v tx规模数为θ，填满一个空簇，继续下一个空簇填入余下的(p
‑
θ)个v2v tx用户的标识，以此类推，直至所有v2v tx用户的标识全部填入，生成个簇；当为整数时，否则其中(p/θ)的运算结果为整数商值，不计余数；
13.步骤1.2、对步骤1.1中生成的簇，计算每个簇中v2v tx车辆彼此之间的距离，并找出各个簇内v2v tx车辆彼此之间距离的最小值，记为tx车辆彼此之间距离的最小值，记为
14.步骤1.3、找出中的最小值，记为定义单个簇内可接受的v2v tx用户之间最小距离为η，若则回到步骤1.1；若则继续步骤1.4；
15.步骤1.4、令t＝2，重复步骤1.1到步骤1.3，若新的v2v tx分簇结果计算得到的大于则选择t＝2时的v2v tx分簇结果作为更新后的分簇结果，然后t＝t 1；
16.步骤1.5、重复步骤1.1到步骤1.4，直至重复连续n0次后，结果不再更新，则把最后一次更新所对应的v2v tx簇作为最终的v2v tx分簇结果。
17.进一步的，步骤2具体包括以下步骤：
18.步骤2.1、计算每个v2i用户与步骤1中生成的每个v2v tx簇之间的距离，计算方法为计算v2i用户与每个簇内v2v tx用户距离的最小值，将其作为v2i用户与该v2v tx簇的距离；
19.步骤2.2、若基站bs覆盖范围内有m个v2i用户，每个v2i用户不与其他v2i用户共享同一子信道，则步骤2.1中得到的v2i用户与v2v tx簇的距离有(其中，为权利要求2中分簇算法生成的v2v tx簇数)个数值，将这些距离数值保存在维度为的矩阵中，记为
20.步骤2.3、若直接对步骤2.2中的距离矩阵的负矩阵，即采用kuhn
‑
munkres算法，获取m个v2i用户与v2v tx簇的子信道共享结果；
21.步骤2.4、若补充个虚拟v2v tx簇，即在的右侧增加个全零列，使得新的矩阵为维度(m
×
m)的方阵，记为对的负矩阵，即采用kuhn
‑
munkres算法，获取m个v2i用户与v2v tx簇的子信道共享结
果；
22.步骤2.5、若补充个v2i用户，在步骤2.2中的距离矩阵的下方补充个全零行，使得新的矩阵为维度的方阵，记为对的负矩阵，即采用kuhn
‑
munkres算法，获取m个v2i用户与v2v tx簇的子信道共享结果。
23.进一步的，步骤3具体包括以下步骤：
24.步骤3.1、对于基站bs覆盖范围内的所有v2v rx用户，计算其与各v2vtx用户之间的距离数值，并将计算出的数值保存在维度为q
×
p的矩阵中，其中q为基站bs覆盖范围内的v2v rx用户数，p为基站bs覆盖范围内的v2vtx用户数记该矩阵为
25.步骤3.2、建立维度为q
×
p的辅助矩阵对于对于(其中，i为v2v rx的标识，j为v2v tx的标识，逻辑符号表示“任意”)，若则否则若则说明v2v rx i接收到v2v tx j的广播数据；
26.步骤3.3、对于基站bs覆盖范围内的每一个v2v rx用户建立noma下的接收各v2v tx广播信息的传输速率计算。
27.进一步的，步骤3.3具体包括以下步骤：
28.步骤3.3.1、记当前所要计算的v2v rx用户i接收v2v tx用户j广播数据的信干噪比为
[0029][0030]
其中i为v2v rx用户的标识，1≤i≤q；j，k为v2v tx用户的标识，1≤j，k≤p；m为v2i用户的标识，1≤m≤m；为v2v tx用户j的发射功率，为v2i用户m的发射功率，为v2v tx用户k的发射功率，σ2为噪声功率；为v2v rx用户i到v2v tx用户j的信道增益系数，为v2v rx用户i到v2i用户m的信道增益系数，为v2v rx用户i到v2v tx用户k的信道增益系数；β
i，m
为二进制变量，当β
i，m
＝1时，表示v2v rx用户i与v2i用户m共享同一子信道；当β
i，m
＝0时，表示v2v rx用户i与v2i用户m占用不同的子信道；τ
i，k
为二进制变量，当时，τ
i，k
＝1，此时v2v rx用户i接收到v2v tx用户j的广播信息；否则当时，τ
i，k
＝0，此时v2v rx用户i不在v2v tx用户j的覆盖范围，无法接收到v2v tx用户j的广播信息；ρ
i，k
为二进制变量，具体定义如下：
[0031][0032]
步骤3.3.2、考虑v2i用户对于同信道下v2v tx用户干扰的承受能力，考虑对于v2i
用户来说，保证v2i用户的信干噪比γ
m
的下限，设定pr{γ
m
≤γ0}≤p0(其中，pr{.}表示概率分布函数，γ0为可接收的v2i用户最小信干噪比，p0为v2i用户信干噪比小于γ0时可承受的概率阈值)；其中，v2i用户的信干噪比γ
m
的信干噪比如下：
[0033][0034]
其中为基站bs对于v2i用户m的发射功率，为v2v tx用户i的发射功率，为基站bs与v2i用户m之间的信道增益常数，为v2vtx用户i与v2i用户m之间的信道增益常数，μ
i，m
为二进制变量，当v2v tx用户i与v2i用户共享同一子信道时，μ
i，m
＝1，否则μ
i，m
＝0；σ2为噪声功率；
[0035]
步骤3.3.3，定义时隙变量slot，令slot＝1，产生的n
s
个维数为p
×
1的向量并计算其相对应的比例公平函数(pf)函数值，将其作为ps0算法的适应度函数值，从而得到第一个时隙的最优pf函数值；
[0036]
步骤3.3.4，令slot＝slot 1，重复步骤3.3.3，直到slot＞n
sps
(其中，n
sps
为sps的数据传输阶段的时隙数)时停止。
[0037]
进一步的，所述步骤3.3.2中v2v tx用户对于v2i用户的影响主要在于距离v2i用户最近的，且与v2i用户共享同一子信道的v2v tx用户的发射功率的影响；在noma应用于v2v tx的情况下，v2v tx与v2v rx之间信道差异较大，且发射功率随距离衰减较快，车辆之间的分布满足泊松分布，因此可对上述γ
m
的信干噪比公式做近似运算，即：
[0038][0039]
其中角标i
*
表示与v2i用户m距离最近的v2v tx用户的标识；考虑到v2i用户m的信干噪比qos条件，即pr{γ
m
≤γ0}≤p0，可将其进行如下转化：
[0040]
记其中α
bm
、为大尺度衰落系数，g
bm
、为小尺度衰落系数，则可进一步转化γ
m
≤γ0(其中，γ0为v2i用户可承受的最小信干噪比)得到：
[0041][0042]
故所以pr{γ
m
≤γ0}可表示如下二重积分：
[0043][0044]
因此，可将pr{γ
m
≤γ0}≤p0表示为：
[0045][0046]
将与v2i用户m距离最近的v2v tx用户i
*
发射功率转换为：
[0047][0048]
进一步的，所述步骤3.3.3具体包括以下步骤：
[0049]
步骤3.3.3.1、在步骤3.3.2中与v2i用户m距离最近且共享同一子信道的v2v tx用户i
*
发射功率约束下，对所有v2v tx用户的传输功率做初始化，即产生n
s
个维数为p
×
1的功率值向量生成方式为：对于接收到v2v tx用户广播信息的v2v rx用户，分担pso算法计算量，不同的v2v rx用户之间通过交叉区域的v2v tx用户交换各自的功率分配、数据接收情况；
[0050]
步骤3.3.3.2，设定一个sps周期中的数据传输阶段有n
sps
个时隙，定义pso算法中的迭代变量itera，令itera＝1时，计算每个v2v rx用户的传输速率，将所有的传输速率值记在维度为q
×
p的raterxtx
q
×
p
矩阵中，raterxtx(i，j)表示v2v rx用户i接收到v2v tx用户j广播信息的传输速率：
[0051][0052]
其中e[.]为期望运算符；表示noma应用于v2v用户场景的情况下，v2v rx用户i接收到v2v tx用户j的广播信息的信干噪比，已在步骤3.3.1中给出；对v2v rx用户i所能接收到的所有v2v tx广播信息的传输速率累加，即对raterxtx
q
×
p
矩阵各行求和，并将求和之后的各v2v rx总传输速率保存在维度为q
×
n
sps
的矩阵中中，初始化为全零矩阵，即：
[0053][0054]
其中p为基站bs覆盖范围内v2v tx用户的总数；
[0055]
步骤3.3.3.3、计算目前在第slot个时隙时的适应度函数fitness值，适应度函数fitness定义如下：
[0056]
当slot＝1时，设维度为1
×
q的权重向量初始值为weight＝(1，1，...，1)
[0057]
fitness＝weight
·
raterxslot(：，slot)
[0058]
＝(1，1，...，1)(raterxslot(1，slot)，...，raterxslot(q，slot))
t
；
[0059]
当1＜slot≤n
sps
时，weight的第i个分量记为weight
i
，其中1≤i≤q，计算方法为：
[0060][0061]
fitness＝weight
·
raterxslot(：，slot)
[0062]
＝(weight1，...，weight
q
)(raterxslot(1，slot)，...，raterxslot(q，slot))
t
[0063]
步骤3.3.3.4、对n
s
个维数为p
×
1的功率值向量分别计算对应的n
s
个fitness值；将最大fitness值对应的功率值向量记为将各个功率向量目前最大fitness值对应的功率值向量记为生成各的更新梯度向量其中：
[0064]
当itera＝1时，
[0065]
v
t
＝rand(1
×
p)；(其中，v
t
为第t个功率值向量对应的更新梯度向量)
[0066]
当itera＞1时，
[0067]
(其中，为上一次pso迭代更新后的第t个功率值向量)
[0068]
其中rand产生(0，1)之间的随机数，c1为各个功率向量自身的更新速度系数，c2为各功率向量的全局更新速度系数；
[0069]
步骤3.3.3.5、利用步骤3.3.3.4的更新梯度向量对n
s
个维数为p
×
1的功率值向量进行更新，
[0070]
步骤3.3.3.6，令itera＝itera 1，重复3.3.3.2至3.3.3.5，直到itera＝maxgen，其中maxgen为最大迭代次数，最终获得maxgen轮迭代的最优fitness值，及其对应的最优功率值向量
[0071]
本发明的一种基于noma的车联网广播通信资源分配方法，具有以下优点：
[0072]
本发明将noma应用于v2v广播通信下行链路场景中，在基站bs为v2v tx用户预留子信道的基础上，基于kuhn
‑
munkres方法，解决了v2v用户与v2i用户子信道共享问题，然后使用基于pf函数的分布式多用户粒子群算法，解决了v2v rx用户吞吐量、传输时延及用户
公平性的问题。
附图说明
[0073]
图1为本发明的平直公路车联网广播通信下行链路场景图；
[0074]
图2为本发明的基于noma的车联网广播通信资源分配方法的流程图；
具体实施方式
[0075]
为了更好地了解本发明的目的、结构及功能，下面结合附图，对本发明一种基于noma的车联网广播通信资源分配方法做进一步详细的描述。
[0076]
如图1所示，其中bs为基站，v2v tx为车辆与车辆通信用户类型中的发送方，v2v rx为车辆与车辆通信用户类型中的接收方接受方，v2i为蜂窝车辆用户，r为基站服务半径，d0为基站到平直公路的距离。
[0077]
一种基于noma的车联网广播通信下行链路资源分配方法，如图1、2所示，包括以下步骤：
[0078]
步骤1、在半静态调度(sps)的开始阶段，基站bs根据覆盖范围内的车辆与车辆通信用户类型中的发送方(v2v tx)用户的地理位置，对v2v tx归并成簇，同簇内的用户数最多为θ，分簇条件：1、各个簇内的v2v tx用户彼此的距离不低于η；2、分簇过程使得各簇内v2v tx的最小间距最大化。同簇内的v2v tx用户可共享同一子信道，且每个v2i用户只与一个v2v tx簇共享子信道，具体步骤如下：
[0079]
步骤1.1、令t＝1，对基站bs覆盖范围内的所有v2v tx建立标识，记为tx1、tx2、
…
、tx
p
(假设基站bs覆盖范围内存在p个v2v tx用户)，采用生成随机数的方法，将p个v2v tx用户的标识自动填入一个空簇中，记单个簇的最大v2v tx规模数为θ，填满一个空簇，继续下一个空簇填入余下的(p
‑
θ)个v2v tx用户的标识，以此类推，直至所有v2v tx用户的标识全部填入，生成个簇。当为整数时，否则其中(p/θ)的运算结果为整数商值，不计余数；
[0080]
步骤1.2、对步骤1.1中生成的簇，计算每个簇中v2v tx车辆彼此之间的距离，并找出各个簇内v2v tx车辆彼此之间距离的最小值，记为tx车辆彼此之间距离的最小值，记为
[0081]
步骤1.3、找出中的最小值，记为定义单个簇内可接受的v2v tx用户之间最小距离为η，若则回到步骤1.1；若则继续步骤1.4；
[0082]
步骤1.4、令t＝2，重复步骤1.1到步骤1.3，如果新的v2v tx分簇结果计算得到的大于则选择t＝2时的v2v tx分簇结果作为更新后的分簇结果，然后t＝t 1；
[0083]
步骤1.5、重复步骤1.1到步骤1.4，直至重复连续n0次后，结果不再更新，则把最后一次更新所对应的v2v tx簇作为最终的v2v tx分簇结果。
[0084]
步骤2、根据步骤1中产生的簇数与基站覆盖范围内的v2i用户数，对kuhn
‑
munkres算法做出调整，将各v2v tx簇与v2i用户的距离矩阵算出，记为d
v2i_cluster
，将调整后的kuhn
‑
munkres算法应用于距离矩阵d
v2i_cluster
的负矩阵，即
‑
d
v2i_cluster
，以获取v2i用户与各v2itx簇子信道共享结果，具体步骤如下：
[0085]
步骤2.1、计算每个v2i用户与步骤1中生成的每个v2v tx簇之间的距离，计算方法为计算v2i用户与每个簇内v2v tx用户距离的最小值，将其作为v2i用户与该v2v tx簇的距离；
[0086]
步骤2.2、若基站bs覆盖范围内有m个v2i用户(假设每个v2i用户不与其他v2i用户共享同一子信道)，则步骤2.1中得到的v2i用户与v2v tx簇的距离有个数值，其中，为权利要求2中分簇算法生成的v2v tx簇数，将这些距离数值保存在维度为的矩阵中，记为
[0087]
步骤2
‑
3、若直接对步骤2.2中的距离矩阵的负矩阵，即采用kuhn
‑
munkres算法，获取m个v2i用户与v2v tx簇的子信道共享结果；
[0088]
步骤2.4、若补充个虚拟v2v tx簇，即在的右侧增加个全零列，使得新的矩阵为维度(m
×
m)的方阵，记为对的负矩阵，即采用kuhn
‑
munkres算法，获取m个v2i用户与v2v tx簇的子信道共享结果；
[0089]
步骤2.5、若补充个v2i用户，在步骤2.2中的距离矩阵的下方补充个全零行，使得新的矩阵为维度的方阵，记为对的负矩阵，即采用kuhn
‑
munkres算法，获取m个v2i用户与v2v tx簇的子信道共享结果；
[0090]
步骤3、在半静态调度(sps)的数据传输阶段，v2v tx用户发送广播信息给d≤r范围内的v2v rx用户，其中r为v2v tx用户的最大发送距离，v2v rx用户为车辆与车辆通信用户类型中的接收方用户；v2v rx用户根据sic解码条件及noma下的接收各v2v tx的信干噪比情况，以及v2v tx对于v2i用户干扰受限的条件下，将pso算法应用于分布式场景中，以pf函数作为适应度函数，设定迭代次数及精确度获取v2v tx用户的发射功率最优解，具体步骤如下：
[0091]
步骤3.1、对于基站bs覆盖范围内的所有v2v rx用户，计算其与各v2vtx用户之间的距离数值，并将计算出的数值保存在维度为q
×
p的矩阵中，其中q为基站bs覆盖范围内的v2v rx用户数，p为基站bs覆盖范围内的v2vtx用户数记该矩阵为
[0092]
步骤3.2、建立维度为q
×
p的辅助矩阵对于对于(其中，i为v2v rx的标识，j为v2v tx的标识，逻辑符号表示“任意”)，若则否则若则说明v2v rx i接收到v2v tx j的广播数据。
[0093]
步骤3.3、对于基站bs覆盖范围内的每一个v2v rx用户建立noma下的接收各v2v tx广播信息的传输速率计算，具体步骤如下：
[0094]
步骤3.3.1、记当前所要计算的v2v rx用户i接收v2v tx用户j广播数据的信干噪比为
[0095][0096]
其中i为v2v rx用户的标识(1≤i≤q)，j，k为v2v tx用户的标识(1≤j，k≤p)，m为v2i用户的标识(1≤m≤m)。为v2v tx用户j的发射功率，为v2i用户m的发射功率，为v2v tx用户k的发射功率，σ2为噪声功率。为v2v rx用户i到v2v tx用户j的信道增益系数，为v2v rx用户i到v2i用户m的信道增益系数，为v2v rx用户i到v2v tx用户k的信道增益系数。β
i，m
为二进制变量，当β
i，m
＝1时，表示v2v rx用户i与v2i用户m共享同一子信道；当β
i，m
＝0时，表示v2v rx用户i与v2i用户m占用不同的子信道。τ
i，k
为二进制变量，当时，τ
i，k
＝1，此时v2v rx用户i可以接收到v2v tx用户j的广播信息；否则当对，τ
i，k
＝0，此时v2v rx用户i不在v2v tx用户j的覆盖范围，无法接收到v2v tx用户j的广播信息。ρ
i，k
为二进制变量，具体定义如下：
[0097][0098]
步骤3.3.2、考虑v2i用户对于同信道下v2v tx用户干扰的承受能力，考虑对于v2i用户来说，保证v2i用户的信干噪比γ
m
的下限，设定pr{γ
m
≤γ0}≤p0(其中，pr{.}表示概率分布函数，γ0为可接收的v2i用户最小信干噪比，p0为v2i用户信干噪比小于γ0时可承受的概率阈值)；其中，v2i用户的信干噪比γ
m
的信干噪比如下：
[0099][0100]
其中为基站bs对于v2i用户m的发射功率，为v2v tx用户i的发射功率，为基站bs与v2i用户m之间的信道增益常数，为v2vtx用户i与v2i用户m之间的信道增益常数，μ
i，m
为二进制变量，当v2v tx用户i与v2i用户共享同一子信道时，μ
i，m
＝1，否则μ
i，m
＝0；σ2为噪声功率。
[0101]
实际场景中，v2v tx用户对于v2i用户的影响主要在于距离v2i用户最近的，且与v2i用户共享同一子信道的v2v tx用户的发射功率的影响。在noma应用于v2v tx的情况下，v2v tx与v2v rx之间信道差异较大，且发射功率随距离衰减较快，车辆之间的分布满足泊松分布，因此可对上述γ
m
的信干噪比公式做近似运算，即：
[0102][0103]
其中角标i
*
表示与v2i用户m距离最近的v2v tx用户的标识。考虑到v2i用户m的信干噪比qos条件，即pr{γ
m
≤γ0}≤p0。可将其进行如下转化：
[0104]
记其中α
bm
、为大尺度衰落系数，g
bm
、为小尺度衰落系数。则可进一步转化γ
m
≤γ0(其中，γ0为v2i用户可承受的最小信干噪比)得到：
[0105][0106]
故所以pr{γ
m
≤γ0}可表示如下二重积分：
[0107][0108]
因此，可将pr{γ
m
≤γ0}≤p0表示为：
[0109][0110]
此时，可将与v2i用户m距离最近的v2v tx用户i
*
发射功率转换为：
[0111][0112]
步骤3.3.3、定义时隙变量slot，令slot＝1，产生的n
s
个维数为p
×
1的向量并计算其相对应的比例公平函数(pf)函数值，将其作为pso算法的适应度函数值，从而得到第一个时隙的最优pf函数值。具体步骤如下：
[0113]
步骤3.3.3.1、在步骤3.3.2中与v2i用户m距离最近且共享同一子信道的v2v tx用户i
*
发射功率约束下，对所有v2v tx用户的传输功率做初始化，即产生ns个维数为p
×
1的功率值向量生成方式为：对于接收到v2v tx用户广播信息的v2v rx用户，分担pso算法计算量，不同的v2v rx用户之间通过交叉区域的v2v tx用户交换各自的功率分配、数据接收情况。
[0114]
步骤3.3.3.2，设定一个sps周期中的数据传输阶段有n
sps
个时隙，定义pso算法中的迭代变量itera，令itera＝1时，计算每个v2v rx用户的传输速率，将所有的传输速率值记在维度为q
×
p的raterxtx
q
×
p
矩阵中，raterxtx(i，j)表示v2v rx用户i接收到v2v tx用户j广播信息的传输速率：
[0115][0116]
其中e[.]为期望运算符；表示noma应用于v2v用户场景的情况下，v2v rx用户i接收到v2v tx用户j的广播信息的信干噪比，已在步骤3.3.1中给出；对v2v rx用户i所能接收到的所有v2v tx广播信息的传输速率累加，即对raterxtx
q
×
p
矩阵各行求和，并将求和之后的各v2v rx总传输速率保存在维度为q
×
n
sps
的矩阵中中，初始化为全零矩阵，即：
[0117][0118]
其中p为基站bs覆盖范围内v2v tx用户的总数。
[0119]
步骤3.3.3.3、计算目前在第slot个时隙时的适应度函数fitness值，适应度函数fitness定义如下：
[0120]
当slot＝1时，设维度为1
×
q的权重向量初始值为weight＝(1，1，...，1)
[0121]
fitness＝weight
·
raterxslot(：，slot)
[0122]
＝(1，1，...，1)(raterxslot(1，slot)，...，raterxslot(q，slot))
t
；
[0123]
当1＜slot≤n
sps
时，weight的第i个分量记为weight
i
，其中1≤i≤q，计算方法为：
[0124][0125]
fitness＝weight
·
raterxslot(：，slot)
[0126]
＝(weight1，...，weight
q
)(raterxslot(1，slot)，...，raterxslot(q，slot))
t
[0127]
步骤3.3.3.4、对n
s
个维数为p
×
1的功率值向量分别计算对应的n
s
个fitness值；将最大fitness值对应的功率值向量记为将各个功率向量目前最大fitness值对应的功率值向量记为生成各的更新梯度向量其中：
[0128]
当itera＝1时，
[0129]
v
t
＝rand(1
×
p)；(其中，v
t
为第t个功率值向量对应的更新梯度向量)
[0130]
当itera＞1时，
[0131]
(其中，为上一次pso迭代更新后的第t个功率值向量)
[0132]
其中rand产生(0，1)之间的随机数，c1为各个功率向量自身的更新速度系数，c2为各功率向量的全局更新速度系数；
[0133]
步骤3.3.3.5、利用步骤3.3.3.4的更新梯度向量对n
s
个维数为p
×
1的功率值向量进行更新，
[0134]
步骤3.3.3.6，令itera＝itera 1，重复3.3.3.2至3.3.3.5，直到itera＝maxgen，其中maxgen为最大迭代次数，最终获得maxgen轮迭代的最优fitness值，及其对应的最优功率值向量
[0135]
步骤3
‑3‑
4，令slot＝slot 1，重复步骤3.3.3，直到slot＞n
sps
(其中，n
sps
为sps的数据传输阶段的时隙数)时停止。最终给出各时隙的功率分配与所对应的最优pf函数值。
[0136]
本发明综合使用noma在v2v广播通信下行链路的应用特征，解决了基站bs在车联网快速时变信道场景下，难获取车辆用户csi以解决资源分配的问题，并保证了v2v rx用户公平性，算法可有效收敛至最优解，有效提升了频谱利用率，扩大了网络接口规模。
[0137]
可以理解，本发明是通过一些实施例进行描述的，本领域技术人员知悉的，在不脱离本发明的精神和范围的情况下，可以对这些特征和实施例进行各种改变或等效替换。另外，在本发明的教导下，可以对这些特征和实施例进行修改以适应具体的情况及材料而不会脱离本发明的精神和范围。因此，本发明不受此处所公开的具体实施例的限制，所有落入本技术的权利要求范围内的实施例都属于本发明所保护的范围内。