
Power System State Estimation and Optimal PMu Placement
441
sum = sum + VM(k)*(Gbus(bus,k)*sin(angle)-...
Bbus(bus,k)*cos(angle));
end
Bbus(bus,bus);
Hmat(nn,jjj) = sum - VM(bus)*Bbus(bus,bus);
end
if jj~=bus
angle = VA(bus)-VA(jj);
Hmat(nn,jjj) = VM(bus)*(Gbus(bus,jj)*sin(angle)-...
Bbus(bus,jj)*cos(angle));
end
end
end
if bus~=1
for jj = 2:Nbus
jjj = jj-1;
if jj==bus
sum = 0;
for k = 1:Nbus
angle = VA(bus)-VA(k);
sum = sum + VM(bus)*VM(k)*(Gbus(bus,k)*...
cos(angle)+Bbus(bus,k)*sin(angle));
end
Hmat(nn,jjj) = sum - VM(bus)*VM(bus)*Gbus(bus,bus);
end
if jj~=bus
angle = VA(bus)-VA(jj);
Hmat(nn,jjj) = VM(bus)*VM(jj)*(-Gbus(bus,jj)*...
cos(angle)-Bbus(bus,jj)*sin(angle)); ...