diff options
Diffstat (limited to 'Antenna_and_Wave_Propagation_by_k.k._sharma')
9 files changed, 3740 insertions, 0 deletions
diff --git a/Antenna_and_Wave_Propagation_by_k.k._sharma/README.txt b/Antenna_and_Wave_Propagation_by_k.k._sharma/README.txt new file mode 100644 index 00000000..62f55e2a --- /dev/null +++ b/Antenna_and_Wave_Propagation_by_k.k._sharma/README.txt @@ -0,0 +1,10 @@ +Contributed By: Ashish Kumar +Course: btech +College/Institute/Organization: COER +Department/Designation: ME +Book Title: Antenna and Wave Propagation +Author: k.k. sharma +Publisher: Shubham Publications, Delhi +Year of publication: 2008 +Isbn: 81-903721-5-7 +Edition: 1
\ No newline at end of file diff --git a/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter1.ipynb b/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter1.ipynb new file mode 100755 index 00000000..01df0fac --- /dev/null +++ b/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter1.ipynb @@ -0,0 +1,745 @@ +{ + "metadata": { + "name": "", + "signature": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "chapter 01 : Antenna Principles" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 1.1 : page 1.42" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#given data :\n", + "E=4.0 #in V/m\n", + "Eta=120*pi #constant\n", + "#Formula : E/H=Eta\n", + "H=E/Eta #in A/m\n", + "print \"Strength of magnetic field in free space = %0.4f A/m \" %H" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Strength of magnetic field in free space = 0.0106 A/m \n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 1.2 : page 1.42" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#given data :\n", + "H=5.2 #in mA/m\n", + "Eta=120*pi #constant\n", + "#Formula : E/H=Eta\n", + "E=H*10**-3*Eta #in V/m\n", + "print \"Strength of Electric field in free space =\",round(E),\"V/m\" " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Strength of Electric field in free space = 2.0 V/m\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 1.3 : page 1.42" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#given data :\n", + "I=20.0 #in A\n", + "Rr=100.0 #in Ohm\n", + "#Formula : Wr=I**2*R\n", + "Wr=I**2*Rr #in W\n", + "print \"Radiated power = %0.f kW \" %(Wr/1000) " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Radiated power = 40 kW \n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 1.4 : page 1.42" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt\n", + "#given data :\n", + "W=625.0 #in KW\n", + "r=30.0 #in Km\n", + "Erms=sqrt(90*W*1000)/(r*1000) #in V/m\n", + "print \"Strength of Electric field at 30Km away = %0.f mV/m \" %(Erms*1000) " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Strength of Electric field at 30Km away = 250 mV/m \n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 1.5 : page 1.43" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#given data :\n", + "le=10.0 #in m\n", + "Irms=450.0 #in A\n", + "f=50.0 #in KHz\n", + "R=1.5 #in Ohm\n", + "lamda=300.0/(f/1000) #in m\n", + "Rr=160*(pi)**2*(le/lamda)**2 #in Ohm\n", + "print \"Radiation resistance = %0.5f ohm\" %Rr\n", + "Wr=Irms**2*Rr #in W\n", + "print \"Radiated power = %0.2f Watts \" %Wr \n", + "Eta=(Rr/(Rr+R))*100 #efficiency in %\n", + "print \"Efficiency of antenna = %0.2f %%\" %Eta\n", + "# Ans in the textbook is not accurate." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Radiation resistance = 0.00439 ohm\n", + "Radiated power = 888.26 Watts \n", + "Efficiency of antenna = 0.29 %\n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 1.6 : page 1.43" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#given data :\n", + "le=50.0 #in m\n", + "f=100.0 #in MHz\n", + "lamda=300.0/(f) #in m\n", + "Rr=(160*(pi)**2)*(le/lamda)**2 #in Ohm\n", + "print \"Radiation Resistance = %0.2f Mohm \" %(Rr/10**6) \n", + "#Note : Answer in the book is wrong" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Radiation Resistance = 0.44 Mohm \n" + ] + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 1.7 : page 1.44" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#given data :\n", + "l=30 #in m\n", + "Irms=20 #in A\n", + "f=1 #in MHz\n", + "r=10 #in Km\n", + "r=r*1000 #in m\n", + "le=2*l/pi #in m\n", + "lamda=300/(f) #in m\n", + "Erms=120*pi*le*Irms/(lamda*r) #in V/m\n", + "print \"Field strength at 10Km distance = %0.2e V/m \" %Erms \n", + "#Note : Answer in the book is wrong" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Field strength at 10Km distance = 4.80e-02 V/m \n" + ] + } + ], + "prompt_number": 14 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 1.8 : page 1.44" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#given data :\n", + "Rl=1.0 #in ohm\n", + "#Formula : Rr=80*pi**2*(l/lamda)**2\n", + "#Given l=lamda/10\n", + "#l/lamda=1/10\n", + "Rr=80*pi**2*(1.0/10)**2 #in Ohm\n", + "print \"Radiation resistance = %0.2f Ohm \" %(Rr) \n", + "Eta=Rr/(Rr+Rl) #Unitless\n", + "print \"Antenna Efficiency = %0.2f %% \" %(Eta*100) " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Radiation resistance = 7.90 Ohm \n", + "Antenna Efficiency = 88.76 % \n" + ] + } + ], + "prompt_number": 15 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 1.9 : page 1.44" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt\n", + "#given data :\n", + "r=100 #in Km\n", + "W=100 #in KW\n", + "Erms=sqrt(90*W*1000)/(r*1000) #in V/m\n", + "print \"Strength of Electric Field = %0.2f V/m \" %Erms " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Strength of Electric Field = 0.03 V/m \n" + ] + } + ], + "prompt_number": 16 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 1.10 : page 1.44" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#given data :\n", + "le=200.0 #in m\n", + "Irms=200 #in A\n", + "f=300 #in KHz\n", + "r=10 #in Km\n", + "c=3*10**8 #speed of light i m/s\n", + "lamda=c/(f*1000) #in m\n", + "Erms=120*pi*le*Irms/(lamda*r*10**3) #in V/m\n", + "print \"Field strength at 10Km distance = %0.4f V/m\" %(Erms) \n", + "Rr=(160*(pi)**2)*(le/lamda)**2 #in Ohm\n", + "W=Irms**2*Rr #in Watts\n", + "print \"Radiated Power = %0.2f MW \" %(W/10**6) \n", + "#Note : Answer is wrong in the book. Unit of answer in the book is written mW instead of MW by mistake." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Field strength at 10Km distance = 1.5080 V/m\n", + "Radiated Power = 2.53 MW \n" + ] + } + ], + "prompt_number": 21 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 1.11 : page 1.45" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#given data :\n", + "#Formula : Rr=80*pi**2*(l/lamda)**2\n", + "#Given l=lamda/60\n", + "#l/lamda=1/60\n", + "Rr=80*pi**2*(1.0/60)**2 #in Ohm\n", + "print \"Radiation resistance = %0.3f Ohm \" %Rr " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Radiation resistance = 0.219 Ohm \n" + ] + } + ], + "prompt_number": 23 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 1.12 : page 1.45" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#given data :\n", + "r=10.0 #in Km\n", + "Erms=10.0 #in mV/m\n", + "r1=20.0 #in Km\n", + "#Formula : Erms=sqrt(90*W)/r #in V/m\n", + "#Let swrt(90*W)=a\n", + "a=Erms*r \n", + "Erms1=a/r1 #in mV/m\n", + "print \"Field strength at 20Km distance = %0.f mV/m \" %Erms1 " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Field strength at 20Km distance = 5 mV/m \n" + ] + } + ], + "prompt_number": 24 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 1.13 : page 1.45" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#given data :\n", + "r=1.0 #in Km\n", + "r=1*10**3 #in m\n", + "l=1.0 #in m\n", + "Irms=10.0 #in A\n", + "f=5.0 #in MHz\n", + "c=3*10**8 #speed of light i m/s\n", + "lamda=c/(f*10**6) #in m\n", + "le=2*l/pi #in m\n", + "Erms=120*pi*le*Irms/(lamda*r) #in V/m\n", + "print \"Field strength at 10Km distance = %0.4f V/m \" %Erms\n", + "#Note : Answer in the book is wrong. Mistake during value putting." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Field strength at 10Km distance = 0.0400 V/m \n" + ] + } + ], + "prompt_number": 26 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 1.14 : page 1.46" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#given data :\n", + "Irms=30.0 #in A\n", + "f=1.0 #in MHz\n", + "Erms=10.0 #in mV/m\n", + "Erms=Erms*10**-3 #in V/m\n", + "r=50.0 #in Km\n", + "r=r*10**3 #in m\n", + "c=3*10**8 #speed of light i m/s\n", + "lamda=c/(f*10**6) #in m\n", + "le=Erms*lamda*r/(120*pi*Irms) #in m\n", + "print \"Effetive height of Antenna = %0.4f meter \" %le " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Effetive height of Antenna = 13.2629 meter \n" + ] + } + ], + "prompt_number": 28 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 1.15 : page 1.46" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from sympy import symbols, solve\n", + "I, r = symbols('I r')\n", + "E = 10*I/r # V/m\n", + "\n", + "#given data :\n", + "Erms_sqr = E**2\n", + "Wt = (Erms_sqr*r**2)/30 \n", + "Rr = Wt/I**2 # ohm\n", + "print \"Radiation resistance = %0.2f Ohm \" %float(Rr) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Radiation resistance = 3.33 Ohm \n" + ] + } + ], + "prompt_number": 41 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 1.16 : page 1.46" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#given data :\n", + "lamda=300/(50*10**-6) #in m\n", + "r=round(lamda)/(2*pi) #in m\n", + "print \"Distance = %0.2e meter \" %r \n", + "#Note : Answer in the book is wrong." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Distance = 9.55e+05 meter \n" + ] + } + ], + "prompt_number": 43 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 1.17 : page 1.47" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt\n", + "#given data :\n", + "r=2 #in Km\n", + "r=r*10**3 #in m\n", + "Wt=1 #in KW\n", + "Wt=Wt*10**3 #in Watt\n", + "Erms=sqrt(30*Wt)/r #in V/m\n", + "print \"Field strength at 2Km distance = %0.3f mV/m \" %(Erms*10**3) " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Field strength at 2Km distance = 86.603 mV/m \n" + ] + } + ], + "prompt_number": 45 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 1.18 : page 1.47" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#given data :\n", + "f=20.0 #in MHz\n", + "f=f*10**6 #in Hz\n", + "le=100.0 #in m\n", + "c=3*10**8 #speed of light in m/s\n", + "lamda=c/f #in m\n", + "Rr=160*(pi*le/lamda)**2 #in ohm\n", + "print \"Radiation Resistance = %0.1f kohm \" %(Rr/1000) " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Radiation Resistance = 70.2 kohm \n" + ] + } + ], + "prompt_number": 47 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 1.19 : page 1.47" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt, pi\n", + "#given data :\n", + "P=10.0 #in W/m**2\n", + "f=40.0 #in MHz\n", + "f=f*10**6 #in Hz\n", + "mu_r=4.0 #constant\n", + "epsilon_r=5 #constant\n", + "#Velocity of propagation\n", + "#formula : v=(1/sqrt(mu_o*epsilon_o))*(1/sqrt(mu_r*epsilon_r)) #in m/s\n", + "#1/sqrt(mu_o*epsilon_o)=c=speed of light=3*10**8 m/s\n", + "c=3*10**8 #speed of light in m/s\n", + "v=c*(1.0/sqrt(mu_r*epsilon_r)) #in m/s\n", + "print \"Velocity of propagation = %0.1e m/s \" %v \n", + "#Wavelength\n", + "lamda=v/f #in meter\n", + "print \"Wavelength = %0.2f m \" %lamda \n", + "#rms electric field\n", + "#Formula : E=P*sqrt(mu_o/epsilon_o)*sqrt(mu_r/epsilon_r) #in V/m\n", + "E=sqrt(1200*pi*sqrt(4.0/5)) #in V/m\n", + "Erms=sqrt(E**2/sqrt(2)) #in V/m\n", + "print \"rms Electric Field = %0.2f V/m\" %Erms \n", + "#Impedence of medium\n", + "Eta=(sqrt(2)*Erms)**2/P #in Ohm\n", + "print \"Impedence of medium = %0.2f ohm \" %Eta " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Velocity of propagation = 6.7e+07 m/s \n", + "Wavelength = 1.68 m \n", + "rms Electric Field = 48.83 V/m\n", + "Impedence of medium = 476.86 ohm \n" + ] + } + ], + "prompt_number": 50 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 1.20 : page 1.48" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "from sympy import symbols, solve, N\n", + "#given data :\n", + "#Hfi = (Im*dlsin(theta)/(4*pi))*[cos(omega*t1)/r-omega*sin(omega*t1)/(c*r)]\n", + "lamda, r = symbols('lamda r')\n", + "#expr = 200.0/r**2-2*pi*f/c\n", + "expr = 200.0/r**2-2*pi/lamda/r # putting f/c = lamda\n", + "r = solve(expr, r)\n", + "print \"r =\",N(r[0],4)" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "r = 31.83*lamda\n" + ] + } + ], + "prompt_number": 3 + } + ], + "metadata": {} + } + ] +} diff --git a/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter10.ipynb b/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter10.ipynb new file mode 100755 index 00000000..4dee9c4a --- /dev/null +++ b/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter10.ipynb @@ -0,0 +1,432 @@ +{ + "metadata": { + "name": "", + "signature": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "chapter 10 : Sky wave propagation - The ionospheric waves" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 10.1 : page 10-19" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt\n", + "#given data :\n", + "H=500 #in km\n", + "n=0.8 #in m\n", + "f_muf=10 #in MHz\n", + "f_muf=f_muf*10**6 #in Hz\n", + "f=10 #in MHz\n", + "f=f*10**6 #in Hz\n", + "# Formula : n=sqrt(1-81*N/f**2)\n", + "Nmax=(1-n**2)*f**2/81 #in Hz \n", + "fc=9*sqrt(Nmax) #in Hz\n", + "Dskip=2*H*sqrt((f_muf/fc)**2-1) #in Km\n", + "print \"Assuming the earth is flat the range = %0.2f km\" %Dskip\n", + "#Note : Answer in the book is wrong." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Assuming the earth is flat the range = 1333.33 km\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 10.2 : page 10-19" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#given data :\n", + "n=0.8 #in m\n", + "H=500 #in km\n", + "a=6370 #in km\n", + "D=1349.07 #in Km\n", + "f_muf=10 #in MHz\n", + "f_muf=f_muf*10**6 #in Hz\n", + "f=10 #in MHz\n", + "f=f*10**6 #in Hz\n", + "# Formula : n=sqrt(1-81*N/f**2)\n", + "Nmax=(1-n**2)*f**2/81 #in Hz \n", + "fc=9*sqrt(Nmax) #in Hz\n", + "# Formula : f_muf/fc=sqrt(D**2/(4*(H+D**2/(8*a))))+1\n", + "D1=2*(H+D**2/(8*a))*sqrt((f_muf/fc)**2-1) #in Km\n", + "Dskip=2*H*sqrt((f_muf/fc)**2-1) #in Km\n", + "print \"Assuming the earth is curved the ground range = %0.2f km\"% D1\n", + "# Answer wrong in the textbook." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Assuming the earth is curved the ground range = 1428.57 km\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 10.3 : page 10-20" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt\n", + "#given data :\n", + "Nmax=2.48*10**6 #in cm**-3\n", + "Nmax=2.48*10**6*10**-6 #in m**-3\n", + "fc=9*sqrt(Nmax) #in MHz\n", + "print \"Critical frequency = %0.2f MHz \" %fc " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Critical frequency = 14.17 MHz \n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 10.4 : page 10-20" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt\n", + "#given data :\n", + "H=200 #in Km\n", + "D=4000 #in Km\n", + "fc=5 #in MHz\n", + "f_muf=fc*sqrt(1+(D/(2*H))**2) #in MHz\n", + "print \"MUF for the given path = %0.2f MHz \" %f_muf\n", + "#Note : Answer in the book is wrong." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "MUF for the given path = 50.25 MHz \n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 10.5 : page 10-20" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt\n", + "#given data :\n", + "#For F1 layer :\n", + "print \"For F1 layer :\" \n", + "Nmax=2.3*10**6 #in cm**3\n", + "Nmax=2.3*10**6*10**-6 #in m**3\n", + "fc=9*sqrt(Nmax) #in MHz\n", + "print \"Critical frequency = %0.2f MHz \" %fc \n", + "\n", + "#For F2 layer :\n", + "print \"For F2 layer :\" \n", + "Nmax=3.5*10**6 #in cm**3\n", + "Nmax=3.5*10**6*10**-6 #in m**3\n", + "fc=9*sqrt(Nmax) #in MHz\n", + "print \"Critical frequency = %0.2f MHz\" %fc\n", + "\n", + "#For F3 layer :\n", + "print \"For F3 layer :\" \n", + "Nmax=1.7*10**6 #in cm**3\n", + "Nmax=1.7*10**6*10**-6 #in m**3\n", + "fc=9*sqrt(Nmax) #in MHz\n", + "print \"Critical frequency = %0.2f MHz \" %fc " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "For F1 layer :\n", + "Critical frequency = 13.65 MHz \n", + "For F2 layer :\n", + "Critical frequency = 16.84 MHz\n", + "For F3 layer :\n", + "Critical frequency = 11.73 MHz \n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 10.6 : page 10-21" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt\n", + "#given data :\n", + "n=0.7 #refractive index\n", + "N=400 #in cm**-3\n", + "#Formula : n=sqrt(1-81*N/f**2)\n", + "f=sqrt(81*N/(1-n**2)) #in KHz\n", + "print \"Frequency of wave propagation = %0.2f kHz\" %f\n", + "#Note : Unit of Answer in the book is MHz. It is written by mistake. It is accurately calculated by scilab in KHz. " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Frequency of wave propagation = 252.05 kHz\n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 10.7 : page 10-21" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt\n", + "#given data :\n", + "HT=169.0 #in meter\n", + "HR=20.0 #in meter\n", + "d=4.12*(sqrt(HT)+sqrt(HR)) #in Km\n", + "print \"Maximum distance = %0.2f km \" %d \n", + "r_dash=(4/3)*6370/1000 #in Km\n", + "RadioHorizon=sqrt(2*r_dash*HT) #in Km\n", + "print \"Radio Horizon = %0.2f km \" %RadioHorizon\n", + "# Answe wrong in thetextbook." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Maximum distance = 71.99 km \n", + "Radio Horizon = 45.03 km \n" + ] + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 10.8 : page 10-21" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import tan , pi, asin, cos\n", + "H=200 #in Km\n", + "Beta=20 #in Degree\n", + "a=6370 #in Km\n", + "D_flat=2*H/tan(Beta*pi/180) #in Km\n", + "print \"If earth assumed to be flat transmission path distance = %0.2f km\" %D_flat\n", + "D_curved=2*a*(90*pi/180-Beta*pi/180)-asin(a*cos(Beta*pi/180)/(a+H))\n", + "print \"If earth assumed to be curved transmission path distance = %0.2f \"%D_curved\n", + "# Answe wrong in thetextbook." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "If earth assumed to be flat transmission path distance = 1098.99 km\n", + "If earth assumed to be curved transmission path distance = 15563.70 \n" + ] + } + ], + "prompt_number": 17 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 10.9 : page 10-22" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import acos\n", + "#given data :\n", + "R=6370 #in Km\n", + "hm=400 #in Km\n", + "#Formula : d=2*R*Q=2*R*acos(R/(R+hm))\n", + "d=2*R*acos(R/(R+hm)) #in Km\n", + "print \"Maximum Range in a single range transmission = %0.2f km \" %d \n", + "# Answe wrong in thetextbook." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Maximum Range in a single range transmission = 20011.95 km \n" + ] + } + ], + "prompt_number": 21 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 10.10 : page 10-22" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt\n", + "#given data :\n", + "n=0.6 #refractive index\n", + "N=4.23*10**4 #in m**-3\n", + "#Formula : n=sqrt(1-81*N/f**2)\n", + "f=sqrt(81*N/(1-n**2)) #in Hz\n", + "print \"Frequency of wave propagation = %0.3f kHz\" %(f/1000)" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Frequency of wave propagation = 2.314 kHz\n" + ] + } + ], + "prompt_number": 23 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 10.11 : page 10-23" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt\n", + "#given data :\n", + "n=0.8 #refractive index\n", + "N=500 #in cm**-3\n", + "#Formula : n=sqrt(1-81*N/f**2)\n", + "f=sqrt(81*N/(1-n**2)) #in KHz\n", + "print \"Frequency of wave propagation = %0.2f kHz\" %f " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Frequency of wave propagation = 335.41 kHz\n" + ] + } + ], + "prompt_number": 24 + } + ], + "metadata": {} + } + ] +} diff --git a/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter3.ipynb b/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter3.ipynb new file mode 100755 index 00000000..2278567d --- /dev/null +++ b/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter3.ipynb @@ -0,0 +1,698 @@ +{ + "metadata": { + "name": "", + "signature": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "chapter 03 : Antenna Terminology" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 3.1 : page 3.42" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#given data :\n", + "E=10.0 #in V/m\n", + "ETA_o=120.0*pi #Constant\n", + "H=E/ETA_o #in A/m\n", + "print \"The Magnetic Field Strength = %0.4f A/m \" %H " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The Magnetic Field Strength = 0.0265 A/m \n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 3.2 : page 3.42" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt\n", + "#given data :\n", + "W=25.0 #in KW\n", + "W=W*10**3 #in W\n", + "r=3 #in Km\n", + "r=r*10**3 #in m\n", + "Erms=sqrt(90*W)/r #in V/m\n", + "print \"Field strength at receiver = %0.2f V/m \" %Erms " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Field strength at receiver = 0.50 V/m \n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 3.3 : page 3.42" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#given data :\n", + "le=125 #in m\n", + "Irms=5 #in A\n", + "lamda=1.25 #in Km\n", + "lamda=lamda*10**3 #in m\n", + "Rl=10 #in Ohm\n", + "#radiation Resistance\n", + "Rr=(80*pi**2)*(le/lamda)**2 #in Ohm\n", + "Rr=round(Rr) #in Ohm : approx\n", + "print \"Radiation resistance = %0.2f Ohm \" %Rr \n", + "#Power radiated\n", + "W=(Irms**2)*Rr #in \n", + "print \"Power radiated = %0.2f W \" %W\n", + "#Antenna efficiency \n", + "ETA=Rr/(Rr+Rl)\n", + "print \"Antenna efficiency = %0.2f %% \" %(ETA*100) " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Radiation resistance = 8.00 Ohm \n", + "Power radiated = 200.00 W \n", + "Antenna efficiency = 44.44 % \n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 3.4 : page 3.43" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import cos, pi, sin\n", + "#given data :\n", + "r=1 #in Km\n", + "r=r*10**3 #in m\n", + "I=0.5 #in A\n", + "#For theta = 45 degree\n", + "theta=45 #in degree\n", + "E=(60*I/r)*((cos(pi*cos(theta*pi/180)/2))/sin(theta*pi/180)) \n", + "print \"E-Field for 45 degree angle = %0.2f mV/m \" %(E*10**3) \n", + "ETA_o=120*pi #constant\n", + "H=E/ETA_o #in A/m\n", + "print \"H-Field for 45 degree angle = %0.5f mV/m \" %(H*10**3) \n", + "\n", + "#For theta = 90 degree\n", + "theta=90 #in degree\n", + "E=(60*I/r)*((cos(pi*cos(theta*pi/180)/2))/sin(theta*pi/180)) \n", + "print \"E-Field for 90 degree angle = %0.2f mV/m \" %(E*10**3) \n", + "ETA_o=120*pi #constant\n", + "H=E/ETA_o #in A/m\n", + "print \"H-Field for 90 degree angle = %0.4f mV/m \" %(H*10**3) " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "E-Field for 45 degree angle = 18.84 mV/m \n", + "H-Field for 45 degree angle = 0.04997 mV/m \n", + "E-Field for 90 degree angle = 30.00 mV/m \n", + "H-Field for 90 degree angle = 0.0796 mV/m \n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 3.5 : page 3.44" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#given data :\n", + "#l=lambda/10 meter\n", + "#Assume pi**2 = 10\n", + "Rl=2.0 #in Ohm\n", + "#Rr=80*pi**2*(dl/lambda)**2\n", + "Rr=80*10*(1.0/10)**2 #in Ohm\n", + "print \"Radiation Resistance = %0.2f Ohm\" %(Rr)\n", + "ETA=Rr/(Rr+Rl) #in Ohm\n", + "print \"Efficiency = %0.2f %%\" %(ETA*100)" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Radiation Resistance = 8.00 Ohm\n", + "Efficiency = 80.00 %\n" + ] + } + ], + "prompt_number": 11 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 3.6 : page 3.44" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#given data :\n", + "#l=lambda/15 meter\n", + "#Assume pi**2 = 10\n", + "Rl=2.0 #in Ohm\n", + "#Gain : \n", + "Gain=5.33/4 #Unitless\n", + "#Directivity\n", + "Rr=80*10*(1.0/15)**2 #in Ohm\n", + "ETA=Rr/(Rr+Rl) #Unitless\n", + "Directivity=Gain/ETA #unitless\n", + "#Beam solid angle \n", + "BSA=4.0*pi/Directivity #in steradian\n", + "print \"Directivity = %0.4f \" %Directivity \n", + "print \"Gain = %0.2f \"%Gain \n", + "#Effective aperture\n", + "print \"Effective aperture = \" ,\n", + "print round((Gain/(4*pi)),3),\"lambda**2\" \n", + "print \"Beam Solid Angle = %0.2f steradian \"%BSA \n", + "Rr=80*10*(1.0/15)**2 #in Ohm\n", + "print \"Radiation Resistance = %0.2f Ohm \" %Rr \n", + "print \"Pt =\",120*10/225,\"I**2\" \n", + "print \"Pr = 4*I**2\" " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Directivity = 2.0820 \n", + "Gain = 1.33 \n", + "Effective aperture = 0.106 lambda**2\n", + "Beam Solid Angle = 6.04 steradian \n", + "Radiation Resistance = 3.56 Ohm \n", + "Pt = 5 I**2\n", + "Pr = 4*I**2\n" + ] + } + ], + "prompt_number": 22 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 3.7 : page 3.45" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#given data :\n", + "D=30.0 #in m\n", + "k=0.55 #illumination efficiency\n", + "f=4.0 #in GHz\n", + "f=f*10**9 #in Hz\n", + "c=3*10**8 #speed of light in m/s\n", + "lamda=c/f #in Meter\n", + "r=D/2 #in m\n", + "A=pi*(r**2) #in m**2\n", + "G=(4*pi/lamda**2)*k*A #Unitless\n", + "print \"Gain = %0.5e\"%G\n", + "HPBW=70*lamda/D #in Degree\n", + "print \"HPBW = %0.3f Degree \" % HPBW\n", + "BWFN=2*70*lamda/D #in Degree\n", + "print \"BWFN = %0.2f Degree \" %BWFN " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Gain = 8.68525e+05\n", + "HPBW = 0.175 Degree \n", + "BWFN = 0.35 Degree \n" + ] + } + ], + "prompt_number": 25 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 3.8 : page 3.46" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#given data :\n", + "Rl=20.0 #in Ohm\n", + "Rr=100.0 #in Ohm\n", + "Gp=25.0 #power gain \n", + "ETA=Rr/(Rr+Rl) #Unitless\n", + "D=Gp/ETA #unitless\n", + "print \"Directivity = %0.2f\" %D" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Directivity = 30.00\n" + ] + } + ], + "prompt_number": 30 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 3.9 : page 3.46" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#given data :\n", + "lamda=10 #in m\n", + "D=80 #unitless\n", + "Aem=D*lamda**2/(4*pi) #in m**2\n", + "print \"Maximum effective aperture = %0.2f m^2\" %Aem" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Maximum effective aperture = 636.62 m^2\n" + ] + } + ], + "prompt_number": 31 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 3.10 : page 3.47" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import log10\n", + "#given data :\n", + "P1=30 #in KW\n", + "P1=P1*1000 #in W\n", + "P2=5000 #in W\n", + "Gdb=10*log10(P1/P2) #unitless\n", + "print \"Front to back ratio, Gdb =\",round(Gdb ,3)" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Front to back ratio, Gdb = 7.782\n" + ] + } + ], + "prompt_number": 34 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 3.11 : page 3.47" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#given data :\n", + "f=10 #in GHz\n", + "f=f*10**9 #in Hz\n", + "Gt=40 #in dB\n", + "Gr=40 #in dB\n", + "print \"Gain = Gt = Gr =\",Gt ,\"dB\"" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Gain = Gt = Gr = 40 dB\n" + ] + } + ], + "prompt_number": 36 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 3.12 : page 3.47" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#given data :\n", + "L=10 #in m\n", + "f=1.5 #in MHz\n", + "f=f*10**6 #in Hz\n", + "X=350 #in Ohm\n", + "Q=100 #Coil parameter\n", + "c=3*10**8 #speed of light in m/s\n", + "lamda=c/f #in Meter\n", + "l_eff=2*L/2 #in m\n", + "Re=2*X/Q #in Ohm\n", + "Rr=40*pi**2*(l_eff/lamda)**2 #in hm\n", + "Gd=(3/2)*(lamda**2/(4*pi)) #unitless\n", + "ETA=Rr/(Rr+Re) #Efficiency unitless\n", + "Gp=Gd*ETA ##unitless\n", + "print \"Antenna Efficiency = %0.1f %%\" %(ETA*100)\n", + "print \"Power gain = %0.2f \" %(Gp)\n", + "#Note : Answer of Gp is wrong in the book." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Antenna Efficiency = 12.4 %\n", + "Power gain = 393.34 \n" + ] + } + ], + "prompt_number": 38 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 3.13 : page 3.48" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#given data :\n", + "delf=600.0 #in KHz\n", + "fr=50 #in MHz\n", + "Q=(fr*10**6)/(delf*10**3) #unitless\n", + "print \"Quality Factor = %0.2f \" %(Q) " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Quality Factor = 83.33 \n" + ] + } + ], + "prompt_number": 41 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 3.14 : page 3.48" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#given data :\n", + "OmegaA=4.0*pi #For isotropic Antenna\n", + "D=4.0*pi/OmegaA #Directivity : Unitless\n", + "print \"Directivity of Isotropic Antenna = %0.2f\" %D " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Directivity of Isotropic Antenna = 1.00\n" + ] + } + ], + "prompt_number": 42 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 3.15 : page 3.48" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from sympy import symbols, N\n", + "lamda = symbols('lamda')\n", + "#given data :\n", + "D=500.0 #Directivity : Unitless\n", + "Aem = D*lamda**2/(4*pi)\n", + "print \"Aem =\",N(Aem,4)" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Aem = 39.79*lamda**2\n" + ] + } + ], + "prompt_number": 50 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 3.16 : page 3.48" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#given data\n", + "Fn_dB=1.1 #in dB\n", + "Fn=10**(Fn_dB/10) #unitless\n", + "To=290 #in Kelvin\n", + "Te=To*(Fn-1) #in Kelvin\n", + "print \"Effective Noise Temperature = %0.2f degree Kelvin \" %Te " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Effective Noise Temperature = 83.59 degree Kelvin \n" + ] + } + ], + "prompt_number": 52 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 3.19 : page 3.50" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi, log10\n", + "#given data\n", + "D=6.0 #in meter\n", + "f=10.0 #in GHz\n", + "f=f*10**9 #in Hz\n", + "Aactual=pi*D**2/4 #in m**2\n", + "Ae=0.6*Aactual #in m**2\n", + "c=3*10**8 #speed of light in m/s\n", + "lamda=c/f #in Meter\n", + "G=4*pi*Ae/lamda**2 #Unitless\n", + "Gdb=10*log10(G) #gain in dB\n", + "BWFN=140*lamda/D #in degree\n", + "print \"Gain = %0.1f \" %G \n", + "print \"Gain = %0.2f dB \" %Gdb \n", + "print \"Beamwidth = %0.2f degree \" %BWFN \n", + "print \"Capture Area = %0.2f m**2 \" %Ae \n", + "#Note : Answer in the book is not accurate." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Gain = 236870.5 \n", + "Gain = 53.75 dB \n", + "Beamwidth = 0.70 degree \n", + "Capture Area = 16.96 m**2 \n" + ] + } + ], + "prompt_number": 54 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 3.20 : page 3.50" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#given data\n", + "Gdb=44 #gain in dB\n", + "G=10**(Gdb/10) #gain unitless\n", + "OmegaB=4*pi/G #n steradian\n", + "THETA3db=sqrt(4*OmegaB/pi) #in Radian\n", + "print \"Beamwidth THETA3db = %0.4f degree \" %THETA3db \n", + "#Note : Answer in the book is not accurate." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Beamwidth THETA3db = 0.0400 degree \n" + ] + } + ], + "prompt_number": 56 + } + ], + "metadata": {} + } + ] +} diff --git a/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter4.ipynb b/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter4.ipynb new file mode 100755 index 00000000..2a68737d --- /dev/null +++ b/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter4.ipynb @@ -0,0 +1,481 @@ +{ + "metadata": { + "name": "", + "signature": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "chapter 04 : Antenna Arrays" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 4.3 : page 4.67" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import acos, pi, cos, sqrt, degrees\n", + "#given data :\n", + "from sympy import symbols\n", + "lamda, Ep = symbols('lamda Ep')\n", + "d = 3.0/2*lamda\n", + "beta = 2*pi/lamda\n", + "delta = 0 # for broad side array\n", + "theta = pi/2 # for maxima\n", + "si = 3*pi/2*cos(theta)\n", + "E0 = Ep/sqrt(2) # at half power beam width\n", + "#Ep = 2*E0*cos(si/2)\n", + "#it leads to cos(3*pi/2*cos(theta))=1/sqrt(2)\n", + "theta=acos(acos(1/sqrt(2))/(3*pi/2)) # radian\n", + "theta = degrees(theta) # degree\n", + "HPBW=2*(90-theta) #in degree\n", + "print \"(i) HPBW = %0.f degree \" %HPBW \n", + "# 2nd part is wrong. Some mistake in question as cos(theta) = 13/12 which >1 not possible" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "(i) HPBW = 19 degree \n" + ] + } + ], + "prompt_number": 18 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 4.4 : page 4.68" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import log10\n", + "from sympy import symbols\n", + "lamda = symbols('lamda')\n", + "#given data :\n", + "n=10 #no. of elements\n", + "d=lamda/4 #separation in meter\n", + "D=2*n/(lamda/d)\n", + "Ddb=10*log10(D) #in db\n", + "print \"For broad side array D = %0.2f db \" %Ddb \n", + "D=4*n/(lamda/d)\n", + "Ddb=10*log10(D) #in db\n", + "print \"For end fire array D = %0.2f db \" %Ddb " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "For broad side array D = 6.99 db \n", + "For end fire array D = 10.00 db \n" + ] + } + ], + "prompt_number": 22 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 4.5 : page 4.68" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt, pi\n", + "#given data :\n", + "from sympy import symbols, N\n", + "lamda = symbols('lamda')\n", + "delta=-90 #in degree\n", + "#Formula : HPBW=57.3/(sqrt(L/(2*lambda))) in Degree\n", + "n=20 #no. of point sources\n", + "d=lamda/4 #in meter\n", + "L=(n-1)*d\n", + "HPBW=57.3/(sqrt(L/lamda/2)) # in Degree\n", + "print \"HPBW = %0.2f Degree \" %HPBW \n", + "D=4*L/lamda #Directivity\n", + "print \"Directivity = %0.2f \" %D \n", + "Ae = D*lamda**2/4/pi\n", + "print \"Effective aperture : Ae =\",N(Ae,3)\n", + "Omega=4*pi/D #in steradian\n", + "print \"Beam Solid Angle : Omega =\",round(Omega,2),\"Steradian\" \n", + "#Note : Answer in the book is not accurate." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "HPBW = 37.18 Degree \n", + "Directivity = 19.00 \n", + "Effective aperture : Ae = 1.51*lamda**2\n", + "Beam Solid Angle : Omega = 0.66 Steradian\n" + ] + } + ], + "prompt_number": 32 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 4.6 : page 4.69" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import degrees\n", + "#given data :\n", + "n=8 #no. of half wave dipoles\n", + "lamda=100 #in cm\n", + "lamda=lamda*10**-2 #in m\n", + "d=50 #in cm\n", + "d=d*10**-2 #in m\n", + "I=0.5 #in A\n", + "Rr=73 #in Ohm\n", + "Pr=n*I**2*Rr #in Watts\n", + "print \"Pr = %0.2f Watts \" %Pr \n", + "BWFN=2*lamda/(n*d) #in radian\n", + "HPBW=BWFN/2 #in radian\n", + "print \"HPBW = %0.2f degree\" % degrees(HPBW)" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Pr = 146.00 Watts \n", + "HPBW = 14.32 degree\n" + ] + } + ], + "prompt_number": 34 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 4.7 : page 4.69" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import log10\n", + "from sympy import symbols, N\n", + "lamda = symbols('lamda')\n", + "#given data :\n", + "n=10 #no. of elements\n", + "d=lamda/4 #separation in meter\n", + "Do=1.789*4*n*d/lamda\n", + "Dodb=10*log10(Do) #in db\n", + "print \"Do = %0.2f db\" %(Dodb) " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Do = 12.53 db\n" + ] + } + ], + "prompt_number": 36 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 4.13 : page 4.74" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sin, pi\n", + "#given data :\n", + "n=8 #no. of elements\n", + "BWFN=45 #in degree\n", + "theta=45 #in degree\n", + "f=40 #in MHz\n", + "f=f*10**6 #in Hz\n", + "#Formula : theta=2*asin(2*pi/(n*dr))\n", + "dr=(2*pi/n)/sin((theta/2)*(pi/180)) #\n", + "c=3*10**8 #speed of light in m/s\n", + "lamda=c/f #in m\n", + "d=dr*lamda/(2*pi) #in m\n", + "print \"Distance = %0.2f m \" %d " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Distance = 2.34 m \n" + ] + } + ], + "prompt_number": 37 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 4.14 : page 4.74" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#given data :\n", + "n=10 #no. of elements\n", + "from sympy import symbols, N\n", + "lamda = symbols('lamda')\n", + "#given : \n", + "d=lamda/4 #in m\n", + "Llambda=n*d/lamda\n", + "D=2*Llambda #in unitless \n", + "print \"Directivity of broadside uniform array = %0.2f \" %D" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Directivity of broadside uniform array = 5.00 \n" + ] + } + ], + "prompt_number": 40 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 4.16 : page 4.75" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "n=4 \n", + "from sympy import symbols, N, cos, solve\n", + "lamda, theta = symbols('lamda theta')\n", + "d=lamda/2\n", + "delta = pi/3\n", + "dr = 2*pi*d/lamda\n", + "# Peaks\n", + "si = pi*cos(theta)+pi/3\n", + "theta = solve(si, theta) # radian\n", + "theta = degrees(theta[0]) # degree\n", + "print \"Peaks : theta =\",round(theta,2),\"degree\"\n", + "# Nulls\n", + "print \"Nulls : \"\n", + "for k in range(0,3):\n", + " theta = degrees(acos(-1.0/3+k/2.0))\n", + " print \"k =\",k,\", theta =\",round(theta,2),\"degree\"\n", + "print \"Side lobes :\"\n", + "for k in range(0,3):\n", + " theta = degrees(acos(-1.0/3+(2*k+1)/4.0))\n", + " print \"k =\",k,\", theta =\",round(theta,2),\"degree\"\n", + "# Ans in the textbook is not accurate." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Peaks : theta = 109.47 degree\n", + "Nulls : \n", + "k = 0 , theta = 109.47 degree\n", + "k = 1 , theta = 80.41 degree\n", + "k = 2 , theta = 48.19 degree\n", + "Side lobes :\n", + "k = 0 , theta = 94.78 degree\n", + "k = 1 , theta = 65.38 degree\n", + "k = 2 , theta = 23.56 degree\n" + ] + } + ], + "prompt_number": 74 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 4.17 : page 4.76" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import cos, sin, pi\n", + "#given data :\n", + "MainBeamwidth=45 #in degree\n", + "thetaN=MainBeamwidth/2 #in degree\n", + "thetaN=thetaN*pi/180 #in radian\n", + "m=5 #no. of elements\n", + "#given : d=lambda/2 in meter\n", + "x=cos(pi/(2*(m-1))) \n", + "xo=x/cos((pi/2)*sin(thetaN)) #unitless\n", + "print \"E5=ao*z+a1*(2*z**2-1)+a2*(8*z**4-8*z**2+1)\" \n", + "print \"We Know that : z=x/xo, E5=T4*xo\" \n", + "print \"ao=a1*(2*(x/xo)**2-1)+a2*[8*(x/xo)**4-8*(x/xo)**2+1]=8*x**4-8*x**2+1\" \n", + "print \"By comparing the term we have : \" \n", + "print \"a2=xo**4 a1=4*a2-4*xo**2 ao=1+a1-a2 \"\n", + "a2=xo**4 \n", + "a1=4*a2-4*xo**2 \n", + "ao=1+a1-a2 \n", + "print \"And therefore the 5 elements array is given by : \" \n", + "print (a2),\" \",(a1),\" \",(2*ao),\" \",(a1),\" \",a2" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "E5=ao*z+a1*(2*z**2-1)+a2*(8*z**4-8*z**2+1)\n", + "We Know that : z=x/xo, E5=T4*xo\n", + "ao=a1*(2*(x/xo)**2-1)+a2*[8*(x/xo)**4-8*(x/xo)**2+1]=8*x**4-8*x**2+1\n", + "By comparing the term we have : \n", + "a2=xo**4 a1=4*a2-4*xo**2 ao=1+a1-a2 \n", + "And therefore the 5 elements array is given by : \n", + "1.5218051188 1.15276185234 1.26191346708 1.15276185234 1.5218051188\n" + ] + } + ], + "prompt_number": 82 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 4.18 : page 4.77" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from __future__ import division\n", + "#given data :\n", + "#Side lobe level below main lobe\n", + "print \"Side lobe level below main lobe : \"\n", + "SideLobe=20 #in dB\n", + "r=10**(SideLobe/20) #\n", + "print \"r=\",r \n", + "#No. of elements are 5, n=5\n", + "print \"No. of elements are 5, n=5 :\" \n", + "print \"Tchebyscheff polynomials of degree (n-1) is\" \n", + "print \"5-1=4\" \n", + "print \"T4(xo)=r\" \n", + "print \"8*xo**4-8*xo**2+1=10\" \n", + "print \"By using alternate formula, we get\" \n", + "m=4 \n", + "r=10 \n", + "xo=(1/2)*((r+sqrt(r**2-1))**(1/m)+(r-sqrt(r**2-1))**(1/m))\n", + "print \"xo=\" ,xo\n", + "print \"E5=T4(xo)\"\n", + "print \"E5=ao*z+a1*(2*z**2-1)+a2*(8*z**4-8*z**2+1)\" \n", + "print \"We Know that : z=x/xo, E5=T4*xo\" \n", + "print \"ao=a1*(2*(x/xo)**2-1)+a2*[8*(x/xo)**4-8*(x/xo)**2+1]=8*x**4-8*x**2+1\" \n", + "print \"By comparing the term we have : \" \n", + "print \"a2=xo**4 a1=4*a2-4*xo**2 ao=1+a1-a2 \"\n", + "a2=xo**4 \n", + "a1=4*a2-4*xo**2 \n", + "ao=1+a1-a2 \n", + "print \"And therefore the 5 elements array is given by : \" \n", + "print round(a2,3),\" \",round(a1,3),\" \",round(2*ao,3),\" \",round(a1,3),\" \",round(a2,3)\n", + "# Ans in the textbook are not accurate." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Side lobe level below main lobe : \n", + "r= 10.0\n", + "No. of elements are 5, n=5 :\n", + "Tchebyscheff polynomials of degree (n-1) is\n", + "5-1=4\n", + "T4(xo)=r\n", + "8*xo**4-8*xo**2+1=10\n", + "By using alternate formula, we get\n", + "xo= 1.29329190052\n", + "E5=T4(xo)\n", + "E5=ao*z+a1*(2*z**2-1)+a2*(8*z**4-8*z**2+1)\n", + "We Know that : z=x/xo, E5=T4*xo\n", + "ao=a1*(2*(x/xo)**2-1)+a2*[8*(x/xo)**4-8*(x/xo)**2+1]=8*x**4-8*x**2+1\n", + "By comparing the term we have : \n", + "a2=xo**4 a1=4*a2-4*xo**2 ao=1+a1-a2 \n", + "And therefore the 5 elements array is given by : \n", + "2.798 4.5 5.405 4.5 2.798\n" + ] + } + ], + "prompt_number": 86 + } + ], + "metadata": {} + } + ] +} diff --git a/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter5.ipynb b/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter5.ipynb new file mode 100755 index 00000000..300f89b1 --- /dev/null +++ b/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter5.ipynb @@ -0,0 +1,225 @@ +{ + "metadata": { + "name": "", + "signature": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "chapter 05 : Practical Antennas-I" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 5.1 : page 5.57" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#For Single Turn:\n", + "from sympy import symbols, sqrt\n", + "lamda = symbols('lamda')\n", + "a=lamda/25\n", + "A=pi*pow(a,2)\n", + "Rr = (A/lamda**2)**2*31171.2\n", + "print \"radiation Resistance =\",round(Rr,4),\"Ohm for single turn \"\n", + "\n", + "#For Eight Turn:\n", + "N=8 #no. of turns\n", + "Rr=Rr*N**2 #in Ohm\n", + "print \"radiation Resistance = %0.2f Ohm for Eight turn \" %Rr" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "radiation Resistance = 0.7876 Ohm for single turn \n", + "radiation Resistance = 50.40 Ohm for Eight turn \n" + ] + } + ], + "prompt_number": 26 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 5.2 : page 5.58" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi, acos, sqrt \n", + "#Given data :\n", + "f=20.0 #in MHz\n", + "N=15.0 #No. of turns\n", + "A=2.0 #in m**2\n", + "Vrms=200.0 #in uV\n", + "theta=acos(1) #in radian\n", + "mu_o=4*pi*10**-7 #in H/m\n", + "#Formula : Vm=2*pi*f*mu_o*H*A*N\n", + "Vm=Vrms*sqrt(2) #in uV\n", + "H=(Vm*10**-6)/(2.0*pi*f*10**6*mu_o*A*N) #in A/m\n", + "print \"Peak Value of magnetic feld intensity = %0.3e mA/m \" %(H*1000) \n", + "#Note : Answer in the book is wrong." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Peak Value of magnetic feld intensity = 5.970e-05 mA/m \n" + ] + } + ], + "prompt_number": 32 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 5.3 : page 5.58" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt\n", + "#Given data :\n", + "f=20 #in MHz\n", + "f=f*10**6 #in Hz\n", + "Wmax=25 #in mW/m**2\n", + "A=10.0 #in m**2\n", + "c=3*10**8 #speed of light in m/s\n", + "lamda=c/f #in meter\n", + "Rr=31171.2*(A/lamda**2)**2 #iin Ohm\n", + "#Formula : Wmax=V**2/(4*Rr)\n", + "V=sqrt(Wmax*10**-3*4*Rr) #in Volts\n", + "print \"Maximum emf in the loop = %0.3f Volts \"%V " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Maximum emf in the loop = 2.481 Volts \n" + ] + } + ], + "prompt_number": 35 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 5.4 : page 5.59" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#Given data :\n", + "N=20.0 #turns\n", + "D=1.0 #in meter\n", + "r=D/2 #in meter\n", + "E=200*10**-6 #in V/m\n", + "L=50*10**-6 #in H\n", + "R=2.0 #in Ohm\n", + "f=1.5 #in MHz\n", + "f=f*10**6 #in Hz\n", + "c=3*10**8 #speed of light in m/s\n", + "lamda=c/f #in meter\n", + "A=pi*r**2 #in m**2\n", + "Vrms=2*pi*E*A*N/lamda #in Volts\n", + "Q=2*pi*f*L/R #unitless\n", + "Vc_rms=Vrms*Q #in Volts\n", + "print \"Voltage across the capacitor = %0.2f mV\" %(Vc_rms*1000) \n", + "#Note : Answer in the book is wrong." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Voltage across the capacitor = 23.25 mV\n" + ] + } + ], + "prompt_number": 36 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 5.5 : page 5.59" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi, cos\n", + "#Given data :\n", + "N=100 #No. of turns\n", + "A=2 #in m**2\n", + "f=10 #in MHz\n", + "f=f*10**6 #in Hz\n", + "Q=150 #Quality factor\n", + "c=3*10**8 #speed of light in m/s\n", + "lamda=c/f #in meter\n", + "Erms=10*10**-6 #in V/m\n", + "theta=60 #in degree\n", + "Vrms=2*pi*Erms*A*N*cos(theta*pi/180)/lamda \n", + "Vin=Vrms*Q #in Volts\n", + "print \"Voltage to the receiver = %0.1f mV \" %(Vin*1000) \n", + "#Note : Answer in the book is wrong." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Voltage to the receiver = 31.4 mV \n" + ] + } + ], + "prompt_number": 40 + } + ], + "metadata": {} + } + ] +} diff --git a/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter6.ipynb b/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter6.ipynb new file mode 100755 index 00000000..62ee756d --- /dev/null +++ b/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter6.ipynb @@ -0,0 +1,524 @@ +{ + "metadata": { + "name": "", + "signature": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "chapter 06 : Practical Antennas - II" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 6.1 : page 6.39" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt\n", + "n=20 #no. of turns\n", + "#Clamda=lamda\n", + "#Slamda=lamda/4\n", + "#HPBW : \n", + "# HPBW=52/(Clamda*sqrt(n*Slamda))\n", + "#Putting values below :\n", + "Clamda=1 #in Meter\n", + "Slamda=1.0/4 #in Meter\n", + "HPBW=52.0/(Clamda*sqrt(n*Slamda)) #in degree\n", + "print \"HPBW = %0.2f degree \" %HPBW \n", + "#Axial Ratio\n", + "Aratio=(2*n+1)/2 #unitless\n", + "print \"Axial Ratio = %0.2f \"%Aratio \n", + "#Gain\n", + "D=12*Clamda**2*n*Slamda #unitless\n", + "print \"Gain = %0.2f \"%D " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "HPBW = 23.26 degree \n", + "Axial Ratio = 20.00 \n", + "Gain = 60.00 \n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 6.2 : page 6.39" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#Part (a): Given data :\n", + "n=20 #no. of turns\n", + "Slamda=0.472 #in meter\n", + "D=12*n*Slamda #in meter\n", + "from sympy import symbols, N, sqrt\n", + "lamda = symbols('lamda', real =True)\n", + "Ae=(lamda**2/(4*pi))*D\n", + "d=(sqrt(Ae))\n", + "print \"Part (a) : d=\",N(d,1)\n", + "print \"Part (b) : With a space of 3*lamda the total effective area : \" \n", + "Ae=9.02*lamda**2*4\n", + "D=4*pi*Ae/lamda**2\n", + "print \"\\t D = %0.2f\" %D" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Part (a) : d= 3.0*Abs(lamda)\n", + "Part (b) : With a space of 3*lamda the total effective area : \n", + "\t D = 453.39\n" + ] + } + ], + "prompt_number": 15 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 6.3 : page 6.40" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import log10, ceil\n", + "#from 7dBi gain graph the data obtained is given below :\n", + "K=1.2 #Scale constant\n", + "alfa=1.5 #Apex angle in degree\n", + "Slamda=0.15 \n", + "print \"K**n=F or n=logF/logK\" \n", + "F=4 \n", + "n=log10(F)/log10(K) \n", + "n=ceil(n) \n", + "nplus1=n+1 \n", + "print \"Apex Angle = %0.2f degree \" %alfa \n", + "print \"Sale constant = %0.2f\" %K \n", + "print \"No. of elements = %d \" %n " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "K**n=F or n=logF/logK\n", + "Apex Angle = 1.50 degree \n", + "Sale constant = 1.20\n", + "No. of elements = 8 \n" + ] + } + ], + "prompt_number": 16 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 6.4 : page 6.40" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import log10\n", + "#Given data :\n", + "#d=10*lamda\n", + "print \"d=10*lamda\" \n", + "print \"Power Gain : G=6*(d/lamda)**2\" \n", + "print \"Putting value of d, we get G=6*10**2\"\n", + "G=6*10**2 #unitless\n", + "print \"Power gain = %0.2f \" %G \n", + "G_dB=10*log10(G) #in dB\n", + "print \"Power Gain = %0.1f dB \" %G_dB " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "d=10*lamda\n", + "Power Gain : G=6*(d/lamda)**2\n", + "Putting value of d, we get G=6*10**2\n", + "Power gain = 600.00 \n", + "Power Gain = 27.8 dB \n" + ] + } + ], + "prompt_number": 18 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 6.5 : page 6.40" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import log10\n", + "#Given Data:\n", + "f=10.0 #in GHz\n", + "f=f*10**9 #in Hz\n", + "BWFN=10 #in degree\n", + "c=3*10**8 #Speed of light in m/s\n", + "lamda=c/f #in meter\n", + "#Part (a):\n", + "d=140*lamda/BWFN #in meter\n", + "print \"Diameter of a parabolic Antenna = %0.2f m\" %d\n", + "#Part (b):\n", + "HPBW=58.0*lamda/d #in degree\n", + "print \"3-dB Beamwidth = %0.2f degree \" %HPBW \n", + "#Part (c):\n", + "Gp=6*(d/lamda)**2 #gain \n", + "Gp_dB=10*log10(Gp) #in dB\n", + "print \"Power Gain = %0.2f dB \" %Gp_dB " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Diameter of a parabolic Antenna = 0.42 m\n", + "3-dB Beamwidth = 4.14 degree \n", + "Power Gain = 30.70 dB \n" + ] + } + ], + "prompt_number": 19 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 6.6 : page 6.41" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import log10\n", + "#Given Data:\n", + "f=1430.0 #in MHz\n", + "f=f*10**6 #in Hz\n", + "d=64 #in meter\n", + "c=3*10**8 #Speed of light in m/s\n", + "lamda=c/f #in meter\n", + "#Part (a):\n", + "HPBW=70*lamda/d #in degree\n", + "print \"HPBW = %0.2f degree \" %HPBW \n", + "#Part (b):\n", + "BWFN=140*lamda/d #in degree\n", + "print \"BWFN = %0.2f degree \" %BWFN \n", + "#Part (c):\n", + "Gp=6*(d/lamda)**2 #gain \n", + "Gp_dB=10*log10(Gp) #in dB\n", + "print \"Power Gain = %0.f dB \" %Gp_dB " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "HPBW = 0.23 degree \n", + "BWFN = 0.46 degree \n", + "Power Gain = 57 dB \n" + ] + } + ], + "prompt_number": 22 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 6.7 : page 6.42" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt\n", + "#Given Data:\n", + "f=15.0 #in GHz\n", + "f=f*10**9 #in Hz\n", + "Gp_dB=75.0 #in dB\n", + "c=3*10**8 #Speed of light in m/s\n", + "lamda=c/f #in meter\n", + "#Formula : Gp=9.87*(d/lamda)**2\n", + "#Formula : Gp_dB=10log10(Gp)\n", + "d=sqrt((10**(Gp_dB/10))*lamda**2/9.87) #in meter\n", + "print \"Diameter of a parabolic reflector = %0.2f m\" %d" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Diameter of a parabolic reflector = 35.80 m\n" + ] + } + ], + "prompt_number": 24 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 6.8 : page 6.42" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Given Data:\n", + "f=5000.0 #in MHz\n", + "f=f*10**6 #in Hz\n", + "d=10 #in feet\n", + "d=d*0.3048 #in meter\n", + "c=3*10**8 #Speed of light in m/s\n", + "lamda=c/f #in meter\n", + "r=2*d**2.0/lamda #in meter\n", + "print \"Minimum distance between primary and secondary antenna = %0.f m\" %r\n", + "# Ans wrong in the textbook" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Minimum distance between primary and secondary antenna = 310 m\n" + ] + } + ], + "prompt_number": 27 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 6.9 : page 6.42" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Given Data:\n", + "K=55.0 #Aperture Efficiency in %\n", + "K=K/100 #Aperture Efficiency\n", + "f=15.0 #in GHz\n", + "f=f*10**9 #in Hz\n", + "c=3*10**8 #Speed of light in m/s\n", + "lamda=c/f #in meter\n", + "G_dB=30 #in dB\n", + "G=10**(G_dB/10) #Gain unitless\n", + "#Formula : G=4*pi*K*A/lamda**2\n", + "A=(G*lamda**2)/(4*pi*K) #in m**2\n", + "print \"Diameter of parabolic reflector = %0.3f m**2\" %A \n", + "#Part (b)\n", + "d=sqrt(4.0*A/pi) #in meter\n", + "HPBW=70*lamda/d #in degree\n", + "print \"HPBW = %0.2f degree \" %HPBW \n", + "#Part (c)\n", + "BWFN=140*lamda/d #in Degree\n", + "print \"BWFN = %0.2f degree \" %BWFN \n", + "#Note : Answer in the book is not accurate." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Diameter of parabolic reflector = 0.058 m**2\n", + "HPBW = 5.16 degree \n", + "BWFN = 10.31 degree \n" + ] + } + ], + "prompt_number": 29 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 6.10 : page 6.43" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Given Data:\n", + "Tau=0.7 #Design Factor\n", + "L1=0.3*2 #in meter\n", + "c=3*10**8 #speednof light in m/s\n", + "f1=(c/(2*L1))/10**6 #in MHz\n", + "#Design factor : L1/L2=L2/L3=L3/L4=.......=0.7\n", + "L2=0.7/L1 #in meter\n", + "f2=f1*0.7 #in MHz\n", + "f3=f2*0.7 #in MHz\n", + "f4=f3*0.7 #in MHz\n", + "f5=f4*0.7 #in MHz\n", + "f6=f5*0.7 #in MHz\n", + "f7=f6*0.7 #in MHz\n", + "f8=f7*0.7 #in MHz\n", + "f9=f8*0.7 #in MHz\n", + "f10=f9*0.7 #in MHz\n", + "print \"Cutoff frequencies in MHz :\"\n", + "print \"f1 = %0.2f MHz \" %f1 \n", + "print \"f2 = %0.2f MHz \" %f2\n", + "print \"f3 = %0.2f MHz \" %f3\n", + "print \"f4 = %0.2f MHz \" %f4\n", + "print \"f5 = %0.2f MHz \" %f5\n", + "print \"f6 = %0.2f MHz \" %f6\n", + "print \"f7 = %0.2f MHz \" %f7\n", + "print \"f8 = %0.2f MHz \" %f8\n", + "print \"f9 = %0.2f MHz \" %f9\n", + "print \"f10 = %0.2f MHz \" %f10\n", + "print \"Passband = %0.2f\"%(f1-f10 )" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Cutoff frequencies in MHz :\n", + "f1 = 250.00 MHz \n", + "f2 = 175.00 MHz \n", + "f3 = 122.50 MHz \n", + "f4 = 85.75 MHz \n", + "f5 = 60.02 MHz \n", + "f6 = 42.02 MHz \n", + "f7 = 29.41 MHz \n", + "f8 = 20.59 MHz \n", + "f9 = 14.41 MHz \n", + "f10 = 10.09 MHz \n", + "Passband = 239.91\n" + ] + } + ], + "prompt_number": 31 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 6.11 : page 6.44" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi, tan, acos\n", + "#Given Data:\n", + "from sympy import symbols, simplify, atan, acos, N\n", + "lamda = symbols('lamda', real = True)\n", + "#Assuming typical values for f \n", + "f1=0.2*lamda #in E-plane \n", + "f2=0.375*lamda # in H-plane\" \n", + "b=10*lamda # mouth height\n", + "delta=0.2*lamda\n", + "print \"Length :\"\n", + "L=pow(b,2)/(8*delta)\n", + "print (L)\n", + "print \"Flare Angle (Theta):\",\n", + "Theta=atan(b/(2*L))*180/pi\n", + "print round(Theta,1),'degree'\n", + "print \"Flare Angle (fi):\",\n", + "delta=0.375*lamda\n", + "fi=acos(L/(L+delta))*180/pi # degree\n", + "print round(fi,1),'degree'\n", + "print \"fi=\",(acos((10**2/(8*0.2))/((10**2/(8*0.2))+0.375))),\" radian\" \n", + "fi=(acos((10**2/(8*0.2))/((10**2/(8*0.2))+0.375))) #in Degree\n", + "print \"Width :\" \n", + "a=2*L*tan(fi)\n", + "print N(a,3)" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Length :\n", + "62.5*lamda\n", + "Flare Angle (Theta): 4.6 degree\n", + "Flare Angle (fi): 6.3 degree\n", + "fi= 0.109271705413178 radian\n", + "Width :\n", + "13.7*lamda\n" + ] + } + ], + "prompt_number": 30 + } + ], + "metadata": {} + } + ] +} diff --git a/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter7.ipynb b/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter7.ipynb new file mode 100755 index 00000000..a5b1cf64 --- /dev/null +++ b/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter7.ipynb @@ -0,0 +1,274 @@ +{ + "metadata": { + "name": "", + "signature": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "chapter 07 : Antenna Measurements" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 7.1 : page 7.28" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#given data :\n", + "f=6.0 #in GHz\n", + "f=f*10**9 #in Hz\n", + "d=10 #in feet\n", + "d=3.048 #in meter\n", + "c=3*10**8 #in m/s\n", + "lamda=c/f #in meters\n", + "rmin=2*d**2/lamda #in meters\n", + "print \"Minimum separation distance = %0.2f m\" %rmin" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Minimum separation distance = 371.61 m\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 7.2 : page 7.28" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#given data :\n", + "GP=12.5 #unitless\n", + "P_dB=23 #in dB\n", + "P=10**(P_dB/10) #unitless\n", + "G=GP*P #unitless\n", + "GdB=GP+P_dB #in dB\n", + "print \"Gain of large antenna = %0.2f \"% GdB\n", + "#Note : Answer in the book is wrong." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Gain of large antenna = 35.50 \n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 7.3 : page 7.28" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import log10\n", + "#given data :\n", + "print \"Open mouth aperture, D = 10*lambda\" \n", + "print \"Power gain : GP = 6*(D/labda)**2\" \n", + "GP=6*10**2 #unitless\n", + "GPdB=10*log10(GP)\n", + "print \"Power gain = %0.1f dB \" %GPdB " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Open mouth aperture, D = 10*lambda\n", + "Power gain : GP = 6*(D/labda)**2\n", + "Power gain = 27.8 dB \n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 7.4 : page 7.28" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#given data :\n", + "f=30000. #in MHz\n", + "f=f*10**6 #in Hz\n", + "d=20 #in feet\n", + "d=20*0.3048 #in meter\n", + "c=3*10**8 #in m/s\n", + "lamda=c/f #in meters\n", + "r=2*d**2/lamda #in meters\n", + "print \"Minimum distance between primary and secondary = %0.2f m\" %r\n", + "# Answe wrong in the textbook." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Minimum distance between primary and secondary = 7432.24 m\n" + ] + } + ], + "prompt_number": 21 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 7.5 : page 7.29" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#given data :\n", + "f=1.2 #in GHz\n", + "f=f*10**9 #in Hz\n", + "BWFN=5 #in degree\n", + "c=3.0*10**8 #in m/s\n", + "lamda=c/f #in meters\n", + "D=140*lamda/BWFN #in meters\n", + "print \"Diameter of a paraboloidal reflector = %0.2f m\" %D" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Diameter of a paraboloidal reflector = 7.00 m\n" + ] + } + ], + "prompt_number": 23 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 7.6 : page 7.29" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi, sqrt, log10\n", + "#given data :\n", + "f=9.0 #in GHz\n", + "f=f*10**9 #in Hz\n", + "c=3*10**8 #in m/s\n", + "lamda=c/f #in meters\n", + "r=35 #in cm\n", + "r=r*10**-2 #in meters\n", + "Attenuation=9.8 #in dB\n", + "#Formula : 10*log10(WT/Wr) = 9.8dB\n", + "WTbyWr=10**(Attenuation/10) #unitless\n", + "D=(4*pi*r/lamda)*(sqrt(1/WTbyWr)) #unitless\n", + "D_dB=10*log10(D) \n", + "print \"Gain of the horn = %0.2f dB \" %D_dB " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Gain of the horn = 16.30 dB \n" + ] + } + ], + "prompt_number": 25 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 7.7 : page 7.29" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi, sqrt, log10\n", + "#given data :\n", + "ratio = 28 # length:diameter\n", + "from sympy import symbols, N\n", + "lamda = symbols('lamda')\n", + "L = 0.925*lamda\n", + "Z = 710+1J*0 # ohm\n", + "Zs = 35476/Z # ohm\n", + "D = L/ratio\n", + "omega = 2*D\n", + "print \"omega =\",N(omega,2)" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "omega = 0.066*lamda\n" + ] + } + ], + "prompt_number": 28 + } + ], + "metadata": {} + } + ] +} diff --git a/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter9.ipynb b/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter9.ipynb new file mode 100755 index 00000000..bac053ad --- /dev/null +++ b/Antenna_and_Wave_Propagation_by_k.k._sharma/chapter9.ipynb @@ -0,0 +1,351 @@ +{ + "metadata": { + "name": "", + "signature": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "chapter 09 : Ground wave propagation" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 9.1 : page 9-23" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt\n", + "#given data :\n", + "HT=50 #in meter\n", + "HR=10 #in meter\n", + "f=60 #in MHz\n", + "P=10 #in KW\n", + "D=10 #in Km\n", + "D=D*10**3 #in m\n", + "c=3*10**8 #speed of light in m/s\n", + "lamda=c/(f*10**6) #in meter\n", + "#Part (i) \n", + "d=3.55*(sqrt(HT)+sqrt(HR)) #in Km\n", + "print \"Maximum line of sight range = %0.2f km \" %d \n", + "#Part (ii)\n", + "Et=88*sqrt(P*1000)*HT*HR/(lamda*D**2)\n", + "print \"The field strength at 10 km = %0.1e V/m\" %Et \n", + "#Part (iii)\n", + "#Formula : Et=88*sqrt(p)*HT*HR/(lambda*D**2)\n", + "Et=1 #in mV/m\n", + "D=sqrt(88*sqrt(P*1000)*HT*HR/(lamda*Et*10**-3)) #in m\n", + "print \"Distance = %0.3f km \" %(D/1000)" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Maximum line of sight range = 36.33 km \n", + "The field strength at 10 km = 8.8e-03 V/m\n", + "Distance = 29.665 km \n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 9.2 : page 9-24" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt\n", + "#given data :\n", + "P=200 #in KW\n", + "D=20 #in Km\n", + "D=D*10**3 #in m\n", + "E=300*sqrt(P)/D #in V/m\n", + "print \"Field Strength at 20 km = %0.2f mV/m \" %(E*10**3)" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Field Strength at 20 km = 212.13 mV/m \n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 9.3 : page 9-24" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#given data :\n", + "HT=10 #in meter\n", + "HR=3 #in meter\n", + "P=200 #in W\n", + "D=50 #in Km\n", + "D=D*10**3 #in Km\n", + "f=150 #in MHz\n", + "c=3*10**8 #speed of light in m/s\n", + "lamda=c/(f*10**6) #in meter\n", + "E=88*sqrt(P)*HT*HR/(lamda*D**2) #in m\n", + "print \"Field Strength at 20 km = %0.2f microV/m \" %(E*10**6)" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Field Strength at 20 km = 7.47 microV/m \n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 9.4 : page 9-25" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt\n", + "#given data :\n", + "HT=100 #in meter\n", + "d=60 #in Km\n", + "#Formula : d=4.12*(sqrt(HT)+sqrt(HR)) #in Km\n", + "HR=(d/4.12-sqrt(HT))**2 #in meter\n", + "print \"Height of receiving antenna = %0.2f m\" %HR" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Height of receiving antenna = 20.82 m\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 9.5 : page 9-25" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt\n", + "#given data :\n", + "HT=3000 #in meter\n", + "HR=6000 #in meter\n", + "d=4.12*(sqrt(HT)+sqrt(HR)) #in Km\n", + "print \"Maximum possible distance = %0.2f km\" %d" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Maximum possible distance = 544.80 km\n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 9.6 : page 9-25" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import log10\n", + "#given data :\n", + "f_MHz=3000 #in MHz\n", + "d_Km=384000 #in Km\n", + "PathLoss=32.45+20*log10(f_MHz)+20*log10(d_Km) #in dB\n", + "print \"Path loss = %0.2f dB \" %PathLoss" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Path loss = 213.68 dB \n" + ] + } + ], + "prompt_number": 9 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 9.7 : page 9-26" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#given data :\n", + "#Part (i)\n", + "D=10 #in Km\n", + "lamda=10000 #in meter\n", + "LP=(4*pi*D*1000/lamda)**2 #in dB\n", + "print \"Path loss = %0.2f dB\" %LP\n", + "#Part (ii)\n", + "D=10**6 #in Km\n", + "lamda=0.3 #in cm\n", + "LP=(4*pi*D*1000/(lamda*10**-2))**2 #in dB\n", + "print \"Path loss = %0.2e dB \" %LP \n", + "#Note : Answer in the book is wrong as value putted in the solution is differ from given in question." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Path loss = 157.91 dB\n", + "Path loss = 1.75e+25 dB \n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 9.8 : page 9-26" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import sqrt\n", + "#given data :\n", + "HT=50 #in meter\n", + "HR=5 #in meter\n", + "d=4.12*(sqrt(HT)+sqrt(HR)) #in Km\n", + "print \"Range of LOS system = %0.2f km\"%d" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Range of LOS system = 38.35 km\n" + ] + } + ], + "prompt_number": 11 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exa 9.9 : page 9-26" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import pi\n", + "#given data :\n", + "PT=5.0 #in KW\n", + "PT=PT*1000 #in W\n", + "D=100.0 #in Km\n", + "D=D*10**3 #in m\n", + "f=300.0 #in MHz\n", + "GT=1.64 #Directivity of transmitter\n", + "GR=1.64 #Directivity of receiver\n", + "c=3*10**8 #speed of light in m/s\n", + "lamda=c/(f*10**6) #in meter\n", + "Pr=PT*GT*GR*(lamda/(4*pi*D))**2\n", + "print \"Maximum power received = %0.3e W\"% Pr\n", + "# Answer wrong in the textbook." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Maximum power received = 8.516e-09 W\n" + ] + } + ], + "prompt_number": 14 + } + ], + "metadata": {} + } + ] +} |