The problem of Inaccurate/Infeasible in SDP


(Ruoyusun320) #1

Hello, everyone.
Here is the Matlab code for state estimation in power system.
I tried to solve the optimization problem through cvx with the solver mosey, but it shows that
Status: Inaccurate/Infeasible
Optimal value (cvx_optval): +Inf

Problem status: The problem is dual infeasible
Optimizer terminated. Time: 14.35

Could you give me some suggestions please?

    clear
    load  Y1.mat
    cvx_begin  sdp
    cvx_solver  mosek
    n = 132 ; Nbus = 44;
    variable  t(n) 
    variable  W(88,88)  semidefinite 

    R = [ones(44,1)*0.015;ones(44,1)*0.015;ones(58,1)*0.02];
    F = sum(t(n));
    subject to


    for ii=1:Nbus
        Pinj(ii,:)=[ii,trace(Yk{ii}*W)];
        Qinj(ii,:)=[ii ,trace(Ykbar{ii}*W)];
        Vbus(ii,:)=[ii, trace(M{ii}*W)];
    end

    Z1=[ones(Nbus,1),(1:Nbus)',zeros(Nbus,1),Pinj(:,2),0.015*ones(Nbus,1)];
    Z2=[2*ones(Nbus,1),(1:Nbus)',zeros(Nbus,1),Qinj(:,2),.015*ones(Nbus,1)];
    Z5=[5*ones(Nbus,1),Vbus(:,1),zeros(Nbus,1),Vbus(:,2),.01*ones(Nbus,1)];

    Z=[Z1;Z2;Z5];
    clear Pinj Qinj Vbus Pline Qline Sline
    for i = 1:n
            
         [R(i)*R(i)*t(i) Znoise(i,4)-Z(i,4);
             Znoise(i,4)-Z(i,4) 1] >= 0;
        
    end 
    W  >=  0;    
    cvx_end


(Ruoyusun320) #2