summaryrefslogtreecommitdiff
path: root/3793/CH6
diff options
context:
space:
mode:
Diffstat (limited to '3793/CH6')
-rw-r--r--3793/CH6/EX6.1/exp_6_1.sce17
-rw-r--r--3793/CH6/EX6.3/exp_6_3.sce31
-rw-r--r--3793/CH6/EX6.4/exp_6_4.sce18
-rw-r--r--3793/CH6/EX6.5/exp_6_5.sce62
-rw-r--r--3793/CH6/EX6.6/exp_6_6.sce52
-rw-r--r--3793/CH6/EX6.7/exp_6_7.sce67
-rw-r--r--3793/CH6/EX6.8/exp_6_8.sce71
7 files changed, 318 insertions, 0 deletions
diff --git a/3793/CH6/EX6.1/exp_6_1.sce b/3793/CH6/EX6.1/exp_6_1.sce
new file mode 100644
index 000000000..8ada421e5
--- /dev/null
+++ b/3793/CH6/EX6.1/exp_6_1.sce
@@ -0,0 +1,17 @@
+clear;
+clc;
+A=[-1 0 0; 0 0 -1; 1 -1 0; 0 -1 1; -1 0 1];
+B=(A');
+z=diag([.05;.10;.5;.40;.25]);
+y=pinv(z);
+Yb=(B*y*A);
+mprintf("Ybus matrix without coupling\n");
+disp(Yb);
+
+//case2
+z1=[.05 0 0 0 0; 0 .10 0 0 0; 0 0 .5 0 0; 0 0 0 .4 .2;0 0 0 .2 .25];
+y1=pinv(z1);
+Y1b=B*y1*A;
+mprintf("\nYbus matrix with coupling\n");
+disp(Y1b);
+Zb=pinv(Y1b);
diff --git a/3793/CH6/EX6.3/exp_6_3.sce b/3793/CH6/EX6.3/exp_6_3.sce
new file mode 100644
index 000000000..52097172a
--- /dev/null
+++ b/3793/CH6/EX6.3/exp_6_3.sce
@@ -0,0 +1,31 @@
+clear;
+clc;
+y12=complex(.0059,-.0235);
+y14=complex(.0055,-.0183);
+y23=complex(.0077,-.0385);
+y24=complex(.0240,-.0320);
+y34=complex(.0100,-.0300);
+y40=complex(.0100,-.0200);
+Y11=y12+y14;
+Y22=y12+y23+y24;
+Y33=y23+y34;
+Y44=y14+y24+y34+y40;
+Y13=0;
+Y31=0;
+Y12=-y12;
+Y21=-y12;
+Y14=-y14;
+Y41=-y14;
+Y23=-y23;
+Y32=-y23;
+Y24=-y24;
+Y42=-y24;
+Y34=-y34;
+Y43=-y34;
+Yb=[Y11 Y12 Y13 Y14; Y21 Y22 Y23 Y24; Y31 Y32 Y33 Y34; Y41 Y42 Y43 Y44];
+mprintf(" Ybus matrix is\n");
+disp(Yb);
+
+
+
+
diff --git a/3793/CH6/EX6.4/exp_6_4.sce b/3793/CH6/EX6.4/exp_6_4.sce
new file mode 100644
index 000000000..82907e511
--- /dev/null
+++ b/3793/CH6/EX6.4/exp_6_4.sce
@@ -0,0 +1,18 @@
+clear;
+clc;
+y=[4 3 6;2 8 5;1 5 9];
+nbus=3;
+l11=4;
+disp(y(1,2));
+u12=(1/l11)*y(1,2);
+u13=(1/l11)*y(1,3);
+l21=2;
+l22=y(2,2)-(l21*u12);
+u23=(1/l22)*(y(2,3)-(l21*u13));
+l31=1;
+l32=y(3,2)-(l31*u12);
+l21=2;
+l33=y(3,3)-(l31*u13)-(l32*u23);
+Yb=[l11 u12 u13; l21 l22 u23;l31 l32 l33];
+mprintf("matrix for table of factor is\n")
+disp(Yb);
diff --git a/3793/CH6/EX6.5/exp_6_5.sce b/3793/CH6/EX6.5/exp_6_5.sce
new file mode 100644
index 000000000..519654d61
--- /dev/null
+++ b/3793/CH6/EX6.5/exp_6_5.sce
@@ -0,0 +1,62 @@
+clear;
+clc;
+
+function[Ybus,I]=fbsub(Ybus,nbus,I);
+ for k=1:nbus;
+ if k==1;
+ for j=2:nbus;
+ Ybus(k,j)=Ybus(k,j)/Ybus(k,k);
+ end
+else
+ for j=2:nbus;
+ if j<=k;
+ for m=1:j-1;
+ Ybus(k,j)=Ybus(k,j)-Ybus(k,m)*Ybus(m,j);
+
+ end
+else
+ for m=1:k-1;
+ Ybus(k,j)=Ybus(k,j)-Ybus(k,m)*Ybus(m,j);
+
+ end
+ Ybus(k,j)=Ybus(k,j)/Ybus(k,k);
+end
+end
+end
+end
+for k=1:nbus;
+ if k==1;
+ I(k)=I(k)/Ybus(k,k);
+ else
+ for j=1:k-1;
+ I(k)=I(k)-Ybus(k,j)*I(j);
+end
+I(k)=I(k)/Ybus(k,k);
+end
+end
+for k=nbus:-1:1;
+ if k==nbus;
+ disp('node voltages');
+ disp(Ybus);
+
+
+ else
+ for j=nbus:-1:k+1;
+ I(k)=I(k)-Ybus(k,j)*I(j);
+end
+end
+end
+endfunction
+Ybus=[4 3 6;2 8 5;1 5 9];
+nbus=3;
+I=[1;1;1];
+[Ybus,I]=fbsub(Ybus,nbus,I);
+V1=1/Ybus(1,1);
+V2=(1/Ybus(2,2))*(1-2*V1);
+V3=(1/Ybus(3,3))*(1-1*V1-4.25*V2);
+VV3=V3;
+VV2=(V2-Ybus(2,3)*V3);
+VV1=(V1-Ybus(1,2)*VV2-Ybus(1,3)*V3);
+V=[VV1 ; VV2 ;VV3]
+disp("V is");
+disp(V);
diff --git a/3793/CH6/EX6.6/exp_6_6.sce b/3793/CH6/EX6.6/exp_6_6.sce
new file mode 100644
index 000000000..5a4c85cb4
--- /dev/null
+++ b/3793/CH6/EX6.6/exp_6_6.sce
@@ -0,0 +1,52 @@
+clear;
+clc;
+Zb=complex(0,.20);
+Z1b=complex(0,.25);
+V1=Z1b;
+
+Z2b=[Z1b 0;0 Z1b];
+//p=2,q=3
+Zbn=complex(0,.10);
+Z22=complex(0,.25);
+Z33=Z22+Zbn;
+Z13=complex(0,0);
+Z31=Z13;
+Z32=complex(0,.25);
+Z23=Z32;
+
+Z2bn=[Z1b 0 Z13;0 Z22 Z23;Z31 Z32 Z33];
+//p=1,q=4
+Z44=Z1b+Zb;
+Z14=complex(0,.25);
+Z41=Z14;
+Z24=complex(0,0);
+Z42=complex(0,0);
+Z34=complex(0,0);
+Z43=Z34;
+Z4b=[Z1b 0 Z13 Z14;0 Z22 Z23 Z24;Z31 Z32 Z33 Z34;Z41 Z42 Z43 Z44];
+//p=3,q=4
+Zbb=complex(0,.3);
+y=1/((Zbb+Z33)-(2*Z34)+Z44);
+Z=[Z13-Z14;Z23-Z24;Z33-Z34;Z43-Z44];
+Z5b=Z4b-((-y)*Z*(Z'));
+//p=0,q=3
+y1=1/(Zbn+Z5b(3,3));
+ZZ=[-Z5b(1,3);-Z5b(2,3);-Z5b(3,3);-Z5b(4,3)];
+Z6b=Z5b-((-y1)*ZZ*(ZZ'));
+mprintf("Z1b = \n");
+disp(Z1b);
+mprintf("Z2b = \n");
+disp(Z2bn);
+
+mprintf("Z4b = \n");
+disp(Z4b);
+mprintf("Z5b = \n");
+disp(Z5b);
+mprintf("Z6b = \n");
+disp(Z6b);
+
+
+
+
+
+
diff --git a/3793/CH6/EX6.7/exp_6_7.sce b/3793/CH6/EX6.7/exp_6_7.sce
new file mode 100644
index 000000000..b97351b7c
--- /dev/null
+++ b/3793/CH6/EX6.7/exp_6_7.sce
@@ -0,0 +1,67 @@
+clear;
+clc;
+function [Zbus]=zeebus(nelemnt,ind,node);
+ if ind==1;
+ Zbus=input('Partial matrix Zbus');
+ ind=0;
+ else
+ end
+ if ind==0;
+ for l=1:nelemnt;
+ p=input('bus number p');
+ q=input('bus number q');
+ zb=input('impedence');
+ typee=input('type of bus');
+ if typee==1;
+ for k=1:q;
+ if k==q;
+ Zbus(k,k)=zb;
+ else
+ Zbus(k,q)=0;
+ end
+ end
+ end
+ if typee==2;
+ for k=1:q;
+ if k==q;
+ Zbus(q,q)=zb+Zbus(p,p);
+ else
+ if k==p;
+ Zbus(p,q)=Zbus(p,p);
+ Zbus(q,p)=Zbus(p,q);
+ else
+ Zbus(k,q)=0;
+ Zbus(q,k)=0;
+ end
+ end
+ end
+ end
+ if typee==3;
+ y=1/(zb+Zbus(p,p)-2*Zbus(p,q)+Zbus(q,q));
+ for k=1:node;
+ X(k,1)=Zbus(k,p)-Zbus(k,q);
+ Xt(1,k)=(k:1);
+ end
+ Zbus=Zbus-(-y)*X*X';
+ end
+ if typee==4;
+ y=1/(zb+Zbus(q,q));
+ for k=1:node;
+ X(k,1)=-Zbus(k,q);
+ Xt(1,k)=(k:1);
+ end
+ Zbus=Zbus-(-y)*X*X';
+ end
+ end
+ else
+
+
+ end
+
+endfunction
+nelemnt=7;
+ind=0;
+node=4;
+[Zbus]=zeebus(nelemnt,ind,node);
+mprintf("Zbus = \n");
+disp(Zbus);
diff --git a/3793/CH6/EX6.8/exp_6_8.sce b/3793/CH6/EX6.8/exp_6_8.sce
new file mode 100644
index 000000000..7bda4e775
--- /dev/null
+++ b/3793/CH6/EX6.8/exp_6_8.sce
@@ -0,0 +1,71 @@
+clear;
+clc;
+function [Zbus]=zeebus(nelemnt,ind,node);
+ if ind==1;
+ Zbus=input('Partial matrix Zbus');
+ ind=0;
+ else
+ end
+ if ind==0;
+ for l=1:nelemnt;
+ p=input('bus number p');
+ q=input('bus number q');
+ zb=input('impedence');
+ typee=input('type of bus');
+ if typee==1;
+ for k=1:q;
+ if k==q;
+ Zbus(k,k)=zb;
+ else
+ Zbus(k,q)=0;
+ end
+ end
+ end
+ if typee==2;
+ for k=1:q;
+ if k==q;
+ Zbus(q,q)=zb+Zbus(p,p);
+ else
+ if k==p;
+ Zbus(p,q)=Zbus(p,p);
+ Zbus(q,p)=Zbus(p,q);
+ else
+ Zbus(k,q)=0;
+ Zbus(q,k)=0;
+ end
+ end
+ end
+ end
+ if typee==3;
+ y=1/(zb+Zbus(p,p)-2*Zbus(p,q)+Zbus(q,q));
+ for k=1:node;
+ X(k,1)=Zbus(k,p)-Zbus(k,q);
+ Xt(1,k)=(k:1);
+ end
+ Zbus=Zbus-(-y)*X*X';
+ end
+ if typee==4;
+ y=1/(zb+Zbus(q,q));
+ for k=1:node;
+ X(k,1)=-Zbus(k,q);
+ Xt(1,k)=(k:1);
+ end
+ Zbus=Zbus-(-y)*X*X';
+ end
+ end
+ else
+
+
+ end
+
+endfunction
+nelemnt=7;
+ind=0;
+node=4;
+[Zbus]=zeebus(nelemnt,ind,node);
+mprintf("Zbus = \n");
+disp(Zbus);
+nelemnt=1;
+ind=1;
+[Zbus]=zeebus(nelemnt,ind,node);
+disp(Zbus);
generated by cgit (git 2.34.1) at 2025-03-08 08:09:25 +0000