{
 "metadata": {
  "name": ""
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "<h1>Chapter 5: Point Source and Their Arrays<h1>"
     ]
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "<h3>Example 5-6.1, Page number: 90<h3>"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "import scipy.integrate\n",
      "\n",
      "#Variable declaration\n",
      "def integrand(theta, phi):\n",
      "    return (math.cos(theta)*math.sin(theta))\n",
      "            #Integrand (unitless)\n",
      "Um = 1      #Maximum radiation intensity (unitless)\n",
      "\n",
      "#Calculation\n",
      "P = scipy.integrate.dblquad(integrand, 0, 2*math.pi,\n",
      "                             lambda x: 0, lambda x: math.pi/2)\n",
      "                        #Total power radiated (relative to Um)\n",
      "D = (4*math.pi)/P[0]    #Directivity (unitless)\n",
      "\n",
      "#Result\n",
      "print \"The directivity is \", round(D)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The directivity is  4.0\n"
       ]
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "<h3>Example 5-6.2, Page number: 91<h3>"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "import scipy.integrate\n",
      "\n",
      "#Variable declaration\n",
      "def integrand(theta, phi):\n",
      "    return (math.cos(theta)*math.sin(theta))\n",
      "            #Integrand (unitless)\n",
      "Um = 1      #Maximum radiation intensity (unitless)\n",
      "\n",
      "#Calculation\n",
      "P = scipy.integrate.dblquad(integrand, 0, 2*math.pi,\n",
      "                             lambda x: 0, lambda x: math.pi/2)\n",
      "                        #Total power radiated (relative to Um)\n",
      "D = (4*math.pi)/(2*P[0])    #Directivity (unitless)\n",
      "\n",
      "#Result\n",
      "print \"The directivity is \", round(D)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The directivity is  2.0\n"
       ]
      }
     ],
     "prompt_number": 2
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "<h3>Example 5-6.3, Page number: 91<h3>"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math, scipy.integrate\n",
      "\n",
      "#Variable declaration\n",
      "def integrand(theta, phi):\n",
      "    return (math.sin(theta)**2)\n",
      "                #Integrand (unitless)\n",
      "Um = 1          #Maximum radiation intensity (unitless)\n",
      "\n",
      "#Calculation\n",
      "P = scipy.integrate.dblquad(integrand, 0, 2*math.pi,\n",
      "                             lambda x: 0, lambda x: math.pi)\n",
      "                    #Total radiated power (relative to Um)\n",
      "D = 4*math.pi/P[0]  #Directivity (unitless)\n",
      "\n",
      "#Result             \n",
      "print \"The directivity is\", round(D,2)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The directivity is 1.27\n"
       ]
      }
     ],
     "prompt_number": 3
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "<h3>Example 5-6.4, Page number: 91<h3>"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math, scipy.integrate\n",
      "\n",
      "#Variable declaration\n",
      "def integrand(theta, phi):\n",
      "    return (math.sin(theta)**3)\n",
      "                #Integrand (unitless)\n",
      "Um = 1          #Maximum radiation intensity (unitless)\n",
      "\n",
      "#Calculation\n",
      "P = scipy.integrate.dblquad(integrand, 0, 2*math.pi,\n",
      "                             lambda x: 0, lambda x: math.pi)\n",
      "                    #Total radiated power (relative to Um)\n",
      "D = 4*math.pi/P[0]  #Directivity (unitless)\n",
      "\n",
      "#Result                \n",
      "print \"The directivity is\", round(D,2)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The directivity is 1.5\n"
       ]
      }
     ],
     "prompt_number": 4
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "<h3>Example 5-6.5, Page number: 92<h3>"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math, scipy.integrate\n",
      "\n",
      "#Variable declaration\n",
      "def integrand(theta, phi):\n",
      "    return (math.sin(theta)*math.cos(theta)**2)\n",
      "                #Integrand (unitless)\n",
      "Um = 1          #Maximum radiation intensity (unitless)\n",
      "\n",
      "#Calculation\n",
      "P = scipy.integrate.dblquad(integrand, 0, 2*math.pi,\n",
      "                             lambda x: 0, lambda x: math.pi/2)\n",
      "                    #Total radiated power (relative to Um)\n",
      "D = 4*math.pi/P[0]  #Directivity (unitless)\n",
      "\n",
      "#Result                \n",
      "print \"The directivity is\", round(D,2)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The directivity is 6.0\n"
       ]
      }
     ],
     "prompt_number": 5
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "<h3>Example 5-6.6, Page number:93<h3>"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "\n",
      "#Variable declaration\n",
      "lobes = [0.25,0.37,0.46,0.12,0.07]  #Normalized power of lobes (unitless)\n",
      "\n",
      "#Calculation\n",
      "ohm_a = 0                   #Beam area (sr)\n",
      "sum_lobes = 0               #Sum of all lobes (unitless)\n",
      "for i in lobes:\n",
      "    ohm_a += 2*math.pi*(math.pi/36)*(i)\n",
      "    sum_lobes += i\n",
      "\n",
      "D = 4*math.pi/ohm_a         #Directivity (unitless)\n",
      "D_db = 10*math.log10(D)     #Directivity (in dBi)\n",
      "e_m = lobes[0]/sum_lobes    #Beam efficiency (unitless)\n",
      "\n",
      "#Result\n",
      "print \"The directivity is\", round(D), \"or\", round(D_db,1), \"dBi\"\n",
      "print \"The beam efficiency is\", round(e_m, 2)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The directivity is 18.0 or 12.6 dBi\n",
        "The beam efficiency is 0.2\n"
       ]
      }
     ],
     "prompt_number": 6
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "<h3>Example 5-21.1, Page number: 146<h3>"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "\n",
      "#Variable declaration\n",
      "a = 25              #Height of vertical conducting wall (m)\n",
      "r = 100             #Distance to the receiver (m)\n",
      "wave_lt = 10e-2     #Transmitter dimension (m)\n",
      "\n",
      "#Calculation\n",
      "k = math.sqrt(2/(r*wave_lt))    #contant (unitless)\n",
      "S_av = (r*wave_lt)/(4*(math.pi**2)*(a**2))  #Relative signal level (unitless)\n",
      "S_av_db = 10*math.log10(S_av)   #Signal level (in db)\n",
      "\n",
      "#Result\n",
      "print \"The signal level at the receiver is\", round(S_av,5), \"or\", round(S_av_db), \"dB\""
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The signal level at the receiver is 0.00041 or -34.0 dB\n"
       ]
      }
     ],
     "prompt_number": 1
    }
   ],
   "metadata": {}
  }
 ]
}