# Sqlp stop: primal problem is suspected of being infeasible

## I have a problem here . How can I solve it

code:

cvx_begin

variable omegai1(M,1) complex
variable omegae1(M,1) complex
variable omegai2(M,1) complex
variable omegae2(M,1) complex
variable pk1 nonnegative
variable pk2 nonnegative
variable pk3 nonnegative
variable x1 nonnegative
variable x2 nonnegative

maximize(min(c1+d1*(2-((norm(hAU1*omegail1))^2)*inv_pos(x1)-(pow_abs(hAU1*omegai2,2)+pow_abs(hAU1*omegae1,2)+pow_abs(hAU1*omegae2,2)+pk1*(norm(hDT1IU1))^2+pk2*(norm(hDT2IU1))^2+pk3*(norm(hDT3IU1))^2+sigmaK2)...
/((norm(hAU1*omegail2))^2+(norm(hAU1*omegael1))^2+(norm(hAU1*omegael2))^2+pkl(1)*(norm(hDT1IU1))^2+pkl(2)*(norm(hDT2IU1))^2+pkl(3)*(norm(hDT3IU1))^2+sigmaK2))-eRT,...
c2+d2*(2-((norm(hAU2*omegail2))^2)*inv_pos(x2)-(pow_abs(hAU2*omegai1,2)+pow_abs(hAU2*omegae1,2)+pow_abs(hAU2*omegae2,2)+pk1*(norm(hDT1IU2))^2+pk2*(norm(hDT2IU2))^2+pk3*(norm(hDT3IU2))^2+sigmaK2)...
/((norm(hAU2*omegail1))^2+(norm(hAU2*omegael1))^2+(norm(hAU2*omegael2))^2+pkl(1)*(norm(hDT1IU2))^2+pkl(2)*(norm(hDT2IU2))^2+pkl(3)*(norm(hDT3IU2))^2+sigmaK2))-eRT))

