Sorry to bother you and this is my code:

% Max-min power control

R_dk_min = min(R_dk);

t_min = 2^R_dk_min-1;

t_max = 2.5;

epsi = max(t_min/5,0.01);

v_k1_phi = zeros(K,K);

I_k = abs(eye(num_UEs)-1);

for k=1:K

for k_dash=1:K

v_k1_phi(k_dash,k)=phi_k(:,k_dash)’*phi_k(:,k);

end

end

v_k1_phi = v_k1_phi.*I_k;

while(t_max-t_min>epsi)

t=(t_max+t_min)/2;

cvx_begin sdp quiet

variables var_sigma(M,K) var_rho(K,K) var_theta(M,1)

minimize(0)

subject to

for k=1:K

norm([v_k1_phi(:,k).*var_rho(:,k);sqrt(beta_mk(:,k)).*var_theta;1/sqrt(rho_d)]) <= 1/sqrt(t)*gamma_mk(:,k)’*var_sigma(:,k);

end

for m=1:M

norm(sqrt(gamma_mk(m,:)).*var_sigma(m,:)) <= var_theta(m);

var_theta(m) <= 1;

var_theta(m) >= 0;

end

for k=1:K

for k_dash=1:K

sum(gamma_mk(:,k_dash).*beta_mk(:,k)./beta_mk(:,k_dash).*var_sigma(:,k_dash)) <= var_rho(k_dash,k);

end

end

for m=1:M

for k=1:K

var_sigma(m,k) >= 0;

end

end

cvx_end

% Feasible bisection

if contains(cvx_status,‘Solved’) % feasible

t_min = t;

else % not feasible

t_max = t;

end

toc

end

R_dk_opt_min(round)=log2(1+t_min);