{
 "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
}