diff options
Diffstat (limited to 'Electrical_Network_by_R._Singh/Chapter4.ipynb')
| -rw-r--r-- | Electrical_Network_by_R._Singh/Chapter4.ipynb | 1846 |
1 files changed, 1846 insertions, 0 deletions
diff --git a/Electrical_Network_by_R._Singh/Chapter4.ipynb b/Electrical_Network_by_R._Singh/Chapter4.ipynb new file mode 100644 index 00000000..da063bd5 --- /dev/null +++ b/Electrical_Network_by_R._Singh/Chapter4.ipynb @@ -0,0 +1,1846 @@ +{
+ "metadata": {
+ "name": "",
+ "signature": "sha256:3e2d2b728039737326c3b0b27fd3059d5efefa657f8300b7cf26aeeb6c7e63cf"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "Chapter4-AC Circuits"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex1-pg4.4"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits:example 4.1:(pg4.4)\n",
+ "import math\n",
+ "i=15.;\n",
+ "\n",
+ "t=3.375*10**-3;\n",
+ "f=40.;\n",
+ "pi=3.14;\n",
+ "Im=(i/math.sin(2.*pi*f*t));\n",
+ "print(\"i=15 Amp\");\n",
+ "print(\"t=3.375 ms\");\n",
+ "print(\"f=40 Hz\");\n",
+ "print(\"i=Im*sin(2*pi*f*t)\");\n",
+ "print'%s %.2f %s'%(\"Im= \",Im, \"Amp\");"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "i=15 Amp\n",
+ "t=3.375 ms\n",
+ "f=40 Hz\n",
+ "i=Im*sin(2*pi*f*t)\n",
+ "Im= 20.00 Amp\n"
+ ]
+ }
+ ],
+ "prompt_number": 1
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex2-pg4.4"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits:example 4.2:(pg4.4)\n",
+ "import math\n",
+ "f=50.;\n",
+ "Im=100.;\n",
+ "i1=86.6;\n",
+ "t=(1/600.);\n",
+ "pi=3.14;\n",
+ "print(\"f=50 c/s\");\n",
+ "print(\"Im=100 A\");\n",
+ "## part(a)\n",
+ "print(\"i=Im*sin(2*pi*f*t)\");\n",
+ "i=Im*math.sin(2*pi*f*t);\n",
+ "print'%s %.2f %s'%(\"i= \",i,\" A\");\n",
+ "## part (b)\n",
+ "print(\"i=Im*sin(2*pi*f*t1)\");\n",
+ "t1=(math.asin(i1/Im)/(2.*pi*f));\n",
+ "print'%s %.2e %s'%(\"t1= \",t1,\" second\");"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "f=50 c/s\n",
+ "Im=100 A\n",
+ "i=Im*sin(2*pi*f*t)\n",
+ "i= 49.98 A\n",
+ "i=Im*sin(2*pi*f*t1)\n",
+ "t1= 3.33e-03 second\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex3-pg4.5"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits:example 4.3:(pg4.5)\n",
+ "f=50.;\n",
+ "import math\n",
+ "I=20.;\n",
+ "t1=0.0025;\n",
+ "t2=0.0125;\n",
+ "I1=14.14;\n",
+ "pi=3.14;\n",
+ "print(\"f=50 c/s\");\n",
+ "print(\"I=20 A\");\n",
+ "print(\"Im=I*sqrt(2)\");\n",
+ "Im=(math.sqrt(2)*I);\n",
+ "print'%s %.2f %s'%(\"\\nIm= \",Im,\" A\");\n",
+ "print(\"\\nEquation of current, \\ni=Im*sin(2*pi*f*t)\");\n",
+ "print(\"=28.28sin(2*pi*f*t)=28.28sin(100*pi*t)\");\n",
+ "print(\"(a)At t=0.0025 seconds\");\n",
+ "i=(Im*math.sin(2.*pi*f*t1));\n",
+ "print'%s %.2f %s'%(\"i= \",i,\" A\"); ##when t=0.0025seconds\n",
+ "print(\"(b)At t=0.0125 seconds\");\n",
+ "i=(Im*math.sin(2*pi*f*t2));\n",
+ "print'%s %.2f %s'%(\"i= \",i,\" A\"); ##when t=0.0125seconds\n",
+ "print(\"(c) i=28.28sin(100*pi*t) \");\n",
+ "t=(math.asin(I1/Im)/(2*math.pi*f));\n",
+ "print'%s %.2e %s'%(\"t= \",t,\" second\");## when I=14.14A"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": []
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex4-pg4.5"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.4 :pg(4.5)\n",
+ "import math\n",
+ "pi=3.14;\n",
+ "Vm=200.;\n",
+ "print(\"v=200sin314t\");\n",
+ "print(\"v=Vmsin(2*pi*f*t)\");\n",
+ "print(\"(2*pi*f)=314\");\n",
+ "f=(314./(2.*pi));\n",
+ "print'%s %.2f %s'%(\"f= \",f,\" Hz\");\n",
+ "Vavg=((2.*Vm)/pi);\n",
+ "Vrms=(Vm/math.sqrt(2.));\n",
+ "print('\\nFor a sinusoidal waveform, \\nVavg=(2*Vm/pi) \\nVrms=(Vm/sqrt(2))');\n",
+ "kf=(Vrms/Vavg);\n",
+ "kc=(Vm/Vrms);\n",
+ "print'%s %.2f %s'%('\\nform fator=',kf,'');\n",
+ "print'%s %.2f %s'%('\\ncrest factor=',kc,'');\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "v=200sin314t\n",
+ "v=Vmsin(2*pi*f*t)\n",
+ "(2*pi*f)=314\n",
+ "f= 50.00 Hz\n",
+ "\n",
+ "For a sinusoidal waveform, \n",
+ "Vavg=(2*Vm/pi) \n",
+ "Vrms=(Vm/sqrt(2))\n",
+ "\n",
+ "form fator= 1.11 \n",
+ "\n",
+ "crest factor= 1.41 \n"
+ ]
+ }
+ ],
+ "prompt_number": 3
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex5-pg4.6"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.5 :(pg 4.6)\n",
+ "kf=1.2;\n",
+ "import math\n",
+ "kp=1.5;\n",
+ "Vavg=10.;\n",
+ "print(\"kf=1.2\");\n",
+ "print(\"kp=1.5\");\n",
+ "print(\"Vavg=10\");\n",
+ "print(\"form factor kf=(Vrms/Vavg)\");\n",
+ "Vrms=(kf*Vavg);\n",
+ "print'%s %.2f %s'%(\"\\nVrms= \",Vrms,\" V\");\n",
+ "print(\"peak factor kp=(Vm/Vrms)\");\n",
+ "Vm=(kp*Vrms);\n",
+ "print'%s %.2f %s'%(\"\\nVm= \",Vm,\" V\");"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "kf=1.2\n",
+ "kp=1.5\n",
+ "Vavg=10\n",
+ "form factor kf=(Vrms/Vavg)\n",
+ "\n",
+ "Vrms= 12.00 V\n",
+ "peak factor kp=(Vm/Vrms)\n",
+ "\n",
+ "Vm= 18.00 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 4
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex14-pg4.11"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits: example 4.14 :(pg 4.11)\n",
+ "v1=0.;\n",
+ "import math\n",
+ "v2=40.;\n",
+ "v3=60.;\n",
+ "v4=80.;\n",
+ "v5=100.;\n",
+ "t=8.;\n",
+ "Vavg=((v1+v2+v3+v4+v5+v4+v3+v2)/t);\n",
+ "Vrms=math.sqrt((v1**2+v2**2+v3**2+v4**2+v5**2+v4**2+v3**2+v2**2)/t);\n",
+ "print(\"Vavg=((0.+40.+60.+80.+100.+80.+60.+40.)/8.)\");\n",
+ "print'%s %.2f %s'%(\"\\nVavg= \",Vavg,\" V\");\n",
+ "print(\"Vrms=sqrt((0+(40)^2+(60)^2+(80)^2+(100)^2+(80)^2+(60)^2+(40)^2)/8)\");\n",
+ "print'%s %.2f %s'%(\"\\nVrms= \",Vrms,\" V\");"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Vavg=((0.+40.+60.+80.+100.+80.+60.+40.)/8.)\n",
+ "\n",
+ "Vavg= 57.50 V\n",
+ "Vrms=sqrt((0+(40)^2+(60)^2+(80)^2+(100)^2+(80)^2+(60)^2+(40)^2)/8)\n",
+ "\n",
+ "Vrms= 64.42 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 5
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex15-pg4.11"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.15 :pg(4.11 & 4.12)\n",
+ "v1=0.;\n",
+ "import math\n",
+ "v2=10.;\n",
+ "v3=20.;\n",
+ "t=3.;\n",
+ "Vavg=((v1+v2+v3)/t);\n",
+ "Vrms=(math.sqrt((v1**2+v2**2+v3**2)/t));\n",
+ "print(\"Vavg=((0+10+20)/3)\");\n",
+ "print'%s %.2f %s'%(\"Vavg= \",Vavg, \"V\");\n",
+ "print(\"Vrms=(((0)^2+(10)^2+(20)^2)/3)\");\n",
+ "print'%s %.2f %s'%(\"Vrms= \",Vrms,\" V\");"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Vavg=((0+10+20)/3)\n",
+ "Vavg= 10.00 V\n",
+ "Vrms=(((0)^2+(10)^2+(20)^2)/3)\n",
+ "Vrms= 12.91 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 6
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex33-pg4.27"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.33 :pg(4.27)\n",
+ "Vm=177.;\n",
+ "import math\n",
+ "Im=14.14;\n",
+ "phi=30.;\n",
+ "V=(Vm/math.sqrt(2));\n",
+ "I=(Im/math.sqrt(2));\n",
+ "pf=math.cos(30/57.3);\n",
+ "P=(V*I*pf);\n",
+ "print(\"v(t)=177sin(314t+10)\");## value of 10 is in degrees\n",
+ "print(\"i(t)=14.14sin(314t-20)\");##value of 20 is in degrees\n",
+ "print(\"\\nCurrent i(t) lags behind voltage v(t) by 30degrees\");\n",
+ "print(\"phi=30degrees\");\n",
+ "print'%s %.2f %s'%(\"Power factor pf=cos(30)= \",pf,\" (lagging)\");\n",
+ "print'%s %.2f %s'%(\"\\nPower consumed P=V*I*cos(phi)= \",P,\" W\");"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "v(t)=177sin(314t+10)\n",
+ "i(t)=14.14sin(314t-20)\n",
+ "\n",
+ "Current i(t) lags behind voltage v(t) by 30degrees\n",
+ "phi=30degrees\n",
+ "Power factor pf=cos(30)= 0.87 (lagging)\n",
+ "\n",
+ "Power consumed P=V*I*cos(phi)= 1083.76 W\n"
+ ]
+ }
+ ],
+ "prompt_number": 7
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex42-pg4.32"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.42 :pg(4.32 & 4.33)\n",
+ "import math\n",
+ "PR=1000.;\n",
+ "VR=200.;\n",
+ "Pcoil=250.;\n",
+ "Vcoil=300.;\n",
+ "R=((VR**2)/PR);\n",
+ "I=(VR/R);\n",
+ "r=((Pcoil/(I**2)));\n",
+ "Zcoil=(Vcoil/I);\n",
+ "XL=math.sqrt((Zcoil**2)-(r**2));\n",
+ "RT=(R+r);\n",
+ "ZT=math.sqrt((RT**2)+(XL**2));\n",
+ "V=(ZT*I);\n",
+ "print(\"\\nPR=1000 W \\nVR=200 V \\nPcoil=250 W \\nVcoil=300 V \\nPR=(VR^2/R)\");\n",
+ "print'%s %.2f %s'%(\"\\nR= \",R,\" Ohms\");\n",
+ "print'%s %.2f %s'%(\"\\nVR=R*I \\nI= \",I,\" A\");\n",
+ "print(\"Pcoil=(I^2)*r\");\n",
+ "print'%s %.2f %s'%(\"\\nResistance of coil r= \",r,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nImpedance of coil Zcoil=(Vcoil/I)= \",Zcoil,\" Ohms\");\n",
+ "print'%s %.2f %s'%(\"\\nReactance of coil XL=sqrt((Zcoil^2)-(r^2)) = \",XL,\" Ohms\");\n",
+ "print'%s %.2f %s'%(\"\\nCombined resistance RT=R+r= \",RT,\"Ohms\");\n",
+ "print'%s %.2f %s'%(\"\\nCombined impedance ZT=sqrt(((R+r)^2)+(XL^2)) = \",ZT,\" Ohms\");\n",
+ "print'%s %.2f %s'%(\"\\nSupply voltage V=ZT*I= \",V,\" V\");"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "PR=1000 W \n",
+ "VR=200 V \n",
+ "Pcoil=250 W \n",
+ "Vcoil=300 V \n",
+ "PR=(VR^2/R)\n",
+ "\n",
+ "R= 40.00 Ohms\n",
+ "\n",
+ "VR=R*I \n",
+ "I= 5.00 A\n",
+ "Pcoil=(I^2)*r\n",
+ "\n",
+ "Resistance of coil r= 10.00 Ohm\n",
+ "\n",
+ "Impedance of coil Zcoil=(Vcoil/I)= 60.00 Ohms\n",
+ "\n",
+ "Reactance of coil XL=sqrt((Zcoil^2)-(r^2)) = 59.16 Ohms\n",
+ "\n",
+ "Combined resistance RT=R+r= 50.00 Ohms\n",
+ "\n",
+ "Combined impedance ZT=sqrt(((R+r)^2)+(XL^2)) = 77.46 Ohms\n",
+ "\n",
+ "Supply voltage V=ZT*I= 387.30 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 8
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex47-pg4.36"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.47 :pg(4.47)\n",
+ "import math\n",
+ "f1=60.;\n",
+ "V=200.;\n",
+ "P=600.;\n",
+ "I=5.;\n",
+ "f=50.;\n",
+ "Z=V/I;\n",
+ "r=(P/(I**2));\n",
+ "XL=math.sqrt((Z**2)-(r**2));\n",
+ "L=(XL/(2.*math.pi*f));\n",
+ "XL1=(2.*math.pi*f1*L);\n",
+ "Z1=math.sqrt((r**2)+(XL1**2));\n",
+ "I=(V/Z1);\n",
+ "print(\"\\nI=5 A \\nV=200 V \\nP=600 W \\nFor f=50 Hz,\");\n",
+ "print'%s %.2f %s'%(\"\\nZ=V/I = \",Z,\" Ohms\");\n",
+ "print'%s %.2f %s'%(\"\\nP=((I^2)*r) \\nr= \",r,\" Ohms\");\n",
+ "print'%s %.2f %s'%(\"\\nXL=sqrt((Z^2)-(r^2)) \\nXL= \",XL,\" Ohms\");\n",
+ "print'%s %.2f %s'%(\"\\nXL=(2*pi*f*L)\\nL= \",L,\" H\");\n",
+ "print'%s %.2f %s'%(\"\\nFor f=60 Hz \\nXL= \",XL1,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nr=24 Ohms \\nZ=sqrt((r^2)+(XL^2))= \",Z1,\" Ohms\");\n",
+ "print'%s %.2f %s'%(\"\\nI=V/Z= \",I,\" A\");"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "I=5 A \n",
+ "V=200 V \n",
+ "P=600 W \n",
+ "For f=50 Hz,\n",
+ "\n",
+ "Z=V/I = 40.00 Ohms\n",
+ "\n",
+ "P=((I^2)*r) \n",
+ "r= 24.00 Ohms\n",
+ "\n",
+ "XL=sqrt((Z^2)-(r^2)) \n",
+ "XL= 32.00 Ohms\n",
+ "\n",
+ "XL=(2*pi*f*L)\n",
+ "L= 0.10 H\n",
+ "\n",
+ "For f=60 Hz \n",
+ "XL= 38.40 Ohm\n",
+ "\n",
+ "r=24 Ohms \n",
+ "Z=sqrt((r^2)+(XL^2))= 45.28 Ohms\n",
+ "\n",
+ "I=V/Z= 4.42 A\n"
+ ]
+ }
+ ],
+ "prompt_number": 9
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex48-pg4.37"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.48 :(pg 4.37)\n",
+ "f=50.;\n",
+ "import math\n",
+ "pi=3.14;\n",
+ "Vdc=12.;\n",
+ "Idc=2.5;\n",
+ "Vac=230.;\n",
+ "Iac=2.;\n",
+ "Pac=50.;\n",
+ "R=(Vdc/Idc);\n",
+ "Z=(Vac/Iac);\n",
+ "Pi=(Pac-((Iac**2)*R));\n",
+ "RT=(Pac/(Iac**2));\n",
+ "XL=math.sqrt((Z**2)-(RT**2));\n",
+ "L=(XL/(2.*pi*f));\n",
+ "pf=(RT/Z);\n",
+ "i=(Pi/(Iac**2));\n",
+ "print(\"\\nFor dc V=12 V, I=2.5 A \\nFor ac V=230 V, I=2 A, P=50 W\");\n",
+ "print(\"\\nIn an iron-cored coil,there are two types of losses \\n(i)Losses in core known as core or iron loss \\n(ii)Losses in winding known as copper loss\");\n",
+ "print(\"\\nP=(I^2)*R+Pi \\nP/(I^2)=R+((Pi)/(I^2)) \\nRT=R+(Pi/(I^2)) \\nwhere R is the resistance of the coil and (Pi/I^2) is the resistance which is equivalent to the effect of iron loss\");\n",
+ "print(\"\\nFor dc supply, f=0 \\nXL=0\");\n",
+ "print'%s %.2f %s'%(\"\\nR= \",R,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nFor ac supply \\nZ= \",Z,\" Ohms\");\n",
+ "print'%s %.2f %s'%(\"\\nIron loss Pi=P-I^2*R= \",Pi,\" W\");\n",
+ "print'%s %.2f %s'%(\"\\nRT=(P/I^2)= \",RT,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nXL=sqrt((Z^2)-(RT^2))= \",XL,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nXL=2*pi*L \\nInductance L= \",L,\" H\");\n",
+ "print'%s %.2f %s'%(\"\\nPower factor =RT/Z= \",pf,\" (lagging)\");\n",
+ "print'%s %.2f %s'%(\"\\nThe series resistance equivalent to the effect of iron loss= Pi/(I^2)= \",i,\" Ohms\");"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "For dc V=12 V, I=2.5 A \n",
+ "For ac V=230 V, I=2 A, P=50 W\n",
+ "\n",
+ "In an iron-cored coil,there are two types of losses \n",
+ "(i)Losses in core known as core or iron loss \n",
+ "(ii)Losses in winding known as copper loss\n",
+ "\n",
+ "P=(I^2)*R+Pi \n",
+ "P/(I^2)=R+((Pi)/(I^2)) \n",
+ "RT=R+(Pi/(I^2)) \n",
+ "where R is the resistance of the coil and (Pi/I^2) is the resistance which is equivalent to the effect of iron loss\n",
+ "\n",
+ "For dc supply, f=0 \n",
+ "XL=0\n",
+ "\n",
+ "R= 4.80 Ohm\n",
+ "\n",
+ "For ac supply \n",
+ "Z= 115.00 Ohms\n",
+ "\n",
+ "Iron loss Pi=P-I^2*R= 30.80 W\n",
+ "\n",
+ "RT=(P/I^2)= 12.50 Ohm\n",
+ "\n",
+ "XL=sqrt((Z^2)-(RT^2))= 114.32 Ohm\n",
+ "\n",
+ "XL=2*pi*L \n",
+ "Inductance L= 0.36 H\n",
+ "\n",
+ "Power factor =RT/Z= 0.11 (lagging)\n",
+ "\n",
+ "The series resistance equivalent to the effect of iron loss= Pi/(I^2)= 7.70 Ohms\n"
+ ]
+ }
+ ],
+ "prompt_number": 10
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex49-pg4.37"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.49 :(pg 4.37 & 4.38)\n",
+ "import math\n",
+ "f=50.;\n",
+ "I1=4.;\n",
+ "pf1=0.5;\n",
+ "V1=200.;\n",
+ "I2=5.;\n",
+ "pf2=0.8;\n",
+ "V2=40.;\n",
+ "Z1=(V2/I2);\n",
+ "R=(Z1*pf2);\n",
+ "XL1=math.sqrt((Z1**2)-(R**2));\n",
+ "L1=(XL1/(2.*math.pi*f));\n",
+ "Z2=(V1/I1);\n",
+ "RT=(Z2*pf1);\n",
+ "XL2=math.sqrt((Z2**2)-(RT**2));\n",
+ "L2=(XL2/(2.*math.pi*f));\n",
+ "Pi=(V1*I1*pf1-(I1**2)*R);\n",
+ "print(\"\\nWith iron core I=4 A pf=0.5, V=200 V \\nWithout iron core I=5 A pf=0.8, V=40 V \\nWhen the iron-core is removed,\");\n",
+ "print'%s %.2f %s'%(\"\\nZ=V/I= \",Z1,\" Ohms\");\n",
+ "print'%s %.2f %s'%(\"\\npf=R/Z \\nR= \",R,\" Ohms\");\n",
+ "print'%s %.2f %s'%(\"\\nXL=sqrt((Z**2)-(RT**2))= \",XL1,\" Ohms\");\n",
+ "print'%s %.2f %s'%(\"\\nXL=(2*pi*f*L) \\nInductance L= \",L1,\" H\");\n",
+ "print'%s %.2f %s'%(\"\\nWith iron core, \\nZ= \",Z2,\" Ohms\");\n",
+ "print'%s %.2f %s'%(\"\\npf=RT/Z \\nRT= \",RT,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nXL=sqrt((Z**2)-(RT**2))= \",XL2,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nXL=(2*pi*f*L) \\nInductance L= \",L2,\" H\");\n",
+ "print'%s %.2f %s'%(\"\\nIron loss Pi=P=(I**2)*R \\n=VIcos(phi)-I**2*R \\n= \",Pi,\" W\");\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "With iron core I=4 A pf=0.5, V=200 V \n",
+ "Without iron core I=5 A pf=0.8, V=40 V \n",
+ "When the iron-core is removed,\n",
+ "\n",
+ "Z=V/I= 8.00 Ohms\n",
+ "\n",
+ "pf=R/Z \n",
+ "R= 6.40 Ohms\n",
+ "\n",
+ "XL=sqrt((Z**2)-(RT**2))= 4.80 Ohms\n",
+ "\n",
+ "XL=(2*pi*f*L) \n",
+ "Inductance L= 0.02 H\n",
+ "\n",
+ "With iron core, \n",
+ "Z= 50.00 Ohms\n",
+ "\n",
+ "pf=RT/Z \n",
+ "RT= 25.00 Ohm\n",
+ "\n",
+ "XL=sqrt((Z**2)-(RT**2))= 43.30 Ohm\n",
+ "\n",
+ "XL=(2*pi*f*L) \n",
+ "Inductance L= 0.14 H\n",
+ "\n",
+ "Iron loss Pi=P=(I**2)*R \n",
+ "=VIcos(phi)-I**2*R \n",
+ "= 297.60 W\n"
+ ]
+ }
+ ],
+ "prompt_number": 11
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex51-pg4.40"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.51 :(pg 4.40 & 4.41)\n",
+ "import math\n",
+ "P=2000.;\n",
+ "pf=0.5;\n",
+ "V=230.;\n",
+ "S=(P/pf);\n",
+ "phi=math.acos(pf)*57.3;\n",
+ "I=(P/(V*pf));\n",
+ "Q=(V*I*math.sin(phi/57.3));\n",
+ "print(\"P=2000 W\");\n",
+ "print(\"pf=0.5 (leading)\");\n",
+ "print(\"V=230 V\");\n",
+ "print(\"P=V*I*cos(phi)\");\n",
+ "print'%s %.2f %s'%(\"\\nI= \",I,\" A\");\n",
+ "print'%s %.2f %s'%(\"\\nS=V*I=P/cos(phi)= \",S,\" VA\");\n",
+ "print'%s %.2f %s'%(\"\\nphi= \",phi,\" degrees\");\n",
+ "print'%s %.2f %s'%(\"\\nQ=V*I*sin(phi)= \",Q,\" VAR\");"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "P=2000 W\n",
+ "pf=0.5 (leading)\n",
+ "V=230 V\n",
+ "P=V*I*cos(phi)\n",
+ "\n",
+ "I= 17.39 A\n",
+ "\n",
+ "S=V*I=P/cos(phi)= 4000.00 VA\n",
+ "\n",
+ "phi= 60.00 degrees\n",
+ "\n",
+ "Q=V*I*sin(phi)= 3464.10 VAR\n"
+ ]
+ }
+ ],
+ "prompt_number": 12
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex52-pg4.41"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.52 :(pg 4.41)\n",
+ "import math\n",
+ "V=240.;\n",
+ "VR=100.;\n",
+ "P=300.;\n",
+ "f=50.;\n",
+ "R=((VR**2)/P);\n",
+ "I=math.sqrt(P/R);\n",
+ "Z=V/I;\n",
+ "XC=math.sqrt((Z**2)-(R**2));\n",
+ "C=(1./(2.*math.pi*f*XC));\n",
+ "VC=math.sqrt((V**2)-(VR**2));\n",
+ "VCmax=(VC*math.sqrt(2.));\n",
+ "Qmax=(C*VCmax);\n",
+ "Emax=((1./2.)*C*(VCmax**2));\n",
+ "print(\"\\nV=240 V \\nVR=100 V \\nP=300 W \\nf=50 Hz\");\n",
+ "print'%s %.2f %s'%(\"\\nP=(VR^2)/R \\nR=((VR^2)/P)= \",R,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nP=(I^2)*R \\nI=sqrt((P/R)) \\nI= \",I,\" A\");\n",
+ "print'%s %.2f %s'%(\"\\nZ=V/I=\",Z,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nXC=sqrt((Z^2)-(R^2))= \",XC,\" Ohm\");\n",
+ "print'%s %.2e %s'%(\"\\nXC=1/2*pi*f*C \\nC= \",C,\" F\");\n",
+ "print'%s %.2f %s'%(\"\\nVoltage across capacitor VC=sqrt((V^2)-(VR^2))= \",VC,\" V\");\n",
+ "print'%s %.2f %s %.2f %s '%(\"\\nMaximum value of max charge \\nVC= \",VCmax,\" V\" and \" \\nQmax=C*VCmax= \",Qmax,\" C\");\n",
+ "print'%s %.2f %s'%(\"\\nMax stored energy Emax=((1/2)*C*(VCmax^2)) \\n= \",Emax,\" J\");\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "V=240 V \n",
+ "VR=100 V \n",
+ "P=300 W \n",
+ "f=50 Hz\n",
+ "\n",
+ "P=(VR^2)/R \n",
+ "R=((VR^2)/P)= 33.33 Ohm\n",
+ "\n",
+ "P=(I^2)*R \n",
+ "I=sqrt((P/R)) \n",
+ "I= 3.00 A\n",
+ "\n",
+ "Z=V/I= 80.00 Ohm\n",
+ "\n",
+ "XC=sqrt((Z^2)-(R^2))= 72.72 Ohm\n",
+ "\n",
+ "XC=1/2*pi*f*C \n",
+ "C= 4.38e-05 F\n",
+ "\n",
+ "Voltage across capacitor VC=sqrt((V^2)-(VR^2))= 218.17 V\n",
+ "\n",
+ "Maximum value of max charge \n",
+ "VC= 308.54 \n",
+ "Qmax=C*VCmax= 0.01 C \n",
+ "\n",
+ "Max stored energy Emax=((1/2)*C*(VCmax^2)) \n",
+ "= 2.08 J\n"
+ ]
+ }
+ ],
+ "prompt_number": 13
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex53-pg4.42"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.53 :(pg 4.42)\n",
+ "import math\n",
+ "C=35.*10**-6;\n",
+ "f=50.;\n",
+ "XC=(1./(2.*math.pi*f*C));\n",
+ "R=math.sqrt(3.*(XC**2));\n",
+ "R2=(3*(XC**2));\n",
+ "print'%s %.2f %s'%(\"\\nC=35*10^-6 F \\nf=50 Hz \\nVC=1/2.V \\nXC=1/(2*pi*f*C)= \",XC,\" Ohm\");\n",
+ "print(\"\\nVC=1/2.V \\nXC.I=1/2.Z.I \\nXC=1/2.Z \\nZ=2.XC \\nZ=sqrt((R^2)+(XC^2)) \\n(2XC)^2=(R^2)+(XC^2) \\n3XC^2=R^2\");\n",
+ "print'%s %.2f %s %.2f %s '%(\"\\nR^2=3*XC^2= \",R2,\" Ohm\" and \" \\nR= \",R,\" Ohm\");\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "C=35*10^-6 F \n",
+ "f=50 Hz \n",
+ "VC=1/2.V \n",
+ "XC=1/(2*pi*f*C)= 90.95 Ohm\n",
+ "\n",
+ "VC=1/2.V \n",
+ "XC.I=1/2.Z.I \n",
+ "XC=1/2.Z \n",
+ "Z=2.XC \n",
+ "Z=sqrt((R^2)+(XC^2)) \n",
+ "(2XC)^2=(R^2)+(XC^2) \n",
+ "3XC^2=R^2\n",
+ "\n",
+ "R^2=3*XC^2= 24813.35 \n",
+ "R= 157.52 Ohm \n"
+ ]
+ }
+ ],
+ "prompt_number": 14
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex54-pg4.42"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.54 :(pg 4.42)\n",
+ "V=125.;\n",
+ "import math\n",
+ "I=2.2;\n",
+ "P=96.8;\n",
+ "f=50.;\n",
+ "Z=V/I;\n",
+ "R=(P/(I**2));\n",
+ "Xc=math.sqrt((Z**2)-(R**2));\n",
+ "C=(1./(2.*math.pi*f*Xc));\n",
+ "print(\"\\nV=125 V \\nP=96.8 W \\nI=2.2 A \\nf=50 Hz\");\n",
+ "print'%s %.2f %s'%(\"\\nZ=V/I= \",Z,\" A\");\n",
+ "print'%s %.2f %s'%(\"\\nP=(I^2)*R \\nR= \",R,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nXc=sqrt((Z^2)-(R^2))= \",Xc,\" Ohm\");\n",
+ "print'%s %.2e %s'%(\"\\nXc=1/(2*pi*f*C) \\n C= \",C,\" F\");"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "V=125 V \n",
+ "P=96.8 W \n",
+ "I=2.2 A \n",
+ "f=50 Hz\n",
+ "\n",
+ "Z=V/I= 56.82 A\n",
+ "\n",
+ "P=(I^2)*R \n",
+ "R= 20.00 Ohm\n",
+ "\n",
+ "Xc=sqrt((Z^2)-(R^2))= 53.18 Ohm\n",
+ "\n",
+ "Xc=1/(2*pi*f*C) \n",
+ " C= 5.99e-05 F\n"
+ ]
+ }
+ ],
+ "prompt_number": 15
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex57-pg4.46"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits :example 4.57 :(pg 4.46)\n",
+ "import math\n",
+ "j=-math.sqrt(1);\n",
+ "f=50.;\n",
+ "L=0.22;\n",
+ "R1=3.;\n",
+ "Z=3.8+j*6.4;\n",
+ "XL=2.*math.pi*f*L;\n",
+ "R2=3.8;\n",
+ "R=R2-R1;\n",
+ "X=6.4;\n",
+ "XC=XL-X;\n",
+ "C=(1./(2.*math.pi*f*XC));\n",
+ "print(\"\\nZ=(3.8+j*6.4) Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nXL=2*pi*f*L= \",XL,\" Ohm\");\n",
+ "print(\"\\nZ=(3+j69.12+R-jXC) \\n=(3+R)+j(69.12-XC)\");\n",
+ "print'%s %.2f %s'%(\"\\n3+R=3.8 \\nR= \",R,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nXC= \",XC,\"Ohm\");\n",
+ "print'%s %.2e %s'%(\"\\nXC=1/2.pi.f.C \\nC= \",C,\" F\");"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "Z=(3.8+j*6.4) Ohm\n",
+ "\n",
+ "XL=2*pi*f*L= 69.12 Ohm\n",
+ "\n",
+ "Z=(3+j69.12+R-jXC) \n",
+ "=(3+R)+j(69.12-XC)\n",
+ "\n",
+ "3+R=3.8 \n",
+ "R= 0.80 Ohm\n",
+ "\n",
+ "XC= 62.72 Ohm\n",
+ "\n",
+ "XC=1/2.pi.f.C \n",
+ "C= 5.08e-05 F\n"
+ ]
+ }
+ ],
+ "prompt_number": 16
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex58-pg4.46"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.58 :(pg 4.46)\n",
+ "import math\n",
+ "R=20.;\n",
+ "phi=45.;\n",
+ "Z=R/math.cos(phi/57.3);\n",
+ "XC=math.sqrt((Z**2)-(R**2));\n",
+ "XL=(2.*XC);\n",
+ "w=1000.;\n",
+ "L=(XL/w);\n",
+ "C=(1./(w*XC));\n",
+ "print(\"\\nvL=300sin(1000t) \\nR=20 Ohm \\nphi=45 \\nVL(max)=2Vcc(max) \\nsqrt(2)*VL=2*sqrt(2)*VC \\nI*XL=2*I*XC \\nXL=2*XC \\ncos(phi)=R/Z\");\n",
+ "print'%s %.2f %s'%(\"\\nZ= \",Z,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nZ=sqrt((R^2)+(XL-XC)^2) \\nXC = \",XC,\" Ohm\"); ##for series R-L-C ckt\n",
+ "print'%s %.2f %s'%(\"\\nXL=2*XC = \",XL,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nXL=w*L \\nL= \",L,\" H\");\n",
+ "print'%s %.2e %s'%(\"\\nXC=1/w*C \\nC= \",C,\" F\");\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "vL=300sin(1000t) \n",
+ "R=20 Ohm \n",
+ "phi=45 \n",
+ "VL(max)=2Vcc(max) \n",
+ "sqrt(2)*VL=2*sqrt(2)*VC \n",
+ "I*XL=2*I*XC \n",
+ "XL=2*XC \n",
+ "cos(phi)=R/Z\n",
+ "\n",
+ "Z= 28.28 Ohm\n",
+ "\n",
+ "Z=sqrt((R^2)+(XL-XC)^2) \n",
+ "XC = 20.00 Ohm\n",
+ "\n",
+ "XL=2*XC = 40.00 Ohm\n",
+ "\n",
+ "XL=w*L \n",
+ "L= 0.04 H\n",
+ "\n",
+ "XC=1/w*C \n",
+ "C= 5.00e-05 F\n"
+ ]
+ }
+ ],
+ "prompt_number": 17
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex59-pg4.47"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.59 :(pg 4.47)\n",
+ "import math\n",
+ "pf=0.5;\n",
+ "C=79.59*10**-6;\n",
+ "f=50.;\n",
+ "XC=(1./(2.*math.pi*f*C));\n",
+ "R=pf*XC;\n",
+ "Zcoil=XC;\n",
+ "XL=math.sqrt((Zcoil**2)-(R**2));\n",
+ "L=(XL/(2.*math.pi*f));\n",
+ "print(\"\\npf=0.5 \\nC=79.57uF \\nf=50 Hz \\nVcoil=VC \");\n",
+ "print'%s %.2f %s'%(\"\\nXC=1/2*pi*f*C = \",XC,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nVcoil=VC \\nZcoil=XC= \",XC,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\npf of coil=cos(phi)=R/Zcoil \\nResistance of coil R= \",R,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nXL=sqrt((Zcoil^2)-(R^2))= \",XL,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nXL=2*pi*f*L \\nInductance of coil= \",L,\" H\");"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "pf=0.5 \n",
+ "C=79.57uF \n",
+ "f=50 Hz \n",
+ "Vcoil=VC \n",
+ "\n",
+ "XC=1/2*pi*f*C = 39.99 Ohm\n",
+ "\n",
+ "Vcoil=VC \n",
+ "Zcoil=XC= 39.99 Ohm\n",
+ "\n",
+ "pf of coil=cos(phi)=R/Zcoil \n",
+ "Resistance of coil R= 20.00 Ohm\n",
+ "\n",
+ "XL=sqrt((Zcoil^2)-(R^2))= 34.64 Ohm\n",
+ "\n",
+ "XL=2*pi*f*L \n",
+ "Inductance of coil= 0.11 H\n"
+ ]
+ }
+ ],
+ "prompt_number": 18
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex60-pg4.48"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.60 :(pg 4.48)\n",
+ "import math\n",
+ "f=50.;\n",
+ "V=250.;\n",
+ "R=5.;\n",
+ "L=9.55;\n",
+ "Vcoil=300.;\n",
+ "XL=2.*math.pi*f*L;\n",
+ "Zcoil=(math.sqrt((R**2)+(XL**2)));\n",
+ "I=Vcoil/Zcoil;\n",
+ "Z=V/I;\n",
+ "XC1=Zcoil-Z;\n",
+ "XC2=Zcoil+Z;\n",
+ "C1=(1./(2.*math.pi*f*XC1));\n",
+ "C2=(1./(2.*math.pi*f*XC2));\n",
+ "print(\"\\nV=250 V \\nR=5 Ohm \\nL=9.55 H \\nVcoil=300 V\");\n",
+ "print'%s %.2f %s'%(\"\\nXL=2*pi*f*L = \",XL,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nZcoil=sqrt(R^2)+(XL^2) = \",Zcoil,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nI=Vcoil/Zcoil = \",I,\" A\");\n",
+ "print'%s %.2f %s'%(\"\\nZ=V/I = \",Z,\" Ohm\");##total impedance\n",
+ "print'%s %.2f %s'%(\"\\nZ=sqrt((R^2)+(XL-XC)^2) \\nXC= \",XC1,\" Ohm\");##when XL>XC\n",
+ "print'%s %.2f %s'%(\"\\nC=1/2*pi*f*XC = \",C1,\" F\");\n",
+ "print'%s %.2f %s'%(\"\\nZ=sqrt((R^2)+(XC-XL)^2) \\nXC= \",XC2,\" Ohm\");##when XC>XL\n",
+ "print'%s %.2e %s'%(\"\\nC= \",C2,\" F\");\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "V=250 V \n",
+ "R=5 Ohm \n",
+ "L=9.55 H \n",
+ "Vcoil=300 V\n",
+ "\n",
+ "XL=2*pi*f*L = 3000.22 Ohm\n",
+ "\n",
+ "Zcoil=sqrt(R^2)+(XL^2) = 3000.23 Ohm\n",
+ "\n",
+ "I=Vcoil/Zcoil = 0.10 A\n",
+ "\n",
+ "Z=V/I = 2500.19 Ohm\n",
+ "\n",
+ "Z=sqrt((R^2)+(XL-XC)^2) \n",
+ "XC= 500.04 Ohm\n",
+ "\n",
+ "C=1/2*pi*f*XC = 0.00 F\n",
+ "\n",
+ "Z=sqrt((R^2)+(XC-XL)^2) \n",
+ "XC= 5500.41 Ohm\n",
+ "\n",
+ "C= 5.79e-07 F\n"
+ ]
+ }
+ ],
+ "prompt_number": 19
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex79-pg4.64"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.79 :(pg 4.64)\n",
+ "import math\n",
+ "R=10.;\n",
+ "L=0.01;\n",
+ "C=100.*10**-6;\n",
+ "f0=(1./(2.*math.pi*math.sqrt(L*C)));\n",
+ "BW=(R/(2.*math.pi*L));\n",
+ "f1=f0-(BW/2.);\n",
+ "f2=f0+(BW/2.);\n",
+ "print(\"\\nR=10 Ohm \\nL=0.01H \\nC=100uF\");\n",
+ "print'%s %.2f %s'%(\"\\nf0=1/2*pi*sqrt(L*C)= \",f0,\" Hz\");##resonant frequency\n",
+ "print'%s %.2f %s'%(\"\\nBW=R/2*pi*L = \",BW,\" Hz\"); ##bandwidth\n",
+ "print'%s %.2f %s'%(\"\\nf1=f0-BW/2 \\n= \",f1,\" Hz\"); ##lower frequency\n",
+ "print'%s %.2f %s'%(\"\\nf2=f0+BW/2 = \",f2,\" Hz\"); ##higher frequency"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "R=10 Ohm \n",
+ "L=0.01H \n",
+ "C=100uF\n",
+ "\n",
+ "f0=1/2*pi*sqrt(L*C)= 159.15 Hz\n",
+ "\n",
+ "BW=R/2*pi*L = 159.15 Hz\n",
+ "\n",
+ "f1=f0-BW/2 \n",
+ "= 79.58 Hz\n",
+ "\n",
+ "f2=f0+BW/2 = 238.73 Hz\n"
+ ]
+ }
+ ],
+ "prompt_number": 28
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex80-pg4.65"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.80 :(pg 4.65)\n",
+ "import math\n",
+ "R=10.;\n",
+ "L=0.2;\n",
+ "C=40.*10**-6;\n",
+ "V=100.;\n",
+ "f0=(1./(2.*math.pi*math.sqrt(L*C)));\n",
+ "I0=(V/R);\n",
+ "P0=((I0**2)*R);\n",
+ "pf=1.;\n",
+ "Vr=(R*I0);\n",
+ "Vl=((2.*math.pi*f0*L)*I0);\n",
+ "Vc=((1/(2.*math.pi*f0*C))*I0);\n",
+ "Q=((1./R)*math.sqrt(L/C));\n",
+ "f1=(f0-(R/(4.*math.pi*L)));\n",
+ "f2=(f0+(R/(4.*math.pi*L)));\n",
+ "print(\"\\nR=10 Ohm \\nL=0.2 H \\nC=40uF \\nV=100 V\");\n",
+ "print'%s %.2f %s'%(\"\\n(i) f0= 1/2*pi*sqrt(LC) = \",f0,\" Hz\"); ##resonant frequency\n",
+ "print'%s %.2f %s'%(\"\\n(ii) I0= V/R = \",I0,\" A\"); ##current\n",
+ "print'%s %.2f %s'%(\"\\n(iii) P0=(I0^2)*R = \",P0,\" W\");##power\n",
+ "print(\"\\n(iv) pf=1\");##power factor\n",
+ "print'%s %.2f %s'%(\"\\n(v) Rv = R.I = \",Vr,\" V\");##voltage across resistor\n",
+ "print'%s %.2f %s'%(\"\\n Lv = XL.I = \",Vl,\" V\");##voltage across inductor\n",
+ "print'%s %.2f %s'%(\"\\n Cv = XC.I = \",Vc,\" V\"); ##voltage across capacitor\n",
+ "print'%s %.2f %s'%(\"\\n(vi) Q =1/R*sqrt(L/C)=\",Q,\"\");##Quality factor\n",
+ "print'%s %.2f %s'%(\"\\n(vii)f1 = f0-R/4.pi.L = \",f1,\" Hz\"); ##half power points\n",
+ "print'%s %.2f %s'%(\"\\nf2=f0+R/4.pi.L = \",f2,\" Hz\");\n",
+ "## x initialisation \n",
+ "import math\n",
+ "%matplotlib inline\n",
+ "import warnings\n",
+ "warnings.filterwarnings('ignore')\n",
+ "from math import log\n",
+ "import numpy\n",
+ "from math import tan\n",
+ "import matplotlib\n",
+ "from matplotlib import pyplot\n",
+ "x=numpy.array([-1,0.1,6.28]);\n",
+ "##simple plot\n",
+ "y=numpy.sin(x)\n",
+ "pyplot.plot(x,y)\n",
+ "\n",
+ "\n",
+ "\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "R=10 Ohm \n",
+ "L=0.2 H \n",
+ "C=40uF \n",
+ "V=100 V\n",
+ "\n",
+ "(i) f0= 1/2*pi*sqrt(LC) = 56.27 Hz\n",
+ "\n",
+ "(ii) I0= V/R = 10.00 A\n",
+ "\n",
+ "(iii) P0=(I0^2)*R = 1000.00 W\n",
+ "\n",
+ "(iv) pf=1\n",
+ "\n",
+ "(v) Rv = R.I = 100.00 V\n",
+ "\n",
+ " Lv = XL.I = 707.11 V\n",
+ "\n",
+ " Cv = XC.I = 707.11 V\n",
+ "\n",
+ "(vi) Q =1/R*sqrt(L/C)= 7.07 \n",
+ "\n",
+ "(vii)f1 = f0-R/4.pi.L = 52.29 Hz\n",
+ "\n",
+ "f2=f0+R/4.pi.L = 60.25 Hz\n"
+ ]
+ },
+ {
+ "metadata": {},
+ "output_type": "pyout",
+ "prompt_number": 37,
+ "text": [
+ "[<matplotlib.lines.Line2D at 0x5aad050>]"
+ ]
+ },
+ {
+ "metadata": {},
+ "output_type": "display_data",
+ "png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAEACAYAAAC9Gb03AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFdpJREFUeJzt3X20XHV97/H3l0SKQRBTbB5MaCKQq1UQ8KFRRAdM7gpc\nSu0ft8oqldtru/yjtN57266mupaevyy61l3l2qerqL3Yax9WKZeFokuSwxkB0fD8IElIUcE83JyQ\nBKJIXU3g2z9mAicnM+fMzJ5zZs+e92utWcye2TP7S7Lymb2/e+/fLzITSVJ1nTDoAiRJc8ugl6SK\nM+glqeIMekmqOINekirOoJekiisc9BGxISK2R8S/RMQft3j/NyLi4Yh4JCK+HRHnFt2mJKlzUeQ6\n+ohYADwOrAN2A/cCV2bmtinrvBPYmpmHImIDMJaZa4uVLUnqVNE9+ncAT2Tmk5l5GPgH4FenrpCZ\n38nMQ83FLcCKgtuUJHWhaNC/Dtg5ZXlX87V2Pgx8veA2JUldWFjw8x33fSLiYuC/AhcW3KYkqQtF\ng343sHLK8koae/XHaJ6AvR7YkJnPtPqiiHDQHUnqQWbGTO8Xbd3cB5wdEasi4kTgA8AtU1eIiDOA\nm4CrMvOJWYot/eOTn/zkwGuoSp3DUKN1WmfZH50otEefmUci4hrgm8AC4IuZuS0iPtJ8/3PAJ4DX\nAH8dEQCHM/MdRbYrSepc0dYNmfkN4BvTXvvclOe/Dfx20e1IknrjnbFdqtVqgy6hI8NQ5zDUCNbZ\nb9Y5/wrdMNVPEZFlqUWShkVEkHN8MlaSVHIGvSRVnEEvSRVn0EtSxRn0klRxBr0kVZxBL0kVZ9BL\nUsUZ9JJUcQa9JFWcQS9JFVd49MpRtmsXPPUU/PzPNx6veQ0s9E9UUsk4qFkBl17aCPrDh+HAATh0\nCE499eXg7/SxaBHEjEMSSVJrnQxqZtD36N/+DU4/HZ58EhYvbrz2wgvwzDON0D94sPHfTh6Zxwb/\n4sWz/zgsXgwLFgz0j0BSCXQS9DYaevTd78KaNS+HPDSC9/TTG49uPP98+x+B3bvhkUeOf92jB0md\nMuh7ND4O73tff75r0aLGY+XK2dc9aurRw/THwYOtfxxaHT20OlLw6EGqFls3PbrwQhgbg/XrB11J\nd2Y6emj38OhBKi979HPkJz+BZctg375GmFXdTEcP048kujl6aHUk4dGD1B179HPkjjvg7W8fjZCH\nuTn3sGsXPPywRw/SfDDoe7B5M6xbN+gqyq/f5x4OHDj+x+HoUUQ3Rw9T73vw6EGjwNZND849Fz7/\neVi7dtCV6Khezj38+McvHz10ckmrRw8qI3v0c2DfvsZllfv3exfssOv03MP0B3j0oPKwRz8Hbr8d\n3vteQ74Kejn3kHns0cP0E9Dtzj1MPXro5sY4jx7UD8ZVlzZv7t/18xo+EXDyyY3HGWd0/rluzz14\n9KB+snXTpdWr4Wtfgze9adCVqOqmHz20O4ro5Oihk5vjPHoYTvbo++wHP2jcKLVnj/8gVF6eexgt\n89Kjj4gNwHXAAuALmfnpFut8FrgUeB74L5n5YNHtDsLRto0hrzLrx7mH6Y+dOzs/99DJkYRHD/Or\nUNBHxALgL4B1wG7g3oi4JTO3TVnnMuCszDw7In4Z+GtgKC9MHB+HDRsGXYXUf3N17mHnTo8eyqBQ\n6yYi3gl8MjM3NJc3AmTmtVPW+d/ARGb+Y3N5O/DezJyc9l2lbt28+CIsWQL339/dPwRJx5rt6KHd\nOYlejh6Onnuosvlo3bwO2DlleRfwyx2sswKYZIg88khjj8KQl4qZ76OHiO4uaa3i0UPRoO90F3z6\nr015d93b6OewxJK6N1fnHh56qPW5h1e/ursjhxUrynveoWjQ7wamjmSyksYe+0zrrGi+dpyxsbGX\nntdqNWq1WsHy+md8HD784UFXIakbvR49HDnSOHpodynr9KOHgwcb04rOx1FAvV6nXq939ZmiPfqF\nwOPA+4A9wD3AlS1Oxl6TmZdFxFrgusw87mRsmXv0raYNlKQymPMefWYeiYhrgG/SuLzyi5m5LSI+\n0nz/c5n59Yi4LCKeAH4K/FaRbQ7Cli1w9tmGvKTh5A1THRgba/T5PvOZQVciScfqZI/+hPkqZpg5\n/rykYeYe/Syeew6WLh2daQMlDRf36PvgjjvgbW8z5CUNL4N+FuPjtm0kDTeDfhaOPy9p2Nmjn8G+\nfY3LKg8ccEYpSeVkj76giQmnDZQ0/Az6Gdi2kVQFBv0MHMhMUhUY9G384AeNu2GdG1bSsDPo2zi6\nN1/WYUclqVMGfRu2bSRVhZdXtvDii41hD+67zxmlJJWbl1f26NFHG7PLGPKSqsCgb8G2jaQqMehb\nMOglVYk9+mkOH25MG/j973c3CbEkDYI9+h5s2QJnnmnIS6oOg34a2zaSqsagn8bx5yVVjT36KY5O\nGzg5CSefPNBSJKkj9ui7dOed8Na3GvKSqsWgn8K2jaQqMuincPx5SVVkj77p6afhrLNg/354xSsG\nVoYkdcUefRcmJuA97zHkJVWPQd9k20ZSVRn0Td4oJamqDHrgyScb19C/+c2DrkSS+q9Q0EfE4ojY\nFBE7IuK2iDitxTorI2IiIh6LiO9FxO8X2eZccNpASVVWdI9+I7ApM9cA483l6Q4D/z0z3wSsBX43\nIt5YcLt9ZX9eUpUVurwyIrYD783MyYhYCtQz8w2zfOZm4M8zc3za6wO5vDKzMezBli2watW8b16S\nCpmPyyuXZOZk8/kksGSWglYB5wNbCm63b773PTjlFENeUnUtnG2FiNgELG3x1senLmRmRkTbXfKI\neBVwI/DRzHyu1TpjY2MvPa/VatRqtdnKK8y2jaRhUq/XqdfrXX2mH62bWmbujYhlwESr1k1EvAL4\nGvCNzLyuzXcNpHVz+eXwoQ/Br//6vG9akgqbj9bNLcDVzedXAze3KCKALwJb24X8oBw+3Bix8pJL\nBl2JJM2dokF/LbA+InYAlzSXiYjlEXFrc50LgauAiyPiweZjQ8Ht9sU998DrX++0gZKqbdYe/Uwy\n8yBw3MC+mbkH+E/N53dR0huzvBtW0igoZQDPF8eflzQKRnaY4p/+FJYscdpAScPNYYpncOedcMEF\nhryk6hvZoLdtI2lUjHTQeyJW0igYyR79/v1w5plOGyhp+Nmjb2NiAi66yJCXNBpGMuht20gaJSMZ\n9A5kJmmUjFzQP/UU/OQnThsoaXSMXNCPjzcGMTth5P7PJY2qkYs72zaSRs1IXV6ZCcuWwXe+A6tX\nz+mmJGleeHnlNI891hjywJCXNEpGKuht20gaRSMV9F4/L2kUjUyP/vDhxkxSTzwBr33tnG1GkuaV\nPfop7r230Zs35CWNmpEJets2kkbVSAW9489LGkUj0aM/Om3g3r3wqlfNySYkaSDs0TfddRecf74h\nL2k0jUTQ27aRNMpGJug9EStpVFW+R3/gQOOyyv374cQT+/71kjRQ9uh5edpAQ17SqKp80Nu2kTTq\nKh/0DmQmadT1HPQRsTgiNkXEjoi4LSJOm2HdBRHxYER8tdft9eJHP4Jnn4VzzpnPrUpSuRTZo98I\nbMrMNcB4c7mdjwJbgXk983u0beO0gZJGWZEIvAK4ofn8BuD9rVaKiBXAZcAXgBnPDPebbRtJKhb0\nSzJzsvl8EljSZr0/A/4IeLHAtrqWCbffbtBL0sKZ3oyITcDSFm99fOpCZmZEHNeWiYjLgX2Z+WBE\n1GYrZmxs7KXntVqNWm3Wj7S1dSu88pXw+tf3/BWSVDr1ep16vd7VZ3q+YSoitgO1zNwbEcuAicx8\nw7R1PgX8JnAEOAk4FfjnzPxQi+/r6w1Tn/0sPPooXH99375Skkpnrm+YugW4uvn8auDm6Stk5scy\nc2VmrgY+CNzeKuTngv15SWooEvTXAusjYgdwSXOZiFgeEbe2+cy8XHVz5AjccQdccsl8bE2Syq2S\nY91897vwkY/Aww/35eskqbRGdqwb2zaS9LJKBr3jz0vSyyrXunn+efiFX3DaQEmjYSRbN3fdBeed\nZ8hL0lGVC3rbNpJ0rEoGvSdiJelllerRHzwIq1Y5baCk0TFyPfqJCXj3uw15SZqqUkFv20aSjmfQ\nS1LFVSbod+5s9OjPPXfQlUhSuVQm6MfHG4OYOW2gJB2rMrFo20aSWqvE5ZWZsHx5467YM8/sc2GS\nVGIjc3nltm1w0klOGyhJrVQi6I+2bWLG3zRJGk2VCHrHn5ek9oa+R3/kCJx+Ojz+OCxZMgeFSVKJ\njUSP/v774YwzDHlJamfog962jSTNbOiD3vHnJWlmQ92j/9d/bUwbuGcPnHLKHBUmSSVW+R79t7/d\nGNvGkJek9oY66Ddvtm0jSbMZ6qB3fBtJmt3Q9uifeQZ+8RedNlDSaKt0j35iAt71LkNekmYztEHv\nZZWS1Jmegz4iFkfEpojYERG3RcRpbdY7LSJujIhtEbE1Itb2Xu7L7M9LUmeK7NFvBDZl5hpgvLnc\nyv8Cvp6ZbwTOBbYV2CYAu3Y1evNveUvRb5Kk6isS9FcANzSf3wC8f/oKEfFq4KLM/BJAZh7JzEMF\ntgk4baAkdaNIVC7JzMnm80mg1bBiq4GnI+JvIuKBiLg+IhYV2CZg20aSurFwpjcjYhOwtMVbH5+6\nkJkZEa2ujVwIXABck5n3RsR1NFo8n2i1vbGxsZee12o1arXacetkNm6U+kTLb5CkaqvX69Tr9a4+\n0/N19BGxHahl5t6IWAZMZOYbpq2zFPhOZq5uLr8b2JiZl7f4vo6uo9+2DS69FH74Q2eUkqS5vo7+\nFuDq5vOrgZunr5CZe4GdEbGm+dI64LEC23TaQEnqUpGgvxZYHxE7gEuay0TE8oi4dcp6vwd8JSIe\npnHVzacKbNP+vCR1aaiGQDhyBF77Wti+3RmlJAkqOATCAw/AihWGvCR1Y6iC3raNJHVvqILe8ecl\nqXtD06M/Om3g7t1w6qnzWJgklVilevR33w3nnGPIS1K3hibobdtIUm+GJug9EStJvRmKHv3RaQOf\nfhp+7ufmuTBJKrHK9OjrdXjnOw15SerFUAS90wZKUu+GJujtz0tSb0of9Lt3w759cN55g65EkoZT\n6YPeaQMlqZjSx6dtG0kqptRBn2nQS1JRpQ76xx+HBQvgrLMGXYkkDa9SB73TBkpScUMR9JKk3pV2\nCIQXXmhMG7h1KyxdOsDCJKnEhnoIhAcegOXLDXlJKqq0QW/bRpL6o7RB7/jzktQfpezR/+xnjf68\n0wZK0syGtkd/993w5jcb8pLUD6UM+s2b7c9LUr+UMugdf16S+qd0Pfpnn4WVK2H/fmeUkqTZDGWP\n/lvfctpASeqnnoM+IhZHxKaI2BERt0XEaW3W+5OIeCwiHo2Iv4uIGSPcyyolqb+K7NFvBDZl5hpg\nvLl8jIhYBfwOcEFmngMsAD4405d6o5Qk9VeRoL8CuKH5/Abg/S3W+TFwGFgUEQuBRcDudl+4Zw9M\nTjptoCT1U5GgX5KZk83nk8CS6Stk5kHgfwI/AvYAz2bm5nZfOD4OF1/cGINektQfC2d6MyI2Aa2G\nFfv41IXMzIg47vKdiDgT+G/AKuAQ8E8R8RuZ+ZVW27NtI0n9N2PQZ+b6du9FxGRELM3MvRGxDNjX\nYrW3AXdn5oHmZ24C3gW0DPqbbhrjlFNgbAxqtRq1Wq3D/w1JGg31ep16vd7VZ3q+jj4iPgMcyMxP\nR8RG4LTM3DhtnbfQCPW3Az8D/g9wT2b+ZYvvy5Urk6eeckYpSerUXF9Hfy2wPiJ2AJc0l4mI5RFx\nK0BmPgx8GbgPeKT5uc+3+0KnDZSk/ivVnbF/+7fJVVcNuhJJGh6d7NGXKuj37EmWLRt0JZI0PIYu\n6MtSiyQNi6Ec60aS1F8GvSRVnEEvSRVn0EtSxRn0klRxBr0kVZxBL0kVZ9BLUsUZ9JJUcQa9JFWc\nQS9JFWfQS1LFGfSSVHEGvSRVnEEvSRVn0EtSxRn0klRxBr0kVZxBL0kVZ9BLUsUZ9JJUcQa9JFWc\nQS9JFWfQS1LFGfSSVHEGvSRVXM9BHxH/OSIei4gXIuKCGdbbEBHbI+JfIuKPe92eJKk3RfboHwV+\nDbij3QoRsQD4C2AD8EvAlRHxxgLbHLh6vT7oEjoyDHUOQ41gnf1mnfOv56DPzO2ZuWOW1d4BPJGZ\nT2bmYeAfgF/tdZtlMCx/+cNQ5zDUCNbZb9Y5/+a6R/86YOeU5V3N1yRJ82ThTG9GxCZgaYu3PpaZ\nX+3g+7OnqiRJfROZxbI4IiaAP8jMB1q8txYYy8wNzeU/AV7MzE+3WNcfBUnqQWbGTO/PuEffhXYb\nuQ84OyJWAXuADwBXtlpxtkIlSb0pcnnlr0XETmAtcGtEfKP5+vKIuBUgM48A1wDfBLYC/5iZ24qX\nLUnqVOHWjSSp3EpzZ2ynN2ANyjDc+BURX4qIyYh4dNC1zCQiVkbERPPv+3sR8fuDrqmViDgpIrZE\nxEMRsTUi/nTQNbUTEQsi4sGI6OQiiYGJiCcj4pFmrfcMup5WIuK0iLgxIrY1/97XDrqm6SLiPzT/\nDI8+Ds3076g0e/QR8QbgReBztDm5OyjNG78eB9YBu4F7gSvL1oaKiIuA54AvZ+Y5g66nnYhYCizN\nzIci4lXA/cD7y/bnCRARizLz+YhYCNwF/GFm3jXouqaLiP8BvBU4JTOvGHQ97UTED4G3ZubBQdfS\nTkTcAHwrM7/U/Hs/OTMPDbqudiLiBBq59I7M3NlqndLs0Xd4A9agDMWNX5l5J/DMoOuYTWbuzcyH\nms+fA7YBywdbVWuZ+Xzz6YnAAqB0ARURK4DLgC/Q/sKIMiltjRHxauCizPwSNM4zljnkm9YB328X\n8lCioC85b/yaI80rss4Htgy2ktYi4oSIeAiYBCYyc+uga2rhz4A/onFEXHYJbI6I+yLidwZdTAur\ngacj4m8i4oGIuD4iFg26qFl8EPi7mVaY16CPiE0R8WiLx6/MZx09KEd/q2KabZsbgY829+xLJzNf\nzMzzgBXAeyKiNuCSjhERlwP7MvNBSrynPMWFmXk+cCnwu812Y5ksBC4A/iozLwB+CmwcbEntRcSJ\nwK8A/zTTev26jr4jmbl+PrfXR7uBlVOWV9LYq1ePIuIVwD8D/zczbx50PbPJzEPNy4bfBtQHXM5U\n7wKuiIjLgJOAUyPiy5n5oQHX1VJm/v/mf5+OiP9Hoy1652CrOsYuYFdm3ttcvpESBz2NH8z7M/Pp\nmVYqa+umbHsmL9341fwF/QBwy4BrGloREcAXga2Zed2g62knIk6PiNOaz18JrAceHGxVx8rMj2Xm\nysxcTeMQ/vayhnxELIqIU5rPTwb+I41RcEsjM/cCOyNiTfOldcBjAyxpNlcCfz/bSqUJ+nY3YJXB\nsNz4FRF/D9wNrImInRHxW4OuqY0LgauAi6dcHrZh0EW1sAy4vdmj3wJ8NTPHB1zTbMrcZlwC3Dnl\nz/NrmXnbgGtq5feAr0TEw8C5wKcGXE9LzR/LdcBNs65blssrJUlzozR79JKkuWHQS1LFGfSSVHEG\nvSRVnEEvSRVn0EtSxRn0klRxBr0kVdy/Axa1VVXVraQUAAAAAElFTkSuQmCC\n",
+ "text": [
+ "<matplotlib.figure.Figure at 0x5932610>"
+ ]
+ }
+ ],
+ "prompt_number": 37
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex81-pg4.66"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.81 :(pg 4.66)\n",
+ "import math\n",
+ "V=200.;\n",
+ "Vc=5000.;\n",
+ "I0=20.;\n",
+ "C=4.*10**-6;\n",
+ "R=V/I0;\n",
+ "Xco=Vc/I0;\n",
+ "f0=(1./(2.*math.pi*Xco*C));\n",
+ "L=(Xco/(2.*math.pi*f0));\n",
+ "print(\"\\nV=200 V \\nI0= 20 A \\nVc=5000 V \\nC=4uF\");\n",
+ "print'%s %.2f %s'%(\"\\nR=V/I0 = \",R,\" Ohm\");##resistance\n",
+ "print'%s %.2f %s'%(\"\\nXco=Vco/Io = \",Xco,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nXco=1/2*pi*f0*C \\nf0=1/2*pi*Xco*C = \",f0,\" Hz\");\n",
+ "print'%s %.2f %s'%(\"\\nat resonance Xco=Xlo \\nXlo= \",Xco,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nXlo=2*pi*f0*L \\nL= \",L,\" H\");\n",
+ "\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "V=200 V \n",
+ "I0= 20 A \n",
+ "Vc=5000 V \n",
+ "C=4uF\n",
+ "\n",
+ "R=V/I0 = 10.00 Ohm\n",
+ "\n",
+ "Xco=Vco/Io = 250.00 Ohm\n",
+ "\n",
+ "Xco=1/2*pi*f0*C \n",
+ "f0=1/2*pi*Xco*C = 159.15 Hz\n",
+ "\n",
+ "at resonance Xco=Xlo \n",
+ "Xlo= 250.00 Ohm\n",
+ "\n",
+ "Xlo=2*pi*f0*L \n",
+ "L= 0.25 H\n"
+ ]
+ }
+ ],
+ "prompt_number": 38
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex82-pg4.66"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.82 :(pg 4.66)\n",
+ "import math\n",
+ "V=230.;\n",
+ "f0=50.;\n",
+ "I0=2.;\n",
+ "Vco=500.;\n",
+ "R=V/I0;\n",
+ "Xco=Vco/I0;\n",
+ "C=(1/(2.*math.pi*f0*Xco));\n",
+ "L=(Xco/(2.*math.pi*f0));\n",
+ "print(\"\\nV = 230 V \\nf0 = 50 Hz \\nI0 = 2A \\nVco = 500 V\");\n",
+ "print'%s %.2f %s'%(\"\\nR=V/I0 = \",R,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nXco=Vco/I0 = \",Xco,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nXco=1/2.pi.f0.C \\nC= \",C,\" F\");##capacitance\n",
+ "print'%s %.2f %s'%(\"\\nXco=Xlo \\nXlo= \",Xco,\" Ohm\");##at resonance\n",
+ "print'%s %.2f %s'%(\"\\nXlo=2.pi.f0.L \\nL= \",L,\" H\");##inductance\n",
+ "\n",
+ "\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "V = 230 V \n",
+ "f0 = 50 Hz \n",
+ "I0 = 2A \n",
+ "Vco = 500 V\n",
+ "\n",
+ "R=V/I0 = 115.00 Ohm\n",
+ "\n",
+ "Xco=Vco/I0 = 250.00 Ohm\n",
+ "\n",
+ "Xco=1/2.pi.f0.C \n",
+ "C= 0.00 F\n",
+ "\n",
+ "Xco=Xlo \n",
+ "Xlo= 250.00 Ohm\n",
+ "\n",
+ "Xlo=2.pi.f0.L \n",
+ "L= 0.80 H\n"
+ ]
+ }
+ ],
+ "prompt_number": 39
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex83-pg4.67"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.82 :(pg 4.66)\n",
+ "import math\n",
+ "R=2.;\n",
+ "L=0.01;\n",
+ "V=200.;\n",
+ "f0=50.;\n",
+ "C=(1./(4.*(math.pi)**2*L*(f0**2)));\n",
+ "I0=V/R;\n",
+ "Vco=I0*(1./(2.*math.pi*f0*C));\n",
+ "print(\"\\nR= 2 Ohm \\nL= 0.01 H \\nV=200 V \\nf0=50 Hz \\nf0=1/(2.pi.sqrt(LC)\");\n",
+ "print'%s %.2f %s'%(\"\\nC = \",C,\" F\");##capacitance\n",
+ "print'%s %.2f %s'%(\"\\nI0= V/R = \",I0,\" A\");##current\n",
+ "print'%s %.2f %s'%(\"\\nVco=I0.Xco \\n= \",Vco,\" V\"); ##voltage across capacitor\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "R= 2 Ohm \n",
+ "L= 0.01 H \n",
+ "V=200 V \n",
+ "f0=50 Hz \n",
+ "f0=1/(2.pi.sqrt(LC)\n",
+ "\n",
+ "C = 0.00 F\n",
+ "\n",
+ "I0= V/R = 100.00 A\n",
+ "\n",
+ "Vco=I0.Xco \n",
+ "= 314.16 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 41
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex84-pg4.67"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.84 :(pg 4.67)\n",
+ "import math\n",
+ "BW=400.;\n",
+ "Vco=500.;\n",
+ "R=100.;\n",
+ "Vm=10.;\n",
+ "V=(Vm/math.sqrt(2.));\n",
+ "I0=V/R;\n",
+ "L=R/BW;\n",
+ "Q0=Vco/V;\n",
+ "C=(L/(Q0*R)**2);\n",
+ "f0=(1/(2.*math.pi*math.sqrt(L*C)));\n",
+ "f1=(f0-(R/(4.*math.pi*L)));##lower cut-off frequency\n",
+ "f2=(f0+(R/(4.*math.pi*L)));##upper cut-off frequency\n",
+ "print(\"\\nv(t)=10sinwt \\nVco=5000V \\nBW=400rad/s \\nR=100 Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nV= \",V,\" V\");\n",
+ "print'%s %.2f %s'%(\"\\nI0=V/R= \",I0,\" A\");\n",
+ "print'%s %.2f %s'%(\"\\nBW=R/L \\nL= \",L,\" H\");\n",
+ "print'%s %.2f %s'%(\"\\nQ0=Vco/V =\",Q0,\"\");\n",
+ "print'%s %.2e %s'%(\"\\nQ0=1/R*sqrt(L/C) \\nC= \",C,\" F\");\n",
+ "print'%s %.2f %s'%(\"\\nf0=1/2.pi.sqrt(LC)= \",f0,\" Hz\");\n",
+ "print'%s %.2f %s'%(\"\\nf1=f0-R/4.pi.L = \",f1,\" Hz\");##lower cut-off frequency\n",
+ "print'%s %.2f %s'%(\"\\nf2=f0+R/4.pi.L = \",f2,\" Hz\"); ##upper cut-off frequency"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "v(t)=10sinwt \n",
+ "Vco=5000V \n",
+ "BW=400rad/s \n",
+ "R=100 Ohm\n",
+ "\n",
+ "V= 7.07 V\n",
+ "\n",
+ "I0=V/R= 0.07 A\n",
+ "\n",
+ "BW=R/L \n",
+ "L= 0.25 H\n",
+ "\n",
+ "Q0=Vco/V = 70.71 \n",
+ "\n",
+ "Q0=1/R*sqrt(L/C) \n",
+ "C= 5.00e-09 F\n",
+ "\n",
+ "f0=1/2.pi.sqrt(LC)= 4501.58 Hz\n",
+ "\n",
+ "f1=f0-R/4.pi.L = 4469.75 Hz\n",
+ "\n",
+ "f2=f0+R/4.pi.L = 4533.41 Hz\n"
+ ]
+ }
+ ],
+ "prompt_number": 42
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex85-pg4.68"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.85 :(pg 4.68)\n",
+ "import math\n",
+ "R=500.;\n",
+ "f1=100.;\n",
+ "f2=10.*10**3;\n",
+ "BW=f2-f1;\n",
+ "f0=((f1+f2)/2.);\n",
+ "L=(R/(2.*math.pi*BW));\n",
+ "XL0=(2.*math.pi*f0*L);\n",
+ "C=(1/(2.*math.pi*f0*XL0));\n",
+ "Q0=((1./R)*(math.sqrt(L/C)));\n",
+ "print'%s %.2f %s'%(\"\\nR= 500 Ohm \\nf1 = 100 Hz \\nf2=10kHz \\nBW= f2-f1 = \",BW,\" Hz\");\n",
+ "print'%s %.2f %s'%(\"\\nf1=f0-BW/2 ------(i) \\nf2=f0+BW/2 ------(ii) \\nf1+f2 =2f0 \\nf0=(f1+f2)/2 = \",f0,\" Hz\");\n",
+ "print'%s %.2f %s'%(\"\\nBW=R/2.pi.f0.L \\nL= \",L,\" H\");\n",
+ "print'%s %.2f %s'%(\"\\nXL0=2.pi.f0.L = \",XL0,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nXL0=XC0 = \",XL0,\" Ohm\");##at resonance\n",
+ "print'%s %.2e %s'%(\"\\nXC0 =1/2.pi.f0.C \\nC= \",C,\" F\");\n",
+ "print'%s %.2f %s'%(\"\\nQ0=(1/R*sqrt(L/C)) =\",Q0,\"\");\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "R= 500 Ohm \n",
+ "f1 = 100 Hz \n",
+ "f2=10kHz \n",
+ "BW= f2-f1 = 9900.00 Hz\n",
+ "\n",
+ "f1=f0-BW/2 ------(i) \n",
+ "f2=f0+BW/2 ------(ii) \n",
+ "f1+f2 =2f0 \n",
+ "f0=(f1+f2)/2 = 5050.00 Hz\n",
+ "\n",
+ "BW=R/2.pi.f0.L \n",
+ "L= 0.01 H\n",
+ "\n",
+ "XL0=2.pi.f0.L = 255.05 Ohm\n",
+ "\n",
+ "XL0=XC0 = 255.05 Ohm\n",
+ "\n",
+ "XC0 =1/2.pi.f0.C \n",
+ "C= 1.24e-07 F\n",
+ "\n",
+ "Q0=(1/R*sqrt(L/C)) = 0.51 \n"
+ ]
+ }
+ ],
+ "prompt_number": 43
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex87-pg4.69"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.87 :(pg 4.69 & 4.70)\n",
+ "import math\n",
+ "f0=10**6;\n",
+ "C1=500.*10**-12;\n",
+ "C2=600.*10**-12;\n",
+ "C=500.*10**-12;\n",
+ "x=((2.*math.pi*f0)**2);\n",
+ "L=(1./(x*C));\n",
+ "XL=(2.*math.pi*f0*L);\n",
+ "y=2.*math.pi*f0*C2;\n",
+ "XC=(1./y);\n",
+ "R=math.sqrt(((XL-XC)**2)/3.);\n",
+ "x=math.sqrt(L/C);\n",
+ "Q0=((1./R)*x);\n",
+ "print(\"\\nf0= 1MHz \\nC1=500pF \\nC2=600pF \\nC=500pF\");##At resonance\n",
+ "print'%s %.2e %s'%(\"\\nf0=1/2.pi.sqrt(LC)\\nL= \",L,\" H\");\n",
+ "print'%s %.2f %s'%(\"\\nXL=2.pi.f0.L = \",XL,\" Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nXC=1/2.pi.f0.C \\nXC= \",XC,\" Ohm\");\n",
+ "print(\"\\nI=1/2.I0 \\nV/Z=1/2.V/R \\nZ=2R\");\n",
+ "print'%s %.2f %s'%(\"\\nsqrt((R^2)-(XL-XC)^2)=2R \\nR= \",R,\" Ohm\");##Resistance of Inductor\n",
+ "print'%s %.2f %s'%(\"\\nQ0=1/R.sqrt(L/C) \\n=\",Q0,\"\");\n",
+ "\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "f0= 1MHz \n",
+ "C1=500pF \n",
+ "C2=600pF \n",
+ "C=500pF\n",
+ "\n",
+ "f0=1/2.pi.sqrt(LC)\n",
+ "L= 5.07e-05 H\n",
+ "\n",
+ "XL=2.pi.f0.L = 318.31 Ohm\n",
+ "\n",
+ "XC=1/2.pi.f0.C \n",
+ "XC= 265.26 Ohm\n",
+ "\n",
+ "I=1/2.I0 \n",
+ "V/Z=1/2.V/R \n",
+ "Z=2R\n",
+ "\n",
+ "sqrt((R^2)-(XL-XC)^2)=2R \n",
+ "R= 30.63 Ohm\n",
+ "\n",
+ "Q0=1/R.sqrt(L/C) \n",
+ "= 10.39 \n"
+ ]
+ }
+ ],
+ "prompt_number": 44
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex88-pg4.72"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.88 :(pg 4.72)\n",
+ "import math\n",
+ "R=20.;\n",
+ "C=100.*10**-6;\n",
+ "L=0.2;\n",
+ "DR=(L/(C*R));\n",
+ "x=(1./(L*C));\n",
+ "y=((R/L)**2);\n",
+ "f0=((1./(2.*math.pi))*math.sqrt(x-y));\n",
+ "DR=(L/(C*R));\n",
+ "print(\"\\nR=20 Ohm \\nL=0.2 H \\nC=100uF\");\n",
+ "print'%s %.2f %s'%(\"\\nf0=1/2.pi.sqrt(1/LC-R^2/L^2) \\n= \",f0,\" Hz\");\n",
+ "print'%s %.2f %s'%(\"\\n dynamic resistance =L/CR \\n= \",DR,\" Ohm\");"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "R=20 Ohm \n",
+ "L=0.2 H \n",
+ "C=100uF\n",
+ "\n",
+ "f0=1/2.pi.sqrt(1/LC-R^2/L^2) \n",
+ "= 31.83 Hz\n",
+ "\n",
+ " dynamic resistance =L/CR \n",
+ "= 100.00 Ohm\n"
+ ]
+ }
+ ],
+ "prompt_number": 45
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex89-pg4.72"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##AC Circuits : example 4.89 :(pg 4.72 & 4.73)\n",
+ "import math\n",
+ "R=20.;\n",
+ "L=200.*10**-6;\n",
+ "f=10**6;\n",
+ "V=230.;\n",
+ "Rs=8000.;\n",
+ "XL=2.*math.pi*f*L;\n",
+ "x=((2.*math.pi*f)**2);\n",
+ "y=((R/L)**2);\n",
+ "C=(1./((x+y)*L));\n",
+ "Q=((2.*math.pi*f*L)/R);\n",
+ "Z=(L/(C*R));\n",
+ "ZT=(Rs+Z);\n",
+ "IT=(V/ZT);\n",
+ "print'%s %.2f %s'%(\"\\nR=20 Ohm \\nL=200uH \\nf=10^6 \\nV=230 V \\nRs=8000 Ohm \\nXL=2.pi.f.L = \",XL,\"Ohm\");\n",
+ "print'%s %.2f %s'%(\"\\nf0=1/2.pi.sqrt(1/LC-R^2/L^2) \\nC= \",C,\" F\");\n",
+ "print'%s %.2f %s'%(\"\\nQ0=2.pi.f.L/R =\",Q,\"\");##quality factor\n",
+ "print'%s %.2f %s'%(\"\\nZ=L/CR \\n \",Z,\" Ohm\");##dynamic impedance\n",
+ "print'%s %.2f %s'%(\"\\nZt= \",ZT,\" Ohm\");##total equivalent Z at resonance\n",
+ "print'%s %.2e %s'%(\"\\nIt= \",IT,\" A\");##total ckt current"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "R=20 Ohm \n",
+ "L=200uH \n",
+ "f=10^6 \n",
+ "V=230 V \n",
+ "Rs=8000 Ohm \n",
+ "XL=2.pi.f.L = 1256.64 Ohm\n",
+ "\n",
+ "f0=1/2.pi.sqrt(1/LC-R^2/L^2) \n",
+ "C= 0.00 F\n",
+ "\n",
+ "Q0=2.pi.f.L/R = 62.83 \n",
+ "\n",
+ "Z=L/CR \n",
+ " 78976.84 Ohm\n",
+ "\n",
+ "Zt= 86976.84 Ohm\n",
+ "\n",
+ "It= 2.64e-03 A\n"
+ ]
+ }
+ ],
+ "prompt_number": 46
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file |