Modern Power System Analysis with MATLAB® Applications by Pearson

Power System State Estimation and Optimal PMu Placement

431

%

Ybus = zeros(Nbus,Nbus);

for k = 1:Nele

p = Edata(k,2);

q = Edata(k,3);

y = 1/Edata(k,4);

shunts = Edata(k,5);

shuntr = Edata(k,6);

offa = Edata(k,7);

offaa = offa*offa;

Ybus(p,p) = Ybus(p,p) + y/offaa + shunts;

Ybus(q,q) = Ybus(q,q) + y + shuntr;

Ybus(p,q) = Ybus(p,q) - y/offa;

Ybus(q,p) = Ybus(q,p) - y/offa;

end

if Nsh~=0

for i = 1:Nsh

q = Shdata(i,1);

y = Shdata(i,2);

Ybus(q,q) = Ybus(q,q) + y;

end

end

Ybus

Gbus = real(Ybus);

Bbus = imag(Ybus);

% Recording element data

for i = 1:Nele

Fbus(i) = Edata(i,2);

Tbus(i) = Edata(i,3);

yele(i) = 1/Edata(i,4);

shunts(i) = Edata(i,5);

shuntr(i) = Edata(i,6); ...