function test05_main2_02
clear;
rand('state',sum(100 * clock));
% Tx and Rx antenna number
n_t = 4;
n_r = 4;
% shadow fading
sigma = 8;
% cell radius
radius = 500;
% save('E:\matlab\02.mimo\mat\radius.mat', 'radius');
% eNB_sig positions
enb_sig(1, :) = [0, 0] * radius;
% eNB_inf positions
intf_num = 6;
enb_inf(1, :) = [3 , -2 * sqrt(3)] * radius;
enb_inf(2, :) = [4.5 , 0.5 * sqrt(3)] * radius;
enb_inf(3, :) = [1.5 , 2.5 * sqrt(3)] * radius;
enb_inf(4, :) = [-3 , 2 * sqrt(3)] * radius;
enb_inf(5, :) = [-5.5, -0.5 * sqrt(3)] * radius;
enb_inf(6, :) = [-1.5, -2.5 * sqrt(3)] * radius;
% SNR
SNR = power(10, 20 / 10) / power(10, -100 / 10);
% users
usr_num = 1
d_sig = zeros(usr_num, 1);
d_inf = zeros(usr_num, intf_num);
for n = 1 : usr_num
% user position
r(n) = (0.1 + rand(1) * 0.9) * radius;
sita = rand(1) * 2 * pi;
x = r(n) * cos(sita);
y = r(n) * sin(sita);
% distances from eNB_sig
d_sig(n) = sqrt((enb_sig(1) - x)^2 + (enb_sig(2) - y)^2);
% distances from eNB_inf
for h = 1 : intf_num
d_inf(n, h) = sqrt((enb_inf(h, 1) - x)^2 + (enb_inf(h, 2) - y)^2);
end
end
% save('E:\matlab\02.mimo\mat\usr_num.mat', 'usr_num');
% save('E:\matlab\02.mimo\mat\r.mat', 'r');
slot_num = 1e5;
avg_rate = ones(1, usr_num) / usr_num;
sch_num = zeros(1, usr_num);
for slot = 1 : slot_num
% 1. user channel condition
req_rate = zeros(1, usr_num);
for n = 1 : usr_num
% H
tmp1 = 35.3 + 37.6 * log10(d_sig(n)) + normrnd(0, sigma);
tmp2 = normrnd(0, 1/sqrt(2), n_t, n_r) + i * normrnd(0, 1/sqrt(2), n_t, n_r);
H = sqrt(power(10, -tmp1 / 10)) * tmp2;
term1 = H * H';
% Rmtx
Rmtx = zeros(n_t, n_t);
for h = 1 : intf_num
tmp1 = 35.3 + 37.6 * log10(d_inf(n, h)) + normrnd(0, sigma);
tmp2 = normrnd(0, 1/sqrt(2), n_t, n_r) + i * normrnd(0, 1/sqrt(2), n_t, n_r);
H = sqrt(power(10, -tmp1 / 10)) * tmp2;
Rmtx = Rmtx + H * H';
end
term2 = inv(Rmtx + eye(n_t) / SNR);
% achievable rate
req_rate(n) = real(log2(det(eye(n_t) + term1 * term2)));
%
end
% 2. scheduling调度
[tmp,sch_usr] = max(req_rate ./ avg_rate);
sch_num(sch_usr) = sch_num(sch_usr) + 1;
avg_rate = 0.99 * avg_rate + 0.01 * ([1 : usr_num] == sch_usr) .* req_rate;
% 平均速率
sys_rate(slot) = sum(avg_rate);
figure;
plot( [tmp,sch_usr]);
grid;
end