subject to
2*real((hAU1*omegai1)*(hAU1*omegail1)')-(norm(hAU1*omegail1))^2>=x1;
2*real((hAU2*omegai2)*(hAU2*omegail2)')-(norm(hAU2*omegail2))^2>=x2;
%(39b)
min((2*real((hAE1*omegai1)*(hAE1*omegail1)')-(norm(hAE1*omegail1))^2+...
2*real((hAE1*omegai2)*(hAE1*omegail2)')-(norm(hAE1*omegail2))^2+...
2*real((hAE1*omegae1)*(hAE1*omegael1)')-(norm(hAE1*omegael1))^2+...
2*real((hAE1*omegae2)*(hAE1*omegael2)')-(norm(hAE1*omegael2))^2+...
pk1*(norm(hDT1EU1))^2+pk2*(norm(hDT2EU1))^2+pk3*(norm(hDT3EU1))^2), ...
(2*real((hAE2*omegai1)*(hAE2*omegail1)')-(norm(hAE2*omegail1))^2+...
2*real((hAE2*omegai2)*(hAE2*omegail2)')-(norm(hAE2*omegail2))^2+...
2*real((hAE2*omegae1)*(hAE2*omegael1)')-(norm(hAE2*omegael1))^2+...
2*real((hAE2*omegae2)*(hAE2*omegael2)')-(norm(hAE2*omegael2))^2+...
pk1*(norm(hDT1EU2))^2+pk2*(norm(hDT2EU2))^2+pk3*(norm(hDT3EU2))^2)) >= ET;
%(39c)
%(RD1)
(a1+b1*(2-pkl(1)*inv_pos(pk1)-(pow_abs(hADR1*omegai1,2)+pow_abs(hADR1*omegae1,2)+pow_abs(hADR1*omegai2,2)+pow_abs(hADR1*omegae2,2)+pk2*(norm(hDT2DR1))^2+pk3*(norm(hDT3DR1))^2+sigmaK2)...
/((norm(hADR1*omegail1))^2+(norm(hADR1*omegael1))^2+(norm(hADR1*omegail2))^2+(norm(hADR1*omegael2))^2+pkl(2)*(norm(hDT2DR1))^2+pkl(3)*(norm(hDT3DR1))^2+sigmaK2))) >= RT;
%(RD2)
(a2+b2*(2-pkl(2)*inv_pos(pk2)-(pow_abs(hADR2*omegai1,2)+pow_abs(hADR2*omegae1,2)+pow_abs(hADR2*omegai2,2)+pow_abs(hADR2*omegae2,2)+pk1*(norm(hDT1DR2))^2+pk3*(norm(hDT3DR2))^2+sigmaK2)...
/((norm(hADR2*omegail1))^2+(norm(hADR2*omegael1))^2+(norm(hADR2*omegail2))^2+(norm(hADR2*omegael2))^2+pkl(1)*(norm(hDT1DR2))^2+pkl(3)*(norm(hDT3DR2))^2+sigmaK2))) >= RT;
%(RD3)
(a3+b3*(2-pkl(3)*inv_pos(pk3)-(pow_abs(hADR3*omegai1,2)+pow_abs(hADR3*omegae1,2)+pow_abs(hADR3*omegai2,2)+pow_abs(hADR3*omegae2,2)+pk1*(norm(hDT1DR3))^2+pk2*(norm(hDT2DR3))^2+sigmaK2)...
/((norm(hADR3*omegail1))^2+(norm(hADR3*omegael1))^2+(norm(hADR3*omegail2))^2+(norm(hADR3*omegael2))^2+pkl(1)*(norm(hDT1DR3))^2+pkl(2)*(norm(hDT2DR3))^2+sigmaK2))) >= RT;
%(33d)
(omegai1'*omegai1+omegae1'*omegae1+omegai2'*omegai2+omegae2'*omegae2) <= P;
%(33e)
pk1 <= pT;
pk2 <= pT;
pk3 <= pT;
%(39f)
(pow_abs(hAE1*omegai1,2)-(exp(eRT)-1)*(2*real((hAE1*omegai2)*(hAE1*omegail2)')-(norm(hAE1*omegail2))^2+...
2*real((hAE1*omegae1)*(hAE1*omegael1)')-(norm(hAE1*omegael1))^2+...
2*real((hAE1*omegae2)*(hAE1*omegael2)')-(norm(hAE1*omegael2))^2+...
pk1*(norm(hDT1EU1))^2+pk2*(norm(hDT2EU1))^2+pk3*(norm(hDT3EU1))^2+sigmaK2)) <= 0;
(pow_abs(hAE1*omegai2,2)-(exp(eRT)-1)*(2*real((hAE1*omegai1)*(hAE1*omegail1)')-(norm(hAE1*omegail1))^2+...
2*real((hAE1*omegae1)*(hAE1*omegael1)')-(norm(hAE1*omegael1))^2+...
2*real((hAE1*omegae2)*(hAE1*omegael2)')-(norm(hAE1*omegael2))^2+...
pk1*(norm(hDT1EU1))^2+pk2*(norm(hDT2EU1))^2+pk3*(norm(hDT3EU1))^2+sigmaK2)) <= 0;
(pow_abs(hAE2*omegai1,2)-(exp(eRT)-1)*(2*real((hAE2*omegai2)*(hAE2*omegail2)')-(norm(hAE2*omegail2))^2+...
2*real((hAE2*omegae1)*(hAE2*omegael1)')-(norm(hAE2*omegael1))^2+...
2*real((hAE2*omegae2)*(hAE2*omegael2)')-(norm(hAE2*omegael2))^2+...
pk1*(norm(hDT1EU2))^2+pk2*(norm(hDT2EU2))^2+pk3*(norm(hDT3EU2))^2+sigmaK2)) <= 0;
(pow_abs(hAE2*omegai2,2)-(exp(eRT)-1)*(2*real((hAE2*omegai1)*(hAE2*omegail1)')-(norm(hAE2*omegail1))^2+...
2*real((hAE2*omegae1)*(hAE2*omegael1)')-(norm(hAE2*omegael1))^2+...
2*real((hAE2*omegae2)*(hAE2*omegael2)')-(norm(hAE2*omegael2))^2+...
pk1*(norm(hDT1EU2))^2+pk2*(norm(hDT2EU2))^2+pk3*(norm(hDT3EU2))^2+sigmaK2)) <= 0;
cvx_end

All but section 1 of Debugging infeasible models - YALMIP also applies to CVX.

I don 't know whether the large magnitude primal and dual objective values are indicative of bad numerical scaling. In any event, you should first make sure there are no non-zero input data that are not within a small number of orders of magnitude of 1. If there are any, you should change the units so that non-zero input data are within a small number of orders of magnitude of 1. The re-run the problem and follow advice in the link as needed.

Thank you very much for your reply, I will try to modify it as you suggested