1 clc;
2 clear all;
3
4 nb=4;
5 nl=6;
6 ld = [1 2 0 0.250
7 2 3 0 0.400
8 2 4 0 0.125
9 3 4 0 0.200
10 0 1 0 1.250
11 0 3 0 1.250
12 ]
13
14 Zbus= zeros(nb);
15
16 for i=1:nl
17 fb(i)= ld(i,1);
18 tb(i)= ld(i,2);
19 r(i)= ld(i,3);
20 x(i)= ld(i,4);
21 z(i)=complex(r(i),x(i));
22 end
23 s=zeros(1,nb);
24 for i=1:nl
25 if fb(i)==0
26 if s(tb(i))==0
27 Zbus(tb(i),tb(i))=Zbus(tb(i),tb(i))+z(i);
28 s(tb(i))=1;
29 else
30 q=tb(i);
31 delz=[Zbus(:,q)*Zbus(q,:)]/(Zbus(q,q)+z(i));
32 Zbus=Zbus-delz;
33 end
34 end
35 end
36 for i=1:nl
37 if fb(i)>0
38 p=fb(i);
39 q=tb(i);
40 if s(p)==1 && s(q)==0
41 Zbus(q,1:(q-1))=Zbus(p,1:(q-1));
42 Zbus(1:(q-1),q)=Zbus(1:(q-1),p);
43 Zbus(q,q)=Zbus(p,p)+z(i);
44 s(q)=1;
45 elseif s(p)==0 && s(q)==1
46 Zbus(p,1:(p-1))=Zbus(q,1:(p-1));
47 Zbus(1:(p-1),p)=Zbus(1:(p-1),q);
48 Zbus(p,p)=Zbus(q,q)+z(i);
49 s(p)=1;
50 elseif s(p)==1 && s(q)==1
51 temp=Zbus(p,p)+Zbus(q,q)-2*Zbus(p,q)+z(i);
52 Zbus=Zbus-[[Zbus(:,p)-Zbus(:,q)]*[Zbus(p,:)-Zbus(q,:)]]/temp;
53 end
54 end
55 end
56
57 display(Zbus);
58
59 Z_addition = Zbus;
60 l_addition = [ 1 4 0.250 0.4]
61 p=l_addition(1);
62 q=l_addition(2);
63 r=l_addition(3);
64 x=l_addition(4);
65 z=complex(r,x);
66 temp=Z_addition(p,p)+Z_addition(q,q)-2*Z_addition(p,q)+z;
67 Z_addition=Z_addition-[[Z_addition(:,p)-Z_addition(:,q)]*[Z_addition(p,:)-Z_addition(q,:)]]/
temp;
68 display(Z_addition);
69
70 Z_subtraction=Zbus;
71 l_subtraction = [ 2 3 0.3 0.750]
72 p=l_subtraction(1);
73 q=l_subtraction(2);
74 r=l_subtraction(3);
75 x=l_subtraction(4);
76 z=complex(r,x);
77 temp=Z_subtraction(p,p)+Z_subtraction(q,q)-2*Z_subtraction(p,q)-z;
78 Z_subtraction=Z_subtraction-[[Z_subtraction(:,p)-Z_subtraction(:,q)]*[Z_subtraction(p,:)-
Z_subtraction(q,:)]]/temp;
79 display(Z_subtraction)
