summaryrefslogtreecommitdiff
path: root/83/CH9/EX9.8/example_9_8.sce
diff options
context:
space:
mode:
Diffstat (limited to '83/CH9/EX9.8/example_9_8.sce')
-rwxr-xr-x83/CH9/EX9.8/example_9_8.sce32
1 files changed, 32 insertions, 0 deletions
diff --git a/83/CH9/EX9.8/example_9_8.sce b/83/CH9/EX9.8/example_9_8.sce
new file mode 100755
index 000000000..c4b0ef7c7
--- /dev/null
+++ b/83/CH9/EX9.8/example_9_8.sce
@@ -0,0 +1,32 @@
+//Chapter 9
+//Example 9.8
+//page 360
+//To build Zbus matrix using Zbus building algorithm
+clear;clc;
+
+disp("Let us go on modifying Zbus by including nodes and the elements connected to it one by one as given in the textbook")
+
+//step-1 type1 modification
+Zbus=[0.25];
+printf('\nstep-1 type1 modification\nZbus=');disp(Zbus);
+
+//step-2 type2 modification
+Zbus=[Zbus,0.25;0.25,0.25+0.1];
+printf('\nstep-2 type2 modification\nZbus=');disp(Zbus);
+
+//step-3 type2 modification
+Zbus=[Zbus [0.25;0.25]; 0.25 0.25 0.35];
+printf('\nstep-3 type2 modification\nZbus=');disp(Zbus);
+
+//step-4 type3 modification
+Zbus=Zbus-(1/(Zbus(3,3)+0.25))*[Zbus(1:3,2:2)]*[Zbus(2:2,1:3)];
+printf('\nstep-4 type3 modification\nZbus=');disp(Zbus);
+
+printf('This is the final Zbus matrix after including all the elements\n');
+//step-5 type4 modification
+Zbus=Zbus-(1/(0.1+Zbus(2,2)+Zbus(3,3)-2*Zbus(2,3)))*[Zbus(1:3,2:2)-Zbus(1:3,3:3)]*[Zbus(2:2,1:3)-Zbus(3:3,1:3)];
+printf('\nstep-5 type4 modification\nZbus=');disp(Zbus);
+
+disp("opening a line between 2-3 is equivalent to connecting (-0.1) between bus3 bus2")
+Zbus=Zbus-(1/(-0.1+Zbus(2,2)+Zbus(3,3)-2*Zbus(2,3)))*[Zbus(1:3,2:2)-Zbus(1:3,3:3)]*[Zbus(2:2,1:3)-Zbus(3:3,1:3)];
+printf('Zbus=');disp(Zbus);printf('(same as in step 4)'); \ No newline at end of file