Sqlp stop: primal problem is suspected of being infeasible

1

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

微信图片_20240205120144
微信图片_202402051201441394×2241 169 KB

2

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.

3

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