1 files changed, 54 insertions, 0 deletions
diff --git a/3640/CH4/EX4.4/Ex4_4.sce b/3640/CH4/EX4.4/Ex4_4.sce
new file mode 100644
index 000000000..4b28c9532
--- /dev/null
+++ b/3640/CH4/EX4.4/Ex4_4.sce
@@ -0,0 +1,54 @@
+clc
+//below is an extension of Ex4_3
+//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
+x1=0.35 //primary reactance in ohms
+r1=0.39 //primary resistance in ohms
+Xm=16 //magnetizing reactance
+r2=0.14 //secondary resistance in ohms 
+x2=0.35 //secondary reactance in ohms
+ws=188.5 //angular frequency in rad/sec
+V=220 //rated voltage in volts
+//part a
+V1m=V/sqrt(3)//ans may vary due to roundoff error
+VTH=V1m*(Xm/(Xm+x2))
+mprintf("VTH=V1m=%fV\n",VTH)//ans may vary due to roundoff error
+X1=x1
+mprintf("X1=%fΩ\n",X1)
+R1=r1*(Xm/(x1+Xm))//ans may vary due to roundoff error
+mprintf("R1=%fΩ\n",R1)
+mprintf("τmax=%fN-m\n",((3/ws)*(VTH^2))/(2*(R1+sqrt((R1^2)+((2*X1)^2)))))//ans may vary due to roundoff error
+//part b
+sM=r2/sqrt((R1^2)+((X1+x1)^2))//ans may vary due to roundoff error
+mprintf("sM=%f\n",sM)
+mprintf("r2/sM=%fΩ\n",r2/sM)//ans may vary due to roundoff error
+Zf=((%i*Xm)*((r2/sM)+(%i*x2)))/((r2/sM)+(%i*(x2+Xm)))//ans may vary due to roundoff error
+disp('Zf='+complexstring(Zf)+'Ω')
+mprintf("Zf=%fohm having a phase angle of %f degrees\n",abs(Zf),(acos(real(Zf)/abs(Zf)))*180/%pi)
+z1=r1+(%i*x1)
+Zin=z1+Zf
+disp('Zin='+complexstring(Zin)+'Ω')//ans may vary due to roundoff error
+mprintf("Zin=%fohm having a phase angle of %f degrees\n",abs(Zin),(acos(real(Zin)/abs(Zin)))*180/%pi)
+I1=V1m/abs(Zin)
+mprintf("I1=%fA\n",I1)//ans may vary due to roundoff error
+Rf=real(Zf) //resistance in ohms
+Pg=3*I1*I1*Rf//ans  in the textbook is wrong
+mprintf("Pg=%fW\n",Pg)
+mprintf("τmax=Pg/ws=%fN-m\n",Pg/ws)//ans may vary due to roundoff error
+
+
+
+
+
+
+