diff options
Diffstat (limited to '3640/CH6/EX6.2/Ex6_2.sce')
-rw-r--r-- | 3640/CH6/EX6.2/Ex6_2.sce | 64 |
1 files changed, 64 insertions, 0 deletions
diff --git a/3640/CH6/EX6.2/Ex6_2.sce b/3640/CH6/EX6.2/Ex6_2.sce new file mode 100644 index 000000000..0ce9922d6 --- /dev/null +++ b/3640/CH6/EX6.2/Ex6_2.sce @@ -0,0 +1,64 @@ +clc
+//Ex6_2 is an extension of Ex6_1
+//code uses userdefined function complexstring
+function s=complexstring(a)
+
+
+ if imag(a)>=0 then
+ s=sprintf('%g+%gi',real(a),imag(a))
+ else
+ s=sprintf('%g%gi',real(a),imag(a))
+
+ end
+ funcprot(0)
+endfunction
+r1a=12//resistance in ohms
+x1a=6.5//reactance in ohms
+Xc=-20 //reactance in ohms
+r1m=1.9 //from E6_1
+x2=2.6 //from Ex6_1
+s=1
+a=1.6 //no unit
+r2=3.6 //resistance in ohms
+x2=2.6 //reactance in ohms
+Xm=56 //magnetizing reactance in ohms
+Vm=115 //applied voltage in volts
+Zf=((%i*Xm)*((r2/s)+(%i*x2)))/((%i*Xm)+(r2/s)+(%i*x2))//from Ex6_1
+Zst=Zf
+Zb=Zf
+z1a=r1a+(%i*x1a)+(%i*Xc)
+disp('z1a='+complexstring(z1a)+'Ω')//ans may vary due to roundoff error
+mprintf("z1a=%fohm havinga phase angle of %f degrees\n",abs(z1a),(acos(real(z1a)/abs(z1a)))*180/%pi)
+Z12=((1/2)*(z1a/(a*a)))-(r1m+(%i*x2))//ans in textbook is wrong
+disp('Z12='+complexstring(Z12)+'Ω')//ans may vary due to roundoff error
+mprintf("Z12=%fohm havinga phase angle of %f degrees\n",abs(Z12),(acos(real(Z12)/abs(Z12)))*180/%pi)
+Vmf=(Vm/2)*(1-(%i/a))
+disp('Vmf='+complexstring(Vmf)+'V')//ans may vary due to roundoff error
+mprintf("Vmf=%fV havinga phase angle of %f degrees\n",abs(Vmf),(-1)*(acos(real(Vmf)/abs(Vmf)))*180/%pi)
+Vmb=(Vm/2)*(1+(%i/a))
+disp('Vmb='+complexstring(Vmb)+'V')//ans may vary due to roundoff error
+mprintf("Vmb=%fV having a phase angle of %f degrees\n",abs(Vmb),(acos(real(Vmb)/abs(Vmb)))*180/%pi)
+Imf=11.77*exp(%i*(-1)*54.93*%pi/180)//textbook doesnt provide any formula or hint for this calculation
+Imb=4.37*exp(%i*(-1)*19.7*%pi/180)//textbook doesnt provide any formula or hint for this calculation
+disp('Imf='+complexstring(Imf)+'A')//ans may vary due to roundoff error
+disp('Imb='+complexstring(Imb)+'A')//ans may vary due to roundoff error
+mprintf("Imf=%fA having a phase angle of %f degrees\n",11.77,-54.93)
+mprintf("Imb=%fA having a phase angle of %f degrees\n",4.37,-19.37)
+mprintf("τst=%fN-m\n",(2*real(Zst)*((abs(Imf)^2)-(abs(Imb)^2)))/(60*%pi))//ans may vary due to roundoff error
+Im=Imf+Imb
+disp('Im='+complexstring(Im)+'A')//ans may vary due to roundoff error
+mprintf("Im=%fA having a phase angle of%f degrees\n",abs(Im),(-1)*(acos(real(Im)/abs(Im)))*180/%pi)
+Ia=(%i*(Imf-Imb))/a
+disp('Ia='+complexstring(Ia)+'A')//ans may vary due to roundoff error
+mprintf("Ia=%fA having a phase angle of %f degrees\n",abs(Ia),(acos(real(Ia)/abs(Ia)))*180/%pi)
+I=Im+Ia
+disp('Line current='+complexstring(I)+'A')//ans may vary due to roundoff error
+mprintf("I=%fA having a phase angle of %f degrees\n",abs(I),(-1)*(acos(real(I)/abs(I)))*180/%pi)
+
+
+
+
+
+
+
+
|