Newtonn Raphson Method Metleb 3 Bus

% Question #5A Newton-Raphson method
clear; clc;
V = [1.0; 1.03; 1.0];
d = [0; 0; 0];
Ps=[4 ; -5];
Qs= -3.5;
YB = [ -j*75 j*50 j*25
j*50 -j*75 j*25
j*25 j*25 -j*50];
Y= abs(YB); t = angle(YB);
iter=0;
while iter < 2
iter = iter +1
P=[V(2)*V(1)*Y(2,1)*cos(t(2,1)-d(2)+d(1))+V(2)^2*Y(2,2)*cos(t(2,2))+V(2)*V(3)*Y(2,3)*cos(t(2,3)-d(2)+d(3));
V(3)*V(1)*Y(3,1)*cos(t(3,1)-d(3)+d(1))+V(3)^2*Y(3,3)*cos(t(3,3))+V(3)*V(2)*Y(3,2)*cos(t(3,2)-d(3)+d(2))];
Q= -V(3)*V(1)*Y(3,1)*sin(t(3,1)-d(3)+d(1))-V(3)^2*Y(3,3)*sin(t(3,3))-V(2)*V(3)*Y(3,2)*sin(t(3,2)-d(3)+d(2));
J(1,1)=V(2)*V(1)*Y(2,1)*sin(t(2,1)-d(2)+d(1))+V(2)*V(3)*Y(2,3)*sin(t(2,3)-d(2)+d(3));
J(1,2)=-V(2)*V(3)*Y(2,3)*sin(t(2,3)-d(2)+d(3));
J(1,3)=V(2)*Y(2,3)*cos(t(2,3)-d(2)+d(3));
J(2,1)=-V(3)*V(2)*Y(3,2)*sin(t(3,2)-d(3)+d(2));
J(2,2)=V(3)*V(1)*Y(3,1)*sin(t(3,1)-d(3)+d(1))+V(3)*V(2)*Y(3,2)*sin(t(3,2)-d(3)+d(2));
J(2,3)=V(1)*Y(3,1)*cos(t(3,1)-d(3)+d(1))+V(2)*Y(3,2)*cos(t(3,2)-d(3)+d(2));
J(3,1)=-V(3)*V(2)*Y(3,2)*cos(t(3,2)-d(3)+d(2));
J(3,2)=V(2)*V(3)*Y(3,2)*cos(t(3,2)-d(3)+d(2))+V(1)*V(3)*Y(3,1)*cos(t(3,1)-d(3)+d(1));
J(3,3)=-V(1)*Y(3,1)*sin(t(3,1)-d(3)+d(1))-2*V(3)*Y(3,3)*sin(t(3,3))-V(2)*Y(3,2)*sin(t(3,2)-d(3)+d(2));
DP = Ps-P;
DQ = Qs-Q;
DC = [DP; DQ]
J
DX = J\DC
d(2) =d(2)+DX(1);
d(3)=d(3) +DX(2);
V(3)= V(3)+DX(3);
V, d, delta =180/pi*d;
end
P1= V(1)^2*Y(1,1)*cos(t(1,1))+V(1)*V(2)*Y(1,2)*cos(t(1,2)-d(1)+d(2))+V(1)*V(3)*Y(1,3)*cos(t(1,3)-d(1)+d(3))
Q1=-V(1)^2*Y(1,1)*sin(t(1,1))-V(1)*V(2)*Y(1,2)*sin(t(1,2)-d(1)+d(2))-V(1)*V(3)*Y(1,3)*sin(t(1,3)-d(1)+d(3))
Q2=-V(2)*V(1)*Y(2,1)*sin(t(2,1)-d(2)+d(1))-V(3)*V(2)*Y(2,3)*sin(t(2,3)-d(2)+d(3))-V(2)^2*Y(2,2)*sin(t(2,2))
P_loss = P1+4-5
Q_loss = Q1+Q2-3.5

Catatan: Gunakan Noteped sebelum dimasukan kedalam Metleb
Sumber: http://imadudd1n.wordpress.com/tag/newton-raphson/