无线携能通信网络在多用户接收用户时的波束成形优化代码

clear all
clc

for j=1:50
q1 = 20;
q2 = 30;
q3 = 15;
q4 = 5;
sigma_sq = 1;

W(:,:,1) = eye(2);
W(:,:,2) = eye(2);
V = eye(2);

h1 = complex(randn(2,1),randn(2,1));
h2 = complex(randn(2,1),randn(2,1));
h3 = complex(randn(2,1),randn(2,1));
h4 = complex(randn(2,1),randn(2,1));

g1 = complex(randn(2,1),randn(2,1));
g2 = complex(randn(2,1),randn(2,1));
g3 = complex(randn(2,1),randn(2,1));
g4 = complex(randn(2,1),randn(2,1));

t = -1;

minus_SINR1 = -h1'*W(:,:,1)*h1/(h1'*W(:,:,2)*h1+g1'*V*g1+sigma_sq);
minus_SINR2 = -h2'*W(:,:,2)*h2/(h2'*W(:,:,1)*h2+g2'*V*g2+sigma_sq);
u1 = max(minus_SINR1,minus_SINR2);
k = 1;
while abs(t) >= 1e-3
cvx_begin sdp
variable W(2,2,2) complex hermitian
variable V(2,2) complex hermitian
variable t
minimize t;
subject to
trace( W(:,:,1)+W(:,:,2) ) <= q1;  % information power constraint
trace(V) <= q2;                           % energy power constraints
W(:,:,1) == hermitian_semidefinite(2);
W(:,:,2) == hermitian_semidefinite(2);
V == hermitian_semidefinite(2);

