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.7 | |
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.7')
-rwxr-xr-x | 14/CH7/EX7.7/example_7_7.sce | 66 | ||||
-rwxr-xr-x | 14/CH7/EX7.7/result_ex_7_7.txt | 26 |
2 files changed, 92 insertions, 0 deletions
diff --git a/14/CH7/EX7.7/example_7_7.sce b/14/CH7/EX7.7/example_7_7.sce new file mode 100755 index 000000000..ff0beda77 --- /dev/null +++ b/14/CH7/EX7.7/example_7_7.sce @@ -0,0 +1,66 @@ +//chapter 7 +//Example 7.7 +//Page 186 +//matrixmodification +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; +//Capacitor +Zb = -%i * 5 +//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; +Z = inv(Y); +disp('Original bus impedance matrix') +disp(Z) +[m,n] = size(Z) +for i = 1:m + for j = 1:n + Z(5,i) = Z(i,j); + Z(i,5) = Z(i,j) + end +end +Z(5,5) = Z(4,4) + Zb; +disp('Modified bus impedance matrix') +disp(Z) +[m1,n1] = size(Z); +Z_new = zeros(m1-1,n1-1); +for c = 1:m1-1 + for d = 1:n1-1 + Z_new(c,d) = Z(c,d) - ((Z(c,5)*Z(5,d)) / Z(5,5)); + end +end +disp('Modified bus impedance matrix after eliminating fifth row and column') +disp(Z_new) +V_4 = Z_new(4,:) * I; +printf("\n V4 = %.2f /_%.2f per unit \n\n",abs(V_4),atan(imag(V_4),real(V_4)) * 180 / %pi) +disp('V4 same as found in Example 7.6')
\ No newline at end of file diff --git a/14/CH7/EX7.7/result_ex_7_7.txt b/14/CH7/EX7.7/result_ex_7_7.txt new file mode 100755 index 000000000..96f733069 --- /dev/null +++ b/14/CH7/EX7.7/result_ex_7_7.txt @@ -0,0 +1,26 @@ + Original bus impedance matrix + + 0.4774405i 0.3705997i 0.4019598i 0.4142155i + 0.3705997i 0.4871730i 0.3922274i 0.4125934i + 0.4019598i 0.3922274i 0.4558128i 0.4231910i + 0.4142155i 0.4125934i 0.4231910i 0.4733096i + + Modified bus impedance matrix + + 0.4774405i 0.3705997i 0.4019598i 0.4142155i 0.4142155i + 0.3705997i 0.4871730i 0.3922274i 0.4125934i 0.4125934i + 0.4019598i 0.3922274i 0.4558128i 0.4231910i 0.4231910i + 0.4142155i 0.4125934i 0.4231910i 0.4733096i 0.4733096i + 0.4142155i 0.4125934i 0.4231910i 0.4733096i - 4.5266904i + + Modified bus impedance matrix after eliminating fifth row and column + + 0.5153434i 0.4083541i 0.4406840i 0.4575258i + 0.4083541i 0.5247796i 0.4307999i 0.4557341i + 0.4406840i 0.4307999i 0.4953761i 0.4674398i + 0.4575258i 0.4557341i 0.4674398i 0.5227987i + + V4 = 1.58 /_-11.97 per unit + + + V4 same as found in Example 7.6 |