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 /1076/CH10/EX10.16/10_16.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 '1076/CH10/EX10.16/10_16.sce')
| -rwxr-xr-x | 1076/CH10/EX10.16/10_16.sce | 58 |
1 files changed, 58 insertions, 0 deletions
diff --git a/1076/CH10/EX10.16/10_16.sce b/1076/CH10/EX10.16/10_16.sce new file mode 100755 index 000000000..1ce736714 --- /dev/null +++ b/1076/CH10/EX10.16/10_16.sce @@ -0,0 +1,58 @@ +clear
+clc
+
+X=[
+ .25 .25 .05
+ .2 .2 .05
+ .06 .06 .06
+ .07 .07 .07
+ .1 .1 .3
+ .1 .1 .3
+]
+
+B=[
+ 100 11
+ 100 11
+ 100 11
+ 100 11
+ 100 220
+ 100 220
+ ]
+V1=11
+V2=220
+S=100
+Xe=3*.03
+//end 9
+X1=(((X(1,1)*B(1,1) *V1/(S*B(1,2))) +(X(3,1)*B(3,1) *V1/(S*B(3,2))) )^-1+((X(2,1)*B(2,1) *V1/(S*B(2,2)))+(X(4,1)*B(4,1) *V1/(S*B(4,2))) +((X(5,1)*B(5,1) *V2/(S*B(5,2)))^-1 +(X(6,1)*B(6,1) *V2/(S*B(6,2)))^-1)^-1)^-1)^-1
+
+X2=(((X(1,2)*B(1,1) *V1/(S*B(1,2))) +(X(3,2)*B(3,1) *V1/(S*B(3,2))) )^-1+((X(2,2)*B(2,1) *V1/(S*B(2,2)))+(X(4,2)*B(4,1) *V1/(S*B(4,2))) +((X(5,2)*B(5,1) *V2/(S*B(5,2)))^-1 +(X(6,2)*B(6,1) *V2/(S*B(6,2)))^-1)^-1)^-1)^-1
+
+X0=((X(3,3)*B(3,1) *V1/(S*B(3,2)))^-1 + ((Xe *B(4,1) *V1/(S*B(4,2))) + (X(2,3)*B(2,1) *V1/(S*B(2,2))) + (X(4,3)*B(4,1) *V1/(S*B(4,2))) +(((X(5,3)*B(5,1) *V2/(S*B(5,2))))^-1 +(X(6,3)*B(6,1) *V2/(S*B(6,2)))^-1)^-1 )^-1 )^-1
+
+Z1=%i * round(X1*1e3)/1e3
+Z2=%i * round(X2*1e3)/1e3
+Z0=%i * round(X0*1e3)/1e3
+
+//end q11
+
+a=exp(%i * 2 * %pi /3)
+A= [1 1 1 ; 1 a^2 a; 1 a a^2]
+
+
+Ia1=1/(Z1 + (Z2*Z0/(Z2+Z0)))
+Ia1=round(Ia1 *1e3)/1e3
+Ia2=(Z0/(Z2+Z0)) * Ia1 *-1
+Ia2=round(Ia2 *1e3)/1e3
+Ia0=(Z2/(Z2+Z0)) * Ia1 *-1
+Ia0=round(Ia0 *1e3)/1e3
+
+mIa1=[ Ia0 Ia1 Ia2]'
+mIa=A*mIa1
+Ib=round(S*1e7/(sqrt(3) * V2 * 1e3))/10
+Ia=round(abs(mIa)) * Ib
+
+
+mprintf("\nline currents are: in Amperes ")
+mprintf("\nIa= %s", string(round(abs(Ia(1))*10)/10) )
+mprintf("\nIb= %s", string(round(abs(Ia(2))*10)/10) +'/_'+ string(round(atand(imag(mIa(2)*Ib)/real(mIa(2)*Ib))*100)/100 +180 ))
+mprintf("\nIc= %s", string(round(abs(Ia(3))*10)/10) +'/_'+ string(round(atand(imag(mIa(3)*Ib)/real(mIa(3)*Ib))*100)/100 ))
|