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 /2252/CH6/EX6.8/Ex6_8.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 '2252/CH6/EX6.8/Ex6_8.sce')
| -rwxr-xr-x | 2252/CH6/EX6.8/Ex6_8.sce | 68 |
1 files changed, 68 insertions, 0 deletions
diff --git a/2252/CH6/EX6.8/Ex6_8.sce b/2252/CH6/EX6.8/Ex6_8.sce new file mode 100755 index 000000000..8a8cd3ddd --- /dev/null +++ b/2252/CH6/EX6.8/Ex6_8.sce @@ -0,0 +1,68 @@ +
+theta=linspace(0,2*%pi,100)
+i=50*sin(theta)
+xset('window',0)
+plot2d(theta,i)
+xtitle("Current wave shape for example 6.8-->(a)","theta","current")
+
+xset('window',1)
+theta=linspace(0,%pi,100)
+i=50*sin(theta)
+plot2d(theta,i)
+theta=linspace(%pi,2*%pi,100)
+i=-50*sin(theta)
+plot2d(theta,i)
+xtitle("Current wave shape for example 6.8-->(b)","theta","current")
+
+xset('window',2)
+theta=linspace(0,0,2)
+i=linspace(0,50,2)
+plot2d(theta,i)
+theta=linspace(0,%pi,2)
+i=linspace(50,50,2)
+plot2d(theta,i)
+theta=linspace(%pi,%pi,2)
+i=linspace(50,-50,2)
+plot2d(theta,i)
+theta=linspace(%pi,2*%pi,2)
+i=linspace(-50,-50,2)
+plot2d(theta,i)
+i=linspace(-50,0,2)
+theta=linspace(2*%pi,2*%pi,2)
+plot2d(theta,i)
+xtitle("Current wave shape for example 6.8-->(c)","theta","current")
+
+xset('window',3)
+theta=linspace(0,%pi/2,2)
+i=linspace(0,50,2)
+plot2d(theta,i)
+theta=linspace(%pi/2,%pi,2)
+i=linspace(50,0,2)
+plot2d(theta,i)
+theta=linspace(%pi,3*%pi/2,2)
+i=linspace(0,-50,2)
+plot2d(theta,i)
+theta=linspace(3*%pi/2,2*%pi,2)
+i=linspace(-50,0,2)
+plot2d(theta,i)
+xtitle("Current wave shape for example 6.8-->(d)","theta","current")
+
+//consider wave shape (a)
+mprintf("For wave shape (a)\nAs the negative and positive parts of the wave are equal, reading of moving coil ammeter is zero\n")
+Irms=sqrt(integrate('(50*sin(theta))^2','theta',0,2*%pi)/(2*%pi))
+mprintf("Reading of moving iron ammeter=%f A\n",Irms)
+
+//consider wave shape (b)
+Iavg=integrate('50*sin(theta)','theta',0,%pi)/%pi
+mprintf("For wave shape (b)\nReading on the moving coil ammeter=%f A\n",Iavg)
+Irms=sqrt(integrate('(50*sin(theta))^2','theta',0,2*%pi)/(2*%pi))
+mprintf("Reading on moving iron ammeter=%f A\n",Irms)
+
+//consider case (c)
+mprintf("For wave shape (c)\nAverage value over one complete period is clearly zero. Thus reading on moving coil ammeter is zero. As the value of current remains constant at 50 A during various intervals, RMS value will be 50 A only. Hence, reading on moving iron ammeter=50 A\n")
+
+//consider case(d)
+Iavg=(0+10+20+30+40+50++40+30+20+10+0+(-10)+(-20)+(-30)+(-40)+(-50)+(-40)+(-30)+(-20)+(-10)+0)/20
+mprintf("For wave shape (d)\nReading on moving coil ammeter=%d\n",Iavg)
+Irms=sqrt((0^2+10^2+20^2+30^2+40^2+50^2+40^2+30^2+20^2+10^2+0^2+(-10)^2+(-20)^2+(-30)^2+(-40)^2+(-50)^2+(-40)^2+(-30)^2+(-20)^2+(-10)^2+0^2)/20)
+mprintf("Reading on moving iron ammeter=%f A\n",Irms)
|