summaryrefslogtreecommitdiff
path: root/14/CH7/EX7.1/example_7_1.sce
diff options
context:
space:
mode:
authorpriyanka2015-06-24 15:03:17 +0530
committerpriyanka2015-06-24 15:03:17 +0530
commitb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch)
treeab291cffc65280e58ac82470ba63fbcca7805165 /14/CH7/EX7.1/example_7_1.sce
downloadScilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip
initial commit / add all books
Diffstat (limited to '14/CH7/EX7.1/example_7_1.sce')
-rwxr-xr-x14/CH7/EX7.1/example_7_1.sce59
1 files changed, 59 insertions, 0 deletions
diff --git a/14/CH7/EX7.1/example_7_1.sce b/14/CH7/EX7.1/example_7_1.sce
new file mode 100755
index 000000000..9657f3cd5
--- /dev/null
+++ b/14/CH7/EX7.1/example_7_1.sce
@@ -0,0 +1,59 @@
+//chapter 7
+//Example 7.1
+//Page 170
+//busadmittancematrix
+clear;clc;
+//Given
+//Voltage Sources
+Ea = 1.5;
+Eb = 1.5*(cos(-36.87 * %pi / 180) + %i * sin(-36.87 * %pi / 180))
+Ec = 1.5;
+//admittances
+Ya = -%i*0.8;
+Yb = Ya;
+Yc= Ya;
+Yd = -%i*5;
+Ye = -%i*8;
+Yf = -%i*4;
+Yg = -%i*2.5;
+Yh = Yd;
+//current sourcs
+I1 = Ea * Ya;
+I2 = Eb * Yb;
+I3 = I1;
+I4 = 0;
+disp('Current Sources are')
+printf("\n I1 = - j%.2f per unit",-imag(I1))
+printf("\n I2 = %.2f - j%.2f per unit",real(I2),-imag(I2))
+printf("\n I3 = - j%.2f per unit \n\n",-imag(I3))
+//Self-admittances
+Y11 = Yd + Yf + Ya;
+Y22 = Yh + Yg + Yb;
+Y33 = Ye + Yc + Yg + Yf;
+Y44 = Yd + Ye + Yh;
+disp('Self-admittances are')
+printf("\n Y11 = - j%.2f per unit",-imag(Y11))
+printf("\n Y22 = - j%.2f per unit",-imag(Y22))
+printf("\n Y33 = - j%.2f per unit",-imag(Y33))
+printf("\n Y44 = - j%.2f per unit\n\n",-imag(Y44))
+//Mutual-admittances
+Y12 = 0;Y21 = Y12;
+Y13 = -Yf;Y31 = Y13;
+Y14 = -Yd;Y41 = Y14;
+Y23 = -Yg;Y32 = Y23;
+Y24 = -Yh;Y42 = Y24;
+Y34 = -Ye;Y43 = Y34;
+disp('Mutual admittances are')
+printf("\n Y12 = Y21 = %.2f per unit",imag(Y12))
+printf("\n Y13 = Y31 = j%.2f per unit",imag(Y13))
+printf("\n Y14 = Y41 = j%.2f per unit",imag(Y14))
+printf("\n Y23 = Y32 = j%.2f per unit",imag(Y23))
+printf("\n Y24 = Y42 = j%.2f per unit",imag(Y24))
+printf("\n Y34 = Y43 = j%.2f per unit \n\n",imag(Y34))
+//Matrix Form
+I = [I1 ; I2 ; I3 ; I4];
+Y = [Y11 Y12 Y13 Y14;Y21 Y22 Y23 Y24;Y31 Y32 Y33 Y34;Y41 Y42 Y43 Y44];
+disp('Current Vector =')
+disp(I)
+disp('Bus admittance matrix =')
+disp(Y) \ No newline at end of file