summaryrefslogtreecommitdiff
path: root/3872/CH2/EX2.1/EX2_1.sce
diff options
context:
space:
mode:
authorprashantsinalkar2017-10-10 12:27:19 +0530
committerprashantsinalkar2017-10-10 12:27:19 +0530
commit7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch)
treedbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3872/CH2/EX2.1/EX2_1.sce
parentb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff)
downloadScilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip
initial commit / add all books
Diffstat (limited to '3872/CH2/EX2.1/EX2_1.sce')
-rw-r--r--3872/CH2/EX2.1/EX2_1.sce47
1 files changed, 47 insertions, 0 deletions
diff --git a/3872/CH2/EX2.1/EX2_1.sce b/3872/CH2/EX2.1/EX2_1.sce
new file mode 100644
index 000000000..ba6cba3b4
--- /dev/null
+++ b/3872/CH2/EX2.1/EX2_1.sce
@@ -0,0 +1,47 @@
+//Book - Power System: Analysis & Design 5th Edition
+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye
+//Chapter - 2 ; Example 2.1
+//Scilab Version - 6.0.0 ; OS - Windows
+
+clc;
+clear;
+Vmax=141.4; //Peak source voltage in Volts
+R=10; //Load resistance in Ohms
+Xl=3.77; //Inductive reactance in Ohms
+Vrms=Vmax/sqrt(2); //RMS value of source voltage in Volts
+Ir=Vrms/(R); //Current through the resistor in Amperes
+Il=Vrms/(%i*Xl); //Current through the inductor in Amperes
+Iload=Ir+Il; //Current through the load in Amperes
+wt=0:0.1:2*%pi;
+v=Vmax*cos(wt); //Instantaneous voltage in Volts
+ir=Vmax*cos(wt)/R; //Instantaneous current through the resistor in Amperes
+il=Vmax*cos(wt+90*%pi/180); //Instantaneous current through the inductor in Amperes
+Pr=Vrms*Ir*(1+cos(2*wt)); //Instantaneous Power absorbed by Resistor in Watts
+Pl=Vrms*abs(Il)*sin(2*wt); //Instantaneous Power absorbed by Inductor in Watts
+del=0;
+bet=atan(imag(Iload),real(Iload));
+P=Vrms*abs(Iload)*cos(del-bet); //Real power absorbed by the load in Watts
+Q=Vrms*abs(Iload)*sin(del-bet); //Reactive power absorbed by the load in VAR
+pf=cos(del-bet); //Power factor
+clf; //To clear figures from previous programs
+subplot(231);
+plot(wt,v);
+xtitle('Input Voltage','Angular displacement','Voltage(Volts)');
+subplot(232);
+plot(wt,ir);
+xtitle('Current through resistor','Angular displacement','Current(Amp.)');
+subplot(233);
+plot(wt,Pr);
+xtitle('Power dissipated in resistor','Angular displacement','Power(Watts)');
+subplot(236);
+xtitle('Power throgh the inductor','Angular displacement','Power(VAR)');
+plot(wt,Pl);
+subplot(234);
+plot(wt,v);
+xtitle('Input Voltage','Angular displacement','Voltage(Volts)');
+subplot(235);
+plot(wt,il);
+xtitle('Current through inductor','Angular displacement','Current(Amp.)');
+printf('\nThe Real power absorbed by the load is %d Watts\n',P);
+printf('The Reactive power absorbed by the load is %d VAR\n',Q);
+printf('The Power factor is %.4f lagging',pf);