ZB & Yb
ZB & Yb
clc;
z1=0.2j;
y1=1/z1;
z2=0.3j;
y2=1/z2;
z3=0.2j;
y3=1/z3;
z4=0.2j;
y4=1/z4;
a=[-1 0 0;
1 -1 0;
0 1 -1;
1 0 -1];
y=[y1 0 0 0;
0 y2 0 0;
0 0 y3 0;
0 0 0 y4];
yb= transpose(a)*y*a
OUTPUT:
PROGRAM:
clear all;
clc;
zb=zeros(branch,branch);
[branch,branch] = size(zb) ;
for i=1:1:branch
zb(i,i)=zb(i,i)+z;
zb(i,i+1)=zb(i,i);
zb(i+1,i)=zb(i,i);
zb(i,i+2)=zb(i,i);
zb(i+2,i)=zb(i,i);
zb(i+1,i+1)=zb(i+1,i+1)+zb(i,i);
end
zb=zb(1:branch, 1:branch);
for l=1:1:link
a=(zb(i:j,i)-zb(i:j,j));
zb=zb-(1/(zb(i,i)+zb(j,j)-(2*zb(i,j))+z))*a*transpose(a);
end
bus_impedance_matrix=zb
OUTPUT: