summaryrefslogtreecommitdiff
path: root/Satellite_Communications/Chapter_4.ipynb
blob: 7f767a844fd3d090265ff50501db90bcd8856160 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
{
 "metadata": {
  "name": "",
  "signature": "sha256:3889f080e541059ccdb3a0ca28e3a8926a614c5d9e6ab3f2de315570240fb716"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "heading",
     "level": 1,
     "metadata": {},
     "source": [
      "Chapter 4: Radio Wave Propagation"
     ]
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 4.1, Page 97"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "#Variable Declaration\n",
      "\n",
      "El=50  #Elevation Angle(degrees)\n",
      "h0=0.6 #Earth station altitude(km)\n",
      "hr=3   #Rain height(km)\n",
      "R01=10 #Point Rain Rate(mm/hr)\n",
      "f=12   #frequency(GHz)\n",
      "ah=0.0188\n",
      "bh=1.217\n",
      "av=0.0168\n",
      "bv=1.2\n",
      "#Calculation\n",
      "Ls=(hr-h0)/math.sin(El*3.142/180) #Slant path length(km)\n",
      "LG=Ls*math.cos(El*3.142/180)      #Horizontal projection(km)\n",
      "r01=90/(90+4*LG)      #Reduction factor\n",
      "L=Ls*r01              #Effective path length(km)\n",
      "alphah=ah*R01**bh     #Specific Attenuation\n",
      "AdBh=round(alphah*L,2)#Rain Attenuation for horizontal polarization\n",
      "alphav=av*R01**bv   #Specific Attenuation\n",
      "AdBv=round(alphav*L,2)#Rain Attenuation for vertical polarization\n",
      "\n",
      "#Results\n",
      "\n",
      "print\"Rain Attenuation for given conditions and horizontal polarization is\",AdBh,\"dB\"\n",
      "print\"Rain Attenuation for given conditions and vertical polarization is\",AdBv,\"dB\"\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Rain Attenuation for given conditions and horizontal polarization is 0.89 dB\n",
        "Rain Attenuation for given conditions and vertical polarization is 0.77 dB\n"
       ]
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 4.2, Page 99"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "#Variable Declaration\n",
      "ah=0.0188\n",
      "bh=1.217\n",
      "av=0.0168\n",
      "bv=1.2\n",
      "R01=10  #Point Rain Rate(mm/hr)\n",
      "L=2.8753812  #Effective path length calculated in Example 4.1(km)\n",
      "\n",
      "#Calculation\n",
      "\n",
      "#Factors depending on frequency and polarization\n",
      "ac=(ah+av)/2 #a for circular polarization\n",
      "bc=(ah*bh+av*bv)/(2*ac) #b for circular polarization\n",
      "\n",
      "alpha=ac*R01**bc  #Specific Attenuation(dB)\n",
      "AdB=round(alpha*L,2)  #Rain Attenuation(dB)\n",
      "\n",
      "\n",
      "#Results\n",
      "\n",
      "print \"The Rain Attenuation for circular polarization is\",AdB,\"dB\""
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The Rain Attenuation for circular polarization is 0.83 dB\n"
       ]
      }
     ],
     "prompt_number": 2
    }
   ],
   "metadata": {}
  }
 ]
}