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
