summaryrefslogtreecommitdiff
path: root/Power_Electronics_Principles_and_Applications_by_Jacob
diff options
context:
space:
mode:
Diffstat (limited to 'Power_Electronics_Principles_and_Applications_by_Jacob')
-rw-r--r--Power_Electronics_Principles_and_Applications_by_Jacob/Chapter1.ipynb558
-rw-r--r--Power_Electronics_Principles_and_Applications_by_Jacob/Chapter2.ipynb65
-rw-r--r--Power_Electronics_Principles_and_Applications_by_Jacob/Chapter3.ipynb481
-rw-r--r--Power_Electronics_Principles_and_Applications_by_Jacob/Chapter4.ipynb733
-rw-r--r--Power_Electronics_Principles_and_Applications_by_Jacob/Chapter5.ipynb647
-rw-r--r--Power_Electronics_Principles_and_Applications_by_Jacob/Chapter6.ipynb528
-rw-r--r--Power_Electronics_Principles_and_Applications_by_Jacob/Chapter7.ipynb505
-rw-r--r--Power_Electronics_Principles_and_Applications_by_Jacob/Chapter8.ipynb326
-rw-r--r--Power_Electronics_Principles_and_Applications_by_Jacob/Chapter9.ipynb405
-rw-r--r--Power_Electronics_Principles_and_Applications_by_Jacob/README.txt10
-rw-r--r--Power_Electronics_Principles_and_Applications_by_Jacob/screenshots/4.pngbin0 -> 15157 bytes
-rw-r--r--Power_Electronics_Principles_and_Applications_by_Jacob/screenshots/5.pngbin0 -> 11087 bytes
-rw-r--r--Power_Electronics_Principles_and_Applications_by_Jacob/screenshots/6.pngbin0 -> 8235 bytes
13 files changed, 4258 insertions, 0 deletions
diff --git a/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter1.ipynb b/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter1.ipynb
new file mode 100644
index 00000000..662bf254
--- /dev/null
+++ b/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter1.ipynb
@@ -0,0 +1,558 @@
+{
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 1: Advanced Operational Amplifier Principles"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 1.1,Page 6"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 1,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "open output voltage is 0.5 V\n",
+ "resistance lower loaded is 333.333 ohm\n",
+ "loaded output voltage is 0.25 V\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding voltage current resistance\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "R1=1000.0;\n",
+ "R2=1000.0;\n",
+ "Rl=500.0#load resistance\n",
+ "V=1.0#input voltage\n",
+ "\n",
+ "#calculation\n",
+ "Vo=(R2/(R1+R2))*V;\n",
+ "Rll=1/((1/R2)+(1/Rl))#lower loaded resistance\n",
+ "Vol=(Rll/(R2+Rll))*V;\n",
+ "\n",
+ "#result\n",
+ "print \"open output voltage is\",round(Vo,3),\"V\"\n",
+ "print \"resistance lower loaded is\",round(Rll,3),\"ohm\"\n",
+ "print \"loaded output voltage is\",round(Vol,3),\"V\"\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 1.2,Page 11"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 1,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "input resistance is 1.01 Kohm\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding voltage current resistance\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Rf=100000.0#resistance\n",
+ "Acl=100.0#amplifier gain\n",
+ "\n",
+ "#calculation\n",
+ "Ri=Rf/(Acl-1);\n",
+ "\n",
+ "#result\n",
+ "print \"input resistance is\",round(Ri/1000,2), \"Kohm\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 1.3,Page 17"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 3,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "current through Ri1 is 178.571 microAmp\n",
+ "current through Ri2 is 31.915 microAmp\n",
+ "current through Ri2 is 31.915 microAmp\n",
+ "current through Rf is 210.486 microAmp\n",
+ "voltage dropped is 2.105 V\n",
+ "output voltage 1 is -2.105 V\n",
+ "output voltage is 2.105 V\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding voltage current resistance\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Vni=0.0#non inverting voltage\n",
+ "Vinv=0.0;#inverting voltage\n",
+ "Vri1=1.0;\n",
+ "Vri2=15.0;\n",
+ "Ri1=5600.0#resistance\n",
+ "Ri2=470000.0;\n",
+ "Rf=10000.0#load resistance\n",
+ "\n",
+ "#calculation\n",
+ "Ir1=Vri1/Ri1;\n",
+ "Ir2=Vri2/Ri2;\n",
+ "Irf=(Vri1/Ri1)+(Vri2/Ri2);\n",
+ "Vr=Irf*Rf;\n",
+ "Vo1=-Vr;\n",
+ "Vo=Irf*Rf;\n",
+ "\n",
+ "#result\n",
+ "print \"current through Ri1 is\",round(Ir1*1e6,3), \"microAmp\"\n",
+ "print \"current through Ri2 is\",round(Ir2*1e6,3), \"microAmp\"\n",
+ "print \"current through Ri2 is\",round(Ir2*1e6,3),\"microAmp\"\n",
+ "print \"current through Rf is\",round(Irf*1e6,3), \"microAmp\"\n",
+ "print \"voltage dropped is\",round(Vr,3), \"V\"\n",
+ "print \"output voltage 1 is\",round(Vo1,3), \"V\"\n",
+ "print \"output voltage is\",round(Vo,3), \"V\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 1.4,Page 25"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 1,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "inverting voltage is 4.955 V\n",
+ "non inverting voltage is 4.955 V\n",
+ "current through Rf2 is 42.698 microA\n",
+ "current through Ri2 is 42.698 microA\n",
+ "voltage dropped is 4.056 V\n",
+ "output voltage is 884.897 mV\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding voltage current resistance\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Ri1=950.00;#ohm\n",
+ "Ri2=1050.00;\n",
+ "Rf1=105000.00;#resistance\n",
+ "Rf2=95000.00;\n",
+ "Vin=5.00;#voltage\n",
+ "\n",
+ "#calculation\n",
+ "Vinv=(Rf1/(Rf1+Ri1))*Vin;\n",
+ "Vni=Vinv;\n",
+ "Irf2=(Vin-Vinv)/Ri2;\n",
+ "Iri2=Irf2;\n",
+ "Vrf2=Irf2*Rf2;\n",
+ "Vo=Vinv-Vrf2-.014;\n",
+ "\n",
+ "#result\n",
+ "print \"inverting voltage is\",round(Vinv,3), \"V\"\n",
+ "print \"non inverting voltage is\",round (Vni,3), \"V\"\n",
+ "print \"current through Rf2 is\",round(Irf2*1e6,3), \"microA\"\n",
+ "print \"current through Ri2 is\",round(Iri2*1e6,3), \"microA\"\n",
+ "print \"voltage dropped is\",round(Vrf2,3), \"V\"\n",
+ "print \"output voltage is\",round(Vo*1000,3), \"mV\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 1.5,Page 27"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 5,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "input resistor current is 272.222 microA\n",
+ "input resistor current is 500.0 microA\n",
+ "feedback resistor current is 227.778 microAmp\n",
+ "resistor voltage is 227.778 mV\n",
+ "1st output voltage is 2.222 V\n",
+ "input resistor current is 327.778 microA\n",
+ "input resistor current is 827.778 microA\n",
+ "feedback resistor voltage is 7.45 V\n",
+ "2nd output voltage is 10.0 V\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding voltage current resistance \n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Vniu1=2.45;#V\n",
+ "Vniu2=2.55;#V\n",
+ "Vinvu1=2.45;\n",
+ "Vinvu2=2.55;\n",
+ "Ri1=9000.0;#ohm\n",
+ "Ri2=1000.0;#ohm\n",
+ "Rf1=1000.0;\n",
+ "Rf2=9000.0;\n",
+ "Rg=200.0;#load resistance\n",
+ "\n",
+ "#calculation\n",
+ "Iri1=Vniu1/Ri1;\n",
+ "Irg=(Vniu2-Vniu1)/Rg;\n",
+ "Irf1=Irg-Iri1;\n",
+ "Vrf1=Irf1*Rf1;\n",
+ "Vou1=Vniu1-Vrf1;\n",
+ "Iri2=(Vniu2-Vou1)/Ri2;\n",
+ "Irf2=Iri2+Irg;\n",
+ "Vrf2=Irf2*Rf2#feedback resistor voltage\n",
+ "Vo=Vrf2+Vniu2;\n",
+ "\n",
+ "#result\n",
+ "print \"input resistor current is\",round(Iri1*1e6,3), \"microA\"\n",
+ "print \"input resistor current is\",round(Irg*1e6,3), \"microA\"\n",
+ "print \"feedback resistor current is\",round(Irf1*1e6,3), \"microAmp\"\n",
+ "print \"resistor voltage is\",round(Vrf1*1000,3), \"mV\"\n",
+ "print \"1st output voltage is\",round(Vou1,3), \"V\"\n",
+ "print \"input resistor current is\",round(Iri2*1e6,3), \"microA\"\n",
+ "print \"input resistor current is\",round(Irf2*1e6,3),\"microA\"\n",
+ "print \"feedback resistor voltage is\",round(Vrf2,3), \"V\"\n",
+ "print \"2nd output voltage is\",round(Vo,3), \"V\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 1.6.a,Page 29"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 2,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "input resistor current is 128.0 microA\n",
+ "feedback resistor current is 128.0 microA\n",
+ "feedback resistor voltage is 5.018 V\n",
+ "output resistor voltage is 5.018 V\n",
+ "output voltage is 3.818 V\n",
+ "load current is 0.5 A\n",
+ "load power is 2.5 W\n",
+ "power dissipated in LM317 is 5.0 W\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding voltage current resistance\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Vniu1=0;#V\n",
+ "Vinvu2=0;#V\n",
+ "Vref=2.56;\n",
+ "Rl=10000.0;#ohm\n",
+ "Rf=39200.0;#ohm\n",
+ "Ro=10.0;#resistance\n",
+ "Vdc1=5.0;\n",
+ "Vdc2=15.0;\n",
+ "Idc=0.5;#current\n",
+ "\n",
+ "#calculation\n",
+ "Iu1=(Vref/Rl)*.5;\n",
+ "Irf=Iu1;\n",
+ "Vrf=Irf*Rf;\n",
+ "Vout=Vrf+Vinvu2;\n",
+ "Eo=Vout-1.2;\n",
+ "Iload=Vdc1/Ro;\n",
+ "Pload=Vdc1**2/Ro;\n",
+ "Plm317=(Vdc2-Vdc1)*Idc;\n",
+ "\n",
+ "#result\n",
+ "print \"input resistor current is\",round(Iu1*1e6,3), \"microA\"\n",
+ "print \"feedback resistor current is\",round(Irf*1e6,3), \"microA\"\n",
+ "print \"feedback resistor voltage is\",round(Vrf,3), \"V\"\n",
+ "print \"output resistor voltage is\",round(Vout,3), \"V\"\n",
+ "print \"output voltage is\",round(Eo,3), \"V\"\n",
+ "print \"load current is\",round(Iload,3), \"A\"\n",
+ "print \"load power is\",round(Pload,3), \"W\"\n",
+ "print \"power dissipated in LM317 is\",round(Plm317,3), \"W\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 1.6.b,Page 31"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 7,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "input resistor current is 360.36 microamp\n",
+ "inverting voltage 1 & 2 is 396.396 mV\n",
+ "current across Rs is 3.964 A\n",
+ "emitter voltage is 8.324 V\n",
+ "output voltage is 10.124 V\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding voltage current resistance\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Vin=4;#V\n",
+ "Vs=1.8;#V\n",
+ "Rf=10000.0;#ohm\n",
+ "Ri=1100.0;#ohm\n",
+ "Rl=2.0;#ohm\n",
+ "Rs=0.1;#ohm\n",
+ "\n",
+ "#calculation\n",
+ "Irf=Vin/(Rf+Ri);\n",
+ "Vni=Irf*Ri;\n",
+ "Ir=Vni/Rs;\n",
+ "Ve=Ir*(Rl+Rs);\n",
+ "Vo=Ve+Vs;\n",
+ "\n",
+ "#result\n",
+ "print \"input resistor current is\",round(Irf*1e6,3),\"microamp\"\n",
+ "print \"inverting voltage 1 & 2 is\",round(Vni*1000,3), \"mV\"\n",
+ "print \"current across Rs is\",round(Ir,3), \"A\"\n",
+ "print \"emitter voltage is\",round(Ve,3), \"V\"\n",
+ "print \"output voltage is\",round(Vo,3), \"V\"\n",
+ "\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 1.7,Page 36"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 8,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "rms voltage is 9.899 V\n",
+ "power delivered is 12.25 W\n",
+ "load voltage is 28.284 V\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding voltage and power\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Vs=18.0;#V\n",
+ "Rl=8.0;#load resistance\n",
+ "Pll=100.0;#power\n",
+ "\n",
+ "#calculation\n",
+ "Vlp=Vs-4;\n",
+ "Vlr=Vlp/(2**(.5));\n",
+ "Pl=(Vlr**2)/Rl;\n",
+ "Vl=(Pll*Rl)**(.5);\n",
+ "\n",
+ "#result\n",
+ "print \"rms voltage is\",round(Vlr,3), \"V\"\n",
+ "print \"power delivered is\",round(Pl,3), \"W\"\n",
+ "print \"load voltage is\",round(Vl,3), \"V\"\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 1.9,Page 44"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 2,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "output voltage is 37.34 V\n",
+ "V+ is 45.34 V ;V- is 29.34 V\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding output volatage and range \n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "import numpy as np\n",
+ "Vp=6.0;#V\n",
+ "Ra=10.0;#Kohm\n",
+ "Rb=1800.0;#ohm\n",
+ "V=8.0;\n",
+ "#solving for Ir & Vo\n",
+ "a=np.array([[1.0,-124.6e-6],[7800.0,-1.0]])\n",
+ "b=np.array([134.6e-6,0.0])\n",
+ "\n",
+ "#calculation\n",
+ "x=np.linalg.solve(a,b);\n",
+ "Vo=x[1];\n",
+ "Va=Vo+V;\n",
+ "Vb=Vo-V;\n",
+ "\n",
+ "#result\n",
+ "print \"output voltage is\",round(Vo,2), \"V\"\n",
+ "print \"V+ is\",round(Va,2), \"V ;V- is\",round(Vb,2), \"V\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 1.11,Page 50"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 1,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "output current is 4.091 mA\n",
+ "output voltage is 45.409 V\n",
+ "gain output voltage 1 is 13.356 V\n",
+ "gain output voltage 2 is 0.38 V\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding output voltage and gain output voltage \n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Vin=4.5;\n",
+ "R1=1100.0;\n",
+ "R2=10000.0;\n",
+ "\n",
+ "G1=3.4#gain 1\n",
+ "G2=120.0#gain 2\n",
+ "\n",
+ "#calculation\n",
+ "Ir=Vin/R1;\n",
+ "Vo=Ir*(R1+R2);\n",
+ "Vuo1=Vo/G1;\n",
+ "Vuo2=Vo/G2;\n",
+ "\n",
+ "#result\n",
+ "print \"output current is\",round(Ir*1000,3),\"mA\"\n",
+ "print \"output voltage is\",round(Vo,3), \"V\"\n",
+ "print \"gain output voltage 1 is\",round(Vuo1,3), \"V\"\n",
+ "print \"gain output voltage 2 is\",round(Vuo2,2),\"V\""
+ ]
+ }
+ ],
+ "metadata": {
+ "kernelspec": {
+ "display_name": "Python 2",
+ "language": "python",
+ "name": "python2"
+ },
+ "language_info": {
+ "codemirror_mode": {
+ "name": "ipython",
+ "version": 2
+ },
+ "file_extension": ".py",
+ "mimetype": "text/x-python",
+ "name": "python",
+ "nbconvert_exporter": "python",
+ "pygments_lexer": "ipython2",
+ "version": "2.7.6"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter2.ipynb b/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter2.ipynb
new file mode 100644
index 00000000..9d88fe28
--- /dev/null
+++ b/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter2.ipynb
@@ -0,0 +1,65 @@
+{
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 2: Power Electronics Circuit Layout"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 2.9,Page 83"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 1,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "load current is 3.75 A\n",
+ "wiring resistance is 26.67 mohm\n",
+ "resistance per inch is 1666.67 microohm/inch\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding an appropriate wire gauge\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "V=15.0;#voltage\n",
+ "R=4.0;#resistance\n",
+ "Vl=.1;\n",
+ "D=8.0;#duty cycle\n",
+ "\n",
+ "#calculation\n",
+ "Il=V/R;\n",
+ "Rw=Vl/Il#wiring resistance\n",
+ "Ri=Rw/(2*D);\n",
+ "\n",
+ "#result\n",
+ "print \"load current is\",round(Il,2), \"A\"\n",
+ "print \"wiring resistance is\",round(Rw*1000,2), \"mohm\"\n",
+ "print \"resistance per inch is\",round(Ri*1e6,2), \"microohm/inch\""
+ ]
+ }
+ ],
+ "metadata": {
+ "kernelspec": {
+ "display_name": "Python 2",
+ "language": "python",
+ "name": "python2"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter3.ipynb b/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter3.ipynb
new file mode 100644
index 00000000..4129f088
--- /dev/null
+++ b/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter3.ipynb
@@ -0,0 +1,481 @@
+{
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 3 : Power Parameter Calculations"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 3.1,Page 109"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 1,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "ramp current is 450.0 kAt/s\n",
+ "current at 5 micro sec is 2.25 A\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding ramp current and current at 5 micro sec\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Ip=3.0;\n",
+ "f=150000.0;\n",
+ "t=5.0e-6;\n",
+ "\n",
+ "#calculation\n",
+ "T=1/f;\n",
+ "It=Ip/T;\n",
+ "I5=It*t;\n",
+ "\n",
+ "#result\n",
+ "print \"ramp current is\",round(It/1000,3), \"kAt/s\"\n",
+ "print \"current at 5 micro sec is\",round(I5,3), \"A\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 3.2,Page 110"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 2,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "current in time 0<=t<800ns is 3.575 A\n",
+ "current in time 800ns<=t<2 microsec is 0.0 A\n",
+ "current in time 400ns is 1.85 A\n",
+ "current in time 1 microsec is 0.0 A\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding current at different time\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Ip=2.0;\n",
+ "f=500000.0;\n",
+ "Ir=.3;\n",
+ "Cd=.4#duty cycle\n",
+ "t1=4.0e-7;\n",
+ "t2=1.0e-6;\n",
+ "I1=0;\n",
+ "\n",
+ "#calculation\n",
+ "T=1/f;\n",
+ "Im=Ip-Ir;\n",
+ "I4=(Ip-Im)*t1/(Cd*T)+Im;\n",
+ "It=(Ip-Im)*t/(Cd*T)+Im;\n",
+ "It1=0\n",
+ "\n",
+ "#resilt\n",
+ "print \"current in time 0<=t<800ns is\",round(It,3),\"A\"\n",
+ "print \"current in time 800ns<=t<2 microsec is\",round(It1,2), \"A\"\n",
+ "print \"current in time 400ns is\",round(I4,2), \"A\"\n",
+ "print \"current in time 1 microsec is\",round(I1,2), \"A\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 3.3,Page 115"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 3,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "average voltage is 54.02 V\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding average voltage\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Vr=120;\n",
+ "\n",
+ "#calculation\n",
+ "V=(Vr*2**.5)/pi;\n",
+ "\n",
+ "#result\n",
+ "print \"average voltage is\",round(V,2), \"V\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 3.4,Page 119"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 5,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "average current is 0.98 A\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding average current\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "f=100000.0;\n",
+ "Cd=.35#duty cycle\n",
+ "Ip=3.0;\n",
+ "Ir=.4;\n",
+ "\n",
+ "#calculation\n",
+ "Im=Ip-Ir;\n",
+ "T=1/f;\n",
+ "I=Cd*((Ip-Im)/2+Im)\n",
+ "\n",
+ "#result\n",
+ "print \"average current is\",round(I,2), \"A\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 3.5,Page 124"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 6,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "rms voltage is 8.87 V\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding rms voltage\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Vp=15.0;\n",
+ "Cd=.35;\n",
+ "f=100000.0;\n",
+ "\n",
+ "#calculation\n",
+ "V=Vp*Cd**.5;\n",
+ "\n",
+ "#result\n",
+ "print \"rms voltage is\",round(V,2), \"V\"\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 3.6,Page 127"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 7,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "rms current is 1.73 A\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding rms current\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Ip=3.0;\n",
+ "f=100000.0;\n",
+ "\n",
+ "#calculation\n",
+ "I=Ip/3**.5;\n",
+ "\n",
+ "#result\n",
+ "print \"rms current is\",round(I,2), \"A\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 3.7,Page 133"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 8,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "rms voltage is 85.0 V\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding rms voltage\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Vp=170.0;\n",
+ "f=60.0;\n",
+ "\n",
+ "#calculation\n",
+ "Vr=Vp/2;\n",
+ "\n",
+ "#result\n",
+ "print \"rms voltage is\",round(Vr,2), \"V\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 3.8,Page 140"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 9,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "power required is 2.42 hp\n",
+ "Pick a 5HP motor\n",
+ "current required is 18.84 amp\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding current and power\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "M=1000.0;\n",
+ "H=40.0;\n",
+ "T=30.0;\n",
+ "E1=.9;\n",
+ "E2=.5;\n",
+ "V=220.0;\n",
+ "P1=5.0;\n",
+ "\n",
+ "#calculation\n",
+ "W=M*H;\n",
+ "P=(W)/(T*550);\n",
+ "Pe=P1/E1;\n",
+ "I=(Pe*746)/V;\n",
+ "\n",
+ "#result\n",
+ "print \"power required is\",round(P,2), \"hp\"\n",
+ "print('Pick a 5HP motor')\n",
+ "print \"current required is\",round(I,2), \"amp\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 3.9,Page 145"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 10,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "power delivered to the load is 6.36 Watt\n",
+ "power provided by each supply is 7.23 Watt\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding power\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Vin=1.0;\n",
+ "Ri=1100.0;\n",
+ "Rf=10000.0;\n",
+ "Rl=8.0;\n",
+ "Vs=18.0;\n",
+ "\n",
+ "#calculation\n",
+ "Ir=Vin/Ri;\n",
+ "Vl=Ir*(Ri+Rf);\n",
+ "Ip=Vl/Rl;\n",
+ "Pl=(Vl*Ip)/2;\n",
+ "Ps=(Vs*Ip)/pi;\n",
+ "\n",
+ "#result\n",
+ "print \"power delivered to the load is\",round(Pl,2),\"Watt\"\n",
+ "print \"power provided by each supply is\",round(Ps,2), \"Watt\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 3.10,Page 149"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 11,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "power delivered is 141.67 Watt\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding power\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "V=170.0;\n",
+ "R=51.0;\n",
+ "\n",
+ "#calculation\n",
+ "I=V/R;\n",
+ "P=(V*I)/4;\n",
+ "\n",
+ "#result\n",
+ "print \"power delivered is\",round(P,2), \"Watt\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 3.11,Page 151"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 12,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "power dissipated is 7.05 watt\n",
+ "power dissipated when transistor resistance is 0.2 hm is 0.35 watt\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding power\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "V=7.2;\n",
+ "Rq=.2;\n",
+ "Rl=4;\n",
+ "D=.6;\n",
+ "\n",
+ "#calculation\n",
+ "Ip=V/(Rq+Rl);\n",
+ "Vl=Ip*Rl;\n",
+ "P=D*Vl*Ip;\n",
+ "Vq=Ip*Rq;\n",
+ "Pq=D*Vq*Ip;\n",
+ "\n",
+ "#result\n",
+ "print \"power dissipated is\",round(P,2), \"watt\"\n",
+ "print \"power dissipated when transistor resistance is 0.2 hm is\",round(Pq,2), \"watt\""
+ ]
+ }
+ ],
+ "metadata": {
+ "kernelspec": {
+ "display_name": "Python 2",
+ "language": "python",
+ "name": "python2"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter4.ipynb b/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter4.ipynb
new file mode 100644
index 00000000..a01cb417
--- /dev/null
+++ b/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter4.ipynb
@@ -0,0 +1,733 @@
+{
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 4: Linear Power Amplifier Integrated Circuits"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 4.1,Page 162"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 1,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "output voltage is 500.0 mV\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding voltage \n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Rf=1;\n",
+ "Ri=10;\n",
+ "Vi=0;\n",
+ "Ip=500;\n",
+ "\n",
+ "#calculation\n",
+ "Vrf=Ip*Rf;\n",
+ "\n",
+ "#result\n",
+ "print \"output voltage is\",round(Vrf,2), \"mV\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 4.2,Page 165"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 2,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "frequency of OPA548 is 67.26 KHz\n",
+ "slew rate of OPA548 is 1.12 Mhz\n",
+ "the OPA548 can be used\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding frequency\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Vi=300.0;\n",
+ "P=35.0;\n",
+ "R=8.0;\n",
+ "S=10000.0;\n",
+ "fh=20.0;\n",
+ "\n",
+ "#calculation\n",
+ "Vl=(P*R)**.5;\n",
+ "Vp=Vl*2**.5;\n",
+ "Il=Vl/R;\n",
+ "f=S/(2*pi*Vp);\n",
+ "Ao=Vl/Vi;\n",
+ "G=Ao*fh;\n",
+ "\n",
+ "#result\n",
+ "print \"frequency of OPA548 is\",round(f,2), \"KHz\"\n",
+ "print \"slew rate of OPA548 is\",round(G,2), \"Mhz\"\n",
+ "print('the OPA548 can be used')"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 4.3,Page 168"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 3,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "power delivered is 3.5 watt\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding power\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Rl=10.0;\n",
+ "V=12.0;\n",
+ "Vl=5.0;\n",
+ "\n",
+ "#calculation\n",
+ "Pl=Vl**2/Rl;\n",
+ "I=Vl/Rl;\n",
+ "Ps=V*I;\n",
+ "Pic=Ps-Pl;\n",
+ "\n",
+ "#result\n",
+ "print \"power delivered is\",round(Pic,2), \"watt\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 4.4,Page 170"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 1,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Vload Iload Pload Ps Pic\n",
+ "0.0 0.0 0.0 0.0 0.0\n",
+ "0.2 0.0 0.0 0.24 0.24\n",
+ "0.4 0.0 0.02 0.48 0.46\n",
+ "0.6 0.1 0.04 0.72 0.68\n",
+ "0.8 0.1 0.06 0.96 0.9\n",
+ "4.8 0.5 2.3 5.76 3.46\n",
+ "5.0 0.5 2.5 6.0 3.5\n",
+ "5.2 0.5 2.7 6.24 3.54\n",
+ "5.4 0.5 2.92 6.48 3.56\n",
+ "5.6 0.6 3.14 6.72 3.58\n",
+ "5.8 0.6 3.36 6.96 3.6\n",
+ "6.0 0.6 3.6 7.2 3.6\n",
+ "6.2 0.6 3.84 7.44 3.6\n",
+ "6.4 0.6 4.1 7.68 3.58\n",
+ "11.4 1.1 13.0 13.68 0.68\n",
+ "11.6 1.2 13.46 13.92 0.46\n",
+ "11.8 1.2 13.92 14.16 0.24\n",
+ "12.0 1.2 14.4 14.4 0.0\n"
+ ]
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAEZCAYAAACXRVJOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcU/W5x/EPi6yiQLGggoKoVL1WKBQBoQRBK4tQsVZF\nrlXqxVuwLrhVkTpVLFdxAVFkcekgICibQBFRJCDKKgNaBEVAZVMpoCDI4kzuH8+JE0JmJjOTk5Nz\n5vt+vfKaLCcnT2bg95zfDiIiIiIiIiIiIiIiIiIiIiIiIiIiIiKl8jnQ0YXzZgEvu3Dez3EnXpEy\nobzXAUhGizg3N87rhth4s3An6aRTHnBGzOOzgdeAncC3wBrgDhL/Pw45798H7AXWAze4F6oEgRKC\niD80BpYBXwD/BdQErgKaAzUKeM8257UTgHuBscA5rkd6rIoefKaIpNhm4GLnfmVgGFbIbAOeAio5\nr9UEZgPfALuBWcCpMedpBCzErlTnASMo+Op9HdA15nFF7Iq4qfO4O7AW2AMsAH6RIN7LgEPAYewK\nOcd5/UbgYyeOjUDfuM++B9gObAVu4ugr9MrA41iB/BXwHFAlQfyVsav382KeOwk4ANRxbrOd+HcB\ni4ByCX8TR3/+eOz3mqwQsCXuuW+AntjfraC/5ULnGICLnBi6OI87kv+7BOiD/T53A3OB0+Ji7wds\nwH7XIuJzsQnhIeB98gu195znAGoDV2AF5PHAq8D0mPMswQrT44B2WIE8roDPHIQVflFdsQQA1mTy\nPVYwVQDuxgqc6BVobLwPJviMLlhyAvgNsB9o5jy+DNiBXUFXdWKILZCfAmZgye94YCbwjwK+wwvA\n4JjH/YE5zv0hWDKp4NwuKuAcxH3+DuCPhRwbL0R+QiiP/X0OAWdR+N/y78DTzv37gc+A/3MeP4T9\nHgB6YL/7Js75BzrniY39Tez3VbkYcYtIhootYD/DCs2oS53XE2mKXTWCXTUewQrZqAkUXENojCWM\nKjHHPuDcHwRMijm2HHY1/5sE8WYV8hlR04FbnfsvAo/ExREtkMthiSi2Pb81sKmA83bEfl9R7wG9\nnft/xxJL4yJig6MTwmHsd56sEJBLfk1kFfAH57XC/pYdsb4JgDeAP2EJHaz28LuY1/rEnKM8lmAb\nxMQeKka8kgHUhyDJOgVrLon60nkOoBowGhvl8x1WcJyIFaSnYIXSDzHvjT1PvI1Ys1F357yXAxOd\n1052Pjcqgl0FxzZPFaYzsBQrIPdgNYafxZw7tolla8z9k5xYPnDetwcrEOsU8Dlh5/iWQEPgAvJr\nTEOxAnke9l3vTTL2XeT/vpO1HaiFfcdfYTU3KPxvuQSrif0cS+zjsEL+Z8CvsSYugNOB4eT/PnY5\nz8f+LeKbrCTDKSFIsrZjhVvUaVj7M8CdWCHSEksE7bFkUA5r6qiFFZBRp1P4SKNXgGuxZomPyb8S\n3+68N6ocVlht41jx568MTAUewwq7WlgzTrT9fgf5V7fE3f8PltDOdd5XC2sKOaGA+HOxwvda5zYL\nu3oGq2nchdUQugMDyK/VFOZt4MokjktGor/lduf+ASzx3Q58hNXu3sf+xp+RX/P7EuuDqRVzq44l\n3Ci3RpOJiAdim2Aexpo+ou3Oi8lvd34UK1wrY/0J07Emg+gFxxLsyvg4oC1WiyioDwGgHlYwLQT+\nEvN8tA/hYudcd2GFVKI+hJuBd8kv8GsAP2LNS+Ww2sL+mO9wGVYo/gJLXtkc3WQzDJiM1RbAroQL\na8JpiSWZj7BaTlRX4Ezyk9l2LIEmEvv5Z2BX4Y8BdZ3nzsSaxU5M8N4QBV+hF/a3BGs6+w7rFwDr\nHN6LDQaI+p3z3c51Hp+IjXpKFLuIBED8KKPhWAG2HSsgoyNTTsZG/OzDxrv3xa6SowmhEdbUsA9r\nKnmawhMC2BXxYexqPtbvsE7mb53PjB1GGRtvbSwh7AZWOs/1w0YI7XE+fyJHF4R/xQrxrcD/YoVa\ntAmkMlZQbsQKy4+BW4r4Dhuw2kXssMvbnTi/xwrsgQneF5XLsfMQXnXO+S2wGusDKWgewpcJnofC\n/5ZgiS4XGwAANsw1l6MLfLB+kQ+x38eXwPOFxC4C2EiKHAoeMvc09h9nDfkjPkS8dg5Wo1CzqkgK\nDcBGisxM8FoX8ofjXcjR7Y8i6XYFdvVcC/v3Os3bcESCpT5W9e9A4hrCKODqmMfryW8fFUm3N7Cm\nmF1YB7T+LUqZ4vaU8qewyUMFjcY4lWOH+tUHvnY5LpFEOnsdgIiX3Gwf7YZNlc+h4Kn5JHhNQ9VE\nRDzgZg2hDTbOugs26/QEbGTH9THHbOPo8d71STCmvHHjxpGNG7UciohIMW3EhidnlPYk7kOI7VRu\nRcGdypEge/DBB70OwVVB/n5B/m6RiL6f31HMFpd0LksbDexm5+doLBl0wSYX7cdWoxQREQ+kKyEs\ndG5giSBWUZN7REQkDTTpJgOEQiGvQ3BVkL9fkL8b6PuVNYWN/skkTnOYiIgkq1y5clCMcl41BBER\nAZQQRETEoYQgIiKAEoKIiDiUEEREBFBCEBERhxKCiIgASggiIuJQQhAREUAJQUREHEoIIiICKCGI\niIhDCUFERAAlBBERcSghiIgIkN4tNEV8JRKBL76AlSshJwe++y759/7859C8ObRoAXXruhejSCpp\ngxwRji78P/gg/1a5shXszZtDnTrJn2vbNjvXqlVQrZq9v3VrGDAAKlVy97uIRBV3gxy3E0IVbC/l\nykAl4HXgvrhjQs7zm5zHU4HBcccoIUjKRCLw+edHF/wffABVquQX/tHbySeX/rM2b7bzjx0L9epB\ndjaU88ulmPhapiUEgGrAAax5ajFwl/MzKgQMALoXcg4lBCmx3bvhnXfcL/yLcuAAhELQtSs8+KC7\nnyUCxU8I6ehDOOD8rARUAHYnOEbXS5JSkQgsXgxjxsCsWdC2Lfz613DbbVb416uX/piqVYOZM6FV\nK2jcGHr3Tn8MIoVJR0IoD6wCGgPPAR/HvR4B2gBrgG1YDSL+GJGk7N4N48ZZIsjLg5tvhmHD4Gc/\n8zoyU68e/Otf0KEDNGgA7dt7HZFIvnRemZ8IvAn8FQjHPF8DyMVqEp2B4cDZce9Vk5EUKL420LWr\nJYJ27TK3rf7tt+G662DRImjSxOtoJKgyscko6jvgX0ALjk4I+2LuvwGMBGoT17SUlZX10/1QKEQo\nFHInSvGN2NpAJAJ9+2ZWbaAwnTrBkCHQpQssXQonneR1RBIE4XCYcDhc4ve7ff1UB/gR+BaoitUQ\n/g7MjzmmLvAN1nTUEngVaBh3HtUQ5CiPPmoFardulggyuTZQmIEDYcECmD8fqlb1OhoJmkwbZXQ+\nkI31I5QHXgaGAjc7r48G+gN/xhLHAWzE0dK48yghyE9efx1uvRWWLIFTTvE6mtLJy7Omo9xcmDQJ\nymvtAEmhTEsIqaKEIABs2AAXXQSzZ0PLll5HkxoHD1oTUrt2VusRSZXiJgRdj4hv7N8PPXvC4MHB\nSQZgcyJmzIApU2zymohXVEMQX4hEbNx+5crwwgv+7C8oyoYNVksYNw4uvdTraCQIVEOQQHrmGVi3\nDp59NpjJAOCss+C11yzxffSR19FIWeSX/1qqIZRh770HV15pnciNGnkdjfsmToT77rPhqG4vpyHB\nlsnzEESK7auv4Oqr4aWXykYyAOjVCzZtgssvh4ULoXp1ryOSskI1BMlYR47Y6JuOHeFvf/M6mvSK\nRKBPH9i1C6ZPhwoVvI5I/EjDTiUw7rwT1q+35SjK4vj8w4ehc2c4/3ybgS1SXOpUlkB49VUbijl+\nfNlMBmAb6UydCvPmwYgRXkcjZYH6ECTjfPwx3HKLFYS1ankdjbdq1rTVUS+6CBo2tH4FEbeU0Wsv\nyVR799rks8cfh6ZNvY4mMzRqZP0IffrY5j4iblEfgmSMSAR+/3vblH7kSK+jyTzTpuWv4dSggdfR\niB9o2Kn41uOP2+b0Eyd6HUlm6tnT9mfu2tX2fzjhBK8jkqBRDUEywjvv2KqfK1ZA/fpeR5O5IhHo\n18/mKcyeDccd53VEksk07FR8Z+tWW6xu/Hi4+GKvo8l8P/4I3btbs9GoUcFdykNKT8NOxVcOHbJ+\ng9tvVzJIVsWKMHmyLW0xdKjX0UiQ+OXaQjWEgOrfH3bssPH2utItnq1boXVrePJJuOoqr6ORTKRO\nZfGNceNss/kVK5QMSqJ+fZg505bKrl/fkoNIafjlv6FqCAGzerUVZAsWwHnneR2Nv82ZA3/6k60K\ne8YZXkcjmUR9CJLx9uyx5axHjFAySIUuXWDQIPu5e7fX0YifuVlDqAIsBCoDlYDXgfsSHPc00Bk4\nANwA5CQ4RjWEgMjLsxEyZ59tbd+SOgMGwKpVtuRHpUpeRyOZIJNqCAeBDkBT4JfO/bZxx3QBzgTO\nAvoCz7kYj2SARx6x5SkefdTrSIJn6FCoXRtuusnmK4gUl9tNRgecn5WACkB8hbY7kO3cXwbUBOq6\nHJN4ZO5cGD3aVjLVhKrUq1DB5nKsXw8PPeR1NOJHbieE8sBq4GtgAfBx3OunAltiHm8FNE81gDZv\nhhtugEmToF49r6MJrmrVbOTRSy9ZchApDreHneZhTUYnAm8CISAcd0x8+1bCym5WVtZP90OhEKFQ\nKDURiut++MEmn91/P7SNbzSUlKtXz5bM7tDBZjO3b+91RJIu4XCYcDhc4venc9jpIOAH4PGY50Zh\nCWKS83g90B6rUcRSp7JPRSI2JPLgQZgwQfMN0untt219qEWLoEkTr6MRL2RSp3IdrE8AoCpwCceO\nIJoJXO/cbwV8y7HJQHzs+edt4tnYsUoG6dapEwwZYsNRd+70OhrxAzebjE7GOozLO7eXgfnAzc7r\no4E52Eijz4D9wI0uxiNptmIFDBxoSzVXr+51NGVTnz6wcSP06AHz50PVql5HJJnML9dsajLymf/8\nB1q0gOHDrTAS7+TlWdNRbq516pfVParLokxqMpIyKjcXrr0WevVSMsgE5cvbqKPt263GJlIQJQRJ\nub/9zX4+/LC3cUi+KlVgxgyYMsX6c0QS0WqnklKvv27j31eutIlSkjnq1LGF8Nq1g9NPt8UFRWKp\nD0FSZsMGuOgi29qxZUuvo5GCvPuuLS44fz6cf77X0Yib1Icgnti/3zaBHzxYySDTtWsHw4ZBt262\nOZFIlGoIUmqRCPTuDZUrwwsvaL6BXwwebP0KCxdqWHBQFbeG4Jf/ukoIGWzECBvF8t57GufuJ5GI\nzVPYtQumT1efTxApIUhavfeetUcvWQKNGnkdjRTX4cPQubP1JQwb5nU0kmrqQ5C0+eoruPpqqx0o\nGfhTpUowdaptqjNihNfRiNc07FRK5MgRSwZ9+9oVpvhXzZq2OupFF0HDhnD55V5HJF5Rk5GUyJ13\n2kYss2ZpKYSgWLbMRh7NnQvNm3sdjaSCmozEda++aqNTxo9XMgiSCy+0He169IAtW4o+XoJHTUZS\nLB9/DLfcYm3OtWp5HY2kWs+etrtd1662Su0JJ3gdkaSTmowkaXv32qSz+++H668v+njxp0gE+vWD\nTZts1rn2v/YvDTsVV0Qitg1m3bowcqTX0YjbfvwRune3LThHjdJkQ79SH4K44vHHYds2eOopryOR\ndKhYESZPto7moUO9jkbSRX0IUqR33oEnn7Qd0CpX9joaSZcaNazJqHVrm2dy1VVeRyRu80tFUE1G\nHtmyxfoNJkyAiy/2OhrxQk6OLZU9c6YlB/EPNRlJyhw6ZFeFd9yhZFCWNWsG2dk2AmnTJq+jETe5\nnRAaAAuAtcC/gVsTHBMCvgNynNsDLsckSbrjDjjlFLj7bq8jEa916QKDBtnP3bu9jkbc4nYfwhHg\nDmA1cDzwAfAWsC7uuIVAd5djkWLIzrYNVFas0AgTMf36wWefWU1h3jxbB0mCxe0awldYMgD4HksE\npyQ4TkVOBlm9Gu66C6ZN08QkOdrQoVC7Ntx0kw1FlmBJZx9CQ6AZsCzu+QjQBlgDzAHOTWNMEmfP\nHlvO+pln4LzzvI5GMk2FCrZkyfr18NBDXkcjqVZUk9FxwKXAb7ACPQJ8ASwC3gR+TPJzjgemALdh\nNYVYq7C+hgNAZ2AGcHb8CbKysn66HwqFCIVCSX60JCsvz3Y+69HDVjIVSaRaNRtx1KoVNG5s/2Yk\nM4TDYcLhcInfX1hTzSDgSmAJsBzYjtUoTgZaAq2wQn5wEZ9xHDAbeANIZguOzUBzILbrSsNO0+Ch\nh+Dtt63vQMsVSFHWroUOHeC116B9e6+jkURSuXRFd6wgzyvg9fJAN2BmEefPBnZhncuJ1AW+wWof\nLYFXsdpILCUEl73xBvzP/8DKlVCvntfRiF/Mnw+9esGiRdCkidfRSLxUJoRmWIdwaUritljz0ocx\n57kfOM25PxroD/wZa346AAwAlsadRwnBRZs3W/V/6lRo29braMRvXnwRHnkEli6Fk07yOhqJlcqE\n8AFwhvPzPee2BNhXivhKSgnBJT/8YDtl3XAD3JpolohIEgYOhAULrMZQtarX0UhUqlc7rY4147QB\nWjv3dwDvY1f16aKE4IJIBPr0sRnJEyZovoGUXF4eXHcd5ObCpEnaOClTuLX89fHAhVgT0PVY/0E6\nt1VXQnDBmDG2sfrSpVC9utfRiN8dPAidOkG7djBkiNfRCKQ2IVyH1QyaAoeAFVjb/hJswlk6KSGk\n2PLltn/u4sVw9jGDfEVK5j//sQXw7rnHBimIt1KZEL4HPgFGYR3Dn5QqstJRQkihnTuhRQt4+mmb\ncyCSShs2WC1h3DhbJVW8k8qEUBG4AOs7aAP8gvz+gyXAOyWOsviUEFIkNxd++1tb0vof//A6Ggmq\nd9+1Ge/z58P553sdTdnl5haadYE/ALdj8wQqFCewUlJCSJH777cF6+bOtWUIRNwycSLcd5/1UZ18\nstfRlE3FTQiFLV1xAVYziN4qYbWDp52f4jMzZthoopUrlQzEfb162f4Jl18OCxdq4IIfFJY5coDF\nWOH/PraGkVdUQyilTz+1SWezZ1tzkUg6RIc279oF06frQiTd3Gwy8pISQins3w8XXmgTz/r29Toa\nKWsOH4bOna0vYVgyq5lJyqRyC81/AVcB1RK8Vh24GluuWjJYJGLD/1q21DBA8UalSrYsyrx5Nu9F\nMldhfQg3ArcAfwdysRFG5YB6zvsmA390O0ApnREjbO36997TTGTxTs2aMGcOtGkDDRtav4JknmSL\niHrA6c79L9DENF9YvNiG/i1dCo3SOa9cpADLl0PXrjbKrXlzr6MJvlT3IVTE9kDuUIqYUkEJoZh2\n7IBf/xrGjrX2W5FMMW2a9WctWQINGngdTbClctgp2JLUeUBN4NuShyXpdOSI7XjWt6+SgWSenj1t\nyfWuXa0Wq327M0cymWMmtjfCW8B+57kIkM7FklVDKIYBA+CTT2DWLK06KZkpEoH+/WHjRhsKrR36\n3OHGsNMbnJ/RErmccz+7OIGVkhJCkiZPttnIK1dCrVpeRyNSsB9/hO7drdlo1CgNenCDW/MQqmG7\nnK0vQUypoISQhLVrIRSCt96Cpk29jkakaPv22UJ4vXrZCqmSWqmchxDVHZu1PNd53IzC91EWD+zd\na22zTzyhZCD+UaOGNRmNGAFTpngdjSSTOVYBFwMLsGQA8G/gv9wKKgHVEAoRidjw0nr1YORIr6MR\nKb6cHFsqe9Ys299bUsONGsIRjh1hlJfk+RtgiWQtlkQK6oh+GtgArCE/6UiShg6F7dvhqae8jkSk\nZJo1g+xsuOIKWxBPvJFMQliL7Z5WETgLGEHyq50eAe4AzgNaAf2Bc+KO6QKc6Zy7L/BckucW4J13\nLBFMmQKVK3sdjUjJdekCgwbZz927vY6mbEomIfwFK9APAa8Ae7E9EZLxFbDauf89sA44Je6Y7uSP\nWFqGzXmom+T5y7QtW2xj8wkToH59r6MRKb1+/Swh9Oxpi+JJeiWTEOoB9wMtnNtA4GAJPqsh1hy0\nLO75U4EtMY+3AireinDoEFx1FdxxB1x8sdfRiKTO0KFQuzbcdJP1j0n6FDVTGeAlrIBege2tvAj4\nqJifczwwBbgNqynEi+/0OOafQVZW1k/3Q6EQoVComCEEyx13wCmnwN13ex2JSGpVqADjx9sQ6oce\nggcf9Doi/wiHw4TD4RK/P9ne58pY7SAE3IwV8LWTfO9xwGzgDSDRauijgDAwyXm8HmgPfB1zjEYZ\nxcjOtv2QV6zQtH8Jrq++gtat4eGHoXdvr6Pxp1SvZQTQFviN87Mmtk/ComTjAV4APiZxMgCb03AL\nlhBaYSOavi7g2DJv9Wq46y4Ih5UMJNjq1bM5Ch062Gzm9u29jij4kskcucAHwBBsQ5xDxTh/Wyx5\nfEh+M9D92KxngNHOz2eAy7C1km7E5j7EUg0B2LMHWrSw2sHVV3sdjUh6zJ9vM5kXLYImTbyOxl/c\nWLqiJlawtwNaYgliKfBACeIrqTKfEPLybFORJk3gySe9jkYkvV58ER55xPb2OOkkr6PxDzeajL4F\nNmEdyw2ANkClkgQnJTd4sK378uijXkcikn59+tjKqD16WI2halWvIwqmZDLHJuAT4F2s+WcFxWs2\nSoUyXUN44w3bD3nlSmtXFSmL8vJs3k1uLkyapKXdk+FGk1EFrJnIS2U2IWzebGu7TJ0Kbdt6HY2I\ntw4ehE6dbIXUIUO8jibzubGW0SnAdGCnc5uKJo6lxQ8/2KJ1AwcqGYgAVKkCM2bYUi1jx3odTfAk\nkzneBiYA453H1zm3S9wKKoEyV0OIRKzd9NAhW5pCm4eI5NuwwWoJ48bZKqmSmBtNRmuAC5J4zk1l\nLiGMGWNrxC9dCtWrex2NSOZ5912rQc+fD+ef73U0mcmNJqNdwH9jfQkVgd7Af0oSnCRn+XJ44AHr\nN1AyEEmsXTsYNgy6dYMdO7yOJhiSSQh9gD9gK5fuAK7CJo+JC3butEXrxo6Fs8/2OhqRzNarl43A\nu/xy2L/f62j8r7CqRDngd9heBR+Rv4WmF8pEk1FuLvz2t9Cypc1GFpGiRfvbdu2C6dNtcTwxqWwy\nGonte1AbeAj4W6kikyINGmSdxw8/7HUkIv5RrhyMHm01hDvv9Doafyssc6wFfonNQagGLAZ+lY6g\nEgh8DWHGDLjtNpt8pqn5IsX37bfQpg38+c/wl794HU1mSOXSFYfJn5B2oDgnleL59FPo29dWdlQy\nECmZmjVhzhxLCg0bWr+CFE9hhfwPwGcxjxsDG537Eaz2kC6BrSHs3w8XXgi33mpJQURKZ/ly6NoV\n5s6F5s29jsZbqZyH0LCI936e7IekQCATQiRia7NUqQIvvKDJZyKpMm2aXWQtWWJ7KZRVqWwy+ry0\nwUjhRoyA9evhvfeUDERSqWdPWwesa1dYvFibSSXLL8VQ4GoIixfbLMulS6FRI6+jEQmeSAT697dl\ns2fPhuOO8zqi9HNjprKk2I4dcM018M9/KhmIuKVcOXj6aZuXcMstliCkcIUlhJ8D5yV4/jxAY2FK\n6MgR2/6yb1/o3NnraESCrWJFmDwZli2DoUO9jibzFZYQRgB1Ejz/M2C4O+EE3733Qo0atlaRiLiv\nRg1rMhoxwpbNloIVlhDOBBYmeH4Rya90+iLwNbb0RSIh4Dsgx7kFupicPBlefx3Gj9duTyLpVL8+\nzJxpk9aWLvU6msxVWLFUo5DXku2eeQm4rIhjFgLNnNvgJM/rO2vXWjvm1KlQq5bX0YiUPc2aQXY2\nXHEFbNrkdTSZqbCE8BnQNcHzXcifoFaUd4E9RRzjl5FOJbZ3rw2De+IJaNrU62hEyq4uXWzNsC5d\nYPdur6PJPIUVxmcDs4H3gQ+cY5sDbYBuwCdJfkZDYBaQaAuL9sA0YCuwDbgL+DjBcb4ddhqJ2PDS\nevVg5EivoxERsEXwPvgA5s2DSpW8jsY9qZyY9im2PEUv4L+w5SoWAv+LLWuRCquABthaSZ2BGVgi\nOkZWVtZP90OhEKFQKEUhuGvoUNi+HV55xetIRCTqscds35GbbrJmpKBMDA2Hw4TD4RK/Px2/hoYU\nXEOItxmrhcRX5nxZQ3jnHVuaYsUK69QSkcxx4ACEQjab+cEHvY7GHamsIXyP1QoSiQCpmAxeF/jG\nOV9LLPBAtOxt2WLJYMIEJQORTFStmo08at0aGjeG3r29jsh7hSWE41Nw/lewfoI6wBbgQfJHKI0G\nfg/8GfgRaza6JgWf6blDh6w6escdcPHFXkcjIgWpV8/mKHToYIvgtW/vdUTe8kvLma+ajPr1g6++\nsiGmQWmbFAmy+fNtf+ZFi6BJE6+jSR2tZeSx7Gz7x/XPfyoZiPhFx44wZIgNR9250+tovOOXIssX\nNYTVq+GSSyAchvMSrQIlIhlt4EBYsMAu6qpW9Tqa0kvlBjmZJOMTwp490KIF/OMftnidiPhPXp4N\nBsnNhUmT/L/EjBKCB/LybP/WJk3gySe9jkZESuPgQejUCdq1s2YkP1MfggcGD4Z9++DRR72ORERK\nq0oVmDHDVkYdO9braNKrsGGnkoQ33oAxY2DlyrK5I5NIENWpA3PmWC3h9NPh0ku9jig91GRUCps3\nQ6tWNry0bVuvoxGRVFu82BamnD8fzk9mrYUMoyajNPnhB1u0buBAJQORoGrbFoYPh27dbOvboFMN\noQQiEejTx2YkT5ig+QYiQTd4sPUrLFwI1at7HU3yNMooDcaMse34li711z8OESmZ6EXgrl0wfTpU\nqOB1RMlRQnDZ8uVWfVy8GM5OuFC3iATR4cPQubP1JQwb5nU0yVEfgot27rRF68aOVTIQKWsqVbIB\nJPPmWQtBEGnYaZJyc+Haa20WY48eXkcjIl6oWdOGo7ZpAw0b2oTUIFGTUZLuv982upk71z/thyLi\njuXLbWOduXOheXOvoymYmoxcMGOGjSaaOFHJQESgZUsYPdpaC7Zs8Tqa1FGTURE+/RT69rVNNE46\nyetoRCRT9Oxpk1O7drVBJiekYg9Jj6nJqBD798OFF8Ktt1pSEBGJFYlA//6wcaNdNGba8jUadpqy\nD7QO5CpeQGJHAAAORklEQVRV4IUXNPlMRBL78Ufo3t224Bw1KrPKCvUhpMiIEbB+PTz7bGb9gUUk\ns1SsCJMnw7JlMHSo19GUjtsJ4UXga+CjQo55GtgArAGauRxPUhYvhkcesTHHQdg1SUTcVaOGNRmN\nGGHLZvuV2wnhJeCyQl7vApwJnAX0BZ5zOZ4i7dgB11xjeyI3auR1NCLiF/Xrw6xZ0K+fLWvjR24n\nhHeBPYW83h3Idu4vA2oCdV2OqUBHjtj2l3372hR1EZHiaNrULiavuAI2bfI6muLzug/hVCB2FO9W\noL5HsXDvvVb1e+ABryIQEb/r0gUGDbKfu3d7HU3xZMI8hPgu24TDibKysn66HwqFCIVCKQ1i8mR4\n/XXb+czvG2uLiLf69bOhqD172tpHlSql53PD4TDhcLjE70/H+JmGwCwg0X5Do4AwMMl5vB5oj3VE\nx3J12OnatRAKwVtvWZVPRKS0cnNtMczjj4fsbG9GK/pt2OlM4HrnfivgW45NBq7au9ey+BNPKBmI\nSOpUqADjx9vw9Yce8jqa5LjdZPQKdsVfB+sreBCIzuUbDczBRhp9BuwHbnQ5nqNEInDDDdCxI1x/\nfZGHi4gUS7VqMHMmtG4NjRtD795eR1Q4v0y5cqXJ6LHHYNo02xavcuWUn15EBLBm6Q4d4LXXoH37\n9H2ulq5I0ocfwqWX2pLWDRqk9NQiIseYPx969YJFi6BJk/R8pt/6EDwzdCgMGKBkICLp0bEjDBli\nw1F37vQ6msTKZA1hyxbrQN640XZAEhFJl4EDYcECqzG4vTSOmoyScNdd1qH8xBMpO6WISFLy8mwl\n5dxcmDTJ3XlPSghF+O47OOMMyMmB005LySlFRIrl4EHo1AnatbNmJLeoD6EIY8fCZZcpGYiId6pU\nsa15p0yxMilTlKkawuHDNhZ45kxolhELbYtIWbZhg9USxo2zUY+pphpCISZPtuFeSgYikgnOOstq\nCb17w0eF7RqTJmUmIUQiNtT07ru9jkREJF/btjB8OHTrZvuxeCkTVjtNi7fesqTgRrVMRKQ0rr3W\nhsFffrmtnFC9ujdxlJk+hEsusWrZH/+YoohERFIoEoE+fWDXLpg+3RbHKy0NO01g9Wqrjm3alL51\nyUVEiuvwYdut8fzzYdiw0p9PncoJPPEE3HqrkoGIZLZKlWDqVNtU5/nn0//5ga8haJkKEfGbnBzr\nTyhtq4ZqCHGGD7c9D5QMRMQvmjWDc86xpS3SKdA1BC1TISJ+9eabNkx+zZqSb7+pGkKMMWO0TIWI\n+FN0iPy8een7zMDWEA4fttrBrFmamSwi/jRuHLz8ss2jKolMrCFcBqwHNgD3Jng9BHwH5Di3B1Lx\noZMnwy9+oWQgIv51zTWwbp01e6eD2zWECsAnQCdgG7ACuBZYF3NMCBgAdC/kPMWqIUQicMEFtmfy\nZZcVN2QRkcwxdKj1I4wfX/z3ZloNoSXwGfA5cASYBPRIcFxKE1N0mYrf/jaVZxURSb++fWHOHPjy\nS/c/y+2EcCqwJebxVue5WBGgDbAGmAOcW9oPffxx2xWtpD3zIiKZ4sQT4cYbbQi929xOCMm086wC\nGgAXACOAGaX5wNWrYe1aWyxKRCQIbr8dXnoJvv3W3c9xe7XTbVhhH9UAqyXE2hdz/w1gJFAb2B17\nUFZW1k/3Q6EQoVAo4QdqmQoRCZoGDaBLFxtKf889BR8XDocJh8Ml/hy3G1UqYp3KHYHtwHKO7VSu\nC3yD1SZaAq8CDePOk1Sn8pYt1pm8aZNmJotIsKxeDV27wubNyV/wZlqn8o/ALcCbwMfAZCwZ3Ozc\nAH4PfASsBoYB15T0w4YPt+WtlQxEJGiaNoVzz3V3OQu/dLsWWUOILlOxahWcfnqaohIRSaPiLmeR\naTWEtBk71uYcKBmISFC5vZxFIGoIhw9D48bw+uvwq1+lMSoRkTQbN85ub79d9LFlsoYweTKcfbaS\ngYgE3zXXwPr17ixn4fuEEInY1O677/Y6EhER91WqZEPrn3gi9ed2ex6C67RMhYiUNTffDI0a2XIW\nqVze3/c1hKFDtUyFiJQtbi1n4ZdiNGGnckkmaoiIBMGXX9rchMIm4papTuXHH9cyFSJSNp12Wv5y\nFqni2xqClqkQkbIuJwe6dSu4laTM1BCGD4cbblAyEJGyq1kzOOec1C1n4csawr59NiM5J0czk0Wk\nbJs7F/76VysP4wfXlIkawvjxEAopGYiIXHop7N8P779f+nP5LiFEIjByJPTv73UkIiLeK18e+vWz\ncrG0fNdktGiR7TG6bp3mHoiIAOzZY6s9r18PdevmPx/4JqNnn7VsqGQgImJq1YIrr4Tnny/defxS\nrEYikQg7dtgGEZ9/bjP1RETE5ORAjx42FL+isyhRoGsIY8fCH/6gZCAiEq9ZM6hfH2bPLvk5fFND\nOHw4QqNGMGcO/PKXXocjIpJ5xo+H7Gxb9BMCXEOYOdNW91MyEBFJ7Kqr4MMP4ZNPSvZ+txPCZcB6\nYANwbwHHPO28vgZoVtCJop3JIiKSWOXK8Kc/wXPPlez9biaECsAzWFI4F7gWOCfumC7AmcBZQF+g\nwK+xbp31ogdROBz2OgRXBfn7Bfm7gb6fH918M7z8sk1WKy43E0JL4DPgc+AIMAnoEXdMdyDbub8M\nqAnUJYGbbgruqqZB/EcZK8jfL8jfDfT9/Oj006FtW5g4sfjvdTMhnApsiXm81XmuqGPqJzpZ374p\njU1EJLD697dm9uJyMyEcu6NNYvE94Anf16BB6YIRESkrOnWCAweK/z43h522ArKwPgSA+4A84NGY\nY0YBYaw5CawDuj3wddy5PgMauxSniEhQbcT6aT1XEQumIVAJWE3iTuU5zv1WwNJ0BSciIunVGfgE\nu8K/z3nuZucW9Yzz+hrgV2mNTkRERERE/CeZyW1+1QBYAKwF/g3c6m04rqgA5ACzvA7EBTWBKcA6\n4GOs2TNI7sP+bX4ETAQqextOqb2I9U9+FPNcbeAt4FNgHvY39atE328o9u9zDTAN8PVKcBWw5qSG\nwHEk7ofws3pAU+f+8VjzWpC+H8AAYAIw0+tAXJAN9HHuV8Tn/9niNAQ2kZ8EJgN/9Cya1GiHrYYQ\nW2A+Btzj3L8X+L90B5VCib7fJeSPJv0//P39aA3MjXn8V+cWVDOAjl4HkUL1gbeBDgSvhnAiVmAG\nVW3sAqUWluxmAZ08jSg1GnJ0gbme/Mmw9ZzHftaQo79frCuA8YW9OdMXt0tmcltQNMSy+zKP40il\np4C7seHGQdMI2Am8BKwCxgLVPI0otXYDTwBfAtuBb7HkHjR1yR/m/jUFrJQQEH3IH9WZUKYnhGQn\nt/nd8Vhb9G3A9x7HkirdgG+w/gO/LLNeHBWxUXEjnZ/7CVbttTFwO3ahcgr2b/Q6LwNKgwjBLXMG\nAoexvqACZXpC2IZ1vEY1wGoJQXIcMBWrys3wOJZUaoOtVbUZeAW4GBjnaUSptdW5rXAeTyFYw6Zb\nAO8Du4AfsQ7JNp5G5I6vsaYigJOxi5iguQGb8+X7hJ7M5DY/K4cVkk95HYjL2hO8PgSARcDZzv0s\njp6F73cXYCPfqmL/TrOB/p5GlBoNObZTOTp68a/4vNOVY7/fZdhIsTqeROOCRJPbgqIt1r6+Gmta\nySF/qY8gaU8wRxldgNUQAjGkL4F7yB92mo3VZv3sFaw/5DDWN3kj1nn+NsEYdhr//fpgw/W/IL98\nGelZdCIiIiIiIiIiIiIiIiIiIiIiIiIiIiLplaolPbKAO1Nwnn8CVzr3b8cmcLnt58C/nM/aBdSI\ne30G8AdshvigNMQjAZPpS1eIRKVqjZlUnid6rttIz8J2t2CJ6AdsFeArYl47EbgImwA4G0tWfp9I\nJmmmhCB+Uw7b9OMj4EPsihhs8bW3gQ+c57vHvGcgNtv9XaBJgnOeCHwe87g6tspnBWy/iqXkz0aO\nnclaDvgLtvjbAmC+8/xz2Azmf2M1kqgu2GYlK4GnyV/Oozq2uckybOXU2Nhj/R6rIYDNSr0m5rUr\nsCRxEJv9vgS4tIDziIj42j7n55XYEgPlsCaUL7DFySqQ34RSB5uyD9AcSxBVnNc3YJv2xJsBhJz7\nVwNjnPsfYhuPAPyd/HWnXgJ6Ovc3Y0sgRNVyflbAEsX5zud/CZzuvDaR/OU8/kH+wmM1seQVX+Oo\nx9Fr1FQCvor5rLlYwom6kWCtrSRpoBqC+E1brDCNYCtTLgR+jSWIIdiV/FvYVXtdrDCfhl0578MK\n4UTLcU/GEgHYlfdkrOZwIlazAFvP5zdJxHg1VlNZBZwHnAv8AttQ5wvnmFdi4rgUW1gtB0sglTl6\nlV+wRLIj5vFh57tchSXApsCbMa9vxxY6E0laRa8DECmmCIkL9N5YwfgrIBe7aq+S4PiC9maYhV2p\n13LO8Q5wQtwxyezr0AjrtG4BfIfVJKJxFHaunuTXagoS/55XsM7jclgNJzfmtfIJPlOkUKohiN+8\ni12BlwdOwq7Yl2GF9zdYodgBu6KOYEtU/478JqNuJC4ov8fa/aNt+xGsQN+D1UoA/hsIJ3jvPvKT\nxwnYZjl7sRpKZ+dcnwBnkN9kdHVMHG8Ct8acr1mCz4g2jcUKY8tv98eSQ6yTya+NiCRFNQTxi2jh\nOR3ba3uN89zdWCKYgBXkH2Kdtuuc43Ow5p81znHLC/mMycCr5PclgG0sPwpr09+Itc3HG4O14W/D\n9sTOwfbm3QIsdo45CPRzjtuPJZ/od3oYGObEXh5rWorvWP4K+/9a3Xk/zvtfw5qNFsYd35Jg7kEh\nIhII1WPuP4sNVy2OLPL7OQpTHttjQxd8IiIZ6nas9rAWeBlrxiqOkyhik3RHd+CBYp5bRERERERE\nREREREREREREREREREREJBn/DwWIrS4GjkBZAAAAAElFTkSuQmCC\n",
+ "text/plain": [
+ "<matplotlib.figure.Figure at 0x7f414b298f90>"
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "#finding Vload vs Pic graph\n",
+ "\n",
+ "#initialisation of variable\n",
+ "%matplotlib inline\n",
+ "import numpy as np\n",
+ "import matplotlib.pyplot as plt\n",
+ "Vload=[0.0, 0.2, 0.4, 0.6, 0.8, 4.8, 5.0, 5.2, 5.4, 5.6, 5.8, 6.0, 6.2, 6.4, 11.4, 11.6, 11.8, 12.0];\n",
+ "Iload=[0.0, 0.0, 0.0, 0.1, 0.1, 0.5, 0.5, 0.5, 0.5, 0.6, 0.6, 0.6, 0.6, 0.6, 1.1, 1.2, 1.2, 1.2];\n",
+ "Pload=[0.00, 0.00, 0.02, 0.04, 0.06, 2.30, 2.50, 2.70, 2.92, 3.14, 3.36, 3.60, 3.84, 4.10, 13.00, 13.46, 13.92, 14.40];\n",
+ "Ps=[0.00, 0.24, 0.48, 0.72, 0.96, 5.76, 6.00, 6.24, 6.48, 6.72, 6.96, 7.20, 7.44, 7.68, 13.68, 13.92, 14.16, 14.40];\n",
+ "Pic=[0.00, 0.24, 0.46, 0.68, 0.90, 3.46, 3.50, 3.54, 3.56, 3.58, 3.60, 3.60, 3.60, 3.58, 0.68, 0.46, 0.24, 0.00];\n",
+ "\n",
+ "#result\n",
+ "print('Vload Iload Pload Ps Pic');\n",
+ "for i in range(0,18):\n",
+ " print Vload[i],\" \",Iload[i],\" \",Pload[i],\" \",Ps[i],\" \", Pic[i]\n",
+ " \n",
+ "plt.plot(Vload,Pic);\n",
+ "plt.xlabel('load voltage (V)')\n",
+ "plt.ylabel('IC Power(W)')\n",
+ "plt.title('load voltage vs IC Power')\n",
+ "plt.show()"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 4.5,Page 173"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 4,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "IC power is 2.57 W\n",
+ "total power is 3.82 W\n",
+ "dc supply current is 159.155 mA\n",
+ "power delivered is 1.25 watt\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding different power and current\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "V=12.0;\n",
+ "Vp=5.0;\n",
+ "R=10.0;\n",
+ "\n",
+ "#calculation\n",
+ "Ip=Vp/R;\n",
+ "Il=Ip/2**.5;\n",
+ "Pl=(Vp*Ip)/2;\n",
+ "Id=Ip/pi;\n",
+ "Pt=2*V*Ip/pi;\n",
+ "Pic=Pt-Pl;\n",
+ "\n",
+ "#result\n",
+ "print \"IC power is\",round(Pic,2), \"W\"\n",
+ "print \"total power is\",round(Pt,2), \"W\"\n",
+ "print \"dc supply current is\",round(Id*1000,3), \"mA\"\n",
+ "print \"power delivered is\",round(Pl,2), \"watt\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 4.6,Page 179"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 5,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "thermal resistance is 24.61 C/W\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding thermal resistance\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Ts=40.0;\n",
+ "P=2.92;\n",
+ "Qj=2.5;\n",
+ "Qc=2.0;\n",
+ "Tj=125.0;\n",
+ "\n",
+ "#calculation\n",
+ "Qs=(Tj-Ts)/P-Qj-Qc;\n",
+ "\n",
+ "#result\n",
+ "print \"thermal resistance is\",round(Qs,2),\"C/W\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 4.7,Page 180"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 48,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Vload Iload Pload Ps Pic\n",
+ "9.4 0.94 4.42 14.36 9.94\n",
+ "9.6 0.96 4.61 14.67 10.06\n",
+ "10.0 power delivered by IC in watt\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding power\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "V=24.0;\n",
+ "R=10.0;\n",
+ "Qs=4.0;\n",
+ "Tj=125.0;\n",
+ "Ta=40.0;\n",
+ "Qj=2.5;\n",
+ "Qc=2.0;\n",
+ "Vload=[9.4, 9.6];\n",
+ "Iload=[.94, .96];\n",
+ "Pload=[4.42, 4.61];\n",
+ "Ps=[14.36, 14.67];\n",
+ "Pic=[9.94, 10.06];\n",
+ "\n",
+ "#calculation\n",
+ "P=(Tj-Ta)/(Qj+Qc+Qs);\n",
+ "\n",
+ "#result\n",
+ "print('Vload Iload Pload Ps Pic');\n",
+ "for i in range(0,2):\n",
+ " print Vload[i],\" \",Iload[i],\" \",Pload[i],\" \",Ps[i],\" \", Pic[i]\n",
+ "print round(P,2),\"power delivered by IC in watt\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 4.8,Page 182"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 7,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "gain is 23.0\n",
+ "limit current is 4.01 A\n",
+ "output voltage is 46.0 V\n",
+ "maximum output voltage is 32.0 V\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding current and voltage\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Rf=22.0;\n",
+ "Ri=1.0;\n",
+ "Rs=15.0;\n",
+ "I=4.75;\n",
+ "Rc=4.0;\n",
+ "Vp=2.0;\n",
+ "Rl=8.0;\n",
+ "Im=4.0;\n",
+ "\n",
+ "#calculation\n",
+ "Av=1+(Rf/Ri);\n",
+ "Il=(Rs*I)/(Rc+13.75);\n",
+ "Vo=Vp*Av;\n",
+ "V=Im*Rl;\n",
+ "\n",
+ "#result\n",
+ "print \"gain is\",round(Av,2)\n",
+ "print \"limit current is\",round(Il,2), \"A\"\n",
+ "print \"output voltage is\",round(Vo,2), \"V\"\n",
+ "print \"maximum output voltage is\",round(V,2), \"V\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 4.9,Page 185"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 8,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "loudness ofsound is 108.06 dB\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding loudness\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan,log\n",
+ "D=8.0;\n",
+ "d=1.0;\n",
+ "I=90.0;\n",
+ "\n",
+ "#calculation\n",
+ "Is=20*log(d/D,10);\n",
+ "Ir=I-Is;\n",
+ "\n",
+ "#result\n",
+ "print \"loudness ofsound is\",round(Ir,2), \"dB\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 4.10,Page 186"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 14,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "19.95 power provided in watt\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding power\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan,log\n",
+ "D=1.0;\n",
+ "I1=108.0;\n",
+ "I2=95.0;\n",
+ "P=1.0;\n",
+ "\n",
+ "#calculation\n",
+ "I=I1-I2;\n",
+ "Pr=P*10**(I/10);\n",
+ "\n",
+ "#result\n",
+ "print \"power provided is\",round(Pr,2), \"watt\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 4.11,Page 188"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 9,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "output voltage is 12.65 V\n",
+ "gain is 10.28\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding output voltage and gain\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan,log\n",
+ "P=20;\n",
+ "R=8;\n",
+ "Vi=1.23;\n",
+ "\n",
+ "#calculation\n",
+ "V=(P*R)**.5;\n",
+ "G=V/Vi;\n",
+ "\n",
+ "#result\n",
+ "print \"output voltage is\",round(V,2), \"V\"\n",
+ "print \"gain is\",round(G,2)\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 4.12,Page 191"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 10,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "resistor b/w pins 1&8 is 600.0 ohm\n",
+ "thus pick a 620 ohm resistor\n",
+ "capacitor b/w pins 1&8 is 22.46 microF\n",
+ "thus pick a 27 microF capacitor\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding resistor and capacitor\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan,log\n",
+ "G=40.0;\n",
+ "f=80.0;\n",
+ "R1=15000.0;\n",
+ "R2=150.0;\n",
+ "\n",
+ "#calculation\n",
+ "R=2*(R1/G)-R2;\n",
+ "R11=620;\n",
+ "C=1/(2*pi*f*R11/7);\n",
+ "\n",
+ "#result\n",
+ "print \"resistor b/w pins 1&8 is\",round(R,2),\"ohm\"\n",
+ "print('thus pick a 620 ohm resistor')\n",
+ "print \"capacitor b/w pins 1&8 is\",round(C*1e6,2), \"microF\"\n",
+ "print('thus pick a 27 microF capacitor')"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 4.13,Page 193"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 11,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "power dissipated is 140.0 mW\n",
+ "thermal resistance is 628.93 degree C/W\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding thermal resistance and power\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan,log\n",
+ "R=8.0#resistance\n",
+ "V=5.0#voltage\n",
+ "Tm=150.0#temperature\n",
+ "Ta=50.0#temperature\n",
+ "Qa=107.0;\n",
+ "Qc=37.0;\n",
+ "Ps=299.0;\n",
+ "\n",
+ "#calculation\n",
+ "Vd=V/2;\n",
+ "Vm=V-1;\n",
+ "Vp=Vm-Vd;\n",
+ "Vr=Vp/2**.5;\n",
+ "Pl=1000*Vr**2/R;\n",
+ "Pl=140;\n",
+ "Pic=Ps-Pl;\n",
+ "Q=(Tm-Ta)/Pic;\n",
+ "\n",
+ "#result\n",
+ "print \"power dissipated is\",round(Pl,2), \"mW\"\n",
+ "print \"thermal resistance is\",round(Q*1000,2),\"degree C/W\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 4.14,Page 197"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 12,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "power deliverd is 562.5 mwatt\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding power delivered\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan,log\n",
+ "R=8.0#resistance\n",
+ "V=5.0#voltage\n",
+ "\n",
+ "#calculation\n",
+ "Vl=V-1;\n",
+ "Vp=Vl-1;\n",
+ "Vr=Vp/2**.5;\n",
+ "P=Vr**2/R;\n",
+ "\n",
+ "#result\n",
+ "print \"power deliverd is\",round(P*1000,2), \"mwatt\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 4.15,Page 201"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 13,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "power delivered is 85.56 watt\n",
+ "thermal resistance is 1.4 degreeC/W\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding thermal resistance and power\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan,log\n",
+ "R=8.0#resistance\n",
+ "Ts=35.0#temperature\n",
+ "Ta=150.0#temperature\n",
+ "Vm=42.0#voltage\n",
+ "\n",
+ "#calcuation\n",
+ "Vp=Vm-5;\n",
+ "Vr=Vp/2**.5;\n",
+ "Pm=Vr**2/R;\n",
+ "P=45;\n",
+ "Qs=(Ta-Ts)/P-1.2;\n",
+ "\n",
+ "#result\n",
+ "print \"power delivered is\",round(Pm,2), \"watt\"\n",
+ "print \"thermal resistance is\",round(round(Qs*10)/10,2), \"degreeC/W\""
+ ]
+ }
+ ],
+ "metadata": {
+ "kernelspec": {
+ "display_name": "Python 2",
+ "language": "python",
+ "name": "python2"
+ },
+ "language_info": {
+ "codemirror_mode": {
+ "name": "ipython",
+ "version": 2
+ },
+ "file_extension": ".py",
+ "mimetype": "text/x-python",
+ "name": "python",
+ "nbconvert_exporter": "python",
+ "pygments_lexer": "ipython2",
+ "version": "2.7.6"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter5.ipynb b/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter5.ipynb
new file mode 100644
index 00000000..58305b46
--- /dev/null
+++ b/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter5.ipynb
@@ -0,0 +1,647 @@
+{
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 5: Discrete Linear Power Amplifier"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 5.1,Page 215"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 1,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Id=0 from 0 to 2 so not shown in the graph\n"
+ ]
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEZCAYAAACTsIJzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4U2Xax/FvW4ogi7K8uIGgiAOIIqKC4EBxFxcUURRE\nQAUGcccNx4UZfd/RUWdwRFFBEBQQFWVEQFywbiOiwyKyuLAJio6gIoso0Lx/3KfTENP2pM3Jc5L8\nPteVi+TkNLmtp7nzPPezgIiIiIiIiIiIiIiIiIiIiIiIiIiIiIgkUROgCMh1HIdIKOkPQzLVauBE\nR+99IbAqzvEqwH+ArqkNRyQxSgySqSLezYUXgb2BzjHHTwN2Aa+kPCKRBCgxSDbIA+4HvgNWAGeU\nce7NwHMxxx70bgD9vNf4CVgJ9IrzGr8AzwKXxBy/BJiEdWMB3AR8DawDLveOH+w91xVY4r3POmBo\nGTGLiIgPq4ATvPt/AJYBBwB1gDexb+7xvhgdCGwFanqP87AP72OBGsAmoJn33D5Ay1Lev4N3bjXv\n8V7ANuAI7/FpwHqgBVAdeJrdE8N6oGPUz7Yp+z9XRETKE50Y5gADo547mbKLz+8AfaLO/cK7XwP4\nAeiOfZiX5zPgIu/+AGBB1HNjgf+NetyU3RPDGi/m2j7eRySp1JUk2WA/YG3U4y/LOX8SJR/ovYCJ\n3v2tQE+sBfI18DLwuzJeZwIl3Ul9vMelxbQu5mfPw7qTVgOFQPtyYhYRkXLEthgGRT1XXovhf7Bu\nnwOwFkK8D/89sLrF22XE0AT4FTgOqzs0iHoutsVwCLu3GIrlAddSfjITSRq1GCQbPAtcTUmN4ZZy\nzv8O+5b+JFZg/tQ73gDohnUp7cBaELvKeJ3VwLvAZOBVbKhqdEz9gebAnsDtUc/lA72x2sIuYHM5\n7yMiIj5EtxjygL8BG7ARRVdQevG52MXYN/jo0UD7YgnjR6wlMQf7YC9LX++9zo/z3C1YkXkd1j1V\nhCWvfGAW8D1WwP4AK2aLpL1G2OiPJcAn2De2eP4BfA4sQiMvJHu1AHaiVrxkuH2BI737NbHmeIuY\nc7oCM7377YC5qQlNJBTOxWoVdYCXgBfchiOSetP47RIFj2KjPIotx8aGi2SDWVi31EZgKrr2Jcs0\nwcZl14w5Pp3d+05fB9qmKCYREYkjFf2ZNYHngWuALXGez4l57Gp9GxERwVZ7DFI+1kR+GutKivUV\nVqQu1tA7tpv69ZtGNmxYEUiAIiIZbAU2RyYhQbYYcoAngKXAiFLOeYmSmaHtsf7Wb2NP2rBhBZFI\nhKVLI9x0U4R9941w3HERHn88wqZNESKR1N/uvPNOJ++rmDI3LsWkmJJ9w5ZaSViQiaEjNha8C7ZG\nzALgdGwGavEs1JnYBKIvgMew8eWlatEC7r0X1q6FW2+FWbPgwAOhb18oLISiorJ+WkRE/AiyK+ld\n/CWeKxN94SpV4Mwz7faf/8DEiXDVVbBtG/TrZ4niwAMTjldERMiAyTQNGsB118HHH8OUKbB+PbRp\nA6eeCs88A9u3B/O+BQUFwbxwJSgm/8IYl2LyRzEFL3ZEUFhFvP4yX37+GaZNg3HjYP586NkTLr0U\njjoKctLlv1hEpJJy7AMv4U+9dPmYTCgxRPvySxg/3pJErVpw2WVw8cVQt26SIxQRCRklhnIUFcFb\nb8Ho0TBzptUnLr8cOndWK0JEMpMSQwI2boSnn7Yk8euvliD69oV9tCCBiGQQJYYKvSjMnWsJ4sUX\n4cQTYcAAOOkkyMtL+tuJiKSUEkMl/fQTTJpkSWLjRitWX3opNGwY6NuKiARGiSGJ5s+HMWNsuGuH\nDtaKOOMMmz8hIpIulBgCsHUrPPectSJWrYL+/W1U08Gxu/KKiIRQRRND2k9wC1KNGjaT+r334LXX\nbGZ1u3ZWg3j+edixw3WEIiLJpxZDgrZvt0L1o4/CZ59ZC2LgQC3BISLhoxZDilSrBhddZHMi3ngD\nNm+2JTjOOgtmzIBdu1xHKCJSOWoxJMG2bVaofvRRW9Rv4EAb0bTvvq4jE5FsphaDQ3vuaYlg3jyY\nOhVWr7Ylwi+4AN580+ZLiIikC7UYArJpk82uHjXKitR/+IPNrtYaTSKSKhquGlKRCPzrX5YgZsyA\nbt0sSbRrpzWaRCRYSgxpYMMGePJJq0XUrAlDhkCvXjYsVkQk2ZQY0khREbz+Ojz8sM2R6NsXrrgC\nmlZod1YRkfhUfE4jublwyinwz3/CRx/ZUhvt29uyG7Nmae9qEXFLLYaQ+PlnG/I6cqQVrocMsVnX\ndeq4jkxE0pW6kjJEJAIffGAJYsYMG/I6ZAgccYTryEQk3agrKUPk5Fi30tNPw/Ll0KgRdO0KnTrB\ns89qfSYRCZ5aDGlgxw6rR4wcCZ9/DoMG2exqzawWkbKoxZDB8vOhRw8oLIRXXoGvv7aZ1b16WbeT\niEgyqcWQpn78EcaNg4cesr2qr70Wune3JCIiAio+Z61du2D6dBgxAlasgCuvtB3ntPSGiKgrKUvl\n5cE551g300svwbJlNlFu8GC7LyKSKCWGDNKmjS25sWyZdS916QKnnw6zZ2uFVxHxT11JGWz7dps0\nN2IE/PILXHMN9OmjtZlEsoVqDFKqSMR2nBsxwtZmuuwymzTXqJHryEQkSKoxSKlycqCgAKZNg7lz\nrSXRujVceKGt1SQiEk0thiz1008wZgw8+CAcdBDccIPNsM7VVwWRjKGuJKmQHTvg+efhvvtsIb+h\nQ+Hii6FaNdeRiUhlKTFIpUQiNuT1vvtg/nybDzF4MNSr5zoyEako1RikUnJybHjrzJm2idDKlXDI\nIZYgVqxwHZ2IpJISg/xGq1YwdiwsWQK1a9v+1Oefr3WZRLKFupKkXFu2WKL4+9+hYUMrVJ91lgrV\nImGnGoMEbudOmDrV6hCbN8NNN1mheo89XEcmIvEoMUjKFBeq77nHupuuv972h6hZ03VkIhJNxWdJ\nmeJC9ezZtoHQ3Lk2F+LOO2HDBtfRiUhlKTFIpbRta1uOvveebSB06KG2N8Tata4jE5GKUmKQpDj0\nUBg9GhYvhipVbMmN/v1t32oRSS9KDJJUBxwA998PX3wBBx8MnTrZznIffug6MhHxS8VnCdTWrfDE\nE5YsmjWDYcPgxBOtTiEiwdKoJAm1HTtg0iT4y1+gTh24/XbbREgJQiQ4SgySFnbtsrkQd98NVavC\nbbfB2WdrspxIEJQYJK0UFdke1XfdZa2J226D886zPaxFJDnCOo9hLPAtsLiU5wuATcAC73ZbwPFI\nSOTmwjnn2EZB99xjy220agVPPWUzrEXEnaBbDL8HtgATgMPjPF8AXA+cXc7rqMWQ4SIReOMNa0F8\n9ZUVqfv0se4mEamYsLYY3gF+KOecdOnOkgDl5MBJJ9ne1GPHwjPP2NyIUaPgl19cRyeSXVyX/CJA\nB2ARMBNo6TYcCYNOneC11yw5vPyy7QvxyCNKECKpUsXx+88HGgHbgNOBacCh8U4cPnz4f+8XFBRQ\nUFAQfHTiVPv2MGOGTY67806491744x+hXz91MYnEU1hYSGFhYaVfJxXdOE2A6cSvMcRaBbQFvo85\nrhqDMHeuJYjPPrN5EH36QH6+66hEwiusNYby7ENJ0Md692OTgghgLYjZs23k0sSJ0KIFTJigUUwi\nyRZ0i2Ey0Bmojw1bvRMo/o73GDAEGAzsxLqTrgfmxnkdtRjkNwoLrQXxzTf2b8+emgchEk0T3CQr\nRSIwZw7ccQf88IMliPPP10xqEVBikCwXidhIpjvusIX7/vQnOPdcrcUk2U2JQQRLEDNnWnE6Lw/+\n7/9sfoQShGQjJQaRKEVF8NxzliAaNrRVXdu1cx2VSGql66gkkUDk5loxeskSuOgi6NHD1mZassR1\nZCLhp8QgGS0/HwYMsLkPv/89dOkCffvC6tWuIxMJLyUGyQrVq8PQobblaJMm0LYtXHUVfPut68hE\nwkeJQbJK7do2YmnZMitOt2xpy2z8+KPryETCQ4lBslKDBjBiBMyfD+vX20quI0ZooT4RUGKQLNe4\nsS3zPWeOzYNo2RKmTLFhryLZSsNVRaLMmQM33mjdTPffb0uAi6QrzWMQSZKiItsL4tZb4YgjbLnv\nFi1cRyWSOM1jEEmS3Fzo1QuWL4fOna3VMGiQLdYnkg2UGERKUa2aDXH99FOoVQsOOwyGD4ctW1xH\nJhIsJQaRctSta/WGjz6yiXK/+x2MH29dTiKZSDUGkQTNnQvXXgu7dtkQ144dXUckEp+KzyIpVFQE\nkyfDLbdYYrj3Xhv6KhImKj6LpFBuLvTubQXq5s3hqKNsJVfVHyQTKDGIVEKNGlaQXrgQVq60JDFh\nguoPkt7UlSSSRO+/b/WHSMTqDx06uI5IsplqDCIhUVQEkybBsGE2B+K++2D//V1HJdlINQaRkMjN\nhYsvtvpD48Y2e/qBB2DHDteRifijFoNIwD77DK6+GtauhZEjbbMgkVRQV5JIiEUi8M9/Wv2hXTtr\nQTRs6DoqyXTqShIJsZwc23N66VKbOX3kkTb34ddfXUcm8ltqMYg4sGIFXHONbTX60ENw8smuI5JM\npK4kkTQ0fboliKOPtuGtGr0kyaSuJJE0dNZZsGSJdS+1bg2PPGJrMIm4pBaDSEgsXWr7Pvz6Kzz2\nmNUhRCpDLQaRNNeyJbz1FgwcCKecYluMbt3qOirJRkoMIiGSmwuXXQaffGI7xh12GLz8suuoJNuo\nK0kkxF5/HQYPtvrDgw/CAQe4jkjSSZBdSXsDpwODgT8ApwF7JfpGIpK4k06CxYutm6m4OK2VWyVo\nZWWS3wM3Ak2ABcDX3vn7AW2A1cBfgXcDjdCoxSBZb+lSuPxyyM+HMWOgWTPXEUnYVbTFUKWM584F\nhgKfl/L8oVgLIhWJQSTrtWwJ77xjE+KOO852j7vuOsjLcx2ZZBrVGETS0IoVMGCA7Rg3diy0auU6\nIgmjIGsM12I1hRzgCaxb6dRE30hEkqdpU3jjDeta6tIF/vxnrbskyeMnMVwKbAJOAeoCfYB7ggxK\nRMqXk2NzHhYsgHnzbFmNf//bdVSSCfwkhuJmyBnAU8AnwYUjIolq2NDWXLr5ZujaFW69Va0HqRw/\nieHfwKtAV2A2UBvQgDmREMnJgd694eOPbe2lY4+1+yIV4acokQe0BlYCPwL1gIbAogDjiqXis4hP\nkQiMH29Lagwdav9q5FJ2CnLZ7RygO3A8EAHeAaZ591NFiUEkQWvWQP/+sH27JQrNe8g+QY5KegQY\nBHyM1RcGAQ8n+kYiklqNG9uSGhdeaPMeHn5Ys6bFHz+ZZDnQkpK6Qi6wFGgeVFBxqMUgUgnLl0Pf\nvlC7Nowbp/2ms0WQLYYvgAOjHh/oHRORNNG8Obz3HnTqBG3bwtSpriOSMPOTSd4GjgHmYXWFY4EP\ngZ+8x2cHFl0JtRhEkuSDD6BXLzjhBNtOtEYN1xFJUIIsPheU8VwEeCvRN60AJQaRJPrpJ7jqKpg7\nFyZPhqOOch2RBCHIxFCsNrsvuvd9om9WCUoMIgF45hm4+uqSoa252rorowSZGAYBfwJ+oaQAHQEO\n9vGzY7EZ0/8BDi/lnH9g+z1sA/phazHFUmIQCciaNTY5rlo1mDAB9t/fdUSSLEEWn28EWgGNgYO8\nm5+kADAO29inNF2BQ4BmwEBglM/XFZEkadwYCgutMH3UUdpKVPwlhpXAzxV8/XeAH8p4/mxgvHf/\nA2y3uH0q+F4iUkFVqsAdd9hopSuugJtugh07XEclrvhJDLcA7wOPAQ95t38k6f0PANZGPV6HLbch\nIg507GirtS5ZAp07w9q15f+MZJ6ydnAr9jjwOrAYqzHkkNzlMGL7v+K+9vDhw/97v6CggIKCgiSG\nICLF6tWz1VofeACOOQaeeALOOMN1VOJHYWEhhYWFlX4dP0WJBdgezxXVBJhO/OLzo0Ah8Iz3eDnQ\nGfg25jwVn0UceO89uOgiu919t+03LekjyOLzLGxk0n7YRj3Ft2R4CbjEu98eW701NimIiCMdO8L8\n+bB4se0U99VXriOSVPCTSVYTv3vnIB8/OxlrAdTHPvDvBIq/czzm/TsSG7m0FegPzI/zOmoxiDhU\nVAT33gsPPWQT4jp3dh2R+JGKCW4uKTGIhMCrr8Ill8Att8A119gGQRJeQSeGVtgKq9Wijk1I9M0q\nQYlBJCRWr4bu3aFFC3j8ca21FGZB1hiGY0NURwJdgL+SmoXzRCSEmjSxonReHnToACtWuI5Iks1P\nYugBnASsx2oArbGJaCKSpapXt13hBg605DBrluuIJJn8JIafgV3ATmAvbN2jRkEGJSLhl5MDQ4bA\nCy/A5ZfDPffYftOS/vwkhg+BOsBo4CNsXsO/ggxKRNJHx44wbx4895ztErd9u+uIpLISLUochC2/\nvSiAWMqi4rNIyG3bZonhq6/gxRdhH6165lwQxeemcY6tYvekEO8cEclCe+4JU6bAySdDu3bw8ceu\nI5KKKiuTTAFqYLOTP8KKzznYDOijsZFJm4ELA44R1GIQSSvFGwCNGQNnawyjM0HNYzgE++DviO3H\nALAGeBeb1bwy0TesICUGkTQzb57Nd7jqKlvGW5PhUk8zn0UkdNatg27dbAOgUaNs3wdJnSASw3mU\nvbz2C4m+WSUoMYikqS1b4PzzbT/pKVOgZk3XEWWPIBLDk1hiaAB0AOZ4x7tgw1XPTPTNKkGJQSSN\n7dgBgwfDwoUwY4ZGLKVKEKOS+mEznati6ySd590O846JiPiSnw+jR1sh+rjj4NNPXUckZfHT49cI\n+Cbq8bfAgcGEIyKZKifH9pVu1MiW7X7hBVtOQ8LHT2J4HZgNTMKaJD2B14IMSkQyV//+sP/+cM45\ntm3oWWe5jkhi+el7ygHOBTphNYe3gReDDCoO1RhEMsyHH1pSeOAB6N3bdTSZScNVRSTtLFkCp50G\nw4bBFVe4jibzVDQxlNWVtIXSh6tGsDWTREQq7LDD4K23bBmNTZtsZzhNhHMvXf4XqMUgksHWr7fk\n0LWr7S2t5JAc6koSkbS2caMlhtatbZZ0Xp7riNJfkFt7iogErl49eP11+PxzuPRS2LXLdUTZS4lB\nREKjVi2bGb12rQ1rVXJwQ4lBREJlzz3h5Zdtwx8lBzeUGEQkdPbcE6ZPh6+/tl3hlBxSS4lBREKp\nODl8+y1ccomSQyopMYhIaFWvDi+9ZMlh4EAoKnIdUXZQYhCRUKteHaZNg2XL4PrrQSPXg6fEICKh\nV7MmzJxps6SHD3cdTebTRnsikhb23htmz4ZOnaB2bRg61HVEmUuJQUTSRoMG8Nprlhxq1bK6gySf\nEoOIpJVGjSw5dO5ss6XPO891RJlHayWJSFpasABOPdV2gjv+eNfRhJPWShKRrNKmDTz9tLUYli1z\nHU1mUWIQkbR1yim2THfXrrZ0tySHagwiktb69YN16+CMM+Dtt21oq1SOagwikvYiEbj8cvj+e5g6\nFXLVFwKoxiAiWSwnxzb32bgRbr/ddTTpT4lBRDJC1arWWpg8GSZOdB1NelNXkohklE8+gRNOsJVZ\n27VzHY1b6koSEQFatYInnoDu3W0/B0mcWgwikpHuugtefRXmzIH8fNfRuFHRFoMSg4hkpKIiOPNM\naNkS7r/fdTRuqCtJRCRKbi489ZQVpF94wXU06UUtBhHJaB99ZDOj330XDj3UdTSppRaDiEgcRx8N\nf/4z9OgBP//sOpr0oBaDiGS8SAR69bJlukeOdB1N6qj4LCJShh9/hCOPtMRw5pmuo0mNsHYlnQYs\nBz4Hbo7zfAGwCVjg3W4LOB4RyVJ7723LdA8YAN984zqacAuyxZAHfAqcBHwFfAhcBESvnF4AXA+c\nXc5rqcUgIklxxx0wbx7MmmVrLGWyMLYYjgW+AFYDO4BngG5xzsvw/zUiEia3326rsI4e7TqS8Aoy\nMRwArI16vM47Fi0CdAAWATOBlgHGIyJCfj6MGwd//CN8+aXraMIpyI16/PT9zAcaAduA04FpQNyR\nxsOHD//v/YKCAgoKCiodoIhkp8MOg+uus3rDK69kTpdSYWEhhYWFlX6dIH8d7YHhWAEaYBhQBNxb\nxs+sAtoC38ccV41BRJJq505o3x4GD4bLLnMdTTDCWGP4CGgGNAGqAj2Bl2LO2YeSoI/17scmBRGR\npKtSxbqUbrlF+0XHCjIx7ASuBGYDS4Ep2IikQd4NoAewGFgIjAAuDDAeEZHdHH64bQl6ww2uIwmX\ndOlZU1eSiARi61arOYwdaxv8ZJIwdiWJiIRejRowYgQMGQK//uo6mnBQYhCRrNetGxx8MPztb64j\nCQd1JYmIACtXwjHHwOLFsP/+rqNJDi2iJyJSSTfdZIvtPf6460iSQ4lBRKSSfvjBNvN56y3bEjTd\nqfgsIlJJderYvIZhw1xH4pZaDCIiUbZvh+bNbYnu4493HU3lqMUgIpIE1arBXXfBjTfazm/ZSIlB\nRCRG796weTPMnu06EjeUGEREYuTm2rLcd92Vna0GJQYRkTguuAC++w6SsIp12lFiEBGJIy8Pbr3V\nWg3ZRqOSRERKsWMHNG0KL74Ibdu6jiZxGpUkIpJk+flw5ZXw4IOuI0kttRhERMrw/fdwyCGwdCns\nu6/raBKjFoOISADq1oWePWHUKNeRpI5aDCIi5Vi2DLp0gS+/hKpVXUfjn1oMIiIBadHCFtebOdN1\nJKmhxCAi4kP//vDkk66jSA11JYmI+LB5MzRqBJ99Bg0auI7GH3UliYgEqFYt2wJ04kTXkQRPiUFE\nxKd+/WDcONdRBE+JQUTEp86dbV7DsmWuIwmWEoOIiE+5udC9O0yd6jqSYCkxiIgkoEcPeP5511EE\nS6OSREQSsGuXLY0xf76NUgozjUoSEUmBvDw49VSYNct1JMFRYhARSdDpp2f2LGh1JYmIJGjDBtun\nYeNGqFLFdTSlU1eSiEiK1K8PjRvDggWuIwmGEoOISAV06gRvv+06imAoMYiIVEAmJwbVGEREKmD9\nemjVyuoNOSH9JFWNQUQkhfbbD/bYA9audR1J8ikxiIhU0JFHwsKFrqNIPiUGEZEKat1aiUFERKIc\neSQsWuQ6iuRTYhARqaBMbTGEtJb+GxqVJCKhs3071K4Nv/wSzpFJGpUkIpJi1arZyKTNm11HklxK\nDCIilVCvns1lyCRKDCIilVC/vi2ml0mUGEREKqF+fbUYREQkihKDiIjsRjUGERHZjWoMiTsNWA58\nDtxcyjn/8J5fBLQJOB4RkaRSV1Ji8oCRWHJoCVwEtIg5pytwCNAMGAiMCjCepCosLHQdwm8oJv/C\nGJdi8idsMdWrB8uWFboOI6mCTAzHAl8Aq4EdwDNAt5hzzgbGe/c/APYG9gkwpqQJ28UJiikRYYxL\nMfkTtpjq14c1awpdh5FUQSaGA4DolcrXecfKO6dhgDGJiCRV/fqwbZvrKJIryMTgd3Gj2HU8tCiS\niKSNTEwMQS771B4YjtUYAIYBRcC9Uec8ChRi3UxgherOwLcxr/UF0DSgOEVEMtUKrI4bGlWwoJoA\nVYGFxC8+z/Tutwfmpio4ERFx43TgU+wb/zDv2CDvVmyk9/wi4KiURiciIiIiIumjEfAmsAT4BLi6\nlPNSOSHOT0y9vVg+Bt4DjghBTMWOAXYC3UMSUwGwwDunMAQx1Qdewbo5PwH6BRwTQDVsaPZCYCnw\nl1LOS+V17iemVF/nfn9PkLrr3G9MBaTuOvcTk4vrPGn2BY707tfEuqDKqkm0I/iahJ+YjgP28u6f\nFpKYwCYYzgFeBs4LQUx7Yx/SxcOR64cgpuGU/CHVBzZitbGg7en9WwW7Xo6PeT7V17mfmFJ9nfuJ\nCVJ7nfuJKdXXuZ+YhpPgdR6mtZK+wTIawBZgGbB/zDmpnhDnJ6b3gU1RMQU9D8NPTABXAc8D3wUc\nj9+YegFTsbkqAEEvIuAnpvVAbe9+bewPZmfAcQEUD26sin2wfR/zvIuJn+XFlOrr3E9MkNrr3E9M\nqb7O/cSU8HUepsQQrQnWfP4g5rjLCXGlxRTtMkq+6aVCE0r/PXWjZImRVM4NKS2mZkBdrHvnI6BP\nCGIaDRwGfI11k1yTonhysaT1Lfb7WBrzvIvrvLyYoqXqOvfze0r1dV5eTC6u8/JicnWdJ1VN7Bd6\nTpznpgMdox6/TmpGMpUVU7Eu2P+QOimIB8qO6TmsCwLgSVLTxC4vppHAv4DqQD3gM+yPyGVMtwEj\nvPtNgZVArRTEVGwvrOlfEHPc1XVeVkzFUn2dQ+kxubrOofSYXF3nZcWU8HUethZDPtYMexqYFuf5\nr7CiYrGG3jGXMYEV4kZjXQA/BByPn5jaYpMGV2F/LI94sbmMaS3wKvAz1pR9G2jtOKYO2IcL2Jyb\nVcDvAo4p2iZgBnB0zHEX13l5MUHqr/PyYnJxnZcXk4vrvLyYXF/nlZIDTAD+XsY5qZ4Q5yemA7F5\nGO0DjqWYn5iijSP40Rp+YmqOffPNw4pli7FVd13G9DfgTu/+PliXTd0AYwIr/u3t3a+OfXCcGHNO\nqq9zPzGl+jr3E1O0VFznfmJK9XXuJyYX13nSHI8tmbEQG+q1AJsg53JCnJ+YxmDfDIqfnxeCmKKl\n4g/Gb0w3YCM2FlP2MNtUxVQf67ZZ5MXUK+CYAA4H5ntxfQzc6B13eZ37iSnV17nf31OxVFznfmNK\n5XXuJyYX17mIiIiIiIiIiIiIiIiIiIiIiIiIiIhItpoDnBJz7FpslmyirqTiSxYXUjIf4VYf5z8L\nHFTB9xIRkTIMAMbGHHuf+Es4lyUHm9hV0WW536QkMWz2cf7J2J4MIiKSZHWxVSeLP9CbAGuwdcIe\nwZbjfhVbY6Z4EbZ7sJmri4D7vGPHA5O9+83ZfaXWJtisU7DlCOZ7j5/AlkIGSwxtvdfeiSWZp7Bl\nE2ZgM1cXAxd45+djM59FRCQA0ylZXO0W4K9AD+wDGWz9mO+xJRXqAcujfrZ21M8NjTq+AEsIADdj\n3UPVgC+BQ7zj4ylZ4ri0FsN5wONx3g/gLeJvwiQSqLCtrioShMnAhd79nt7jjlg/PpSsYw/wI7Ad\n+7Z/LrbnORmYAAABXElEQVRKJtgicuujXvNZ77XAvuVPwVasXEXJN/3xQKdyYvsY6za6B2uV/BT1\n3NeUJB+RlFFikGzwEtbF0wbrulngHc+Jc+4u4FhsV7Azsb1yiXP+FCwhNMM2iFkR57XivX6sz724\nFgN3A7fH/HyRj9cQSSolBskGW7AWwThgknfsPawbJwfrSirwjtfAljGeBVxPyVr6a7B9pIutxJLI\n7dieAGD7SjfBNkMB272rME48OyipeeyHtVAmAvez+0qq+3nvKyIiAeiGfZAf6j3OwbaELC4+v4a1\nKvbFCsuLsG6e4q0ZO1JSfC421HvNA6OOnUBJ8XkMVkSG3WsM92C7oD2FDaVdhLViPog6J5/4rRAR\nEQlQDe/felhdoEEZ5xYPV61axjnJdArwYIreS0REPG9iH/ZLgEt8nH8F0D/QiEo8iwrPIiIiIiIi\nIiIiIiIiIiIiIiIiIiIiknn+HwC9g+9EmLPsAAAAAElFTkSuQmCC\n",
+ "text/plain": [
+ "<matplotlib.figure.Figure at 0x7fb0a512a8d0>"
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "#finding Id vs Vgs\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "%matplotlib inline\n",
+ "import numpy as np\n",
+ "import matplotlib.pyplot as plt\n",
+ "Vth=3.6;\n",
+ "Vgs=4;#voltage\n",
+ "#volt change beyond 3.6 causes a major increase in Id as it is cut off voltage\n",
+ "\n",
+ "#result\n",
+ "print('Id=0 from 0 to 2 so not shown in the graph')\n",
+ "x=np.linspace(2,3.6,300);\n",
+ "y=(-2.5*(x-3.6))**.5;\n",
+ "plt.plot(x,y)\n",
+ "plt.xlabel('Vgs(volts)');\n",
+ "plt.ylabel('Id(amps)');\n",
+ "plt.title('Id vs Vgs');\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 5.2,Page 217"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 1,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "drain current is 3.8 A\n",
+ "Vth=4V is assumed\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding drain current\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "V=4.5;#voltage\n",
+ "T=25;#degreeC\n",
+ "Id=3.8;\n",
+ "\n",
+ "#result\n",
+ "print \"drain current is\",round(Id,2), \"A\"\n",
+ "print('Vth=4V is assumed')"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 5.3,Page 219"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 2,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "MOSFET is IRF530N\n",
+ "lower limit of Vth is -4.0 V\n",
+ "upper limit of Vth is -2.0 V\n",
+ "current is 2.3 A\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding drain current of IRF530N\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Vgs=-5;#voltage\n",
+ "Vthl=-4;\n",
+ "Vthu=-2;\n",
+ "Id=2.3;#current\n",
+ "\n",
+ "#result\n",
+ "print('MOSFET is IRF530N')\n",
+ "print \"lower limit of Vth is\",round(Vthl,2), \"V\"\n",
+ "print \"upper limit of Vth is\",round(Vthu,2), \"V\"\n",
+ "print \"current is\",round(Id,2), \"A\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example5.5,Page 225"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 3,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "resistance is 1.51 kohm\n",
+ "load voltage is 36.68 V\n",
+ "Pq is 40.02 watt\n",
+ "Ps is 82.0 watt\n",
+ "Pl is 41.97 watt\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding Pq,Pl,Ps,resistance,load voltage\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "R1=22.0;#resistance\n",
+ "Vg=3.6;#voltage\n",
+ "Vd=56.0;\n",
+ "G=.98;#gain\n",
+ "Vi=40.0;\n",
+ "Rl=8.0;#load resistance\n",
+ "Vp=36.5;\n",
+ "\n",
+ "#calculation\n",
+ "Vr=Vd-Vg;\n",
+ "Ir=Vr/R1;\n",
+ "R2=Vg/Ir;\n",
+ "Va=(R1/(R1+R2))*Vi;\n",
+ "Vl=G*Va;\n",
+ "Il=Vp/Rl;\n",
+ "Pl=Vp*4.6/4;\n",
+ "Ps=Vd*4.6/pi;\n",
+ "Pq=Ps-Pl;\n",
+ "\n",
+ "#result\n",
+ "print \"resistance is\",round(R2,2), \"kohm\"\n",
+ "print \"load voltage is\",round(Vl,2), \"V\"\n",
+ "print \"Pq is\",round(Pq,2), \"watt\"\n",
+ "print \"Ps is\",round(Ps,2), \"watt\"\n",
+ "print \"Pl is\",round(Pl,2), \"watt\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 5.6,Page 232"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 4,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "current is 2.45 mA\n",
+ "resistance2 is 814.815 ohm\n",
+ "pick R2=R3=820ohm R1=R4=22 kohm\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding resistance and current \n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "R1=22.0;#resistance\n",
+ "V1=56.0;#voltage\n",
+ "V2=2.0;#voltage\n",
+ "\n",
+ "#calculation\n",
+ "I=(V1-V2)/R1;\n",
+ "R2=V2/I;\n",
+ "\n",
+ "#result\n",
+ "print \"current is\",round(I,2), \"mA\"\n",
+ "print \"resistance2 is\",round(R2*1000,3), \"ohm\"\n",
+ "print('pick R2=R3=820ohm R1=R4=22 kohm')"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 5.7,Page 234"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 5,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "load voltage is 10.01 V\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding load voltage \n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Vi=350.0;#voltage\n",
+ "f=100.0;#frequency\n",
+ "Rf=10000.0;#resistance\n",
+ "Ri=520.0;\n",
+ "\n",
+ "#calculation\n",
+ "Vp=(1+(Rf/Ri))*Vi*2**.5;\n",
+ "\n",
+ "#result\n",
+ "print \"load voltage is\",round(Vp/1000,2), \"V\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 5.8,Page 238"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 6,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "load rms voltage is 20.0 V\n",
+ "resistance is 54.04 kohm\n",
+ "current is 1.18 mA\n",
+ "load current is 4.41 A\n",
+ "supply power is 39.3 watt\n",
+ "load power is 38.9 W\n",
+ "power is 19.552 W\n",
+ "thermal resistance is 3.01 degreC/W\n"
+ ]
+ }
+ ],
+ "source": [
+ "#designing amplifier\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "P=50.0;#power\n",
+ "Z=4.7#impedence\n",
+ "R=4.0;#resistance\n",
+ "Ta=40.0;#degreeC\n",
+ "Tj=140.0;#degreeC\n",
+ "Vd=28.0;\n",
+ "R2=22.0;\n",
+ "\n",
+ "#calculation\n",
+ "Vr=(P*R)**.5;\n",
+ "Vp=Vr*2**.5;\n",
+ "Av=-Vr/1.23;\n",
+ "Rf=-Av*Z;\n",
+ "I=(Vd-2)/R2;\n",
+ "Vm=.63*Vd;\n",
+ "Ip=Vm/R;\n",
+ "Ps=Vd*Ip/pi;\n",
+ "Pl=Ip**2/2*R;\n",
+ "Pq=round(Ps)-Pl/2;\n",
+ "Qs=(Tj-Ta)/Pq-2.1;\n",
+ "\n",
+ "#result\n",
+ "print \"load rms voltage is\",round(Vp,2), \"V\"\n",
+ "print \"resistance is\",round(Rf,2), \"kohm\"\n",
+ "print \"current is\",round(I,2), \"mA\"\n",
+ "print \"load current is\",round(Ip,2), \"A\"\n",
+ "print \"supply power is\",round(Ps,2), \"watt\"\n",
+ "print \"load power is\",round(Pl,2), \"W\"\n",
+ "print \"power is\",round(Pq,3), \"W\"\n",
+ "print \"thermal resistance is\",round(Qs,2), \"degreC/W\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 5.9,Page 243"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 7,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "output voltage is 129.37 mV\n",
+ "load current is 32.34 mA\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding load current,output voltage \n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Vi=7.5e-3;#voltage\n",
+ "Ib=800e-9;#current\n",
+ "R=53.9e3;#resistance\n",
+ "\n",
+ "#calculation\n",
+ "Vo=11.5*Vi+Ib*R;\n",
+ "Id=Vo/4;\n",
+ "\n",
+ "#result\n",
+ "print \"output voltage is\",round(Vo*1000,2), \"mV\"\n",
+ "print \"load current is\",round(Id*1000,2), \"mA\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 5.10,Page 249"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 8,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "resistance is 0.31 ohm\n",
+ "thus pick a .33ohm rsistance\n",
+ "voltage is 0.55 V\n",
+ "power is 0.23 W\n",
+ "thermal resistance is 8.1 degreeC/W\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding resistance,voltage,power \n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "G=6.4;#A/V\n",
+ "I=5.0;#current\n",
+ "Pq=9.8;#W\n",
+ "Tj=140.0;\n",
+ "Ta=40.0;\n",
+ "R1=.33;\n",
+ "\n",
+ "#calculation\n",
+ "R=2/G;\n",
+ "Im=I/3;\n",
+ "Vr=Im*R1;\n",
+ "P=Vr*Im/4;\n",
+ "Qs=(Tj-Ta)/Pq-2.1;\n",
+ "\n",
+ "#result\n",
+ "print \"resistance is\",round(R,2), \"ohm\"\n",
+ "print('thus pick a .33ohm rsistance')\n",
+ "print \"voltage is\",round(Vr,2), \"V\"\n",
+ "print \"power is\",round(P,2), \"W\"\n",
+ "print \"thermal resistance is\",round(Qs,2), \"degreeC/W\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 5.11,Page 251"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 9,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "limit level current is 8.49 A\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding limit current \n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "P=200;#power\n",
+ "R=8;#ohm\n",
+ "\n",
+ "#calculation\n",
+ "Il=(P/R)**.5*2**.5;\n",
+ "Ilm=1.2*Il;\n",
+ "\n",
+ "#result\n",
+ "print \"limit level current is\",round(Ilm,2), \"A\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 5.12,Page 253"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 10,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "resistance is 0.1 ohm\n",
+ "power is 1.8 W\n",
+ "MOSFET power is 84.0 W\n",
+ "temperature is 468.4 degreeC\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding resistance,power,temperature \n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "I=6;#current\n",
+ "V=.6;#voltage\n",
+ "D=.5;#duty cycle\n",
+ "T=40;#temperature\n",
+ "\n",
+ "#calculation\n",
+ "Rs=V/I;\n",
+ "Pr=D*V*I;\n",
+ "Vp=28;\n",
+ "Pm=D*Vp*I;\n",
+ "Tj=T+Pm*5.1;\n",
+ "\n",
+ "#result\n",
+ "print \"resistance is\",round(Rs,2), \"ohm\"\n",
+ "print \"power is\",round(Pr,2), \"W\"\n",
+ "print \"MOSFET power is\",round(Pm,2), \"W\"\n",
+ "print \"temperature is\",round(Tj,2), \"degreeC\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 5.13,Page 255"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 11,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "maximum safe temperature is 89.05 degreeC\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding maximum safe temperature \n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "T=130;#temperature\n",
+ "P=19.5;#power\n",
+ "\n",
+ "#calculation\n",
+ "Ts=T-P*2.1;\n",
+ "\n",
+ "#result\n",
+ "print \"maximum safe temperature is\",round(Ts,2), \"degreeC\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 5.14,Page 257"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 12,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "reactance is 8.86 ohm\n",
+ "voltage across resistor is 10.03 V\n",
+ "-48 is the angle of the voltage in degrees\n",
+ "power dissipated by load is 12.5 watts\n",
+ "current across the resistance is 1.77 A\n",
+ "power supply is 15.8 W\n",
+ "power dissipated by transistor is 9.55 watts\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding 3 powers and current across resistance\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "V=15.0;#voltage\n",
+ "f=300.0;#frequency\n",
+ "L=4.7;#inductance\n",
+ "Vdc=28;#V\n",
+ "Pr=12.5;\n",
+ "\n",
+ "#calculation\n",
+ "Xl=2*pi*f*L;\n",
+ "Zload=sqrt(8**2+8.9**2);#magnitude of Zload\n",
+ "Vload=15.0;#msgnitude of Vload\n",
+ "Vr=Vload*8/Zload;\n",
+ "I=Vr/8*sqrt(2);\n",
+ "Psupply=Vdc*I/pi;\n",
+ "Pq=Psupply-Pr/2;\n",
+ "\n",
+ "#result\n",
+ "print \"reactance is\",round(Xl/1000,2), \"ohm\"\n",
+ "print \"voltage across resistor is\",round(Vr,2), \"V\"\n",
+ "print \"-48 is the angle of the voltage in degrees\";\n",
+ "print \"power dissipated by load is\",round(Pr,2), \"watts\"\n",
+ "print \"current across the resistance is\",round(I,2), \"A\"\n",
+ "print \"power supply is\",round(Psupply,2), \"W\"\n",
+ "print \"power dissipated by transistor is\",round(Pq,2), \"watts\""
+ ]
+ }
+ ],
+ "metadata": {
+ "kernelspec": {
+ "display_name": "Python 2",
+ "language": "python",
+ "name": "python2"
+ },
+ "language_info": {
+ "codemirror_mode": {
+ "name": "ipython",
+ "version": 2
+ },
+ "file_extension": ".py",
+ "mimetype": "text/x-python",
+ "name": "python",
+ "nbconvert_exporter": "python",
+ "pygments_lexer": "ipython2",
+ "version": "2.7.6"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter6.ipynb b/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter6.ipynb
new file mode 100644
index 00000000..0c7858a4
--- /dev/null
+++ b/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter6.ipynb
@@ -0,0 +1,528 @@
+{
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 6: Power Switches"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 6.1,Page 274"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 2,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "collector current is 1.81 A\n",
+ "load power is 49.32 W\n",
+ "transistor power is 1.45 W\n",
+ "least value of base current is 90.67 mA\n",
+ "max value of base resistance is 4.85 ohm\n",
+ "thus pick Rb=33ohm\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding Ic,Pload,Pq,Ibase,Rbase\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Vs=28.0;#V\n",
+ "Vi=5.0;#V\n",
+ "Rl=15.0;#ohm\n",
+ "Vc=.8;#V\n",
+ "b=20.0;\n",
+ "\n",
+ "#calculation\n",
+ "Ic=(Vs-Vc)/Rl;\n",
+ "Pl=Ic**2*Rl;\n",
+ "Pq=Ic*Vc;\n",
+ "Ib=Ic/b*1000;\n",
+ "Rb=(Vi-.6)/Ib;\n",
+ "\n",
+ "#result\n",
+ "print \"collector current is\",round(Ic,2), \"A\"\n",
+ "print \"load power is\",round(Pl,2), \"W\"\n",
+ "print \"transistor power is\",round(Pq,2), \"W\"\n",
+ "print \"least value of base current is\",round(Ib,2), \"mA\"\n",
+ "print \"max value of base resistance is\",round(Rb*100,2), \"ohm\"\n",
+ "print ('thus pick Rb=33ohm')"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 6.4,Page 282"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 3,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "load resistance is 554.0 ohm\n",
+ "thus pick Rl=560ohm\n",
+ "max value of Rb is 3.0 kohm\n",
+ "thus pick Rb=2.2kohm\n",
+ "load current is 49.46 mA\n",
+ "load power is 685.08 mW\n",
+ "power delivered is 7.42 mW\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding Pload,Pq,Iload,resistances\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Vd=28.0;#V\n",
+ "f=100.0;#frequency\n",
+ "I=50.0;#current\n",
+ "Rl1=560.0;\n",
+ "Vp=2.4;\n",
+ "Ib=500.0;#microAmp\n",
+ "D=.5;#duty cycle\n",
+ "\n",
+ "#calculation\n",
+ "Rl=(Vd-.3)/I;\n",
+ "Rb=(Vp-.9)/Ib;\n",
+ "Vl=Vd-.3;\n",
+ "Ip=Vl/Rl1;\n",
+ "Pl=D*Vl*Ip;\n",
+ "Pq=D*Ip*.3;\n",
+ "\n",
+ "#result\n",
+ "print \"load resistance is\",round(Rl*1000), \"ohm\"\n",
+ "print('thus pick Rl=560ohm')\n",
+ "print \"max value of Rb is\",round(Rb*1000,2),\"kohm\"\n",
+ "print('thus pick Rb=2.2kohm')\n",
+ "print \"load current is\",round(Ip*1000,2), \"mA\"\n",
+ "print \"load power is\",round(Pl*1000,2), \"mW\"\n",
+ "print \"power delivered is\",round(Pq*1000,2), \"mW\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 6.5,Page 286"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 4,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "time of rise is 788.48 ns\n",
+ "time of rise is 4.65 microsec\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding time of rise\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "C=640.0;#capacitor\n",
+ "R1=560.0;#load resistance\n",
+ "R2=3.3;#kohm\n",
+ "\n",
+ "#calculation\n",
+ "t1=2.2*R1*C;\n",
+ "t2=2.2*R2*C;\n",
+ "\n",
+ "#result\n",
+ "print \"time of rise is\",round(t1/1000,2), \"ns\"\n",
+ "print \"time of rise is\",round(t2/1000,2), \"microsec\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 6.6,Page 287"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 2,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "resistance is 682.5 ohm\n",
+ "pick up resistance=680 ohm\n",
+ "rise time is 957.44 ns\n"
+ ]
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEZCAYAAACervI0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGQJJREFUeJzt3XuQXGWdh/FnkgkxkISEiyEEcCQsiBguLuAFLy0qxVYE\nZNeyQNwFRbC2RNDVFXBdmdWqJUKJWrulW6uAQTQqoEi0QIhLLxfxgiQkhITInUAyYU0gIQESkt4/\n3tNMz6RnpufS/b7n9POp6prTZ87p/tFkzve873nf0yBJkiRJkiRJkiRJkiRJkiRJasADwLtiFyFJ\nKbsF+Lc6608B1gDjRvHa3we+Oor9U3sfadhG8wckNdv3gY/WWf/3wLXAjpZWI0lqmUnAc8A7a9ZN\nB14E5gATgW8CT2ePbwC7ZNudBdzZ7/V2ALOBc4GtwMvAJuAXdd77O8Dl/db9AvhMtnwhsBrYCKwE\njq/zGgO9z+M123cD1wE/yF5rKfBXwMVAD/AE8P6a19wduBJ4Jnv/r+KJnKSC+m/guzXPPwncly1/\nBfgtsFf2uDtbBwMHwIHZ8tU129bzTuDJmufTgS3APsAh2e/2yX53QM3r9lfvfR6jbwC8SDjIjwfm\nEwLi4uz5J4BHa/b9OSGcJgF7A78nBI0kFc5xwAZ6z+zvBi7Ilh8GTqzZ9gTCwRUaC4DB+uY7CGff\n1dbHOcCibPkgwtn5e4EJQ9Rf7336B8Cva353EqG10JE9n5LVPRWYAbwEvKZm+9OB/xmiBqkum45K\n3d3A/wGnErpvjgF+lP1uX8JBuurJbN1YqAA/JhxgAT4C/DBbfpjQFdRNCIIFwMxRvNe6muUXCf+9\nlZrnAJOB1xECZw0hFDcA/0VoCUjDZgAoD64B/oFwQfgW4Nls/TNAV812B2TrADYDu9b8bh/6qjC0\nBcCHCAfeY4Eb+v3undnvKsDXBniNRt6nUU8RrifsSeiSmk64JjBnDN9DbcQAUB5cQ+gj/wShj7xq\nAfAleq8BfJlwMRXgfuAw4AhCl0l3v9fsYeB++6olhLPx7xGCZ2O2/mBCF85EwgH5JWD7AK/RyPs0\nag1wK3AFoWtoHKFV5JwCSYV2O/AX+va5TwS+RTjrf4YwImiXmt9/kdBaeAI4g3CQrh6MDwIWE7pR\nfjbI+34p2+/vatbNIVx83ZjVdBM7tzCq6r1P7TWASwgBV/U++l707czev9q1NRX4NqE18BzhgviH\nB6lfiuI1hD+SJcCDwKXZ+j2A24BVhLOZaVGqkyQ1VbUPthP4HfAO4DLgC9n6C4F5EeqSJLXIrsAf\nCX2yKwnD2SA0m1fGKkqS1DzjCF1Amwhn/hD6Qqs6+j2XJBXM7oQuoPew8wF/fevLkSR1tuh9ngd+\nBfw1YVjcPsBawuSZdf03nj17duWRRx5pUWmSVBiPEEaeNaSZ8wD2oneEzyTCOO7FhCFzZ2brzwRu\n7L/jI488QqVS8VGpcMkll0SvIZWHn4WfhZ/F4A/CvJCGNbMFMJMwaWdc9vgB8BtCCPwUOJtw0yvH\nMEtSBM0MgGXAm+usX0+Y7CJJishbQSSuVCrFLiEZfha9/Cx6+VmMXMfQm0RRyfqzJEkN6ujogGEc\n120BSFKbMgAkqU0ZAJLUpgwASWpTBoAktSkDQJLalAEgSW3KAJCkNmUASFKbatXtoCWpuCoVePll\n2Lix72PTpp3XnXYazJkTu2LAAJDUzrZv3/kgPdTzgdaNHw9TpsDUqb2P/s+nToUJE2L/V7/KewFJ\nypdKBV58cfgH6HrrXnoJJk/e+SBd78A92MF9yhSYODH2JzPsewGlGwDbt8M4L1FIhbFtWzgIj+aA\nXV23yy6NnW0PtW7XXQt1nClOAGzdmlRTSWpLlQps3jw2B+2tW4d/Zj3Q2bbHhrqKEwCbN4d0ljR8\nW7eOvC+7dt2mTTBp0sjOtvs/nzQJOlI95BTDcAMg3YvAr7wSuwIpTS+8AF/+MvT0DHwQ37595wN0\nvYPyAQcMvs3kydCZ7mFCo5Pu/1kDQKrvppvgnnvgvPMGPrBPnOjZtoZkAEh5s3AhnH02nHFG7EqU\nc6meIlQqTz8N++4buw4pLdu2wYwZsHw5zJwZuxolpjhfCWkLQNrZXXfB7Nke/DUmDAApTxYuhJNO\nil2FCiLdANi2LXYFUnp++UsDQGMm3QCwBSD19dBDsGULHHlk7EpUEAaAlBcLF8IHPuDwTo0ZA0DK\nC/v/NcYMACkP1q+HJUvg+ONjV6ICaWYA7A/cDiwHHgDOz9Z3A6uBxdnjxLp7GwBSr1tugVIp3E9H\nGiPNnAm8DfgssASYDPwJuA2oAFdkj0H2dhSQ9Cq7f9QEzWwBrCUc/AFeAFYAs7LnQ1/FsgUgBdu2\nwa9/DXPnxq5EBdOqawBdwFHA77LnnwbuB64EptXdwwCQAmf/qklaEQCTgeuBCwgtge8ArweOBNYA\nX6+7lwEgBU7+UpM0+26gE4AbgGuBG7N162p+/z1gYb0du6+9Fu69F4BSqUSpVGpelVLKFi6En/wk\ndhVKULlcplwuj3j/Zs4o6QDmA38hXAyumkk48ydbfwzwkX77VirXXQcf+lATy5Ny4KGH4L3vhaee\ncgKYhpTSN4IdB3wUWEoY7gnwReB0QvdPBXgM+GTdvR0FJDn7V03VzAC4i/rXGG5uaG+vAUghAL7w\nhdhVqKCcCSylytm/ajIDQEqVs3/VZAaAlKpq/7/UJAaAlKLq7F8DQE1kAEgpcvavWiDdAHAYqNqZ\ns3/VAukGgC0AtTPv/qkWMACk1Pjdv2oRA0BKjbN/1SIGgJQau3/UIgaAlBJn/6qF0g0ARwGpHTn7\nVy2UbgDYAlA7cvavWsgAkFLh7F+1mAEgpeLuu539q5YyAKRUOPpHLWYASKkwANRi6QaAo4DUTpz9\nqwjSDQBbAGonzv5VBAaAlALv/qkIDAAptg0b4L77nP2rljMApNhuvtnZv4rCAJBic/SPIjEApJic\n/auI0g0Ah4GqHTj7VxGlGwC2ANQObr4Z5s6NXYXalAEgxbR4MRxzTOwq1KYMACmmpUvhiCNiV6E2\nZQBIsfT0wNatMGtW7ErUppoZAPsDtwPLgQeA87P1ewC3AauAW4Fpdfc2AFR0y5bB4Yd7+wdF08wA\n2AZ8FjgMeCvwKeBQ4CJCABwM/CZ7XmdvRwGp4JYuDQEgRdLMAFgLLMmWXwBWALOAk4H52fr5wAfr\n7m0LQEVnACiyVl0D6AKOAn4PzAB6svU92fOdGQAqOgNAkXW24D0mAzcAFwCb+v2ukj120v3889Dd\nDUCpVKJUKjWvQqnVXnkFVq6Eww6LXYlyrFwuUy6XR7x/s68+TQB+CdwMfDNbtxIoEbqIZhIuFL+h\n336VyvTpsH59k8uTIlm+HE49FVatil2JCqQjDCho+LjezC6gDuBK4EF6D/4ANwFnZstnAjfW3dsu\nIBWZ3T9KQDO7gI4DPgosBRZn6y4G5gE/Bc4GHgc+XHdvRwGpyAwAJaCZAXAXA7cw3jfk3rYAVGRL\nl8I558SuQm0u7ZnAlbrXh6X8swWgBKQbAOPGwY4dsauQxt769fD889DVFbsStbl0A6Cz024gFdOy\nZTBnTjjJkSJK91+gAaCisvtHiUg7ABwJpCIyAJSIdANgwgRbAComA0CJSDcA7AJSEW3fHmYBv+lN\nsSuRDACppR59FPbeG3bfPXYlkgEgtZTdP0qIASC10v33GwBKhgEgtZItACUk7QBwGKiKxgBQQtIN\nAIeBqmg2boSeHjjooNiVSEDKAWAXkIrmgQfgjW+E8eNjVyIBBoDUOnb/KDEGgNQqS5fCEUfErkJ6\nlQEgtYotACUm7QBwFJCKolLpvQ20lIh0A8BRQCqSJ56AKVNgzz1jVyK9Kt0AsAtIRWL3jxJkAEit\nYAAoQQaA1AreA0gJMgCkVrAFoAQ1EgAHNrhubDkKSEWxZQs8+SQcckjsSqQ+GgmAG+qsu26sC9mJ\no4BUFMuXh4P/hAmxK5H66Bzkd4cCbwR2B/4W6AAqwFTgNc2vzC4gFYTdP0rUYAFwMHASIQBOqlm/\nCTinmUUBBoCKwwBQogYLgF9kj7cB97SmnBoGgIpi6VKYOzd2FdJOBguAqnOzR1Ul+/nxBva9CpgL\nrAOqc+C7gU8Az2bPLwZu2bkyA0AFUKl4Ezglq5EA+BW9B/1JwKnAMw2+/tXAfwDX1KyrAFdkj0Eq\nMwBUAM88E/4tz5gRuxJpJ40EwPX9nv8IuLvB178T6KqzvmPIPR0GqiKw/18JG8lEsIOBvUf5vp8G\n7geuBKbV3cJhoCoCA0AJa6QF8AK9XUAVoAe4cBTv+R3gK9nyV4GvA2f336j7jjteDYBSqUSpVBrF\nW0qRLF0KJ5wQuwoVVLlcplwuj3j/obtiRq8LWEjvReBGflepXHopPPcczJvX3OqkZpozB665Bo46\nKnYlagMdHR0wjON6Iy0AgFOAdxFaAP9LOGiP1ExgTbZ8KrCsfmVeBFbOvfwyPPwwHHpo7EqkuhoJ\ngHnAMcAPCclyPvB2wvDNoSwA3g3sBTwFXAKUgCMJYfIY8Mn6lRkAyrkVK+DAA+E1zZ84L41EIwEw\nl3DA3p49/z6whMYC4PQ6665qrDJHASnnvACsxDUyCqhC35E60+i9KNw8jgJS3hkAStxgAfBt4B3A\nvwP3Ec785wN/ytY1l11AyjsDQIkbrAtoFXA5sC+wCHiC0PVzIbC2+ZUZAMo5A0CJG6wF8E3CjeDe\nDfyZcEvoywkXbQ9uemUGgPKspyeMAtpvv9iVSANq5BrA44SRQEcCpxGGbq5oYk2BAaA8W7Ys3ACu\noxVTbaSRaSQAOoGTCfcAugVYSWgNNJejgJRndv8oBwa7BnAC4Yx/LvAHwpj+cwm3hmg+WwDKswcf\nhKOPjl2FNKjBWgAXEb4I5lDCN4L9iFYd/MFhoMq31ath//1jVyENarAWwPEtq6IeWwDKs6efhlmz\nYlchDWokt4NuDQNAeWYAKAcMAGmsvfgibN4Me+0VuxJpUAaANNaqZ/8OAVXi0g4Ah4Eqj+z+UU6k\nGwCOAlJeGQDKiXQDwC4g5dXq1QaAcsEAkMaaLQDlhAEgjTUDQDlhAEhjzQBQTqQdAI4CUh4ZAMqJ\ndAPAUUDKox07YO1a2Hff2JVIQ0o3AOwCUh6tWwfTpsHEibErkYZkAEhjye4f5YgBII0lA0A5YgBI\nY8lJYMqRtAPAUUDKG1sAypF0A2D8+DCiYseO2JVIjTMAlCPpBkBHR2gFbN8euxKpcQaAciTdAACv\nAyh/DADlSLMD4CqgB1hWs24P4DZgFXArMG3AvQ0A5Y0BoBxpdgBcDZzYb91FhAA4GPhN9rw+A0B5\n8sILYeDC9OmxK5Ea0uwAuBPY0G/dycD8bHk+8MEB9zYAlCd+FaRyJsY1gBmEbiGynzMG3NKhoMoT\nu3+UM52R37+SPXbS3d0NW7bA5ZdTOuUUSqVSSwuThs1JYGqxcrlMuVwe8f6taKt2AQuBOdnzlUAJ\nWAvMBG4H3tBvn0qlUoEDD4RFi8JPKXWXXgobNsBll8WuRG2qI3Q/Nnxcj9EFdBNwZrZ8JnDjgFt6\nDUB5YheQcqbZAbAA+C1wCPAU8DFgHvB+wjDQ47Pn9RkAyhMDQDnT7GsApw+w/n0N7W0AKE8MAOVM\n+jOBHQWkvDAAlDNpB4BfC6m8eOUVePZZmDkzdiVSw9IOALuAlBc9PbDnnuGkRcoJA0AaC3b/KIcM\nAGksGADKIQNAGgvOAlYOpR8AjgJSHtgCUA6lHwC2AJQHBoByKO0AcBio8sIAUA6lHQC2AJQXBoBy\nyACQRqtSCQGw336xK5GGxQCQRmvjxvAtYFOnxq5EGhYDQBotu3+UU+kHgMNAlToDQDmVdgA4Ckh5\n4CQw5VTaAWAXkPLAFoByygCQRssAUE4ZANJoGQDKKQNAGi0DQDmVfgA4CkipcxKYcirtAHAUkFK3\nbRusXw8zZsSuRBq2tAPALiClbs0aeO1rYfz42JVIw2YASKNh/79yzACQRsNJYMoxA0AaDVsAyrH0\nA8BRQEqZAaAcSz8AbAEoZQaAciztAHAYqFJnACjHOiO+9+PARmA7sA04dqctbAEodU4CU47FDIAK\nUALWD7iFAaCUVb8K0haAcip2F1DHoL81AJSyDRtg4kTYbbfYlUgjEjMAKsAi4F7gnLpbGABKmWf/\nyrmYXUDHAWuAvYHbgJXAnX22cBioUmYAKOdiBsCa7OezwM8JF4FfDYDu7m549FFYtYpSuUypVGp9\nhdJg1q6FffaJXYXaWLlcplwuj3j/wfvgm2dXYDywCdgNuBX4t+wnQKVSqUC5DN3d4aeUmnnzwp1A\nL7ssdiUSAB0dHTCM43qsFsAMwll/tYYf0nvw7+U1AKWsp8choMq1WAHwGHDkkFsZAErZ2rVw9NGx\nq5BGLPYw0MEZAEpZT49fBKNcSz8AHAWkVHkRWDmXdgB4LyClzBaAci7tALALSKnauhU2boQ994xd\niTRiBoA0EuvWwd57w7i0/4SkwaT9r9cAUKrs/1cBGADSSNj/rwJIPwAcBaQU2QJQAaQfALYAlCJb\nACqAtAPAYaBKlS0AFUDaAWALQKmyBaACMACkkbAFoAIwAKSRsAWgAkg7AKqTbLZvj1uH1J8tABVA\n2gEAtgKUnpdegi1bYPr02JVIo5J+ADgSSKmpdv90xPpCPWlspB8AtgCUGvv/VRAGgDRc9v+rIAwA\nabhsAaggDABpuHp6bAGoEPIRAN4QTimxC0gFkX4AOApIqbELSAWRfgDYBaTU2AJQQRgA0nDZAlBB\nGADScNkCUEEYANJwbNkCW7fC1KmxK5FGLR8B4CggpaI6BNTbQKgA8hEAtgCUCvv/VSCxAuBEYCXw\nZ+DCQbd0GKhSYv+/CiRGAIwH/pMQAm8ETgcOHXDrNm8BlMvl2CUkI4nPIpEWQBKfRSL8LEYuRgAc\nCzwMPA5sA34MnDLg1gZA7BKSkcRnkUgLIInPIhF+FiPXGeE9ZwFP1TxfDbxlwK07O+H662HFimbX\nlaZ77oErrohdRRpS+CwWLYLTTotbgzRGYgRAZVhbn3023HEHrF7dpHISt3Fj+/6395fCZ3HMMXDi\niXFrkMZIjLFsbwW6CdcAAC4GdgBfq9nmYWB2a8uSpNx7BDgodhGD6SQU2QXsAixhsIvAkqRC+Rvg\nIcKZ/sWRa5EkSZIUU+OTxIptf+B2YDnwAHB+3HKSMB5YDCyMXUhk04DrgRXAg4Trau3qYsLfyDLg\nR8DEuOW01FVAD+G/vWoP4DZgFXAr4d9KbowndAt1ARNo7+sD+wBHZsuTCV1m7fpZVP0T8EPgptiF\nRDYf+Hi23AnsHrGWmLqAR+k96P8EODNaNa33TuAo+gbAZcAXsuULgXmtLmo03gbcUvP8ouwhuBF4\nb+wiItoPWAS8h/ZuAexOOOgpnO0+BEwnBOFC4H1RK2q9LvoGwEqgOlV9n+z5gFK7GVy9SWKzItWS\nki5C0v8+ch0xfQP4Z8KQ4Xb2euBZ4GrgPuC7wK5RK4pnPfB14EngGeA5wklCO5tB6BYi+znofUtS\nC4DhTRJrD5MJ/b0XAC9EriWWDwDrCP3/7X4f5k7gzcC3s5+bad9W8mzgM4QTpH0JfytnxCwoMRWG\nOKamFgBPEy5+Vu1PaAW0qwnADcC1hC6gdvV24GTgMWABcDxwTdSK4lmdPf6YPb+eEATt6Gjgt8Bf\ngFeAnxH+rbSzHkLXD8BMwolTbjhJrFcH4SD3jdiFJObdtPc1AIA7gIOz5W76zqJvJ0cQRshNIvy9\nzAc+FbWi1uti54vA1dGTF5Gzi8DgJLGqdxD6u5cQuj4W03v7jHb2bhwFdAShBXA/4ay3XUcBQRjx\nUh0GOp/Qam4XCwjXPrYSrp1+jHBhfBE5HQYqSZIkSZIkSZIkSZIkSZIkqS3tDvxjtjwTuG4MX/s8\n4KwxfL3+Tgb+tYmvL0mF1kXfmZFjpYMwGa+zCa9d+x5LaK/JTIootXsBSaM1j3CTsMXAT+kNg7MI\n91O6lXBPofOAzxPuqHkP4ZbCZPveDNxLuOXCIdn64wi31n0le14GriDMyF0BHAP8nDAD86vZNl30\nDaPPA5dky+cTZrDeT5jRCeHGXfcAJ4zgv1uS2t7r6D3o1i6fRfiWud2AvYDngXOz311BuNsqwG+A\ng7Llt2TPIdxX5XM173M7cGm2fD5hSv4Mwj2sniIEShd9A+BzwJez5afpPdOfWrPNx2jfe/uoxZrZ\nnJVi6BhgGcJBe3P2eI7em8otAw4nhMPb6XvdYJfs5wHAXf1er3pPogeyR/U+7I8S7mS7cZD6lhK+\nwvBG+t7p9Rm855NaxABQO3m5ZnlHzfMdhL+FccAGwpfv1NM/UGr37//anYTuotpu1kk1y3OBdwEn\nAf8CvCnbbxx+L4ZaxGsAKppNwJRh7lM9sG8iXB/4UM36w7PlJ+i9z3qjeoDXEu7QOJHwxTaV7HUP\nIFxHuIgwcmlyts/M7L2kprMFoKL5C3A3oVtnBb1n0/2/Han/cvX5GcB3gC8R+ugXELpr7iJcOK5n\noG9e2gZ8BfgDoc//wWz9eOAHhAN/B/AteruLjsXvO5CkpFSHge4y1IajMI4wDNQTM7XE+NgFSDky\nnvANdUua9PonEVow5Sa9viRJkiRJkiRJkiRJkiRJkqT8+3+DM0mrOvS/VQAAAABJRU5ErkJggg==\n",
+ "text/plain": [
+ "<matplotlib.figure.Figure at 0x7f4e9da3f2d0>"
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "#finding Rpick up,time of rise\n",
+ "\n",
+ "#initialisation of variable\n",
+ "%matplotlib inline\n",
+ "import numpy as np\n",
+ "import matplotlib.pyplot as plt\n",
+ "Vol=0.7;\n",
+ "Iol=40.0/1000;#current\n",
+ "Rpullup1=680.0;\n",
+ "C=640.0;\n",
+ "Epullup=28.0;\n",
+ "#for plotting\n",
+ "x=[0, .1, 1.9, 4.1, 5, 5.1, 5.3, 5.6, 6.0, 9.3];\n",
+ "y=[27.8, .1, .1, .1, .1, 5, 13.5, 21.0, 27.0, 27.8];\n",
+ "\n",
+ "#calculation\n",
+ "Rpullup=(Epullup-Vol)/Iol;\n",
+ "trise=2.2*Rpullup1*C;\n",
+ "plt.plot(x,y,'r');\n",
+ "plt.xlabel ('time(mus)')\n",
+ "plt.ylabel ('Vout')\n",
+ "plt.title ('Vout vs time')\n",
+ "\n",
+ "#result\n",
+ "print \"resistance is\",round(Rpullup,2), \"ohm\"\n",
+ "print('pick up resistance=680 ohm');\n",
+ "print \"rise time is\",round(trise/1000,2), \"ns\"\n",
+ "\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 6.7,Page 289"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 5,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "worst case resistance is 0.286 ohm\n",
+ "load current is 4.45 A\n",
+ "load voltage is 26.73 V\n",
+ "load power is 47.62 W\n",
+ "drop voltage is 1.27 V\n",
+ "power is 2.27 W\n",
+ "temperature is 182.6 deg.C\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding worst case resistance and power characteristics\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "R1=.11;#resistance\n",
+ "Vd=28.0;#voltage\n",
+ "R2=6.0;#ohm\n",
+ "D=.4;#duty cycle\n",
+ "Q=62.0;\n",
+ "\n",
+ "#calculation\n",
+ "Ro=2.6*R1;\n",
+ "Ip=Vd/(R2+Ro);\n",
+ "Vl=Ip*R2;\n",
+ "Pl=D*Vl*Ip;\n",
+ "Vq=Ip*Ro;\n",
+ "Pq=D*Vq*Ip;\n",
+ "T=40+round(Pq*10)/10*Q;\n",
+ "\n",
+ "#result\n",
+ "print \"worst case resistance is\",round(Ro,3), \"ohm\"\n",
+ "print \"load current is\",round(Ip,2),\"A\"\n",
+ "print \"load voltage is\",round(Vl,2), \"V\"\n",
+ "print \"load power is\",round(Pl,2), \"W\"\n",
+ "print \"drop voltage is\",round(Vq,2), \"V\"\n",
+ "print \"power is\",round(Pq*10,2)/10, \"W\"\n",
+ "print \"temperature is\",round(T,2), \"deg.C\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 6.8,Page 292"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 6,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "voltage is 150.0 V\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding voltage\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "L=10.0;#inductor\n",
+ "I=4.5;#current\n",
+ "t=300.0#time\n",
+ "\n",
+ "#calculation\n",
+ "V=L*I/t;\n",
+ "\n",
+ "#result\n",
+ "print \"voltage is\",round(V*1000,2), \"V\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 6.9,Page 298"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 1,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "R1+R2 is 682.5 ohm\n",
+ "pick R1=330ohm & R2=360ohm as they divide Vd setting 8V<Vg<18V\n",
+ "node voltage for V1 is 28.0 V\n",
+ "node voltage for V2 is 0.7 V\n",
+ "gate voltage is 15.31 V\n",
+ "gate & source diff is -12.69 V\n",
+ "load voltage is 26.73 V\n",
+ "load current is 2.23 A\n",
+ "load power is 47.63 W\n",
+ "Pq is 2.26 W\n",
+ "thermal resistance is 44.92 degreeC/W\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding resistance and power characteristics\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Rl=12.0;#load resistance\n",
+ "V1=.8;#voltage\n",
+ "V2=2.4;#voltage\n",
+ "D=.8;#duty cycle\n",
+ "Tj=150.0;#degreeC\n",
+ "Ta=40.0;#degreeC\n",
+ "Vd=28.0;\n",
+ "Vo=.7;\n",
+ "I=40.0;#mA;\n",
+ "R1=330;\n",
+ "R2=360;\n",
+ "Vn1=28;\n",
+ "Vn2=.7;\n",
+ "\n",
+ "#calculation\n",
+ "k=(Vd-Vo)/I;\n",
+ "Vg=R2*Vd/(R1+R2)+Vn2;\n",
+ "Vgs=Vg-Vd;\n",
+ "Vl=Vd*Rl/(Rl+.57);\n",
+ "Il=Vl/Rl;\n",
+ "Pl=D*Vl*Il;\n",
+ "Vq=Il*.57;\n",
+ "Pq=D*Vq*Il;\n",
+ "Q=(Tj-Ta)/Pq-3.7;\n",
+ "\n",
+ "#result\n",
+ "print \"R1+R2 is\",round(k*1000,2), \"ohm\"\n",
+ "print('pick R1=330ohm & R2=360ohm as they divide Vd setting 8V<Vg<18V')\n",
+ "print \"node voltage for V1 is\",round(Vn1,2),\"V\"\n",
+ "print \"node voltage for V2 is\",round(Vn2,2), \"V\"\n",
+ "print \"gate voltage is\",round(Vg,2), \"V\"\n",
+ "print \"gate & source diff is\",round(Vgs,2), \"V\"\n",
+ "print \"load voltage is\",round(Vl,3), \"V\"\n",
+ "print \"load current is\",round(Il,2), \"A\"\n",
+ "print \"load power is\",round(Pl,2), \"W\"\n",
+ "print \"Pq is\",round(Pq,2), \"W\"\n",
+ "print \"thermal resistance is\",round(Q,2), \"degreeC/W\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 6.10,Page 305"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 8,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "time is 1.1 micro s\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding time\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "I=40.0;#current\n",
+ "Q=44.0;#nC\n",
+ "\n",
+ "#calculation\n",
+ "t=Q/I;\n",
+ "\n",
+ "#result\n",
+ "print \"time is\",round(t,2), \"micro s\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 6.11,Page 313"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 9,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "current is 1.8 A\n",
+ "load voltage is 26.97 V\n",
+ "power is 40.85 W\n",
+ "high side voltage is 0.67 V\n",
+ "high side power is 1.03 W\n",
+ "low side voltage is 0.36 V\n",
+ "low side power is 0.55 W\n",
+ "IC power is 1.56 W\n",
+ "thermal resistance is 55.49 degreeC/W\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding different voltages and power\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Rl=15.0;#load resistance\n",
+ "D=.85;#duty cycle\n",
+ "Ts=60.0;#degreeC\n",
+ "Vd=28.0;#voltage\n",
+ "R1=.375;\n",
+ "R2=.2;\n",
+ "\n",
+ "#calculation\n",
+ "I=Vd/(R1+R2+Rl);\n",
+ "Vl1=I*Rl;\n",
+ "P=D*Vl*I;\n",
+ "Vh=I*R1;\n",
+ "Ph=D*Vh*I;\n",
+ "Vl=I*R2;\n",
+ "Pl=D*Vl*I;\n",
+ "Pic=Ph+Pl;\n",
+ "Pic=1.56;\n",
+ "Tj=150;\n",
+ "Ta=60;\n",
+ "Q=(Tj-Ta)/Pic-2.2;\n",
+ "\n",
+ "#result\n",
+ "print \"current is\",round(I,2), \"A\"\n",
+ "print \"load voltage is\",round(Vl1,2), \"V\"\n",
+ "print \"power is\",round(P,2), \"W\"\n",
+ "print \"high side voltage is\",round(Vh,2), \"V\"\n",
+ "print \"high side power is\",round(Ph,2), \"W\"\n",
+ "print \"low side voltage is\",round(Vl,2), \"V\"\n",
+ "print \"low side power is\",round(Pl,2), \"W\"\n",
+ "print \"IC power is\",round(Pic,2), \"W\"\n",
+ "print \"thermal resistance is\",round(Q,2), \"degreeC/W\""
+ ]
+ }
+ ],
+ "metadata": {
+ "kernelspec": {
+ "display_name": "Python 2",
+ "language": "python",
+ "name": "python2"
+ },
+ "language_info": {
+ "codemirror_mode": {
+ "name": "ipython",
+ "version": 2
+ },
+ "file_extension": ".py",
+ "mimetype": "text/x-python",
+ "name": "python",
+ "nbconvert_exporter": "python",
+ "pygments_lexer": "ipython2",
+ "version": "2.7.6"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter7.ipynb b/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter7.ipynb
new file mode 100644
index 00000000..45af7d98
--- /dev/null
+++ b/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter7.ipynb
@@ -0,0 +1,505 @@
+{
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 7: Switiching Power Supplies"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 7.1,Page 326"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 2,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "duty cycle is 0.25\n",
+ "average voltage is 3.0 V\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding duty cycle and average voltage\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "T=20.0;#time\n",
+ "Vp=12.0;#voltage\n",
+ "t=5.0;\n",
+ "\n",
+ "#calculation\n",
+ "D=t/T;\n",
+ "Vd=(D*Vp);\n",
+ "\n",
+ "#result\n",
+ "print \"duty cycle is\",round(D,3)\n",
+ "print \"average voltage is\",round(Vd,3), \"V\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 7.2,Page 238"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 3,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "duty cycle is 42.0 %\n",
+ "time is 10.0 microsec\n",
+ "on time is 4.167 microsec\n",
+ "ripple current is 133.636 mA\n",
+ "load current is 500.0 mA\n",
+ "peak inductor current is 566.818 mA\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding on time ripple,load,peak inductor current\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Vd=12.0;#voltage\n",
+ "Vl=5.0;#load voltage\n",
+ "Rl=10.0;#load resistance\n",
+ "f=100.0;#frequency\n",
+ "L=220.0;#inductor\n",
+ "\n",
+ "#calculation\n",
+ "D=Vl/Vd;\n",
+ "T=1/f;\n",
+ "t=D*T;\n",
+ "Vr=Vd-Vl;\n",
+ "I=Vr*round(t*10000)/10/L;\n",
+ "Il=Vl/Rl;\n",
+ "Ip=Il+I/2;\n",
+ "\n",
+ "#result\n",
+ "print \"duty cycle is\",round(D*100), \"%\"\n",
+ "print \"time is\",round(T*1000,3), \"microsec\"\n",
+ "print \"on time is\",round(t*10000,2)/10, \"microsec\"\n",
+ "print \"ripple current is\",round(I*1000,3),\"mA\"\n",
+ "print \"load current is\",round(Il*1000,3), \"mA\"\n",
+ "print \"peak inductor current is\",round(Ip*1000,3), \"mA\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 7.3,Page 335"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 4,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "rms current is 325.01 mA\n",
+ "by trapezium method\n",
+ "rms current is 324.04 mA\n",
+ "by rectangle method\n",
+ "\n",
+ " rectangle method gives good result than trapezium method\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding rms current\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Id=500.0;#load current\n",
+ "i=134;#mA\n",
+ "D=.42;#duty cycle\n",
+ "\n",
+ "#calculation\n",
+ "Ip=Id+i/2;\n",
+ "Im=Id-i/2;\n",
+ "I1=((D/3)*(Ip**2+Im*Ip+Im**2))**.5;\n",
+ "I2=D**.5*Id;\n",
+ "\n",
+ "#result\n",
+ "print \"rms current is\",round(I1,2), \"mA\"\n",
+ "print('by trapezium method')\n",
+ "print \"rms current is\",round(I2,2), \"mA\"\n",
+ "print('by rectangle method')\n",
+ "print '\\n rectangle method gives good result than trapezium method'"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 7.4,Page 336"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 5,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "voltage is 0.3 V\n",
+ "dissipated power is 63.0 mW\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding voltage and power\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "Vp=.3;#voltage\n",
+ "I=500.0;#current\n",
+ "D=.42;#duty cycle\n",
+ "T=150.0;#temperature\n",
+ "R=.6;#ohm\n",
+ "\n",
+ "#calculation\n",
+ "Vq=I*R;\n",
+ "Pq=D*Vq*I;\n",
+ "\n",
+ "#result\n",
+ "print \"voltage is\",round(Vq/1000,2), \"V\"\n",
+ "print \"dissipated power is\",round(Pq/1000,2), \"mW\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 7.5,Page 341"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 7,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "duty cycle is 42.0 %\n",
+ "time period is 6.667 microsec\n",
+ "on time is 2.778 microsec\n",
+ "load current is 500.0 mA\n",
+ "ripple current is 125.0 mA\n",
+ "inductor voltage is 7.0 V\n",
+ "inductor is 155.556 microH\n",
+ "inductor current is 562.5 mA\n",
+ "minimum capacitor current is 250.0 mA\n",
+ "minimum capacitor voltage is 18.0 V\n",
+ "Rf/Ri is 3.07\n",
+ "power of LM2595 is 0.33 W\n",
+ "thermal resistance is 210.998 degreeC/W\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding all componenets\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "R=10.0;#resistance\n",
+ "V1=5.0;#V\n",
+ "V2=12.0;#V\n",
+ "Ta=80.0;#degreeC\n",
+ "Tb=150.0;\n",
+ "f=150.0;#frequency\n",
+ "\n",
+ "#calculation\n",
+ "D=V1/V2;\n",
+ "T=1/f;\n",
+ "t=D*T;\n",
+ "Id=V1/R;\n",
+ "i=.25*Id;\n",
+ "Vl=V2-V1;\n",
+ "L=Vl*t/i;\n",
+ "Ip=Id+i/2;\n",
+ "Ic=Id/2;\n",
+ "Vc=1.5*V2;\n",
+ "K=V1/1.23-1;\n",
+ "P=.01*V2+D*Id*1;\n",
+ "Q=(Tb-Ta)/P-2.2;\n",
+ "\n",
+ "#result\n",
+ "print \"duty cycle is\",round(D*100), \"%\"\n",
+ "print \"time period is\",round(T*1000,3), \"microsec\"\n",
+ "print \"on time is\",round(t*1000,3), \"microsec\"\n",
+ "print \"load current is\",round(Id*1000,3), \"mA\"\n",
+ "print \"ripple current is\",round(i*1000,3), \"mA\"\n",
+ "print \"inductor voltage is\",round(Vl,2), \"V\"\n",
+ "print \"inductor is\",round(L*1000,3), \"microH\"\n",
+ "print \"inductor current is\",round(Ip*1000,2), \"mA\"\n",
+ "print \"minimum capacitor current is\",round(Ic*1000,2), \"mA\"\n",
+ "print \"minimum capacitor voltage is\",round(Vc,3), \"V\"\n",
+ "print \"Rf/Ri is\",round(K,2)\n",
+ "print \"power of LM2595 is\",round(P,2), \"W\"\n",
+ "print \"thermal resistance is\",round(Q,3), \"degreeC/W\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 7.6,Page 349"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 9,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "load power is 15.4 W\n",
+ "supply power is 17.11 W\n",
+ "dc current is 1.4 A\n",
+ "inductor current is 1.57 A\n",
+ "duty cycle is 0.45\n",
+ "inductor is 154.29 microH\n",
+ "transistor power is 352.8 mW\n",
+ "diode power is 385.0 mW\n",
+ "capacitor is 157.5 microF\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding different power,inductor current,inductor value\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "V1=12.0#V\n",
+ "V2=22.0;#V\n",
+ "I=.7;#A\n",
+ "f=100.0;#kHz\n",
+ "R=.4;#ohm\n",
+ "Vd=.5;\n",
+ "\n",
+ "#calculation\n",
+ "Pl=V2*I;\n",
+ "Ps=Pl/.9;\n",
+ "Id=round(Ps/V1*10)/10;\n",
+ "i=.25*Id;\n",
+ "Ip=Id+i/2;\n",
+ "D=round((1-V1/V2)*100)/100;\n",
+ "t=D/f;\n",
+ "L=V1*t/i;\n",
+ "Vp=Id*R;\n",
+ "Pq=D*Vp*Id;\n",
+ "Pd=(1-D)*.5*Id;\n",
+ "C=Id*t/2/20;\n",
+ "\n",
+ "#result\n",
+ "print \"load power is\",round(Pl,2), \"W\"\n",
+ "print \"supply power is\",round(Ps,2), \"W\"\n",
+ "print \"dc current is\",round(Id,2), \"A\"\n",
+ "print \"inductor current is\",round(Ip,2), \"A\"\n",
+ "print \"duty cycle is\",round(D,2)\n",
+ "print \"inductor is\",round(L*1000,2), \"microH\"\n",
+ "print \"transistor power is\",round(Pq*1000,2), \"mW\"\n",
+ "print \"diode power is\",round(Pd*100,2)*10, \"mW\"\n",
+ "print \"capacitor is\",round(C*1e6,2), \"microF\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 7.7,Page 355"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 11,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Rf/Ri is 16.886\n",
+ "pick Rf=22; Ri=1.3;\n",
+ "rms current is 1.4 A\n",
+ "switch power is 132.3 mW\n",
+ "IC power is 151.2 mW\n",
+ "total power is 283.5 mW\n",
+ "IC temperature is 98.43 degreeC\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding feedback resistor,power,current and temperature\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "V1=12.0;#V\n",
+ "V2=22.0;#V\n",
+ "I=.7;#A\n",
+ "Ta=80.0;#degreeC\n",
+ "Ps=17.1#supply power\n",
+ "\n",
+ "#calculation\n",
+ "K=V2/1.23-1;\n",
+ "Id=round(Ps/V1*10)/10;\n",
+ "D=round((1-(V1/V2))*100)/100;\n",
+ "Ir=D**.5*Id;\n",
+ "Ps=Ir**2*.15;\n",
+ "Pi=D*Id*V1/50;\n",
+ "P=Ps+Pi;\n",
+ "T=Ta+P*65;\n",
+ "\n",
+ "#result\n",
+ "print \"Rf/Ri is\",round(K,3)\n",
+ "print('pick Rf=22; Ri=1.3;')\n",
+ "print \"rms current is\",round(Id,2), \"A\"\n",
+ "print \"switch power is\",round(Ps*1000,2), \"mW\"\n",
+ "print \"IC power is\",round(Pi*1000,2), \"mW\"\n",
+ "print \"total power is\",round(P*1000,2), \"mW\"\n",
+ "print \"IC temperature is\",round(T,2), \"degreeC\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 7.8,Page 359"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 13,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "maximum voltage is 18.25 V\n",
+ "diode voltage is 20.0 V\n",
+ "duty cycle is 0.34\n",
+ "power delivered is 5.0 W\n",
+ "average current is 466.67 mA\n",
+ "mid primary current is 1.37 A\n",
+ "rms current is 800.33 mA\n",
+ "ramp current is 480.0 mA\n",
+ "maximum transistor current is 1.61 A\n",
+ "minimum transistor current is 1.13 A\n",
+ "diode peak current is 2.02 A\n",
+ "secondary rms current is 1.23 A\n",
+ "capacitor is 170.0 microF\n"
+ ]
+ }
+ ],
+ "source": [
+ "#designing circuit and finding circuit parameter\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "V1=12.0;#V\n",
+ "V2=5.0;#V\n",
+ "Il=1.0;#load current\n",
+ "T=10.0;#microsec\n",
+ "K=1.25;#Npri/Nsec\n",
+ "L=85.0;#microH\n",
+ "\n",
+ "#calculation\n",
+ "Vq=V1+V2*K;\n",
+ "Vd=V1*K+V2;\n",
+ "D=round((K*V2)*100/(V1+K*V2))/100;\n",
+ "Po=V2*Il;\n",
+ "Pi=round(Po/.09)/10;\n",
+ "Id=Pi/V1;\n",
+ "Im=Id/D;\n",
+ "Ir=(Im*D**.5);\n",
+ "i=V1*D*T/L;\n",
+ "IM=Im+.24;\n",
+ "Imin=Im-.24;\n",
+ "Ip=K*IM;\n",
+ "Imid=Il/(1-D);\n",
+ "Irms=Imid*(1-D)**.5;\n",
+ "C=D*Il*T/20;\n",
+ "\n",
+ "#result\n",
+ "print \"maximum voltage is\",round(Vq,2), \"V\"\n",
+ "print \"diode voltage is\",round(Vd,2), \"V\"\n",
+ "print \"duty cycle is\",round(D,2)\n",
+ "print \"power delivered is\",round(Po,2), \"W\"\n",
+ "print \"average current is\",round(Id*1000,2), \"mA\"\n",
+ "print \"mid primary current is\",round(Im,2), \"A\"\n",
+ "print \"rms current is\",round(Ir*1000,2),\"mA\"\n",
+ "print \"ramp current is\",round(i*1000,2), \"mA\"\n",
+ "print \"maximum transistor current is\",round(IM,2),\"A\"\n",
+ "print \"minimum transistor current is\",round(Imin,2),\"A\"\n",
+ "print \"diode peak current is\",round(Ip,2), \"A\"\n",
+ "print \"secondary rms current is\",round(Irms,2),\"A\"\n",
+ "print \"capacitor is\",round(C*1000,2), \"microF\"\n"
+ ]
+ }
+ ],
+ "metadata": {
+ "kernelspec": {
+ "display_name": "Python 2",
+ "language": "python",
+ "name": "python2"
+ },
+ "language_info": {
+ "codemirror_mode": {
+ "name": "ipython",
+ "version": 2
+ },
+ "file_extension": ".py",
+ "mimetype": "text/x-python",
+ "name": "python",
+ "nbconvert_exporter": "python",
+ "pygments_lexer": "ipython2",
+ "version": "2.7.6"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter8.ipynb b/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter8.ipynb
new file mode 100644
index 00000000..7cf6b405
--- /dev/null
+++ b/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter8.ipynb
@@ -0,0 +1,326 @@
+{
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 8: Thyristors"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 8.3,Page 397"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 1,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "inductance is 6.22 microH\n",
+ "load impedence at angle 90 degree is 0.00195 ohm\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding inductance,load impedence\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "V=220.0;#line voltage\n",
+ "f=50.0;#hertz\n",
+ "R=80.0;#load resistance\n",
+ "K=50.0;#di/dt\n",
+ "\n",
+ "#calculation\n",
+ "L=V*2**.5/K;\n",
+ "Z=2*pi*f*L;\n",
+ "\n",
+ "#result\n",
+ "print \"inductance is\",round(L,2),\"microH\"\n",
+ "print \"load impedence at angle 90 degree is\",round(Z*1e-6,5), \"ohm\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 8.4,Page 400"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 3,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "minimum value of capacitor is 0.067 micfoF\n",
+ "\n",
+ "choose C=.1 micoF\n",
+ "capacitor impedence at angle -90degree is 31.83 ohm\n",
+ "Load current in mA at an angle 90 degrees is 6.91\n",
+ "Potential drop in V at an angle 90 degrees is 0.55\n",
+ "Power dissipated is 3 mW\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding capacitor,current\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "V=220.0;#line voltage\n",
+ "f=50.0;#hertz\n",
+ "R=80.0;#load resistance\n",
+ "K=75.0;#dv/dt\n",
+ "Vd=400.0;#DRM voltage\n",
+ "\n",
+ "\n",
+ "#calculation\n",
+ "C=Vd/R/K;\n",
+ "C1=.1;\n",
+ "Z=1/(2*pi*f*C1);\n",
+ "Iload=V/1000/(-Z*cos(180*pi/180)+R*round(cos(90*pi/180)));\n",
+ "Vload=Iload/1000*R;\n",
+ "P=Vload*Iload;\n",
+ "\n",
+ "#result\n",
+ "print \"minimum value of capacitor is\",round(C,3), \"micfoF\"\n",
+ "print('\\nchoose C=.1 micoF')\n",
+ "print \"capacitor impedence at angle -90degree is\",round(Z*1000,2), \"ohm\"\n",
+ "print \"Load current in mA at an angle 90 degrees is\",round(Iload,2)\n",
+ "print \"Potential drop in V at an angle 90 degrees is\",round(Vload,2)\n",
+ "print \"Power dissipated is\",int(P), \"mW\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 8.5,Page 402"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 4,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "snubbing resistnce is 7.39 ohm\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding snubbing resistance\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "V=220;#line voltage\n",
+ "f=50;#hertz\n",
+ "R=80;#load resistance\n",
+ "I=46;#TSM current\n",
+ "\n",
+ "#calculation\n",
+ "Rs=V*2**.5/(I-V*2**.5/R);\n",
+ "\n",
+ "#result\n",
+ "print \"snubbing resistnce is\",round(Rs,2), \"ohm\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 8.6,Page 414"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 5,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "line period is 16.67 ms\n",
+ "half-cycle time is 8.333 ms\n",
+ "no. of cycles is 10.0\n",
+ "voltage for t1 is 54.0 V\n",
+ "power for t1 is 291.6 W\n",
+ "voltage for t2 is 100.0 V\n",
+ "voltage for t2 is 1000.0 W\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding voltage , power and cycles\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "R=10.0;#load\n",
+ "V=120.0;#rms voltage\n",
+ "f=60.0;#hertz\n",
+ "T=83.3;#ms\n",
+ "t1=15;#ms\n",
+ "t2=55;#ms\n",
+ "\n",
+ "#calculation\n",
+ "Tl=1/f;\n",
+ "Th=Tl/2;\n",
+ "C=round(T/Th/100)*100;\n",
+ "D1=.2;\n",
+ "V1=round(V*D1**.5);\n",
+ "P1=V1**2/R;\n",
+ "D2=.7;\n",
+ "V2=round(V*D2**.5);\n",
+ "P2=V2**2/R;\n",
+ "\n",
+ "#result\n",
+ "print \"line period is\",round(Tl*1000,2), \"ms\"\n",
+ "print \"half-cycle time is\",round(Th*1000,3), \"ms\"\n",
+ "print \"no. of cycles is\",C/1000\n",
+ "print \"voltage for t1 is\",round(V1,3), \"V\"\n",
+ "print \"power for t1 is\",round(P1,3), \"W\"\n",
+ "print \"voltage for t2 is\",round(V2,3), \"V\"\n",
+ "print \"voltage for t2 is\",round(P2,3), \"W\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 8.8,Page 427"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 7,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "average voltage is 42.0 V\n",
+ "dc voltage is 41.0 V\n",
+ "\n",
+ "the markers indicae Vp=163V Vave=41\n",
+ "full-wave rms voltage is 108.0 V\n",
+ "rms voltage is 108.0 V\n",
+ "\n",
+ "the markers indicate Vp=169V ;Vave=106V\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding dc volatge,average voltage,rms voltage\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "V=120.0;#line voltage\n",
+ "A=60.0;#degree\n",
+ "D=0.35;\n",
+ "\n",
+ "#calculation\n",
+ "Va=D*V;\n",
+ "Vd=V*2**.5*(cos(A*pi/180)+1)/2/pi;\n",
+ "Vr=.9*V;\n",
+ "Vrms=V*(2**.5)*(.5*(pi-1.047)+.25*sin(2*A*pi/180))**.5/pi**.5;\n",
+ "\n",
+ "#result\n",
+ "print \"average voltage is\",round(Va,3), \"V\"\n",
+ "print \"dc voltage is\",round(Vd), \"V\"\n",
+ "print('\\nthe markers indicae Vp=163V Vave=41')\n",
+ "print \"full-wave rms voltage is\",round(Vr), \"V\"\n",
+ "print \"rms voltage is\",round(Vrms), \"V\"\n",
+ "print('\\nthe markers indicate Vp=169V ;Vave=106V')\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 8.9,Page 430"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 8,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "rms voltage is 141.18 V\n",
+ "double checked value of rms voltage is 141.18 V\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding rms voltage and double checked rms voltage\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "V=220.0;#line voltage\n",
+ "P=1.3;#kW\n",
+ "R=15.0;#ohm\n",
+ "\n",
+ "#calculation\n",
+ "Vr=round((P*1000*R)**.5);\n",
+ "D=Vr/V;\n",
+ "Vr=V*2**.5*(.5*(pi-1.710)+sin(196*pi/180)/4)**.5/pi**.5;\n",
+ "\n",
+ "#result\n",
+ "print \"rms voltage is\",round(Vr,2), \"V\"\n",
+ "print \"double checked value of rms voltage is\",round(Vr,2), \"V\""
+ ]
+ }
+ ],
+ "metadata": {
+ "kernelspec": {
+ "display_name": "Python 2",
+ "language": "python",
+ "name": "python2"
+ },
+ "language_info": {
+ "codemirror_mode": {
+ "name": "ipython",
+ "version": 2
+ },
+ "file_extension": ".py",
+ "mimetype": "text/x-python",
+ "name": "python",
+ "nbconvert_exporter": "python",
+ "pygments_lexer": "ipython2",
+ "version": "2.7.6"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter9.ipynb b/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter9.ipynb
new file mode 100644
index 00000000..d50d131a
--- /dev/null
+++ b/Power_Electronics_Principles_and_Applications_by_Jacob/Chapter9.ipynb
@@ -0,0 +1,405 @@
+{
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 9: Power Conversation and Motor Drive Operations"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 9.1,Page 457"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 1,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "peak voltage is 37.6 V\n",
+ "load voltage is 35.7 V\n",
+ "ripple voltage is 3.96 V\n",
+ "approx. load voltage is 35.62 V\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding peak,load,ripple voltages\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "V=28.0;#V\n",
+ "C=4700.0;#microF\n",
+ "R=16.0;#load\n",
+ "f=120.0;#hertz\n",
+ "\n",
+ "#calculation\n",
+ "Vp=V*2**.5-2;\n",
+ "Vd=0.95*Vp;\n",
+ "Id=Vd/R;\n",
+ "v=Id/f/C;\n",
+ "#approximation\n",
+ "Vd1=Vp-v*1e6/2;\n",
+ "\n",
+ "#result\n",
+ "print \"peak voltage is\",round(Vp,2), \"V\"\n",
+ "print \"load voltage is\",round(Vd,1), \"V\"\n",
+ "print \"ripple voltage is\",round(v*1e6,2), \"V\"\n",
+ "print \"approx. load voltage is\",round(Vd1,2), \"V\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 9.2,Page 459"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 2,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Zth is 1.0 + 1.0 in ohm\n",
+ "inductor is 2.65 mH\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding inductor,Zth\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "V1=120.0;#pri voltage\n",
+ "V2=28.0;#sec voltage\n",
+ "I=2.0;#pri current\n",
+ "f=60.0;#Hz\n",
+ "Vth=28.8;#open voltage\n",
+ "V3=12.1;#pri-short voltage\n",
+ "Is=2.0;#short current at 45 degree\n",
+ "\n",
+ "#calculation\n",
+ "Zi=(V2*V3)/V1/Is*cos(45*pi/180);\n",
+ "Zj=(V2*V3)/V1/Is*sin(45*pi/180);\n",
+ "L=Zi/(2*pi*f);\n",
+ "\n",
+ "#result\n",
+ "print'Zth is',round(Zi),'+',round(Zj),'in ohm'\n",
+ "print \"inductor is\",round(L*1000,2), \"mH\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 9.4,Page 463"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 3,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "power factor is 0.32\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding power factor\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "I1=1.8;#current\n",
+ "R=16.0;#resistance\n",
+ "I2=5.7;#A\n",
+ "V=28.8;#Voltage\n",
+ "\n",
+ "#calculation\n",
+ "P=I1**2*R;\n",
+ "S=I2*V;\n",
+ "Pf=P/S;\n",
+ "\n",
+ "#result\n",
+ "print \"power factor is\",round(Pf,2)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 9.5, Page 468"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 2,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "aparrent power is 8.14 kVA\n",
+ "dissipated power is 7.84 kW\n",
+ "power factor is 0.96\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding power factor\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "I=22.6;#current\n",
+ "I2=28.00;\n",
+ "V=120.0;#Voltage\n",
+ "V2=280.0;\n",
+ "\n",
+ "#calculation\n",
+ "Pt=3*I*V;\n",
+ "Pl=I2*V2;\n",
+ "Pf=Pl/Pt;\n",
+ "\n",
+ "#result\n",
+ "print \"aparrent power is\",round(Pt/1000,2),\"kVA\"\n",
+ "print \"dissipated power is\",round(Pl/1000,2),\"kW\"\n",
+ "print \"power factor is\",round(Pl/Pt,2)\n",
+ "\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 9.6,Page 474"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 4,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "ratio is 0.72\n",
+ "firing angle is 58 degrees\n",
+ "dc voltage is 148.85 V\n",
+ "time delay is 2.69 ms\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding firing angle, time delay,Vd\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "V=208.0;#voltage\n",
+ "R=100.0;#load\n",
+ "Vd=150.0;#V\n",
+ "\n",
+ "#calculation\n",
+ "r=Vd/V;\n",
+ "a=58;#degree\n",
+ "Vd=3*2**.5*208*(cos(pi/3+a*pi/180)-cos(2*pi/3+a*pi/180))/pi;\n",
+ "t=a*16.7/360;\n",
+ "\n",
+ "#result\n",
+ "print \"ratio is\",round(r,2)\n",
+ "print('firing angle is 58 degrees');\n",
+ "print \"dc voltage is\",round(Vd,2), \"V\"\n",
+ "print \"time delay is\",round(t,2), \"ms\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 9.7,Page 480"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 5,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "max. current is 41.67 A\n",
+ "dissipated power is 8.68 W\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding maximum current and power dissipated\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "P=150.0;#power\n",
+ "V=8.0;#voltage\n",
+ "R=.01;#resistance\n",
+ "D=.5;#duty cycle\n",
+ "\n",
+ "#calculation\n",
+ "I=P/.9/D/V;\n",
+ "Ir=I*D**.5;\n",
+ "Pq=Ir**2*R;\n",
+ "\n",
+ "#result\n",
+ "print \"max. current is\",round(I,2), \"A\"\n",
+ "print \"dissipated power is\",round(Pq,2),\"W\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 9.8,Page 489"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 6,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "pwm fundamental frequency is 30.72 kHz\n",
+ "output voltage is 9.46 V\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding fundamental frequency and output voltage\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "f1=60.0;#frequency\n",
+ "V=150.0;#voltage\n",
+ "f2=31.0;#kHz\n",
+ "\n",
+ "#calculation\n",
+ "f3=f1*4;\n",
+ "Vo=V*10**(-4.2);\n",
+ "\n",
+ "#result\n",
+ "print \"pwm fundamental frequency is\",round(f3*2**7/1000,2), \"kHz\"\n",
+ "print \"output voltage is\",round(Vo*1000,2), \"V\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Example 9.9,Page 491"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 7,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "average voltage is 127.32 V\n",
+ "\n",
+ "Va-d @ 200Vin=4.2V\n",
+ "\n",
+ "\n",
+ "pick R1=47kohm\n",
+ "current through dividers is 2.62 mA\n",
+ "R2 is 1.6 kohm\n",
+ "capacitor is 27.01 microF\n"
+ ]
+ }
+ ],
+ "source": [
+ "#finding resistances,capacitor,average voltage\n",
+ "\n",
+ "#initialisation of variable\n",
+ "from math import pi,tan,sqrt,sin,cos,acos,atan\n",
+ "V=120.0;#load voltage\n",
+ "f=60.0;#Hz\n",
+ "Vp=200.0;#V\n",
+ "Vd=5.0;#V\n",
+ "\n",
+ "\n",
+ "#calculation\n",
+ "Vdc=2*Vp/pi;\n",
+ "Va=4.2;\n",
+ "R1=47.0;\n",
+ "I=(Vdc-Va)/R1;\n",
+ "R2=Va/I;\n",
+ "K=1.0/(1/R1+1/R2)# R1 \\\\ R2\n",
+ "C=1.0/2/pi/3.8/K;\n",
+ "\n",
+ "#result\n",
+ "print \"average voltage is\",round(Vdc,2), \"V\"\n",
+ "print('\\nVa-d @ 200Vin=4.2V')\n",
+ "print('\\n\\npick R1=47kohm')\n",
+ "print \"current through dividers is\",round(I,2), \"mA\"\n",
+ "print \"R2 is\",round(R2,2), \"kohm\"\n",
+ "print \"capacitor is\",round(C*1000,2), \"microF\""
+ ]
+ }
+ ],
+ "metadata": {
+ "kernelspec": {
+ "display_name": "Python 2",
+ "language": "python",
+ "name": "python2"
+ },
+ "language_info": {
+ "codemirror_mode": {
+ "name": "ipython",
+ "version": 2
+ },
+ "file_extension": ".py",
+ "mimetype": "text/x-python",
+ "name": "python",
+ "nbconvert_exporter": "python",
+ "pygments_lexer": "ipython2",
+ "version": "2.7.6"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Power_Electronics_Principles_and_Applications_by_Jacob/README.txt b/Power_Electronics_Principles_and_Applications_by_Jacob/README.txt
new file mode 100644
index 00000000..274ef8b6
--- /dev/null
+++ b/Power_Electronics_Principles_and_Applications_by_Jacob/README.txt
@@ -0,0 +1,10 @@
+Contributed By: Amit Kumar Das
+Course: btech
+College/Institute/Organization: ISM Dhanbad
+Department/Designation: Mining Machinery Engineering
+Book Title: Power Electronics Principles and Applications
+Author: Jacob
+Publisher: Cengage Learning, New Delhi
+Year of publication: 2009
+Isbn: 9780766823327
+Edition: 4 \ No newline at end of file
diff --git a/Power_Electronics_Principles_and_Applications_by_Jacob/screenshots/4.png b/Power_Electronics_Principles_and_Applications_by_Jacob/screenshots/4.png
new file mode 100644
index 00000000..5199042e
--- /dev/null
+++ b/Power_Electronics_Principles_and_Applications_by_Jacob/screenshots/4.png
Binary files differ
diff --git a/Power_Electronics_Principles_and_Applications_by_Jacob/screenshots/5.png b/Power_Electronics_Principles_and_Applications_by_Jacob/screenshots/5.png
new file mode 100644
index 00000000..76933e63
--- /dev/null
+++ b/Power_Electronics_Principles_and_Applications_by_Jacob/screenshots/5.png
Binary files differ
diff --git a/Power_Electronics_Principles_and_Applications_by_Jacob/screenshots/6.png b/Power_Electronics_Principles_and_Applications_by_Jacob/screenshots/6.png
new file mode 100644
index 00000000..607dd0a1
--- /dev/null
+++ b/Power_Electronics_Principles_and_Applications_by_Jacob/screenshots/6.png
Binary files differ