summaryrefslogtreecommitdiff
path: root/Working_Examples/215/CH6/EX6.19/ex6_19.sce
blob: d7e6f08ba95b6c4e12bfea3d530914c12b8f5af7 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
clc
//Example 6.19
//Calculate total current through load
//On applying source transformation 
//From figure 6.65
i=%i
V1=10;V2mag=5;V2ph=90;V3mag=14.4;V3ph=225;
x=V2mag * cos (( V2ph * %pi ) /180) ;
y=V2mag * sin (( V2ph * %pi ) /180) ;
V2= complex (x,y)
a=V3mag * cos (( V3ph * %pi ) /180) ;
b=V3mag * sin (( V3ph * %pi ) /180) ;
V3= complex (a,b)
G1=1/2;G2=1/(2+i*3);G3=1/(2-i*2);
//By applying Millman Theorem
disp('V=((V1*G1)+(V2*G2)+(V3*G3))/(G1+G2+G3)')
V=((V1*G1)+(V2*G2)+(V3*G3))/(G1+G2+G3)
[Vmag Vang]=polar(V)
R=1/(G1+G2+G3)
printf("V=%3.2f(%3.2f deg)V",Vmag,(Vang*180)/%pi)
disp(R,'R=')
//Consider the resultant circuit from figure 6.66
disp('Let the total current through 3+i4 be I')
//Applying KVL to the circuit
I=V/(3+i*4+R)
[Imag Iang]=polar(I)
printf("I=%3.2f(%3.2f deg)V",Imag,(Iang*180)/%pi)