中继辅助的设备对设备覆盖/底层通信(D2D通信技术)的传输容量分析(Matlab代码实现)

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

---

💥1 概述

文献来源：

摘要—设备对设备（D2D）通信通过在蜂窝网络中为移动用户提供直接链接，能够有效地满足高数据速率的需求。本文分析了在覆盖和底层模式下与蜂窝网络共存的中继辅助的D2D通信的传输容量。在覆盖模式下，D2D用户可以使用精细的频谱资源，而在底层模式下则重用蜂窝资源。基于随机几何，网络中的蜂窝用户、D2D发送器和中继节点（RNs）都被建模为泊松点过程。然后，我们计算中继节点存在概率和中继链路距离的期望，以获得D2D通信的成功传输概率。根据两种增强D2D传输距离的中继机制，我们进一步获得了在两种模式下中继辅助的D2D通信的传输容量，这反映了D2D密度和功率的影响。此外，还分析了两种模式下具有可变D2D链路距离的D2D传输容量。模拟结果验证了通过中继传输可以增强D2D传输容量，并且受到多种因素的影响，包括用户密度、功率、D2D链路距离以及使用中继节点的方式。

关键词—设备对设备（D2D）通信，中继，随机几何，传输容量。

无线通信的快速发展导致了频谱资源的短缺，成为进一步推动无线服务爆炸式发展的瓶颈【1】。第五代（5G）网络引入了许多令人期待的新技术，如非正交多址（NOMA）【2】、大规模MIMO【3】、设备对设备（D2D）通信等。作为5G关键技术之一，D2D通信通过直接链路实现用户之间的无线传输，带来了许多优势，如增强传输容量、降低传输功率、提高频谱效率等【4】。然而，由于D2D通信会重用频率资源，因此可能对现有网络造成有害干扰，可能降低系统的预期性能【5】–【7】。

作为通信系统的基本性能指标之一，D2D传输容量已经在一些早期文献中进行了分析【8】–【11】。为了提高蜂窝和D2D系统的容量，文献【8】提出了一种基于干扰控制的干扰受限区域管理方案，以管理网络中的干扰。进一步，Yu等人【9】研究了蜂窝和D2D用户之间的资源分配和功率控制，并提出了一种优化方法以提高系统容量。考虑到网络的复杂性，文献【10】研究了双向数据传输所带来的网络容量。此外，在考虑了D2D中的合作通信的情况下，文献【11】提出了一种优化的D2D传输容量，其中D2D传输可以由网络中的蜂窝用户辅助。

 在本文中，我们分析了在覆盖和底层模式下，中继节点（RNs）辅助的设备对设备（D2D）通信的传输容量。基于随机几何，网络通过异构泊松点过程（PPP）进行建模。我们以封闭形式获得了中继节点用于传输距离延伸和容量提升的D2D传输容量表达式，这导致了两种模式下的传输容量增益。此外，还考虑了两种模式中可变D2D直接链路的影响，这可以用于评估实际网络设计。进一步的模拟结果表明，传输容量受到如何使用中继节点的方式的影响。此外，结果还表明D2D传输容量受到D2D直接链路距离以及蜂窝、D2D和RN的密度的影响，这些因素以不同程度引入了对D2D系统的干扰。

📚2 运行结果

 

部分代码：

%% Program Intitialization
clear;
close all;
clc;
format long g;
%% 
% This MATLAB script is to simulate and plot the relationship between relay
% node density (# of relay nodes / sq. m) and the transmission capacity of
% the D2D network (in overlay mode)

% This is the code for CASE 1, where we consider that D2D comms can only
% take place, provided a relay node is present in the "shadow region"

lambda_1 = 0.0003; % D2D density in the area
lambda_2 = 0.0004; % RN density in the area
P_1 = 0.031622; % in Watts ,15 dbm, D2D transmission power
P_2 = 0.031622; % in Watts ,15 dbm, relay node transmission power

T_1 = 1.25893; % 2 db, D2D SIR threshold
T_2 = 1.25893; % 2 db, RN SIR threshold

alpha = 4; % Path loss exponent
R = 35; % D2D link distance

% Values defined in the paper
S = (2*pi/3 - sqrt(3)/2)*R^2;
Pr_e = 1 - exp(-1*lambda_2*S);
A = (4*pi - 36*sqrt(3) + 64)/(12*pi - 9*sqrt(3));
E = A*R;
C_alpha = pi*gamma(1 + 2/alpha)*gamma(1-(2/alpha));

% SIR threshold probability
sir_d_to_rn(lambda_1, Pr_e, C_alpha, T_2, alpha, E , P_2, P_1);
sir_rn_to_d(lambda_1, Pr_e, C_alpha, T_2, alpha, E , P_2, P_1);

transmission_capacity(lambda_1, Pr_e, C_alpha, T_2, alpha, E , P_2, P_1)

% What we want
transmission_capacity(lambda_1, Pr_e, C_alpha, T_2, alpha, E , P_2, P_1); % transmission capacity in bps/Hz

rn_density = [0:0.00001:0.002];
tcs = zeros(1, length(rn_density));
for i = 1: length(rn_density)
    x = rn_density(i);
    Pr_e = 1 - exp(-1*x*S);
    tcs(i) = transmission_capacity(lambda_1, Pr_e, C_alpha, T_2, alpha, E , P_2, P_1);
end

plot(rn_density, tcs);
%% Function definitions
function [prob_d_to_rn] = sir_d_to_rn(lambda_1, Pr_e, C_alpha, T_2, alpha, E , P_2, P_1)
    prob_d_to_rn= exp(-0.5*lambda_1*Pr_e*C_alpha*(T_2^(2/alpha))*(E^2)*(((P_2/P_1)^(2/alpha)) + 1));
end

function [prob_rn_to_d] = sir_rn_to_d(lambda_1, Pr_e, C_alpha, T_2, alpha, E , P_2, P_1)
    prob_rn_to_d= exp(-0.5*lambda_1*Pr_e*C_alpha*(T_2^(2/alpha))*(E^2)*(((P_1/P_2)^(2/alpha)) + 1));
end

function [tc] = transmission_capacity(lambda_1, Pr_e, C_alpha, T_2, alpha, E , P_2, P_1)
    tc = 0.5*lambda_1*Pr_e*sir_d_to_rn(lambda_1, Pr_e, C_alpha, T_2, alpha, E , P_2, P_1)*sir_rn_to_d(lambda_1, Pr_e, C_alpha, T_2, alpha, E , P_2, P_1);

end

🎉3 参考文献

部分理论来源于网络，如有侵权请联系删除。

🌈4 Matlab代码实现