✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,代码获取、论文复现及科研仿真合作可私信。
🍎个人主页:Matlab科研工作室
🍊个人信条:格物致知。
更多Matlab完整代码及仿真定制内容点击👇
🔥 内容介绍
一个大规模的多输入多输出(MIMO)系统利用大量基站(BS)天线为一组用户服务,由于细胞间干扰(ICI)而受到飞行员污染。在这封信中,提出了一个基于图形着色的试点分配(GC-PA)方案,以减轻多单元大规模MIMO系统的试点污染。具体而言,利用衰落通道的大规模特征,首先构建了一个干涉图来描述所有用户的潜在ICI关系。然后,在有限的试点资源下,拟议的GC-PA方案旨在通过在干涉图中有效地向用户分配试点来减少潜在的ICI。通过仿真验证了该方案的性能增益.
📣 部分代码
clear;
clc;
tic;
% system parameters
M = 512; % number of BS antennas M
L = 7; % cells
K = 10; % users
S = 15; % number of pilots
lambda=0.1;
R_Cell = 1000;
r_Min = 100;
alpha = 3.8;
sigma_shadow = 8;
% q = 15;
rho_ul = 10^(15/10); % 15 dB
% Test_num = 300;
Test_num = 15;
UL_SINR_LS_sw = zeros(K, L, Test_num);
UL_SINR_MMSE_sw = zeros(K, L, Test_num);
UL_SINR_MLE_sw = zeros(K, L, Test_num);
UL_SINR_LS_wg = zeros(K, L, Test_num);
UL_SINR_MMSE_wg = zeros(K, L, Test_num);
UL_SINR_MLE_wg = zeros(K, L, Test_num);
UL_SINR_LS_sp = zeros(K, L, Test_num);
UL_SINR_MMSE_sp = zeros(K, L, Test_num);
UL_SINR_MLE_sp = zeros(K, L, Test_num);
UL_SINR_LS_ra = zeros(K, L, Test_num);
UL_SINR_MMSE_ra = zeros(K, L, Test_num);
UL_SINR_MLE_ra = zeros(K, L, Test_num);
% Simulation
for i_test = 1:Test_num
display(i_test);
% Generate Channel vector
[H, Beta] = F_H_Generate(M, L, K, R_Cell, r_Min, sigma_shadow, alpha);
%% random
N = S; % the length of the pilot sequence (N>=K)
temp = zeros(L,S);
for i=1:L
temp(i,:) = randperm(S);
end
pilot = F_ZC(N,S);
pilots = zeros(N,K,L);
for l=1:L
for k=1:K
pilots(:,k,l) = pilot(:,temp(i,k));
end
end
%--- Received Pilot Signal ---
% Y = zeros(M,N,L); % received uplink training sequences at i-th BS
Y = received_signal(M,N,L,rho_ul,H,pilots);
H_est = CE_aid(M, K, L, H, pilots, rho_ul);
UL_SINR_LS_ra(:, :, i_test) = F_SINR(L,K,H,H_est,rho_ul,pilots);
% % ===== LS --> SINR =====%
% H_est = CE_LS(M, K, L, Y, pilots, rho_ul, N);
%
% UL_SINR_LS_ra(:, :, i_test) = F_SINR(L,K,H,H_est,rho_ul,pilots);
%
% % ===== MMSE --> SINR =====%
% H_est = CE_ideal_MMSE(M, K, L, Y, pilots, rho_ul, N, Beta);
%
% UL_SINR_MMSE_ra(:, :, i_test) = F_SINR(L,K,H,H_est,rho_ul,pilots);
% ===== MLE --> SINR =====%
H_est = CE_MLE(M, K, L, Y, pilots, rho_ul, N, Beta);
UL_SINR_MLE_ra(:, :, i_test) = F_SINR(L,K,H,H_est,rho_ul,pilots);
%% sprs
[pilots, N] = F_SPRS2(L,K,Beta,lambda);
Y = received_signal(M,N,L,rho_ul,H,pilots);
H_est = CE_aid(M, K, L, H, pilots, rho_ul);
UL_SINR_LS_sp(:, :, i_test) = F_SINR(L,K,H,H_est,rho_ul,pilots);
% % ===== LS --> SINR =====%
% H_est = CE_LS(M, K, L, Y, pilots, rho_ul, N);
%
% UL_SINR_LS_sp(:, :, i_test) = F_SINR(L,K,H,H_est,rho_ul,pilots);
%
% % ===== MMSE --> SINR =====%
% H_est = CE_ideal_MMSE(M, K, L, Y, pilots, rho_ul, N, Beta);
%
% UL_SINR_MMSE_sp(:, :, i_test) = F_SINR(L,K,H,H_est,rho_ul,pilots);
% ===== MLE --> SINR =====%
H_est = CE_MLE(M, K, L, Y, pilots, rho_ul, N, Beta);
UL_SINR_MLE_sp(:, :, i_test) = F_SINR(L,K,H,H_est,rho_ul,pilots);
%% wgcpa
[P, N] = F_WGCPA(L,K,S,Beta);
Y = received_signal(M,N,L,rho_ul,H,P);
H_est = CE_aid(M, K, L, H, P, rho_ul);
UL_SINR_LS_wg(:, :, i_test) = F_SINR(L,K,H,H_est,rho_ul,P);
% % ===== LS --> SINR =====%
% H_est = CE_LS(M, K, L, Y, P, rho_ul, N);
%
% UL_SINR_LS_wg(:, :, i_test) = F_SINR(L,K,H,H_est,rho_ul,P);
%
% % ===== MMSE --> SINR =====%
% H_est = CE_ideal_MMSE(M, K, L, Y, P, rho_ul, N, Beta);
%
% UL_SINR_MMSE_wg(:, :, i_test) = F_SINR(L,K,H,H_est,rho_ul,P);
% ===== MLE --> SINR =====%
H_est = CE_MLE(M, K, L, Y, P, rho_ul, N, Beta);
UL_SINR_MLE_wg(:, :, i_test) = F_SINR(L,K,H,H_est,rho_ul,P);
%% SPRS+WGC-PA
[pilots, N] = F_SPRS_WGCPA2(L,K,Beta,lambda,S);
%===== Received Pilot Signal =====%
% Y = zeros(M,N,L); % received uplink training sequences at i-th BS
Y = received_signal(M,N,L,rho_ul,H,pilots);
H_est = CE_aid(M, K, L, H, pilots, rho_ul);
UL_SINR_LS_sw(:, :, i_test) = F_SINR(L,K,H,H_est,rho_ul,pilots);
% % ===== LS --> SINR =====%
% H_est = CE_LS(M, K, L, Y, pilots, rho_ul, N);
%
% UL_SINR_LS_sw(:, :, i_test) = F_SINR(L,K,H,H_est,rho_ul,pilots);
%
% % ===== MMSE --> SINR =====%
% H_est = CE_ideal_MMSE(M, K, L, Y, pilots, rho_ul, N, Beta);
%
% UL_SINR_MMSE_sw(:, :, i_test) = F_SINR(L,K,H,H_est,rho_ul,pilots);
% ===== MLE --> SINR =====%
H_est = CE_MLE(M, K, L, Y, pilots, rho_ul, N, Beta);
UL_SINR_MLE_sw(:, :, i_test) = F_SINR(L,K,H,H_est,rho_ul,pilots);
end
%% UL SINR
figure;
Tau = 20;
y = 10*log10(reshape(UL_SINR_LS_ra,1,L*K*Test_num));
ymin=min(y);
ymax=max(y);
x=linspace(ymin,ymax,Tau);
yy=hist(y,x);
yy=yy/length(y)/(x(2)-x(1));
s=0;
for i=2:length(x)
s=[s,trapz(x([1:i]),yy([1:i]))];
end
plot(x,s,'cx-');
hold on;
% y = 10*log10(reshape(UL_SINR_MMSE_ra,1,L*K*Test_num));
% ymin=min(y);
% ymax=max(y);
% x=linspace(ymin,ymax,Tau);
% yy=hist(y,x);
% yy=yy/length(y)/(x(2)-x(1));
% s=0;
% for i=2:length(x)
% s=[s,trapz(x([1:i]),yy([1:i]))];
% end
% plot(x,s,'c*-');
% hold on;
y = 10*log10(reshape(UL_SINR_MLE_ra,1,L*K*Test_num));
ymin=min(y);
ymax=max(y);
x=linspace(ymin,ymax,Tau);
yy=hist(y,x);
yy=yy/length(y)/(x(2)-x(1));
s=0;
for i=2:length(x)
s=[s,trapz(x([1:i]),yy([1:i]))];
end
plot(x,s,'co-');
hold on;
y = 10*log10(reshape(UL_SINR_LS_sp,1,L*K*Test_num));
ymin=min(y);
ymax=max(y);
x=linspace(ymin,ymax,Tau);
yy=hist(y,x);
yy=yy/length(y)/(x(2)-x(1));
s=0;
for i=2:length(x)
s=[s,trapz(x([1:i]),yy([1:i]))];
end
plot(x,s,'bx-');
hold on;
%
% y = 10*log10(reshape(UL_SINR_MMSE_sp,1,L*K*Test_num));
% ymin=min(y);
% ymax=max(y);
% x=linspace(ymin,ymax,Tau);
% yy=hist(y,x);
% yy=yy/length(y)/(x(2)-x(1));
% s=0;
% for i=2:length(x)
% s=[s,trapz(x([1:i]),yy([1:i]))];
% end
% plot(x,s,'b*-');
% hold on;
y = 10*log10(reshape(UL_SINR_MLE_sp,1,L*K*Test_num));
ymin=min(y);
ymax=max(y);
x=linspace(ymin,ymax,Tau);
yy=hist(y,x);
yy=yy/length(y)/(x(2)-x(1));
s=0;
for i=2:length(x)
s=[s,trapz(x([1:i]),yy([1:i]))];
end
plot(x,s,'bo-');
hold on;
y = 10*log10(reshape(UL_SINR_LS_wg,1,L*K*Test_num));
ymin=min(y);
ymax=max(y);
x=linspace(ymin,ymax,Tau);
yy=hist(y,x);
yy=yy/length(y)/(x(2)-x(1));
s=0;
for i=2:length(x)
s=[s,trapz(x([1:i]),yy([1:i]))];
end
plot(x,s,'gx-');
hold on;
%
% y = 10*log10(reshape(UL_SINR_MMSE_wg,1,L*K*Test_num));
% ymin=min(y);
% ymax=max(y);
% x=linspace(ymin,ymax,Tau);
% yy=hist(y,x);
% yy=yy/length(y)/(x(2)-x(1));
% s=0;
% for i=2:length(x)
% s=[s,trapz(x([1:i]),yy([1:i]))];
% end
% plot(x,s,'g*-');
% hold on;
y = 10*log10(reshape(UL_SINR_MLE_wg,1,L*K*Test_num));
ymin=min(y);
ymax=max(y);
x=linspace(ymin,ymax,Tau);
yy=hist(y,x);
yy=yy/length(y)/(x(2)-x(1));
s=0;
for i=2:length(x)
s=[s,trapz(x([1:i]),yy([1:i]))];
end
plot(x,s,'go-');
hold on;
y = 10*log10(reshape(UL_SINR_LS_sw,1,L*K*Test_num));
ymin=min(y);
ymax=max(y);
x=linspace(ymin,ymax,Tau);
yy=hist(y,x);
yy=yy/length(y)/(x(2)-x(1));
s=0;
for i=2:length(x)
s=[s,trapz(x([1:i]),yy([1:i]))];
end
plot(x,s,'rx-');
hold on;
%
% y = 10*log10(reshape(UL_SINR_MMSE_sw,1,L*K*Test_num));
% ymin=min(y);
% ymax=max(y);
% x=linspace(ymin,ymax,Tau);
% yy=hist(y,x);
% yy=yy/length(y)/(x(2)-x(1));
% s=0;
% for i=2:length(x)
% s=[s,trapz(x([1:i]),yy([1:i]))];
% end
% plot(x,s,'r*-');
% hold on;
y = 10*log10(reshape(UL_SINR_MLE_sw,1,L*K*Test_num));
ymin=min(y);
ymax=max(y);
x=linspace(ymin,ymax,Tau);
yy=hist(y,x);
yy=yy/length(y)/(x(2)-x(1));
s=0;
for i=2:length(x)
s=[s,trapz(x([1:i]),yy([1:i]))];
end
plot(x,s,'ro-');
hold on;
grid on;
%axis([-80 40 0 1]);
title('CDF of uplink SINR');
legend('Random','MLE:Random','SPRS','MLE:SPRS','WGC-PA','MLE:WGC-PA','SPRS+WGC-PA','MLE:SPRS+WGC-PA');
xlabel('SINR(dB)');
ylabel('CDF');
hold off;
toc;
⛳️ 运行结果
🔗 参考文献
X. Zhu, L. Dai and Z. Wang, "Graph Coloring Based Pilot Allocation to Mitigate Pilot Contamination for Multi-Cell Massive MIMO Systems," in IEEE Communications Letters, vol. 19, no. 10, pp. 1842-1845, Oct. 2015, doi: 10.1109/LCOMM.2015.2471304.
🎈 部分理论引用网络文献,若有侵权联系博主删除
🎁 关注我领取海量matlab电子书和数学建模资料
👇 私信完整代码和数据获取及论文数模仿真定制
1 各类智能优化算法改进及应用
生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化、背包问题、 风电场布局、时隙分配优化、 最佳分布式发电单元分配、多阶段管道维修、 工厂-中心-需求点三级选址问题、 应急生活物质配送中心选址、 基站选址、 道路灯柱布置、 枢纽节点部署、 输电线路台风监测装置、 集装箱船配载优化、 机组优化、 投资优化组合、云服务器组合优化、 天线线性阵列分布优化
2 机器学习和深度学习方面
2.1 bp时序、回归预测和分类
2.2 ENS声神经网络时序、回归预测和分类
2.3 SVM/CNN-SVM/LSSVM/RVM支持向量机系列时序、回归预测和分类
2.4 CNN/TCN卷积神经网络系列时序、回归预测和分类
2.5 ELM/KELM/RELM/DELM极限学习机系列时序、回归预测和分类
2.6 GRU/Bi-GRU/CNN-GRU/CNN-BiGRU门控神经网络时序、回归预测和分类
2.7 ELMAN递归神经网络时序、回归\预测和分类
2.8 LSTM/BiLSTM/CNN-LSTM/CNN-BiLSTM/长短记忆神经网络系列时序、回归预测和分类
2.9 RBF径向基神经网络时序、回归预测和分类