real( h3'*( W(:,:,1)+W(:,:,2) )*h3 + g3'*V*g3 ) >= q3;  % EH constraint
imag( h3'*( W(:,:,1)+W(:,:,2) )*h3 + g3'*V*g3 ) == 0;
real( h4'*( W(:,:,1)+W(:,:,2) )*h4 + g4'*V*g4 ) <= q4; % Interference constraint
imag( h4'*( W(:,:,1)+W(:,:,2) )*h4 + g4'*V*g4 ) == 0;
real( - h1'*W(:,:,1)*h1 - u1*(h1'*W(:,:,2)*h1+g1'*V*g1+sigma_sq) ) <= t; % mimus_SINR for user 1
real( - h2'*W(:,:,2)*h2 - u1*(h2'*W(:,:,1)*h2+g2'*V*g2+sigma_sq) ) <= t; % minus_SINR for user 2
imag( - h1'*W(:,:,1)*h1 - u1*(h1'*W(:,:,2)*h1+g1'*V*g1+sigma_sq) ) == 0;
imag( - h2'*W(:,:,2)*h2 - u1*(h2'*W(:,:,1)*h2+g2'*V*g2+sigma_sq) ) == 0;
cvx_end
imag(u1)==0
q(k)=real(u1)

p(k) = t;

if p(k) == 0
W_opt = W;
V_opt = V;
break;
else
minus_SINR1 = -h1'*W(:,:,1)*h1/(h1'*W(:,:,2)*h1+g1'*V*g1+sigma_sq);
minus_SINR2 = -h2'*W(:,:,2)*h2/(h2'*W(:,:,1)*h2+g2'*V*g2+sigma_sq);
u1 = max(minus_SINR1,minus_SINR2);
k = k+1;
end
end
z1(j)=u1

W(:,:,1) = eye(2);
W(:,:,2) = eye(2);
V = eye(2);
t = -1;

minus_SINR1 = -h1'*W(:,:,1)*h1/(h1'*W(:,:,2)*h1+g1'*V*g1+sigma_sq);
minus_SINR2 = -h2'*W(:,:,2)*h2/(h2'*W(:,:,1)*h2+g2'*V*g2+sigma_sq);
u2 = max(minus_SINR1,minus_SINR2);
k = 1;
g_value1 = h1'*W(:,:,2)*h1+g1'*V*g1+sigma_sq
g_value2 = h2'*W(:,:,1)*h2+g2'*V*g2+sigma_sq

while abs(t) >= 1e-3
cvx_begin sdp
variable W(2,2,2) complex hermitian
variable V(2,2) complex hermitian
variable t
minimize t;
subject to
trace( W(:,:,1)+W(:,:,2) ) <= q1;  % information power constraint
trace(V) <= q2;                           % energy power constraints
W(:,:,1) == hermitian_semidefinite(2);
W(:,:,2) == hermitian_semidefinite(2);
V == hermitian_semidefinite(2);

real( h3'*( W(:,:,1)+W(:,:,2) )*h3 + g3'*V*g3 ) >= q3;  % EH constraint
imag( h3'*( W(:,:,1)+W(:,:,2) )*h3 + g3'*V*g3 ) == 0;
real( h4'*( W(:,:,1)+W(:,:,2) )*h4 + g4'*V*g4 ) <= q4; % Interference constraint
imag( h4'*( W(:,:,1)+W(:,:,2) )*h4 + g4'*V*g4 ) == 0;
real( ( - h1'*W(:,:,1)*h1 - u2*(h1'*W(:,:,2)*h1+g1'*V*g1+sigma_sq))/g_value1 ) <= t; % mimus_SINR for user 1
real( (- h2'*W(:,:,2)*h2 - u2*(h2'*W(:,:,1)*h2+g2'*V*g2+sigma_sq))/ g_value2 ) <= t; % minus_SINR for user 2
imag( ( - h1'*W(:,:,1)*h1 - u2*(h1'*W(:,:,2)*h1+g1'*V*g1+sigma_sq))/g_value1 ) == 0;
imag( (- h2'*W(:,:,2)*h2 - u2*(h2'*W(:,:,1)*h2+g2'*V*g2+sigma_sq))/ g_value2 ) == 0;
cvx_end
imag(u2)==0
q(k)=real(u2)
p(k) = t;

if p(k) == 0
W_opt = W;
V_opt = V;
break;
else
minus_SINR1 = -h1'*W(:,:,1)*h1/(h1'*W(:,:,2)*h1+g1'*V*g1+sigma_sq);
minus_SINR2 = -h2'*W(:,:,2)*h2/(h2'*W(:,:,1)*h2+g2'*V*g2+sigma_sq);
u2 = max(minus_SINR1,minus_SINR2);
k = k+1;
end
g_value1 = h1'*W(:,:,2)*h1+g1'*V*g1+sigma_sq;
g_value2 = h2'*W(:,:,1)*h2+g2'*V*g2+sigma_sq;
end
z2(j)=u2

W(:,:,1) = eye(2);
W(:,:,2) = eye(2);

V = eye(2);

t = -1;

minus_SINR1 = -h1'*W(:,:,1)*h1/(h1'*W(:,:,2)*h1+sigma_sq);
minus_SINR2 = -h2'*W(:,:,2)*h2/(h2'*W(:,:,1)*h2+sigma_sq);
u3 = max(minus_SINR1,minus_SINR2);
k = 1;

while abs(t) >= 1e-3
cvx_begin sdp
variable W(2,2,3) complex hermitian
variable V(2,2) complex hermitian
variable t
minimize t;
subject to
trace( W(:,:,1)+W(:,:,2)+W(:,:,3) ) <= q1;  % information power constraint
trace(V) <= q2;                           % energy power constraints
W(:,:,1) == hermitian_semidefinite(2);
W(:,:,2) == hermitian_semidefinite(2);
W(:,:,3) == hermitian_semidefinite(2);
V == hermitian_semidefinite(2);

real( h3'*( W(:,:,1)+W(:,:,2) + W(:,:,3))*h3 + g3'*V*g3 ) >= q3;  % EH constraint
imag( h3'*( W(:,:,1)+W(:,:,2) + W(:,:,3) )*h3 + g3'*V*g3 ) == 0;
% real( h4'*( W(:,:,1)+W(:,:,2) )*h4 + g4'*V*g4 ) <= q4; % Interference constraint
%imag( h4'*( W(:,:,1)+W(:,:,2) )*h4 + g4'*V*g4 ) == 0;
real( - h1'*W(:,:,1)*h1 - u3*(h1'*W(:,:,2)*h1+sigma_sq) ) <= t; % mimus_SINR for user 1
real( - h2'*W(:,:,2)*h2 - u3*(h2'*W(:,:,1)*h2+sigma_sq) ) <= t; % minus_SINR for user 2
imag( - h1'*W(:,:,1)*h1 - u3*(h1'*W(:,:,2)*h1+sigma_sq) ) == 0;
imag( - h2'*W(:,:,2)*h2 - u3*(h2'*W(:,:,1)*h2+sigma_sq) ) == 0;
cvx_end

imag(u3)==0
z(k)=real(u3)

p(k) = t;

if p(k) == 0
W_opt = W;
V_opt = V;
break;
else
minus_SINR1 = -h1'*W(:,:,1)*h1/(h1'*W(:,:,2)*h1+sigma_sq);
minus_SINR2 = -h2'*W(:,:,2)*h2/(h2'*W(:,:,1)*h2+sigma_sq);
u3 = max(minus_SINR1,minus_SINR2);
k = k+1;
end
end
z3(j)=u3

W(:,:,1) = eye(2);
W(:,:,2) = eye(2);

V = eye(2);

t = -1;

minus_SINR1 = -h1'*W(:,:,1)*h1/(h1'*W(:,:,2)*h1+sigma_sq);
minus_SINR2 = -h2'*W(:,:,2)*h2/(h2'*W(:,:,1)*h2+sigma_sq);
u4 = max(minus_SINR1,minus_SINR2);
k = 1;
g_value1 = h1'*W(:,:,2)*h1+sigma_sq
g_value2 = h2'*W(:,:,1)*h2+sigma_sq

while abs(t) >= 1e-3
cvx_begin sdp
variable W(2,2,3) complex hermitian
variable V(2,2) complex hermitian
variable t
minimize t;
subject to
trace( W(:,:,1)+W(:,:,2)+W(:,:,3) ) <= q1;  % information power constraint
trace(V) <= q2;                           % energy power constraints
W(:,:,1) == hermitian_semidefinite(2);
W(:,:,2) == hermitian_semidefinite(2);
W(:,:,3) == hermitian_semidefinite(2);
V == hermitian_semidefinite(2);

real( h3'*( W(:,:,1)+W(:,:,2) + W(:,:,3))*h3 + g3'*V*g3 ) >= q3;  % EH constraint
imag( h3'*( W(:,:,1)+W(:,:,2) + W(:,:,3) )*h3 + g3'*V*g3 ) == 0;
% real( h4'*( W(:,:,1)+W(:,:,2) )*h4 + g4'*V*g4 ) <= q4; % Interference constraint
%imag( h4'*( W(:,:,1)+W(:,:,2) )*h4 + g4'*V*g4 ) == 0;
real( ( - h1'*W(:,:,1)*h1 - u4*(h1'*W(:,:,2)*h1+sigma_sq))/g_value1 ) <= t; % mimus_SINR for user 1
real( (- h2'*W(:,:,2)*h2 - u4*(h2'*W(:,:,1)*h2+sigma_sq))/ g_value2 ) <= t; % minus_SINR for user 2
imag( ( - h1'*W(:,:,1)*h1 - u4*(h1'*W(:,:,2)*h1+sigma_sq))/g_value1 ) == 0;
imag( (- h2'*W(:,:,2)*h2 - u4*(h2'*W(:,:,1)*h2+sigma_sq))/ g_value2 ) == 0;
cvx_end
imag(u4)==0
l(k)=real(u4)
p(k) = t;

if p(k) == 0
W_opt = W;
V_opt = V;
break;
else
minus_SINR1 = -h1'*W(:,:,1)*h1/(h1'*W(:,:,2)*h1+sigma_sq);
minus_SINR2 = -h2'*W(:,:,2)*h2/(h2'*W(:,:,1)*h2+sigma_sq);
u4 = max(minus_SINR1,minus_SINR2);
k = k+1;
end
g_value1 = h1'*W(:,:,2)*h1+g1'*V*g1+sigma_sq;
g_value2 = h2'*W(:,:,1)*h2+g2'*V*g2+sigma_sq;
end
z4(j)=u4
end
figure;
plot(1:length(z1),real(-z1),'-ro');
grid on;
hold on;
plot(1:length(z2),real(-z2),'-bo');
hold off;
hold on;
plot(1:length(z3),real(-z3),'-go');
hold off;
hold on;
plot(1:length(z4),real(-z4),'-mo');
hold off;