diff options
author | priyanka | 2015-06-24 15:03:17 +0530 |
---|---|---|
committer | priyanka | 2015-06-24 15:03:17 +0530 |
commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
tree | ab291cffc65280e58ac82470ba63fbcca7805165 /14/CH7/EX7.5/example_7_5.sce | |
download | Scilab-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.5/example_7_5.sce')
-rwxr-xr-x | 14/CH7/EX7.5/example_7_5.sce | 52 |
1 files changed, 52 insertions, 0 deletions
diff --git a/14/CH7/EX7.5/example_7_5.sce b/14/CH7/EX7.5/example_7_5.sce new file mode 100755 index 000000000..fca56d5d3 --- /dev/null +++ b/14/CH7/EX7.5/example_7_5.sce @@ -0,0 +1,52 @@ +//chapter 7 +//Example 7.5 +//Page 181 +//introcapacitor +clear;clc; +//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; +//Value of capacitor introduced in node 4 +C = 5.0 ;//in per unit +Xc = %i*C; +//current sourcs +I1 = Ea * Ya; +I2 = Eb * Yb; +I3 = I1; +I4 = 0; +//Self-admittances +Y11 = Yd + Yf + Ya; +Y22 = Yh + Yg + Yb; +Y33 = Ye + Yc + Yg + Yf; +Y44 = Yd + Ye + Yh; +//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; +//Matrix Form +I = [I1 ; I2 ; I3 ; I4]; +Y = [Y11 Y12 Y13 Y14;Y21 Y22 Y23 Y24;Y31 Y32 Y33 Y34;Y41 Y42 Y43 Y44]; +V = Y\I; +E_th = V(4,1); +Z = inv(Y); +Z_th = Z(4,4); +I_c = E_th / (Z_th - Xc); +disp('Thevenin equivalent of the circuit behind node four') +printf("\n Eth = %.2f /_%.2f per unit \n\n",abs(E_th),atan(imag(E_th),real(E_th)) * 180 / %pi) +disp('Thevenin equivalent impedance') +printf("\n Z_th = j%.2f per unit \n\n",imag(Z_th)) +disp('Current drawn by the capacitor') +printf("\n Ic = %.2f /_%.2f per unit \n\n",abs(I_c),atan(imag(I_c),real(I_c)) * 180 / %pi)
\ No newline at end of file |