This is my code:

variable Psi_hat(M+1, M+1) hermitian

for k=1 : K

for j=1 : K

V1(k, j) =real(trace(V(:, :, k, j)* Psi_hat)) + (h_d(k, j))^2;

d1(k, j)= trans_power_ini(j)*V1(k, j);

d1_all = d1_all +d1(k, j);

end

end

for k=1 : K

D1(k) = (log(d1_all + (noise_maxpower^2)*(norm(W_ini(:, k)’, 2))^2))/(log(2)) ;

end

The problem I detected may be d1_ all is an affine function, while log () should convert the parentheses into concave functions, that is, V1 into concave functions, but I will not convert it. I have also referred to this article, Log( {complex affine} ) but I still haven’t solved it.

Please help me