Microwave Engineering - G. S. Raghuvanshi

Scilab Textbook Companion for
Microwave Engineering
by G. S. Raghuvanshi1
Created by
Ankit Suri
B. Tech.
Computer Engineering
ABES Engineering College
College Teacher
None
Cross-Checked by
Spandana
June 6, 2016
1 Funded by a grant from the National Mission on Education through ICT,

http://spoken-tutorial.org/NMEICT-Intro. This Textbook Companion and Scilab
codes written in it can be downloaded from the ”Textbook Companion Project”
section at the website http://scilab.in
Book Description
Title: Microwave Engineering
Author: G. S. Raghuvanshi
Publisher: Cengage Learning, New Delhi
Edition: 1
Year: 2012
ISBN: 978-81-315-1721-5
1
Scilab numbering policy used in this document and the relation to the
above book.
Exa Example (Solved example)
Eqn Equation (Particular equation of the above book)
AP Appendix to Example(Scilab Code that is an Appednix to a particular

Example of the above book)
For example, Exa 3.51 means solved example 3.51 of this book. Sec 2.3 means
a scilab code whose theory is explained in Section 2.3 of the book.
2
Contents
List of Scilab Codes 4
1 Microwaves 6
2 Waveguides 35
3 Microwave Network Analysis 83
4 Microwave Resonators and Waveguide Components 92
5 Microwave Tubes Klystrons 129
6 Microwave Travelling Wave Tubes O type 166
7 Cross Field Microwave Tubes M Type 187
8 Microwave Solid State Control Devices 211
9 Microwave Solid State Generators and Amplifiers 230
10 Striplines and Microstrip Lines 267
11 Microwave Integrated Circuits 290
12 Microwave Measurements 294
3
List of Scilab Codes
Exa 1.2 Lossless line . . . . . . . . . . . . . . . . . . . . . . . . 6

Exa 1.3 Microwave line . . . . . . . . . . . . . . . . . . . . . . 9
Exa 1.4 Quater wave transformer . . . . . . . . . . . . . . . . 10
Exa 1.5 Parallel resonant circuit . . . . . . . . . . . . . . . . . 10
Exa 1.6 Lossless line . . . . . . . . . . . . . . . . . . . . . . . . 13
Exa 1.7 Lossy cable . . . . . . . . . . . . . . . . . . . . . . . . 15
Exa 1.8 Transmission Line . . . . . . . . . . . . . . . . . . . . 17
Exa 1.12 Active Device . . . . . . . . . . . . . . . . . . . . . . . 23
Exa 1.13 Transmission line . . . . . . . . . . . . . . . . . . . . . 24
Exa 1.16 EM Plane . . . . . . . . . . . . . . . . . . . . . . . . . 28
Exa 1.17 Polyethylene . . . . . . . . . . . . . . . . . . . . . . . 29
Exa 1.18 Electromagnetic wave . . . . . . . . . . . . . . . . . . 31
Exa 1.19 Range in sea water . . . . . . . . . . . . . . . . . . . . 32
Exa 2.1 Dominant mode . . . . . . . . . . . . . . . . . . . . . 35
Exa 2.2 Rectangular Waveguide . . . . . . . . . . . . . . . . . 37
Exa 2.4 Wider dimension . . . . . . . . . . . . . . . . . . . . . 40
Exa 2.5 Rectangular waveguide . . . . . . . . . . . . . . . . . 43
Exa 2.8 Lossless Rectangular Waveguide . . . . . . . . . . . . 46
Exa 2.9 Dimensions . . . . . . . . . . . . . . . . . . . . . . . . 48
Exa 2.14 Waveguide . . . . . . . . . . . . . . . . . . . . . . . . 52
Exa 2.17 Air filled Rectangular Waveguide . . . . . . . . . . . . 54
4
Exa 2.20 Ratio of cross section . . . . . . . . . . . . . . . . . . 60
Exa 2.23 Circular waveguide . . . . . . . . . . . . . . . . . . . . 65
Exa 2.24 Air filled circular waveguide . . . . . . . . . . . . . . . 67
Exa 2.25 Air filled rectangular waveguide . . . . . . . . . . . . . 69
Exa 2.27 Circular Waveguide . . . . . . . . . . . . . . . . . . . 72
Exa 2.29 Circular waveguide . . . . . . . . . . . . . . . . . . . . 77
Exa 2.30 Cutoff frequencies . . . . . . . . . . . . . . . . . . . . 77
Exa 2.31 Dielectric constant . . . . . . . . . . . . . . . . . . . . 80
Exa 3.4 Scattering matrix . . . . . . . . . . . . . . . . . . . . 83
Exa 3.5 Network . . . . . . . . . . . . . . . . . . . . . . . . . . 83
Exa 3.6 Microwave network . . . . . . . . . . . . . . . . . . . . 86
Exa 3.12 Transistor amplifier circuit . . . . . . . . . . . . . . . 88
Exa 4.1 Rectangular cavity resonator . . . . . . . . . . . . . . 92
Exa 4.2 Rectangulr resonator . . . . . . . . . . . . . . . . . . . 94
Exa 4.3 X band resonator . . . . . . . . . . . . . . . . . . . . . 95
Exa 4.4 Rectangular resonator . . . . . . . . . . . . . . . . . . 97
Exa 4.5 Cylindrical resonator . . . . . . . . . . . . . . . . . . . 99
Exa 4.6 Resonator comparison . . . . . . . . . . . . . . . . . . 102
Exa 4.7 Cubical Resonator . . . . . . . . . . . . . . . . . . . . 103
Exa 4.8 Rectangular Resonant Cavity . . . . . . . . . . . . . . 104
Exa 4.9 Rectangular resonator . . . . . . . . . . . . . . . . . . 106
Exa 4.10 Cylindrical resonator . . . . . . . . . . . . . . . . . . . 108
Exa 4.11 Cylindrical resonantor . . . . . . . . . . . . . . . . . . 111
Exa 4.12 Cylindrical Resonator . . . . . . . . . . . . . . . . . . 111
Exa 4.13 Lossless plane H tee . . . . . . . . . . . . . . . . . . . 115
Exa 4.14 E plane tee . . . . . . . . . . . . . . . . . . . . . . . . 115
Exa 4.15 Magic Tee . . . . . . . . . . . . . . . . . . . . . . . . . 117
Exa 4.18 Directional Coupler . . . . . . . . . . . . . . . . . . . 119
Exa 4.20 Directional coupler . . . . . . . . . . . . . . . . . . . . 122
Exa 4.21 Isolator Matrix . . . . . . . . . . . . . . . . . . . . . . 122
Exa 4.22 Circulator Matrix . . . . . . . . . . . . . . . . . . . . 124
5
Exa 4.24 Attenuator matrix . . . . . . . . . . . . . . . . . . . . 127
Exa 5.1 Two Cavity Klystron . . . . . . . . . . . . . . . . . . 129
Exa 5.2 Two cavity Klystron . . . . . . . . . . . . . . . . . . . 131
Exa 5.9 Four cavity Klystron . . . . . . . . . . . . . . . . . . . 146
Exa 5.10 Four cavity Klystron . . . . . . . . . . . . . . . . . . . 148
Exa 5.11 Reflex Klystron . . . . . . . . . . . . . . . . . . . . . . 150
Exa 5.17 Electron Gun . . . . . . . . . . . . . . . . . . . . . . . 161
Exa 5.18 Re entrant Coaxial Cavity . . . . . . . . . . . . . . . . 163
Exa 6.1 TWT . . . . . . . . . . . . . . . . . . . . . . . . . . . 166
Exa 6.2 Helix TWT . . . . . . . . . . . . . . . . . . . . . . . . 168
Exa 6.3 Helical TWT . . . . . . . . . . . . . . . . . . . . . . . 170
Exa 6.4 O type TWT . . . . . . . . . . . . . . . . . . . . . . . 170
Exa 6.5 Cavity coupled . . . . . . . . . . . . . . . . . . . . . . 173
Exa 6.6 O Type Backward Wave amplifier . . . . . . . . . . . 175
Exa 6.7 Multicavity TWT . . . . . . . . . . . . . . . . . . . . 176
Exa 6.8 Gridded TWT . . . . . . . . . . . . . . . . . . . . . . 177
Exa 6.9 Helix TWT . . . . . . . . . . . . . . . . . . . . . . . . 179
Exa 6.10 Low Power TWT . . . . . . . . . . . . . . . . . . . . . 180
Exa 6.11 TWT . . . . . . . . . . . . . . . . . . . . . . . . . . . 182
Exa 6.12 TWT . . . . . . . . . . . . . . . . . . . . . . . . . . . 184
Exa 7.1 X band Magnetron . . . . . . . . . . . . . . . . . . . . 187
Exa 7.2 Cylindrical Magnetron . . . . . . . . . . . . . . . . . . 189
Exa 7.3 Cylindrical Magnetron . . . . . . . . . . . . . . . . . . 191
Exa 7.4 Conventional Magnetron . . . . . . . . . . . . . . . . . 192
Exa 7.5 Conventional Magnetron . . . . . . . . . . . . . . . . . 194
Exa 7.6 Carcinotron . . . . . . . . . . . . . . . . . . . . . . . . 196
Exa 7.7 Frequency Aglile Magnetron . . . . . . . . . . . . . . . 198
6
Exa 7.8 Cross field amplifier . . . . . . . . . . . . . . . . . . . 200
Exa 7.9 Inverted coaxial Magnetron . . . . . . . . . . . . . . . 202
Exa 7.10 Inverted coaxial Magnetron . . . . . . . . . . . . . . . 203
Exa 7.11 Linear Magnetron . . . . . . . . . . . . . . . . . . . . 205
Exa 7.12 Inverted Coaxial Magnetron . . . . . . . . . . . . . . . 206
Exa 7.13 Agile coaxial Magnetron . . . . . . . . . . . . . . . . . 208
Exa 8.1 Single pole Switch . . . . . . . . . . . . . . . . . . . . 211
Exa 8.2 Pin diode switches . . . . . . . . . . . . . . . . . . . . 213
Exa 8.3 Silicon switching diode . . . . . . . . . . . . . . . . . . 216
Exa 8.6 Parametric upconverter . . . . . . . . . . . . . . . . . 217
Exa 8.7 Parametric amplifier . . . . . . . . . . . . . . . . . . . 219
Exa 8.8 Negative resistance parametric amplifier . . . . . . . . 221
Exa 8.9 330 stage CCD . . . . . . . . . . . . . . . . . . . . . . 223
Exa 8.10 3 phase CCD . . . . . . . . . . . . . . . . . . . . . . . 223
Exa 8.11 Surface channel CCD . . . . . . . . . . . . . . . . . . 225
Exa 8.12 3 phase CCD . . . . . . . . . . . . . . . . . . . . . . . 226
Exa 8.13 Surface channel CCD . . . . . . . . . . . . . . . . . . 228
Exa 9.2 Bipolar transistor . . . . . . . . . . . . . . . . . . . . 230
Exa 9.3 Heterojunction transistor . . . . . . . . . . . . . . . . 233
Exa 9.4 GaAs FET . . . . . . . . . . . . . . . . . . . . . . . . 233
Exa 9.5 Microwave transistor . . . . . . . . . . . . . . . . . . . 235
Exa 9.6 Transistor Amplifier . . . . . . . . . . . . . . . . . . . 237
Exa 9.7 Microwave transistor . . . . . . . . . . . . . . . . . . . 241
Exa 9.8 JEFT . . . . . . . . . . . . . . . . . . . . . . . . . . . 241
Exa 9.9 MESFET . . . . . . . . . . . . . . . . . . . . . . . . . 243
Exa 9.10 Gunn device . . . . . . . . . . . . . . . . . . . . . . . 245
Exa 9.11 Gunn oscillator . . . . . . . . . . . . . . . . . . . . . . 247
Exa 9.12 Tunnel diode . . . . . . . . . . . . . . . . . . . . . . . 248
Exa 9.13 Gunn diode . . . . . . . . . . . . . . . . . . . . . . . . 248
Exa 9.14 Gunn diode . . . . . . . . . . . . . . . . . . . . . . . . 251
Exa 9.15 Gunn diode . . . . . . . . . . . . . . . . . . . . . . . . 252
Exa 9.17 IMPATT diode . . . . . . . . . . . . . . . . . . . . . . 253
Exa 9.18 IMPATT diode . . . . . . . . . . . . . . . . . . . . . . 255
Exa 9.19 IMPATT diode . . . . . . . . . . . . . . . . . . . . . . 257
Exa 9.20 IMPATT diode . . . . . . . . . . . . . . . . . . . . . . 258
Exa 9.21 TRAPATT diode . . . . . . . . . . . . . . . . . . . . . 260
Exa 9.22 BARITT diode . . . . . . . . . . . . . . . . . . . . . . 260
Exa 9.23 Laser . . . . . . . . . . . . . . . . . . . . . . . . . . . 263
7
Exa 9.24 Laser . . . . . . . . . . . . . . . . . . . . . . . . . . . 264
Exa 9.25 Heterojunction laser . . . . . . . . . . . . . . . . . . . 265
Exa 10.1 Copper stripline . . . . . . . . . . . . . . . . . . . . . 267
Exa 10.2 Microstrip line . . . . . . . . . . . . . . . . . . . . . . 269
Exa 10.3 Microstrip . . . . . . . . . . . . . . . . . . . . . . . . . 271
Exa 10.4 Stripline . . . . . . . . . . . . . . . . . . . . . . . . . . 274
Exa 10.7 Parallel stripline . . . . . . . . . . . . . . . . . . . . . 279
Exa 10.8 Stripline coupler . . . . . . . . . . . . . . . . . . . . . 279
Exa 10.9 Branch coupler . . . . . . . . . . . . . . . . . . . . . . 281
Exa 10.10 Broadside stripline . . . . . . . . . . . . . . . . . . . . 283
Exa 10.11 Paralle stripline . . . . . . . . . . . . . . . . . . . . . 285
Exa 10.12 Shielded stripline . . . . . . . . . . . . . . . . . . . . . 286
Exa 10.13 Lossless stripline . . . . . . . . . . . . . . . . . . . . . 289
Exa 11.1 Costs . . . . . . . . . . . . . . . . . . . . . . . . . . . 290
Exa 11.2 Yield . . . . . . . . . . . . . . . . . . . . . . . . . . . 293
Exa 12.1 Microwave diode . . . . . . . . . . . . . . . . . . . . . 294
Exa 12.2 Detector mismatch . . . . . . . . . . . . . . . . . . . . 296
Exa 12.3 Transmission waveguide . . . . . . . . . . . . . . . . . 297
Exa 12.4 VSWR of waveguide . . . . . . . . . . . . . . . . . . . 298
Exa 12.5 Directional couplers . . . . . . . . . . . . . . . . . . . 299
Exa 12.9 Three port circulator . . . . . . . . . . . . . . . . . . . 304
Exa 12.10 Air filled cavity . . . . . . . . . . . . . . . . . . . . . . 306
8
List of Figures
1.1 Lossless line . . . . . . . . . . . . . . . . . . . . . . . . . . . 7

1.2 Microwave line . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.3 Quater wave transformer . . . . . . . . . . . . . . . . . . . . 11
1.4 Parallel resonant circuit . . . . . . . . . . . . . . . . . . . . 12
1.5 Lossless line . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
1.6 Lossy cable . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.7 Transmission Line . . . . . . . . . . . . . . . . . . . . . . . . 18
1.10 Microwave line . . . . . . . . . . . . . . . . . . . . . . . . . 23
1.11 Active Device . . . . . . . . . . . . . . . . . . . . . . . . . . 24
1.12 Transmission line . . . . . . . . . . . . . . . . . . . . . . . . 25
1.13 Microwave line . . . . . . . . . . . . . . . . . . . . . . . . . 27
1.14 EM Plane . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
1.15 Polyethylene . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
1.16 Electromagnetic wave . . . . . . . . . . . . . . . . . . . . . . 31
1.17 Range in sea water . . . . . . . . . . . . . . . . . . . . . . . 33
2.1 Dominant mode . . . . . . . . . . . . . . . . . . . . . . . . . 36

2.2 Rectangular Waveguide . . . . . . . . . . . . . . . . . . . . . 37
2.4 Wider dimension . . . . . . . . . . . . . . . . . . . . . . . . 41
2.5 Rectangular waveguide . . . . . . . . . . . . . . . . . . . . . 42
2.7 Lossless Rectangular Waveguide . . . . . . . . . . . . . . . . 47
2.8 Dimensions . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
2.10 Waveguide . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
2.11 Air filled Rectangular Waveguide . . . . . . . . . . . . . . . 55
9
2.14 Ratio of cross section . . . . . . . . . . . . . . . . . . . . . . 61
2.17 Circular waveguide . . . . . . . . . . . . . . . . . . . . . . . 66
2.18 Air filled circular waveguide . . . . . . . . . . . . . . . . . . 68
2.19 Air filled rectangular waveguide . . . . . . . . . . . . . . . . 70
2.21 Circular Waveguide . . . . . . . . . . . . . . . . . . . . . . . 73
2.23 Circular waveguide . . . . . . . . . . . . . . . . . . . . . . . 78
2.24 Cutoff frequencies . . . . . . . . . . . . . . . . . . . . . . . . 79
2.25 Dielectric constant . . . . . . . . . . . . . . . . . . . . . . . 81
3.1 Scattering matrix . . . . . . . . . . . . . . . . . . . . . . . . 84

3.2 Network . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
3.3 Microwave network . . . . . . . . . . . . . . . . . . . . . . . 87
3.4 Transistor amplifier circuit . . . . . . . . . . . . . . . . . . . 89
4.1 Rectangular cavity resonator . . . . . . . . . . . . . . . . . . 93

4.2 Rectangulr resonator . . . . . . . . . . . . . . . . . . . . . . 94
4.3 X band resonator . . . . . . . . . . . . . . . . . . . . . . . . 96
4.4 Rectangular resonator . . . . . . . . . . . . . . . . . . . . . 98
4.5 Cylindrical resonator . . . . . . . . . . . . . . . . . . . . . . 100
4.6 Resonator comparison . . . . . . . . . . . . . . . . . . . . . 102
4.7 Cubical Resonator . . . . . . . . . . . . . . . . . . . . . . . 103
4.8 Rectangular Resonant Cavity . . . . . . . . . . . . . . . . . 105
4.9 Rectangular resonator . . . . . . . . . . . . . . . . . . . . . 107
4.10 Cylindrical resonator . . . . . . . . . . . . . . . . . . . . . . 109
4.11 Cylindrical resonantor . . . . . . . . . . . . . . . . . . . . . 110
4.12 Cylindrical Resonator . . . . . . . . . . . . . . . . . . . . . . 112
4.13 Lossless plane H tee . . . . . . . . . . . . . . . . . . . . . . . 114
4.14 E plane tee . . . . . . . . . . . . . . . . . . . . . . . . . . . 116
4.15 Magic Tee . . . . . . . . . . . . . . . . . . . . . . . . . . . . 118
4.16 Directional Coupler . . . . . . . . . . . . . . . . . . . . . . . 120
4.17 Directional coupler . . . . . . . . . . . . . . . . . . . . . . . 121
4.18 Isolator Matrix . . . . . . . . . . . . . . . . . . . . . . . . . 123
10
4.19 Circulator Matrix . . . . . . . . . . . . . . . . . . . . . . . . 124
4.21 Attenuator matrix . . . . . . . . . . . . . . . . . . . . . . . 128
5.1 Two Cavity Klystron . . . . . . . . . . . . . . . . . . . . . . 130

5.2 Two cavity Klystron . . . . . . . . . . . . . . . . . . . . . . 132
5.9 Four cavity Klystron . . . . . . . . . . . . . . . . . . . . . . 147
5.10 Four cavity Klystron . . . . . . . . . . . . . . . . . . . . . . 149
5.11 Reflex Klystron . . . . . . . . . . . . . . . . . . . . . . . . . 151
5.12 Reflex Klystron . . . . . . . . . . . . . . . . . . . . . . . . . 153
5.13 Reflex Klystron . . . . . . . . . . . . . . . . . . . . . . . . . 155
5.14 Reflex Klystron . . . . . . . . . . . . . . . . . . . . . . . . . 157
5.15 Reflex Klystron . . . . . . . . . . . . . . . . . . . . . . . . . 158
5.16 Reflex Klystron . . . . . . . . . . . . . . . . . . . . . . . . . 160
5.17 Electron Gun . . . . . . . . . . . . . . . . . . . . . . . . . . 162
5.18 Re entrant Coaxial Cavity . . . . . . . . . . . . . . . . . . . 164
6.1 TWT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 167

6.2 Helix TWT . . . . . . . . . . . . . . . . . . . . . . . . . . . 169
6.3 Helical TWT . . . . . . . . . . . . . . . . . . . . . . . . . . 171
6.4 O type TWT . . . . . . . . . . . . . . . . . . . . . . . . . . 172
6.5 Cavity coupled . . . . . . . . . . . . . . . . . . . . . . . . . 174
6.6 O Type Backward Wave amplifier . . . . . . . . . . . . . . . 175
6.7 Multicavity TWT . . . . . . . . . . . . . . . . . . . . . . . . 176
6.8 Gridded TWT . . . . . . . . . . . . . . . . . . . . . . . . . . 178
6.9 Helix TWT . . . . . . . . . . . . . . . . . . . . . . . . . . . 179
6.10 Low Power TWT . . . . . . . . . . . . . . . . . . . . . . . . 181
6.11 TWT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 183
6.12 TWT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 185
7.1 X band Magnetron . . . . . . . . . . . . . . . . . . . . . . . 188

7.2 Cylindrical Magnetron . . . . . . . . . . . . . . . . . . . . . 190
11
7.3 Cylindrical Magnetron . . . . . . . . . . . . . . . . . . . . . 191
7.4 Conventional Magnetron . . . . . . . . . . . . . . . . . . . . 193
7.5 Conventional Magnetron . . . . . . . . . . . . . . . . . . . . 195
7.6 Carcinotron . . . . . . . . . . . . . . . . . . . . . . . . . . . 197
7.7 Frequency Aglile Magnetron . . . . . . . . . . . . . . . . . . 199
7.8 Cross field amplifier . . . . . . . . . . . . . . . . . . . . . . . 201
7.9 Inverted coaxial Magnetron . . . . . . . . . . . . . . . . . . 202
7.10 Inverted coaxial Magnetron . . . . . . . . . . . . . . . . . . 204
7.11 Linear Magnetron . . . . . . . . . . . . . . . . . . . . . . . . 205
7.12 Inverted Coaxial Magnetron . . . . . . . . . . . . . . . . . . 207
7.13 Agile coaxial Magnetron . . . . . . . . . . . . . . . . . . . . 209
8.1 Single pole Switch . . . . . . . . . . . . . . . . . . . . . . . . 212

8.2 Pin diode switches . . . . . . . . . . . . . . . . . . . . . . . 214
8.3 Silicon switching diode . . . . . . . . . . . . . . . . . . . . . 216
8.4 Parametric upconverter . . . . . . . . . . . . . . . . . . . . . 218
8.5 Parametric amplifier . . . . . . . . . . . . . . . . . . . . . . 220
8.6 Negative resistance parametric amplifier . . . . . . . . . . . 221
8.7 330 stage CCD . . . . . . . . . . . . . . . . . . . . . . . . . 222
8.8 3 phase CCD . . . . . . . . . . . . . . . . . . . . . . . . . . 224
8.9 Surface channel CCD . . . . . . . . . . . . . . . . . . . . . . 225
8.10 3 phase CCD . . . . . . . . . . . . . . . . . . . . . . . . . . 227
8.11 Surface channel CCD . . . . . . . . . . . . . . . . . . . . . . 229
9.1 Bipolar transistor . . . . . . . . . . . . . . . . . . . . . . . . 231

9.2 Heterojunction transistor . . . . . . . . . . . . . . . . . . . . 232
9.3 GaAs FET . . . . . . . . . . . . . . . . . . . . . . . . . . . . 234
9.4 Microwave transistor . . . . . . . . . . . . . . . . . . . . . . 236
9.5 Transistor Amplifier . . . . . . . . . . . . . . . . . . . . . . 238
9.6 Microwave transistor . . . . . . . . . . . . . . . . . . . . . . 240
9.7 JEFT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 242
9.8 MESFET . . . . . . . . . . . . . . . . . . . . . . . . . . . . 244
9.9 Gunn device . . . . . . . . . . . . . . . . . . . . . . . . . . . 246
9.10 Gunn oscillator . . . . . . . . . . . . . . . . . . . . . . . . . 247
9.11 Tunnel diode . . . . . . . . . . . . . . . . . . . . . . . . . . 249
9.12 Gunn diode . . . . . . . . . . . . . . . . . . . . . . . . . . . 250
9.13 Gunn diode . . . . . . . . . . . . . . . . . . . . . . . . . . . 251
9.14 Gunn diode . . . . . . . . . . . . . . . . . . . . . . . . . . . 252
12
9.15 IMPATT diode . . . . . . . . . . . . . . . . . . . . . . . . . 254
9.16 IMPATT diode . . . . . . . . . . . . . . . . . . . . . . . . . 256
9.17 IMPATT diode . . . . . . . . . . . . . . . . . . . . . . . . . 257
9.18 IMPATT diode . . . . . . . . . . . . . . . . . . . . . . . . . 259
9.19 TRAPATT diode . . . . . . . . . . . . . . . . . . . . . . . . 261
9.20 BARITT diode . . . . . . . . . . . . . . . . . . . . . . . . . 262
9.21 Laser . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 263
9.22 Laser . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 265
9.23 Heterojunction laser . . . . . . . . . . . . . . . . . . . . . . 266
10.1 Copper stripline . . . . . . . . . . . . . . . . . . . . . . . . . 268

10.2 Microstrip line . . . . . . . . . . . . . . . . . . . . . . . . . . 270
10.3 Microstrip . . . . . . . . . . . . . . . . . . . . . . . . . . . . 272
10.4 Stripline . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 273
10.5 Microstrip line . . . . . . . . . . . . . . . . . . . . . . . . . . 275
10.6 Microstrip line . . . . . . . . . . . . . . . . . . . . . . . . . . 277
10.7 Parallel stripline . . . . . . . . . . . . . . . . . . . . . . . . . 278
10.8 Stripline coupler . . . . . . . . . . . . . . . . . . . . . . . . . 280
10.9 Branch coupler . . . . . . . . . . . . . . . . . . . . . . . . . 282
10.10Broadside stripline . . . . . . . . . . . . . . . . . . . . . . . 284
10.11Paralle stripline . . . . . . . . . . . . . . . . . . . . . . . . . 285
10.12Shielded stripline . . . . . . . . . . . . . . . . . . . . . . . . 287
10.13Lossless stripline . . . . . . . . . . . . . . . . . . . . . . . . 288
11.1 Costs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 291

11.2 Yield . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 292
12.1 Microwave diode . . . . . . . . . . . . . . . . . . . . . . . . 295

12.2 Detector mismatch . . . . . . . . . . . . . . . . . . . . . . . 296
12.3 Transmission waveguide . . . . . . . . . . . . . . . . . . . . 297
12.4 VSWR of waveguide . . . . . . . . . . . . . . . . . . . . . . 298
12.5 Directional couplers . . . . . . . . . . . . . . . . . . . . . . . 299
12.6 Microwave line . . . . . . . . . . . . . . . . . . . . . . . . . 300
12.7 Microwave line . . . . . . . . . . . . . . . . . . . . . . . . . 301
12.9 Three port circulator . . . . . . . . . . . . . . . . . . . . . . 304
12.10Air filled cavity . . . . . . . . . . . . . . . . . . . . . . . . . 305
12.11Rectangular Waveguide . . . . . . . . . . . . . . . . . . . . . 307
13
Chapter 1
Microwaves
Scilab code Exa 1.2 Lossless line
1 // Page Number : 12
2 // Example 1 . 2
3 clc ;
4 // Given
5 z0 =50; //ohm
6 zg =50; //ohm
7 l =0.25; //m
8 f =4 D +9; // hz
9 zl =100; //ohm
10 vg =10; //V
11 w =2* %pi * f ; // r a d / s e c
12 c =3 D +8; //m/ s
13
14 // ( i ) V o l t a g e and c u r r e n t a t any p o i n t
15 tg =( zg - z0 ) /( zg + z0 ) ;
16 tl =( zl - z0 ) /( zl + z0 ) ;
17 vi = z0 * vg /( z0 + zg ) ; //V
18 disp ( ’V ’ ,vi , ’ V o l t a g e a t any p o i n t : ’ ) ;
19 ii = vg /(2* z0 ) ; //A
14
15
Figure 1.1: Lossless line

20 disp ( ’A ’ ,ii , ’ C u r r e n t a t any p o i n t : ’ ) ;
21
22 // ( i i ) V o l t a g e a t g e n e r a t o r end
23 // Taking z=1
24 z =1;
25 bet = w / c ;
26 vz =( vg /2) * exp ( - %i * bet *( z + l ) ) *(1+( tl * exp (2* %i * bet * z ) )
) ; //V
27 disp ( ’V ’ ,vz , ’ V o l t a g e a t g e n e r a t o r end : ’ ) ;
28 iz = ii * exp ( - %i * bet *( z + l ) ) *(1 -( tl * exp (2* %i * bet * z ) ) ) ; //
A
29 vz1 =( vg /2) * exp ( - %i * bet *( z + l ) ) *(1+( tl * exp (2* %i * bet * z )
) ) ; //V
30
31 // V o l t a g e a t l o a d end , z=0
32 z11 =0;
33 vl =( vg /2) * exp ( - %i * bet * l ) *(1+( tl * exp (2* %i * bet * z11 ) ) ) ;
//V
34 disp ( ’V ’ ,vl , ’ V o l t a g e a t l o a d end : ’ ) ;
35
36 // ( i i i ) R e f l e c t i o n c o e f f i c i e n t
37 zx =0.25;
38 tz = tl * exp ( %i *2* bet * zx ) ;
39 disp ( tz , ’ R e f l e c t i o n c o e f f i c i e n t : ’ ) ;
40
41 // ( i v ) VSWR
42 p =(1+ tl ) /(1 - tl ) ;
43 disp (p , ’VSWR: ’ ) ;
44
45 // ( v ) A v e r a g e power d e l i v e r e d t o t h e l o a d
46 vl =20/3;
47 pl0 = vl ^2/(2* zl ) ; //W
48 disp ( ’W’ ,pl0 , ’ A v e r a g e power d e l i v e r e d t o t h e l o a d : ’ )
;
16
Figure 1.2: Microwave line
Scilab code Exa 1.3 Microwave line
2 // Example 1 . 3
3 clc ;
4 // Given
5 pm =3;
6 pl =4;
7 l =24; //cm
8 l1 = l /100; //m
9
10 // A t t e n u a t i o n
11 tin =( pm -1) /( pm +1) ;
12 tl =( pl -1) /( pl +1) ;
13 alp =(1/(2* l1 ) ) * log ( tl / tin ) ; //Np/m
17
14 disp ( ’ Np/m ’ ,alp , ’ A t t e n u a t i o n i n t h e l i n e ’ ) ;
Scilab code Exa 1.4 Quater wave transformer
2 // Example 1 . 4
3 clc ;
4 // Given
5 c =3 D +8; //m/ s
6 z0 =200; //ohm
7 zl =800; //ohm
8 f =30 D +6; // hz
9
10 // C h a r a c t e r s t i c i m p e d a n c e
11 z00 = sqrt ( z0 * zl ) ; //ohm
12 disp ( ’ ohm ’ ,z00 , ’ C h a r a c t e r s t i c i m p e d a n c e : ’ ) ;
13
14 // Length o f l i n e
15 lam = c / f ; //m
16 l = lam /4; //m
17 disp ( ’m ’ ,l , ’ Length o f l i n e : ’ ) ;
Scilab code Exa 1.5 Parallel resonant circuit
2 // Example 1 . 5
3 clc ;
4 // Given
5 l =1.2; //mH
18
Figure 1.3: Quater wave transformer
19
Figure 1.4: Parallel resonant circuit
20
6 r =8; //ohm
7 c =200 D -12; //F
8
9 // ( i ) R e s o n a n t f r e q u e n c y
10 f0 =(1/(2* %pi ) ) * sqrt (1/( l * c ) ) ; // hz
11 disp ( ’ hz ’ ,f0 , ’ R e s o n a n t f r e q u e n c y : ’ ) ;
12
13 // ( i i ) Impedance o f c i r c u i t
14 disp ( ’ ohm ’ ,r , ’ Impedance o f c i r c u i t : ’ ) ;
15
16 // ( i i i )Q f a c t o r o f t h e c i r c u i t
17 q =1/(2* %pi * f0 * c * r ) ;
18 disp (q , ’Q f a c t o r o f t h e c i r c u i t : ’ ) ;
19
20 // ( i v ) Bandwidth
21 df = f0 / q ; // hz
22 disp ( ’ hz ’ ,df , ’ Bandwidth : ’ ) ;
23
24 // The v a l u e o f r e s o n a n t f r e q u e n c y i s c a l c u l a t e d
wrong i n book
25 // Hence Q f a c t o r and bandwidth , a l l t h e s e a n s w e r s
d o n t match
Scilab code Exa 1.6 Lossless line
1 // Page Number :
2 // Example 1 . 6
3 clc ;
4 // Given
5 c =3 D +8; //m/ s
6 le =25; //m
7 zl =40+( %i *30) ; //ohm
8 f =10 D +6; // hz
21
Figure 1.5: Lossless line
22
9 cap =40 D -12; //F
10 l =300 D -9; //H/m
11
12 // I n p u t i m p e d a n c e
13 z0 = sqrt ( l / cap ) ; //ohm
14 zl1 = zl / z0 ;
15 lam = c / f ; //m
16 bet =(2* %pi * le ) / lam ; // r a d
17 zin =(( zl1 * cos ( bet ) ) +( %i * sin ( bet ) ) ) /( cos ( bet ) +( %i * zl1
* sin ( bet ) ) ) ; //ohm
18 disp ( ’ ohm ’ ,zin , ’ I n p u t i m p e d a n c e : ’ ) ;
19
20 // R e f l e c t i o n c o e f f i c i e n t
21 t =( zl1 -1) /( zl1 +1) ;
22 disp (t , ’ R e f l e c t i o n c o e f f i c i e n t : ’ ) ;
Scilab code Exa 1.7 Lossy cable
2 // Example 1 . 7
3 clc ;
4 // Given
5 c =3 D +8; //m/ s
6 R =2.25; //ohm
7 L =1 D -9; //H/m
8 C =1 D -12; //F/m
9 f =0.5 D +9; // hz
10 G =0;
11 w =2* %pi * f ; // r a d / s e c
12
13 // C h a r a c t e r s t i c i m p e d a n c e
14 z0 = sqrt (( R +( %i * w * L ) ) /( G +( %i * w * C ) ) ) ; //ohm
15 disp ( ’ ohm ’ ,z0 , ’ C h a r a c t e r s t i c i m p e d a n c e : ’ ) ;
23
Figure 1.6: Lossy cable
24
16
17 // P r o p a g a t i o n c o n s t a n t
18 gam = sqrt (( R +( %i * w * L ) ) *( G +( %i * w * C ) ) ) ;
19 disp ( gam , ’ P r o p a g a t i o n c o n s t a n t : ’ ) ;
Scilab code Exa 1.8 Transmission Line
2 // Example 1 . 8
3 clc ;
4 // Given
5 c =3 D +8; //m/ s
6 f =3 D +9; // Hz
7 ZL =50 -( %i *100) ; // ohms
8 Z0 =50; //ohm
9 // Wavelength
10 lam = c / f ;
11 disp ( ’ cm ’ , lam *100 , ’ Wavelength : ’ ) ;
12
13 // N o r m a l i z e d l o a d i m p e d a n c e
14 z = ZL / Z0 ;
15 disp (z , ’ N o r m a l i z e d l o a d i m p e d a n c e : ’ ) ;
16
17 // From c h a r t
18 zin =0.45+( %i *1.2) ;
19 yin =0.27 -( %i *0.73) ;
20 ZINN = Z0 * zin ;
21 disp ( ’ ohm ’ , ZINN , ’ L i n e i m p e d a n c e : ’ ) ;
22 YINN = yin / Z0 ;
23 disp ( ’ mho ’ , YINN , ’ L i n e a d m i t t a n c e : ’ ) ;
25
Figure 1.7: Transmission Line
26
27
2 // Example 1 . 9
3 clc ;
4 // Given
5 ZL =75+( %i *100) ; // ohms
6 Z0 =50; //ohm
7
9 z = ZL / Z0 ;
11
12 // ( i ) 0 . 0 5 1 ∗ lam
14 r =4.6;
15 Zi1 = r * Z0 ;
16 disp ( ’ ohm ’ ,Zi1 , ’ I n p u t i m p e d a n c e a t 0 . 0 5 1 lam : ’ ) ;
17
18 // ( i i ) 0 . 1 0 2 ∗ lam
19 r1 =1.5 -( %i *2) ;
20 Zi2 = r1 * Z0 ;
22
23 // ( i i i ) 0 . 3 0 1 ∗ lam
24 r2 =0.22;
25 Zi3 = r2 * Z0 ;
28
29
2 // Example 1 . 1 0
3 clc ;
4 // Given
5 ZL =15+( %i *20) ; // ohms
6 Z0 =50; //ohm
7
9 z = ZL / Z0 ;
11
13 T =0.6;
14 disp (T , ’ R e f l e c t i o n c o e f f i c i e n t : ’ ) ;
15
16 //VSWR
17 p =4;
18 disp (p , ’VSWR: ’ ) ;
2 // Example 1 . 1 1
3 clc ;
4 // Given
5 Z0 =50; //ohm
6 p =2.4;
7
8 // From c h a r t
9 zl =1.4+ %i ;
10 L = Z0 * zl ;
11 disp ( ’ ohm ’ ,L , ’ Load : ’ ) ;
30
Scilab code Exa 1.12 Active Device
2 // Example 1 . 1 2
3 clc ;
4 // Given
5 Z0 =50; //ohm
6 T =2.23;
7
8 // From c h a r t
9 zl =2+ %i ;
10 ZLd = Z0 * zl ;
11 disp ( ’ ohm ’ ,ZLd , ’ N o r m a l i z e d i m p e d a n c e : ’ ) ;
31
Figure 1.11: Active Device
12
13 // Impedance o f d e v i c e i s by n e g a t i n g t h e r e a l part
14 imp = - real ( ZLd ) +( imag ( ZLd ) * %i ) ;
15 disp ( ’ ohm ’ ,imp , ’ Impedance o f d e v i c e : ’ ) ;
Scilab code Exa 1.13 Transmission line
2 // Example 1 . 1 3
3 clc ;
4 // Given
5 p =3;
6 m1 =54; //cm
32
Figure 1.12: Transmission line
33
7 m2 =204; //cm
8
9 // P o i n t A
10 disp ( ’ P o i n t A ’ ) ;
11 lam =4*( m2 - m1 ) ;
12 dA =0.083* lam ;
13 L = m1 - dA ;
14 disp ( ’ cm ’ ,L , ’ L o c a t i o n o f s t u b : ’ ) ;
15 IA =0.114* lam ;
16 disp ( ’ cm ’ ,IA , ’ Length : ’ ) ;
17
18 // P o i n t B
19 disp ( ’ P o i n t B ’ ) ;
20 dB =0.083* lam ;
21 IB =0.386* lam ;
22 Lb = dB + m1 ;
23 disp ( ’ cm ’ ,Lb , ’ L o c a t i o n o f s t u b : ’ ) ;
2 // Example 1 . 1 5
3 clc ;
4 // Given
5 Z0 =50; //ohm
6 ZL =100; // ohms
7 f =10 D +9; // Hz
8 c =0.159 D -12; //F
9
11 z = ZL / Z0 ;
13
34
35
Figure 1.14: EM Plane
15 zin =0.4+( %i *0.55) ;
16 ZINN = Z0 * zin ;
17 disp ( ’ ohm ’ , ZINN , ’ N o r m a l i z e d i m p e d a n c e : ’ ) ;
Scilab code Exa 1.16 EM Plane
2 // Example 1 . 1 6
36
3 clc ;
4 // From g i v e n wave e q u a t i o n we can s e e
5 w =1 D +9; // r a d / s e c
6 bet =30; // r a d /m
7 c =3 D +8; //m/ s
8 u0 =1; // l e t
9 e0 =1/(9 D +16) ;
10
11 vp = w / bet ; //m/ s e c
12 disp ( ’m/ s ’ ,vp , ’ Phase v e l o c i t y : ’ ) ;
13
14 e =1/( vp ^2* u0 ) ;
15 er = e /( e0 * u0 ) ;
16 disp ( er , ’ D i e l e c t r i c c o n s t a n t : ’ ) ;
Scilab code Exa 1.17 Polyethylene
2 // Example 1 . 1 7
3 clc ;
4 // Given
5 c =3 D +8; //m/ s
6 f =10 D +9; // hz
7 er =6;
8 tandel =2 D -4;
9
10 vp = c / er ; //m/ s e c
11 disp ( ’m/ s e c ’ ,vp , ’ Phase v e l o c i t y : ’ ) ;
12 al =( %pi * f * tandel ) / vp ; //Np/m
13 disp ( ’ Np/m ’ ,al , ’ A t t e n u a t i o n c o n s t a n t : ’ ) ;
14
15 // Answer f o r v e l o c i t y i s c a l c u l a t e d wrong i n book ,
h e n c e a n s w e r s d o n t match f o r b o t h
37
Figure 1.15: Polyethylene
38
Figure 1.16: Electromagnetic wave
Scilab code Exa 1.18 Electromagnetic wave
2 // Example 1 . 1 8
3 clc ;
4 // Given
5 er =2.2;
6 n0 =377; //ohm
7 n2 = n0 / sqrt ( er ) ; //ohm
8 n1 =377; //ohm
9
39
10 // R e f l e c t i o n c o e f f i c i e n t
11 t =( n2 - n1 ) /( n2 + n1 ) ;
12 disp (t , ’ R e f l e c t i o n c o e f f i c i e n t : ’ ) ;
13
14 // Vswr
15 // Taking mod o f r e f l e c t i o n c o e f f i c i e n t
16 t1 = - t ;
17 p =(1+ t1 ) /(1 - t1 ) ;
18 disp (p , ’VSWR: ’ ) ;
Scilab code Exa 1.19 Range in sea water
2 // Example 1 . 1 9
3 clc ;
4 // Given
5 sig =5; //mohm/m
6 er =80*8.85 D -12;
7 eaz =0.1;
8 u =1.26 D -6;
9
10 az = - log (0.1) ;
11 // ( i ) Range a t 25 Khz
12 f =25 D +3; // Khz
13 w =2* %pi * f ; // r a d / s e c
14 a = w *( sqrt (( u * er /2) *( sqrt ( sig ^2/( w ^2* er ^2) +1) -1) ) ) ;
15 z = az / a ; //m
16 disp ( ’m ’ ,z , ’ Range a t 25 khz : ’ ) ;
17
18 // ( i i ) Range a t 25Mhz
19 f1 =25 D +6; //Mhz
20 w1 =2* %pi * f1 ; // r a d / s e c
21 a1 = w1 *( sqrt (( u * er /2) *( sqrt ( sig ^2/( w1 ^2* er ^2) +1) -1) ) )
40
Figure 1.17: Range in sea water
41
;
22 z1 = az / a1 ; //m
23 disp ( ’m ’ ,z1 , ’ Range a t 25Mhz : ’ ) ;
42
Chapter 2
Waveguides
Scilab code Exa 2.1 Dominant mode
1 // Page Number : 9 1
2 // Example 2 . 1
3 clc ;
4 // Given ,
5
6 a =6; //cm
7 b =4; //cm
8 d =4.47; //cm
9 c =3 D +8; //m/ s
10 lamc =2* a ;
11 lamg =2* d ;
12
13 // S i g n a l w a v e l e n g t h
14 lam = lamg * lamc /( sqrt ( lamg ^2+ lamc ^2) ) ;
15 lam = lam /100; //m
16 f = c / lam ;
17 disp ( ’ Ghz ’ ,f /1 D +9 , ’ S i g n a l f r e q u e n c y o f dominant mode
: ’ );
43
Figure 2.1: Dominant mode
44
Figure 2.2: Rectangular Waveguide
Scilab code Exa 2.2 Rectangular Waveguide
2 // Example 2 . 2
3 clc ;
4 // Given ,
45
5 c =3 D +8; //m/ s
6 a =2.5; //cm
7 b =5; //cm
8 lam =4.5; //cm
9
10 lamc =2* b ;
11
12 // Guide w a v e l e n g t h
13 lamg = lam /( sqrt (1 -(( lam / lamc ) ^2) ) ) ;
14 disp ( ’ cm ’ , lamg , ’ Guide w a v e l e n g t h : ’ ) ;
15
16 // Phase c o n s t a n t
17 bet =(2* %pi ) / lamg ;
18 bet = bet *100; // r a d /m
19 disp ( ’ r a d /m ’ ,bet , ’ Phase c o n s t a n t : ’ ) ;
20
21 // Phase v e l o c i t y
22 w =(2* %pi * c ) / lam ;
23 vp = w / bet ;
2 // Example 2 . 3
3 clc ;
4 // Given ,
5
6 c =3 D +8; //m/ s
7 a =4; //cm
8 b =2; //cm
9 f =10 D +9; // Hz
10 m =1;
46
47
11 n =1;
12
13
14 // C u t o f f w a v e l e n g t h
15 lamc =2/ sqrt (( m / a ) ^2+( n / b ) ^2) ;
16 disp ( ’ cm ’ , lamc , ’ Cut− o f f w a v e l e n g t h : ’ ) ;
17
18 //Wave i m p e d a n c e
19 lam = c / f ; //m
20 lam = lam *100; //cm
21 eeta =120* %pi ;
22 z0 = eeta * sqrt (1 -( lam / lamc ) ^2) ;
23 disp ( ’ ohm ’ ,z0 , ’ Wave i m p e d a n c e : ’ ) ;
Scilab code Exa 2.4 Wider dimension
2 // Example 2 . 4
3 clc ;
4 // Given ,
5 c =3 D +8; //m/ s
6 f =10 D +9; // Hz
7 zte =410; //ohm
8
9 // Wider d i m e n s i o n
10 lam = c / f ; //m
11 lam = lam *100; //cm
12 a =3/(2*( sqrt (1 -(120* %pi / zte ) ^2) ) ) ;
13 disp ( ’ cm ’ ,a , ’ Wider d i m e n s i o n : ’ ) ;
48
Figure 2.4: Wider dimension
49
50
Scilab code Exa 2.5 Rectangular waveguide
2 // Example 2 . 5
3 clc ;
4 // Given ,
5 c =3 D +8; //m/ s
6 a =3.0; //cm
7 b =1.5; //cm
8 mur =1;
9 er =2.25;
10 x = mur * er ;
11
12 // ( i ) C u t o f f w a v e l e n g t h and f r e q u e n c u y
13 disp ( ’ TE10 mode ’ ) ;
14 m1 =1;
15 n1 =0;
16 lamc10 =2/ sqrt (( m1 / a ) ^2+( n1 / b ) ^2) ;
17 disp ( ’ cm ’ , lamc10 , ’ Cut− o f f w a v e l e n g t h : ’ ) ;
18 lamc10 = lamc10 /100;
19 f10 = c /( lamc10 * sqrt ( x ) ) ;
20 disp ( ’ Ghz ’ , f10 /1 D +9 , ’ C u t o f f f r e q u e n c y : ’ ) ;
21
22 disp ( ’ TE20 mode ’ ) ;
23 m2 =2;
24 n2 =0;
25 lamc20 =2/ sqrt (( m2 / a ) ^2+( n2 / b ) ^2) ;
27 lamc20 = lamc20 /100;
28 f20 = c /( lamc20 * sqrt ( x ) ) ;
30
31 disp ( ’ TE11 mode ’ ) ;
32 m3 =1;
51
33 n3 =1;
34 lamc11 =2/ sqrt (( m3 / a ) ^2+( n3 / b ) ^2) ;
36 lamc11 = lamc11 /100;
37 f11 = c /( lamc11 * sqrt ( x ) ) ;
39
40 // ( i i ) lambg and Z0
41 f =4 D +9; // Hz
42 lam = c / f ;
43 lamg = lam /( sqrt (x -(( lam / lamc10 ) ^2) ) ) ;
44 disp ( ’ cm ’ , lamg *100 , ’ Guide w a v e l e n g t h : ’ ) ;
45
46 fc =3.33 D +9; // Hz
47 Z0 =(120* %pi *(1/ sqrt ( x ) ) *( b / a ) ) / sqrt (1 -(( fc / f ) ^2) ) ;
48 disp ( ’ ohm ’ , round ( Z0 ) , ’ Impedance : ’ ) ;
2 // Example 2 . 5
3 clc ;
4 // Given ,
5 c =3 D +8; //m/ s
6 a =4; //cm
7 b =2; //cm
8
9 // ( i ) Mode
10 lamc =2* a ; //cm
11 lamcm = lamc /100; //m
12 fc = c / lamcm ;
13 // 20% a b o v e f c
14 f =1.2* fc ; // Hz
52
Figure 2.6: Rectangular waveguide
53
15
16 // O p e r a t i n g w a v e l e n g t h
17 lam1 = c / f ; //cm
18
19 // For TE10 mode
20 lamc10 =2* b ; //cm
21 lamcm10 = lamc10 /100; //m
22 fc10 = c / lamcm10 ;
23 disp ( ’ Hence mode o f o p e r a t i o n i s TE10 ’ , ’ Hz ’ ,fc , ’
S i n c e gu i d e i s o p e r a t i n g at ’ );
24
25 // ( i i ) Guide w a v e l e n g t h
26 lamm1 = lam1 *100; //cm
27 lamg = lamm1 /( sqrt (1 -( lamm1 / lamc ) ^2) ) ;
29
30 // ( i i i ) Phase v e l o c i t y
31 vp = f * lamg ;
32 disp ( ’m/ s ’ , vp /100 , ’ Phase v e l o c i t y : ’ ) ;
33
34 // ( i i i ) Group v e l o c i t y
35 vg = c ^2/ vp ;
36 disp ( ’m/ s ’ ,vg , ’ Group v e l o c i t y : ’ ) ;
Scilab code Exa 2.8 Lossless Rectangular Waveguide
2 // Example 2 . 8
3 clc ;
4 // Given ,
5 c =3 D +8; //m/ s
6 a =7; //cm
7 b =3.5; //cm
54
Figure 2.7: Lossless Rectangular Waveguide
55
8 f =3 D +9; // Hz
9 h0 =10; //amp/m
10
11 //Wave i m p e d a n c e
12 lamc =2* a ;
13 lam = c / f ; //m
14 lam = lam *100; //cm
15 lamg = lam / sqrt (1 -( lam / lamc ) ^2) ; //cm
16 z0 =377* lamg / h0 ; //ohm
17
18 a1 = a /100; //m
19 b1 = b /100; //m
20 // A v e r a g e power t r a n s m i t t e d
21 p =( z0 * h0 * h0 * a1 * b1 ) /4;
22 disp ( ’W’ ,p , ’ A v e r a g e power t r a n s m i t t e d : ’ ) ;
23
24 // Peak e l e c t r i c f i e l d
25 e0 = z0 * h0 ;
26 disp ( ’ kV/m ’ , e0 /1000 , ’ Peak e l e c t r i c f i e l d : ’ );
27
28 // Answer f o r p i s g i v e n a s 2 8 . 3 W but i t s h o u l d be
3 2 . 9 9W
Scilab code Exa 2.9 Dimensions
2 // Example 2 . 9
3 clc ;
4 // Given ,
5 c =3 D +8; //m/ s
6 fc =3 D +9; // Hz
7
8 // C u t o f f w a v e l e n g t h
56
Figure 2.8: Dimensions
57
9 lamc = c / fc ;
10 a = lamc /2; //m
11 a = a *100; //cm
12 disp ( ’ D i m e n s i o n s : ’ ) ;
13 disp ( ’ cm ’ ,a , ’ a : ’ ) ;
14 b = a /2; //cm
15 disp ( ’ cm ’ ,b , ’ b : ’ ) ;
1 // Page Number :
2 // Example 2 . 1 0
3 clc ;
4 // Given ,
5
6 c =3 D +8; //m/ s
7 a =3; //cm
8 a1 = a /100; //m
9 b =2; //cm
10 b1 = b /100; //m
11 f =7.5 D +9; //HZ
12 p =5 D +3; //W
13
14 mu = %pi *4 D -7;
15 w =2* %pi * f ;
16 bet = sqrt ((( w / c ) ^2) -(( %pi / a1 ) ^2) ) ;
17 // C h a r e c t e r i s t i c i m p e d a n c e
18 z0 = w * mu *2* b /( bet * a ) ;
19 disp ( ’ ohm ’ ,z0 , ’ C h a r e c t e r i s t i c i m p e d a n c e : ’ ) ;
20
21 // Peak e l e c t r i c f i e l d
22 e0 =4* w * mu * p /( bet * a * b ) ;
23 disp ( ’V/m ’ ,e0 , ’ Peak e l e c t r i c f i e l d : ’ );
58
59
24
25 //Maximum v o l t a g e
26 v0 = e0 * b1 ;
27 disp ( ’ kV ’ , v0 /1000 , ’ Maximum v o l t a g e : ’ ) ;
28
29 // Answer f o r v0 i s g i v e n a s 3 . 1 7 2 kV i t s h o u l d be
3 3 . 7 1 kV
Scilab code Exa 2.14 Waveguide
2 // Example 2 . 1 4
3 clc ;
4 // Given ,
5 c =3 D +8; //m/ s
6 a =1.5; //cm
7 a1 = a /100; //m
8 b =0.8; //cm
9 b1 = b /100; //m
10 mu =1/ c * c ;
11 e =4;
12 w = %pi *1 D +11;
13 n =377;
14
15 // ( i ) F r e q u e n c y o f o p e r a t i o n
16 f = w /(2* %pi ) ;
17 f1 = f /1 D +9; // ghz
18 disp ( ’ Ghz ’ ,f1 , ’ F r e q u e n c y o f o p e r a t i o n : ’ ) ;
19
20 // ( i i ) Cutt o f f f r e q u e n c y
21 fc =( c * sqrt ((1/ a1 ) ^2+(3/ b1 ) ^2) ) /(2* sqrt ( e ) ) ;
22 fc1 = fc /1 D +9; // ghz
23 disp ( ’ Ghz ’ ,fc1 , ’ Cut o f f f r e q u e n c y : ’ ) ;
60
Figure 2.10: Waveguide
61
24
25 // ( i i i ) Phase c o n s t a n t
26 bet =( w * sqrt ( e ) * sqrt (1 -( fc / f ) ^2) ) /( c ) ;
27 disp ( ’ r a d /m ’ ,bet , ’ Phase c o n s t a n t : ’ ) ;
28
29 // ( i v ) P r o p o g a t i o n c o n s t a n t
30 gam = %i * bet ;
31 disp ( ’ r a d / s ’ ,gam , ’ P r o p o g a t i o n c o n s t a n t : ’ ) ;
32
33 // ( v ) I n t r e n s i c wave i m p e d a n c e
34 zte =( n / sqrt ( e ) ) / sqrt (1 -( fc / f ) ^2) ;
35 ztm =( n / sqrt ( e ) ) * sqrt (1 -( fc / f ) ^2) ;
36 disp ( ’Ohm ’ ,ztm , ’ZTM13 ’ , ’Ohm ’ ,zte , ’ ZTE13 ’ , ’ I n t r i n s i c
wave i m p e d a n c e : ’ ) ;
Scilab code Exa 2.17 Air filled Rectangular Waveguide

2 // Example 2 . 1 7
3 clc ;
4 // Given
5 a =2; //cm
6 a1 =1/100; //m
7 b =1; //cm
8 b1 = b /100; //m
9 p =10 D -3; //W
10 c =3 D +8; //m/ s
11 f0 =10 D +9; // Hz
12
13 // Peak v a l u e o f e l e c t r i c f i e l d
14 fc = c /(2* a ) ;
15 E02 =(4* p *377) /( a1 * b1 * sqrt (1 -( fc / f0 ) ^2) ) ;
16 E0 = sqrt ( E02 ) ;
62
Figure 2.11: Air filled Rectangular Waveguide
63
17 disp ( ’V/m ’ ,E0 , ’ Peak v a l u e o f e l e c t r i c f i e l d : ’ ) ;
18
19 //Maximum power t r a n s m i t t e d
20 Ed =3 D +6; //V/m
21 Pt =2.6 D +13*( Ed / f0 ) ^2;
22 disp ( ’kW ’ , Pt /1000 , ’ Maximum power t r a n s m i t t e d : ’ ) ;
23
24 // Answer i s g i v e n a s 2 3 0 0kW but i t i s 2 3 4 0kW

2 // Example 2 . 1 8
3 clc ;
4 // Given
5 f =5 D +9; // Hz
6 c =3 D +8; //m/ s
7 a =7.5; //cm
8 a1 = a /100; //m
9 b =3.5; //cm
10 b1 = b /100; //m
11 lam = c / f ;
12 lamm = lam *100; //m
13
14 disp ( ’ TE10 mode ’ ) ;
15 lamc10 =2* a ;
16 bet10 =(2* %pi * sqrt ((( lamc10 / lamm ) ^2) -1) ) / lamc10 ;
17 disp ( ’ r a d /cm ’ , bet10 , ’ P r o p o g a t i o n c o n s t a n t : ’ ) ;
18 vp10 =(2* %pi * f ) / bet10 ;
19 disp ( ’m/ s ’ , vp10 /100 , ’ Phase v e l o c i t y : ’ ) ;
20
21 disp ( ’ TE01 mode ’ ) ;
22 lamc01 =2* b ;
64
65
25 vp01 =(2* %pi * f ) / bet01 ;
27
28 disp ( ’ TE11 mode ’ ) ;
29 lamc11 =(2* a * b ) / sqrt (( a * a ) +( b * b ) ) ;
32 vp11 =(2* %pi * f ) / bet11 ;
34
35 disp ( ’ TE02 mode ’ ) ;
36 lamc02 = b ;
39 disp ( ’ As b e t a i s i m a g i n a r y , mode g e t s a t t e n u a t e d ’ ) ;
40 alp =(2* %pi * sqrt (1 -(( lamc02 / lamm ) ^2) ) ) / lamc02 ;
41 disp ( ’ Np/m ’ ,alp , ’ P r o p o g a t i o n c o n s t a n t a l p h a : ’ ) ;

2 // Example 2 . 1 9
3 clc ;
4 // Given
5 c =3 D +8; //m/ s
6 a =2.29; //cm
7 b =1.02; //cm
8 a1 = a /100 ; //m
9 b1 = b /100; //m
10 f =6 D +9; // Hz
11 e =1;
66
67
12 mu =1/( c ^2) ;
13
14 // Cut o f f f r e q u e n c y
15 lamc =2* a1 ;
16 fc = c / lamc ;
17 w =2* %pi * fc ;
18
19 // A t t e n u a t i o n c o n s t a n t
20 a =( w * sqrt (1 -(( f / fc ) ^2) ) ) / c ;;
21 adb = -20* log10 ( exp ( - a ) ) ;
22 disp ( ’ dB/m ’ ,adb , ’ A t t e n u a t i o n c o n s t a n t : ’ ) ;
Scilab code Exa 2.20 Ratio of cross section

2 // Example 2 . 2 0
3 clc ;
4 // Given ,
5 a1 =1.84;
6 a2 = %pi ;
7
8 r =2* %pi *( a1 / a2 ) ^2;
9 disp (r , ’ C r o s s s e c t i o n r a t i o : ’ ) ;

2 // Example 2 . 2 1
3 clc ;
68
Figure 2.14: Ratio of cross section
69
70
4 // Given
5 c =3 D +8; //m/ s
6 f =15 D +9; // hz
7 a =1.07; //cm
8 a1 = a /100; //m
9 b =0.43; //cm
10 b1 = b /100; //m
11 er =2.08;
12 tandel =0.0004;
13 lam = c / f ;
14
15
16 // ( i ) Cut o f f f r e q u e n c y
17 m1 =1;
18 n1 =0;
19 fc10 =( c /(2* %pi * sqrt ( er ) ) * sqrt (( m1 * %pi / a1 ) ^2+( n1 * %pi /
b1 ) ^2) ) ;
20 disp ( ’GHz ’ , fc10 /10^9 , ’ Cut o f f f r e q u e n c y f o r mode
TE10 : ’ ) ;
21
22 m2 =2;
23 n2 =0;
b1 ) ^2) ) ;
25 disp ( ’ Ghz ’ , fc20 /10^9 , ’ Cut o f f f r e q u e n c y a t mode TE20
: ’ );
26
27 m3 =0;
28 n3 =1;
b1 ) ^2) ) ;
30 disp ( ’ Ghz ’ , fc01 /10^9 , ’ Cut o f f f r e q u e n c y a t mode TE01
: ’ );
31
32 // D i e l e c t r i c a t t e n u a t i o n c o n s t a n t
33 ad =( %pi * tandel ) /( lam * sqrt (1 -( fc10 / f ) ^2) ) ;
34 adb = -20* log10 ( exp ( - ad ) ) ;
35 disp ( ’ dB/m ’ ,adb , ’ A t t e n u a t i o n c o n s t a n t : ’ ) ;
71

2 // Example 2 . 2 2
3 clc ;
4 // Given
72
5 c =3 D +8; //m/ s
6 a =2.286; //cm
7 a1 = a /100; //m
8 b =1.016; //cm
9 b1 = b /100; //m
10 sig =5.8 D +7; // s /m
11 f =9.6 D +9; // Hz
12
13 w =2* %pi * f ;
14 mu = %pi *4 D -7;
15 et =377;
16
17 lam = c / f ;
18 lamc =2* a1 ;
19 r = lam / lamc ;
20
21 Rs = sqrt (( w * mu ) /(2* sig ) ) ;
22 ac =( Rs *(1+(2*( b1 / a1 ) * r * r ) ) ) /( et * b1 * sqrt (1 -( r ^2) ) ) ;
23 adb = -20* log10 ( exp ( - ac ) ) ;
24 disp ( ’ dB/m ’ ,adb , ’ C o n d u c t o r a t t e n u a t i o n c o n s t a n t : ’ ) ;
Scilab code Exa 2.23 Circular waveguide

2 // Example 2 . 2 3
3 clc ;
4 // Given
5 c =3 D +8; //m/ s
6 f =9 D +9; // hz
7 a =5; //cm
8 a1 = a /100; //m
9 e =1;
10 mu =1/( c * c ) ;
73
Figure 2.17: Circular waveguide
74
11 p11 =1.841;
12
13 fc =( p11 * c ) /(2* %pi * a1 ) ;
14 //Maximum power t r a n s m i t t e d
15 pmax =1790*( a1 * a1 ) * sqrt (1 -(( fc / f ) ^2) ) ;
16 disp ( ’kW ’ , pmax , ’ Maximum power t r a n s m i t t e d : ’ ) ;
Scilab code Exa 2.24 Air filled circular waveguide

2 // Example 2 . 2 6
3 clc ;
4 // Given
5 c =3 D +8; //m/ s
6 a =5; //cm
7 a1 = a /100; //m
8 f =3 D +9; // hz
9 p11 =1.841;
10 e =1;
11 w =2* %pi * f ;
12
14 fc =( p11 * c ) /(2* %pi * a1 ) ;
15 disp ( ’ Ghz ’ , fc /10^9 , ’ Cut o f f f r e q u e n c y : ’ ) ;
16
17 // ( i i ) Guide w a v e l e n g t h
18 bet = sqrt ((( w * w ) /( c * c ) ) -(( p11 / a1 ) ^2) ) ;
19 lamg =(2* %pi ) / bet ;
20 lamg1 = lamg *100; //cm
21 disp ( ’ cm ’ , lamg1 , ’ Guide w a v e l e n g t h : ’ ) ;
22
23 // ( i i i ) Wave i m p e d a n c e
24 zte =( w * %pi *4 D -7) / bet ;
75
Figure 2.18: Air filled circular waveguide
76
25 disp ( ’ ohm ’ , round ( zte ) , ’ Wave i m p e d a n c e : ’ ) ;
Scilab code Exa 2.25 Air filled rectangular waveguide
1 // Page Number : 1 0 8
2 // Example 2 . 2 5
3 clc ;
4 // Given
5 c =3 D +8; //m/ s
6 p01 =2.405;
7 a =1/100;; //cm
8 p11 =1.841;
9
10 fc01 =(( c * p01 ) /(2* %pi * a ) ) ;
11 fc11 =(( c * p11 ) /(2* %pi * a ) ) ;
12 bw = fc01 - fc11 ;
13 disp ( ’ Ghz ’ , bw /10^9 , ’ Bandwidth : ’ ) ;

2 // Example 2 . 2 6
3 clc ;
4 // Given
5 c =3 D +8; //m/ s
6 a =2.286; //cm
7 f =5 D +9; // Hz
8 er =2.25;
9 tandel =1 D -3;
77
Figure 2.19: Air filled rectangular waveguide
78
79

10 w =2* %pi * f ;
11 mu =4 D -7;
12 sig =5.8 D +7; // s /m
13
14 lamc =2* a ;
15 lamm = c / f ; //m
16 lam = lamm *100; //cm
17
18 ermax =( lam / a ) ^2;
19 disp ( ermax , ’ Maximum v a l u e o f d i e l e c t r i c c o n s t a n t : ’ ) ;
20 ermin =( lam /(2* a ) ) ^2;
21 disp ( ermin , ’ Minimum v a l u e o f d i e l e c t r i c c o n s t a n t : ’ ) ;
22
23 // Guide w a v e l e n g t h
24 lam1 = lam / sqrt ( er ) ; //cm
25 lamg = lam1 / sqrt (1 -( lam1 / lamc ) ^2) ;
27
28 lamm1 = lam1 /100;
29 ad =( %pi / lamm1 ) *( tandel / sqrt (1 -( lam1 / lamc ) ^2) ) ;
30 disp ( ’ Np/m ’ ,ad , ’ ad : ’ ) ;
31 bet =2* %pi / lamg ;
32 disp ( ’ r a d /cm ’ ,bet , ’ Beta : ’ ) ;
33 vp = w /( bet *100) ;
Scilab code Exa 2.27 Circular Waveguide

2 // Example 2 . 2 7
3 clc ;
4 // Given
5 c =3 D +8; //m/ s
80
81
Figure 2.21: Circular Waveguide

6 a =0.5; //cm
7 a1 = a /100; //m
8 f =14 D +9; // Hz
9 er =2.08;
10 p11 =1.841;
11 p01 =2.405;
12 tandel =4 D -4;
13 w =2* %pi * f ;
14 u = %pi *4 D -7;
15 sig =4.1 D +7;
16 et =377;
17
18 // ( i ) Cut o f f f r e q u e n c i e s
19 fcte11 = p11 * c /(2* %pi * a1 * sqrt ( er ) ) ;
20 fctm01 = p01 * c /(2* %pi * a1 * sqrt ( er ) ) ;
21 disp ( ’ Ghz ’ , fcte11 /10^9 , ’ Cut o f f f r e q u e n c i e s f o r TE11
mode : ’ ) ;
22 disp ( ’ Ghz ’ , fctm01 /10^9 , ’ Cut o f f f r e q u e n c i e s f o r TM01
mode : ’ ) ;
23
24 // ( i i ) O v e r a l l n o i s e
25 // D i e l e c t r i c a t t e n u a t i o n
26 ad =( %pi * sqrt ( er ) * tandel * f ) /( c * sqrt (1 -(( fcte11 / f ) ^2) )
);
27 disp ( ’ dB/m ’ , ad *8.686 , ’ D i e l e c t r i c a t t e n u a t i o n : ’ ) ;
28
29 // C o n d u c t o r a t t e n u a t i o n
30 k =(2* %pi * f * sqrt ( er ) ) / c ;
31 bet = sqrt (( k * k ) -(( p11 / a1 ) ^2) ) ;
32 // S u r f a c e r e s i s t a n c e
33 rs = sqrt (( w * u ) /(2* sig ) ) ;
34 kc2 =( p11 / a1 ) ^2;
35
36 ac =( rs *( kc2 -(( k ^2) /(( p11 ^2) -1) ) ) ) /( a1 * k * et * bet ) ;
37 disp ( ’ dB/m ’ , ac *8.686 , ’ C o n d u c t o r a t t e n u a t i o n : ’ ) ;
38
39 // T o t a l a t t e n u a t i o n
40 a =( ac + ad ) *8.686;
82
41 disp ( ’ dB/m ’ ,a , ’ T o t a l a t t e n u a t i o n : ’ ) ;
42 ta = a *0.3;
43 disp ( ’ dB ’ ,ta , ’ T o t a l a t t e n u a t i o n i n 30 cm l i n e : ’ ) ;
44
45 // Answer f o r c o n d c u t o r a t t e n u a t i o n i s wrong i n book ,
hence answer f o r t o t a l l o s s i s different

2 // Example 2 . 2 8
3 clc ;
4 // Given
5 c =3 D +8; //m/ s
6 er =9;
7 a =7; //cm
8 a1 = a /100; //m
9 b =3.5; //cm
10 b1 = b /100; //m
11 ur =1;
12 f1 =2 D +9; // Hz
13
15 lamc =2* a1 ;
16 fc = c /( lamc * sqrt ( ur * er ) ) ;
17 disp ( ’ Ghz ’ , fc /10^9 , ’ Cut o f f f r e q u e n c y : ’ ) ;
18
19 // ( i i ) Phase v e l o c i t y
20 lam = c / f1 ; //m
21 lam1 = lam *100; //cm
22 lamc1 = lamc *100; //cm
23 lamg = lam1 /( sqrt (( ur * er ) -(( lamc1 / lam1 ) ^2) ) ) ; //cm
24 lamg1 = lamg /100; //m
83
84
25 vp = f1 * lamg1 ;
27
28 // / ( i i i ) Guide w a v e l e n g t h
Scilab code Exa 2.29 Circular waveguide

2 // Example 2 . 2 9
3 clc ;
4 // Given
5 c =3 D +8; //m/ s
6 fc =9 D +9; // Hz
7 er =1;
8 er1 =4;
9 p11 =1.841;
10
11 // ( i ) a i r f i l l e d
12 a =( p11 * c ) /(2* %pi * fc * sqrt ( er ) ) ;
13 disp ( ’ cm ’ ,a *100 , ’ I n s i d e d i a m e t e r i f a i r f i l l e d : ’ ) ;
14 // ( i i ) d i e l e c t r i c f i e l d
15 a1 =( p11 * c ) /(2* %pi * fc * sqrt ( er1 ) ) ;
16 disp ( ’ cm ’ , a1 *100 , ’ I n s i d e d i a m e t e r i f d i e l e c t r i c
f i l l e d : ’ );
17
18 // Answers a r e c a l c u l a t e d wrong i n book
Scilab code Exa 2.30 Cutoff frequencies
85
Figure 2.23: Circular waveguide
86
Figure 2.24: Cutoff frequencies
87
2 // Example 2 . 3 0
3 clc ;
4 // Given
5 c =3 D +8; //m/ s
6 er =2.55;
7 d =1; //mm
8 d1 = d /1000; //m
9
10 // Cut o f f f r e q u e n c i e s
11 fctm0 =0;
12 disp ( ’ Ghz ’ , fctm0 , ’ Cut o f f f r e q u e n c y f o r mode TM0 : ’ ) ;
13
14 fcte1 = c /(4* d1 * sqrt ( er -1) ) ;
15 disp ( ’ Ghz ’ , fcte1 /10^9 , ’ Cut o f f f r e q u e n c y a t mode TE1
: ’ );
16
17 fctm1 = c /(2* d1 * sqrt ( er -1) ) ;
18 disp ( ’ Ghz ’ , fctm1 /10^9 , ’ Cut o f f f r e q u e n c y a t mode TM1
: ’ );
19
20
21 // Answers a r e c a l c u l a t e d wrong i n book
Scilab code Exa 2.31 Dielectric constant

2 // Example 2 . 3 1
3 clc ;
4 // Given ,
5 c =3 D +8; //m/ s
6 f =15 D +9; // hz
7 d =5; //mm
88
Figure 2.25: Dielectric constant
89
8 d1 = d /1000; //m
9
10 // Cut o f f f r e q u e n c y
11 fc =0.8* f ;
12 // D i e l c t r i c c o n s t a n t
13 er =( c /(2* d1 * fc ) ) ^2+1;
14 disp ( er , ’ D i e l e c t r i c c o n s t a n t : ’ ) ;
90
Chapter 3
Microwave Network Analysis
Scilab code Exa 3.4 Scattering matrix

2 // Example 3 . 4
3 clc ;
4 // Given
5
6 [ z ]=[4 2;2 4];
7 [ I ]=[1 0;0 1];
8
9 // S c a t t e r i n g m a t r i x
10 [ s ]={[ z ] -[ I ]}* inv ({[ z ]+[ I ]}) ;
11 disp ([ s ] , ’ S c a t t e r i n g M a t r i x : ’ ) ;
Scilab code Exa 3.5 Network
91
Figure 3.1: Scattering matrix
92
Figure 3.2: Network
93
2 // Example 3 . 5
3 clc ;
4 // Given
5 P =12.8 D -3; //W
6 l =3; //cm
7 lamb =4.2; //cm
8 vswr =2.2;
9 jfi = %i *4.49;
10
11 // ap
12 ap = sqrt (2* P ) ;
13
14 // Phase s h i f t
15 bl =(2* %pi * l ) / lamb ;
16 // bp
17 bp =( ap *( vswr -1) ) /( vswr +1) ;
18
19 a = ap * exp ( jfi ) ;
20 b = bp * exp ( jfi ) ;
21 disp (a ,b , ’ R e q u i r e d Waves : ’ ) ;
Scilab code Exa 3.6 Microwave network

2 // Example 3 . 6
3 clc ;
4 // Given
5 S11 =0.10;
6 S12 =0.90;
7 A12 = -45;
8 S21 =0.90;
9 A21 =45;
10 S22 =0.3;
94
Figure 3.3: Microwave network
95
11
12 // ( i ) Network i s r e c i p r o c a l
13 if ( A12 == A21 )
14 disp ( ’ Network i s r e c i p r o c a l ’ ) ;
15 else
16
17 disp ( ’ Network i s n o t r e c i p r o c a l ’ ) ;
18 end
19
20 // ( i i ) Network i s l o s s l e s
21 x =( S11 ^2) +( S12 ^2) ;
22 if ( x ==1)
23 disp ( ’ Network i s l o s s l e s s ’ ) ;
24 else
25
26 disp ( ’ Network i s n o t l o s s l e s s ’ ) ;
27 end
28
29 // ( i i i ) Return l o s s
30 T = S11 -(( S12 * S21 ) /(1+ S22 ) ) ;
31 Tm = - T ; //mod o f T
32 L = -20* log10 ( Tm ) ;
33 disp ( ’ dB ’ ,L , ’ Return L o s s : ’ ) ;
Scilab code Exa 3.12 Transistor amplifier circuit

2 // Example 3 . 1 2
3 clc ;
4 // Given
5 S11 =0.6;
6 S12 =0.045;
7 S21 =2.5;
96
Figure 3.4: Transistor amplifier circuit
97
8 S22 =0.50;
9 TS =0.5;
10 TL =0.4;
11 Z0 =50; //ohm
12 Vrms =10; //V
13
14 // ( i ) Gain P a r a m e t e r s
15 // ( i ) R e f l e c t i o n c o e f f i c i e n t s o f i n p u t and o u t p u t
16 Tin = S11 +(( S12 * S21 * TL ) /(1 -( S22 * TL ) ) ) ;
17 Tout = S22 +(( S12 * S21 * TS ) /(1 -( S22 * TS ) ) ) ;
18
19 // T r a n s d u c e r Gain
20 x =(1 -( TS ) ^2) /((1 -( S11 * TS ) ) ^2) ;
21 y =( S21 * S21 ) ;
22 z =(1 -( TL ) ^2) /((1 -( Tout * TL ) ) ^2) ;
23 GT = x * y * z ;
24 disp ( GT , ’ T r a n s d u c e r Gain : ’ ) ;
25
26 // A v a i l a b l e Power Gain
27 z1 =1 -( Tout ) ^2;
28 GA =( x * y ) / z1 ;
29 disp ( GA , ’ A v a i l a b l e power Gain : ’ ) ;
30
31 // Power Gain
32 z2 =1 -( Tin ) ^2;
33 GP =( x * y ) / z2 ;
34 disp ( GP , ’ Power Gain : ’ ) ;
35
36 // ( i i ) Power l e v e l s
37 // Power a v a i l a b l e a t s o u r c e
38 Pavs =( sqrt (2) * Vrms ) ^2/(8* Z0 ) ;
39 disp ( ’W’ , Pavs , ’ Power a v a i l a b l e a t s o u r c e : ’ ) ;
40
41 Pl =9.4* Pavs ;
42 // Power a v a i l a b l e a t i n p u t
43 Pin = Pl /13.5;
44 disp ( ’W’ ,Pin , ’ Power a v a i l a b l e a t i n p u t : ’ ) ;
45
98
46 // ( i i i ) VSWRs
47 M1 = Pin / Pavs ;
48 M2 = Pl /(9.6* Pavs ) ;
49
50 Tin1 = sqrt (1 - M1 ) ;
51 Tout1 = sqrt (1 - M2 ) ;
52
53 vswrin =(1+ Tin1 ) /(1 - Tin1 ) ;
54 disp ( vswrin , ’ I n p u t VSWR: ’ ) ;
55 vswrout =(1+ Tout1 ) /(1 - Tout1 ) ;
56 disp ( vswrout , ’ Output VSWR: ’ ) ;
57
58 // C a l c u l a t i o n s f o r g a i n a r e done wrong i n book ,
h e n c e a n s w e r s d o n t match
99
Chapter 4
Microwave Resonators and

Waveguide Components
Scilab code Exa 4.1 Rectangular cavity resonator

2 // Example 4 . 1
3 clc ;
4 // Given
5 a =5; //cm
6 a1 = a /100; //m
7 b =2; //cm
8 b1 = b /100; //m
9 c =15; //cm
10 c1 = c /100; //m
11
12 // ( i ) A i r f i l l e d c a v i t y
13 m =1;
14 n =0;
15 p =1;
16 c =3 D +8; // f o r a i r
17 fr =(1/2) * c * sqrt (( m / a1 ) ^2+( n / b1 ) ^2+( p / c1 ) ^2) ; // hz
100
Figure 4.1: Rectangular cavity resonator
101
Figure 4.2: Rectangulr resonator
18 disp ( ’ Ghz ’ , fr /10^9 , ’ R e s o n a n t f r e q u e n c y f o r an a i r

f i l l e d c a v i t y : ’ );
19
20 // ( i i ) D i e l c t r i c f i l l e d c a v i t y
21 er =2.56;
22 fr1 =(1/2) *( c / sqrt ( er ) ) * sqrt (( m / a1 ) ^2+( n / b1 ) ^2+( p / c1 )
^2) ; // hz
23 disp ( ’ Ghz ’ , fr1 /10^9 , ’ R e s o n a n t f r e q u e n c y f o r
d i e l e c t r i c c a v i t y : ’ );
Scilab code Exa 4.2 Rectangulr resonator
102
2 // Example 4 . 2
3 clc ;
4 // Given
5 a =0.38; //cm
6 a1 = a /100; //m
7 b =0.76; //cm
8 b1 = b /100; //m
9 f =50 D +9;
10 c =3 D +8;
11
12 // Length f o r TE102
13 m =1;
14 n =0;
15 p =2;
16 l =1/ sqrt (( f / c ) ^2 -(1/(4* b1 ^2) ) ) ; //m
17 disp ( ’ cm ’ ,l *100 , ’ Length c : ’ ) ;
Scilab code Exa 4.3 X band resonator

2 // Example 4 . 3
3 clc ;
4 // Given
5 c =3 D +8; //m/ s
6 a =2.286; //cm
7 a1 = a /100; //m
8 b =1.024; //cm
9 b1 = b /100; //m
10 f =10 D +9; // hz
11 sig =6 D +7;
12 u =4 D -7* %pi ;
13 w =2* %pi * f ;
103
Figure 4.3: X band resonator
104
14 eet =377;
15
16 // S h o r t e s t c a v i t y l e n g t h
17 lamc =2* a1 ; //m
18 fc = c / lamc ; // hz
19 lam = c / f ; //m
20 lamg = lam / sqrt (1 -( fc / f ) ^2) ; //m
21 sc = lamg /2; //m
22 disp ( ’ cm ’ , sc *100 , ’ S h o r t e s t c a v i t y l e n g t h : ’ ) ;
23
24 //Qw o f t h e r e s o n a t o r o p e r a t i n g i n TE101 mode
25 rs = sqrt (( w * u ) /(2* sig ) ) ; //ohm
26 lamr = c / f ;
27 x =((( a1 * b1 ) /( sc ^2) ) +(( sc ^2+ a1 ^2) /(2* sc * a1 ) ) +( b1 * sc /
a1 ^2) ) ;
28 qw =(2* %pi * eet * a1 * b1 * sc ) /( rs *( lamr ^3) * x ) ;
29 disp ( qw , ’Qw o f t h e r e s o n a t o r o p e r a t i n g i n TE101 mode
’ );
Scilab code Exa 4.4 Rectangular resonator

2 // Example 4 . 4
3 clc ;
4 // Given
5 c =3 D +8; //m/ s
6 a =4.8; //cm
7 a1 = a /100; //m
8 b =2.2; //cm
9 b1 = b /100; //m
10 f =5 D +9; // hz
11 er =2.25;
12 tandel =4 D -4;
105
Figure 4.4: Rectangular resonator
106
13 sig =5.813 D +7;
14 oneby =3 D +8;
15 u =4 D -7* %pi ;
16 w =2* %pi * f ;
17 eet =377;
18
19 // Length a t p=1
20 m =1;
21 n =0;
22 p =1;
23 z =( f *2* sqrt ( er ) ) / c ;
24 cp1 = p / sqrt (( z ^2) -(( m / a1 ) ^2) -(( n / b1 ) ^2) ) ;
25 disp ( ’ cm ’ , cp1 *100 , ’ Length o f r e s o n a t o r a t p =1: ’ ) ;
26
27 // At p=2
28 cp2 = cp1 *2;
29 disp ( ’ cm ’ , cp2 *100 , ’ Length o f r e s o n a t o r a t p =2: ’ ) ;
30
31 //Qw
32 rs = sqrt (( w * u ) /(2* sig ) ) ; //ohm
33 lamr = c /( f * sqrt ( er ) ) ;
34 x =((( a1 * b1 ) /( cp1 ^2) ) +(( cp1 ^2+ a1 ^2) /(2* cp1 * a1 ) ) +( b1 *
cp1 / a1 ^2) ) ;
35 qw =(2* %pi *( eet / sqrt ( er ) ) * a1 * b1 * cp1 ) /( rs *( lamr ^3) * x ) ;
36 qd =1/ tandel ;
37 q =( qw * qd ) /( qw + qd ) ;
38 disp (q , ’Q f o r TE101 mode : ’ ) ;
Scilab code Exa 4.5 Cylindrical resonator

2 // Example 4 . 5
3 clc ;
107
Figure 4.5: Cylindrical resonator
108
4 // Given
5 c =3 D +8; //m/ s
6 a =2; //cm
7 a1 = a /100; //m
8 b =2.5; //cm
9 b1 = b /100; //m
10
11 disp ( ’TE modes ’ ) ;
12 h01 =3.832;
13 fr =( c /(2* %pi ) ) * sqrt (( h01 / a1 ) ^2+( %pi / b1 ) ^2) ; // hz
14 disp ( ’ Ghz ’ , fr /10^9 , ’ R e s o n a n t f r e q u e n c y f o r mode
TE010 : ’ ) ;
15
16 h11 =1.841;
17 fr1 =( c /(2* %pi ) ) * sqrt (( h11 / a1 ) ^2+( %pi / b1 ) ^2) ; // hz
18 disp ( ’ Ghz ’ , fr1 /10^9 , ’ R e s o n a n t f r e q u e n c y f o r mode
TE111 : ’ ) ;
19
20 h21 =3.054;
21 fr2 =( c /(2* %pi ) ) * sqrt (( h21 / a1 ) ^2+( %pi / b1 ) ^2) ; // hz
TE211 : ’ ) ;
23
24 disp ( ’TM modes : ’ ) ;
25 l1 =0;
26 h011 =2.405;
27 fr3 =( c /(2* %pi ) ) * sqrt (( h011 / a1 ) ^2+( %pi * l1 / b1 ) ^2) ; // hz
TM010 ’ ) ;
29
30 l2 =1;
32 disp ( ’ Ghz ’ , fr4 /10^9 , ’ r e s o n a n t f r e q u e n c y f o r mode
TM011 : ’ ) ;
33
34 l3 =1;
35 h111 =3.832;
109
Figure 4.6: Resonator comparison

TM111 : ’ ) ;
Scilab code Exa 4.6 Resonator comparison

2 // Example 4 . 6
3 clc ;
4 // Given
5 QTM010 =1.202;
6 QTE101 =1.11;
7
8 r = QTM010 / QTE101 ;
9 disp (r , ’ R a t i o o f Qs o f c y l i n d r i c a l and r e c t a n g u l a r
r e s o n a t o r s : ’ );
110
Figure 4.7: Cubical Resonator
Scilab code Exa 4.7 Cubical Resonator

2 // Example 4 . 7
3 clc ;
4 // Given
5 f =7.07 D +9; // hz
6 a =3; //cm
111
7 a1 = a /100; //m
8 sig =5.8 D +7;
9 er =2.25;
10 tandel =4 D -4;
11 ur =1;
12 n =377;
13 w =2* %pi * f ;
14 u =4 D -7* %pi ;
15
16 //Q o f r e s o n a n t o r
17 rs = sqrt ( w * u /(2* sig ) ) ; //ohm
18 qw =(0.7419* n ) /( rs * sqrt (2.25) ) ;
19 qd =1/ tandel ;
20 q =( qw * qd ) /( qw + qd ) ;
21 disp (q , ’Q o f r e s o n a t o r : ’ ) ;
Scilab code Exa 4.8 Rectangular Resonant Cavity

2 // Example 4 . 8
3 clc ;
4 // Given
5 a =5; //cm
6 a1 = a /100; //m
7 b =4; //cm
8 b1 = b /100; //m
9 c =10; //cm
10 c1 = c /100; //m
11 sig =5.8 D +7;
12 u0 =4 D -7* %pi ;
13 er =3;
14 eet =377;
15
112
Figure 4.8: Rectangular Resonant Cavity
113
16 ur =1;
17 spl =3 D +8;
18 tandel =2.5 D -4;
19
20 // TE101 mode
21 m =1;
22 n =0;
23 p =1;
24 fr =( spl /(2* sqrt ( er * ur ) ) ) * sqrt (( m / a1 ) ^2+( n / b1 ) ^2+( p /
c1 ) ^2) ; // hz
25 disp ( ’ Ghz ’ , fr /10^9 , ’ R e s o n a n t f r e q u e n c y : ’ ) ;
26
27 w =2* %pi * fr ;
28 rs = sqrt (( w * u0 ) /(2* sig ) ) ; //ohm
29 lamr = spl /( fr * sqrt ( er ) ) ;
30 x =((( a1 * b1 ) /( c1 ^2) ) +(( c1 ^2+ a1 ^2) /(2* c1 * a1 ) ) +(( b1 * c1 )
/ a1 ^2) ) ;
31 qw =(2* %pi *( eet / sqrt ( er ) ) * a1 * b1 * c1 ) /( rs *( lamr ^3) * x ) ;
32 disp ( qw , ’Q f o r TE101 mode : ’ ) ;
33
34 qd =1/ tandel ;
35 q =( qw * qd ) /( qw + qd ) ;
36 disp (q , ’Q f o r l o s s y d i e l e c t r i c : ’ ) ;
37
38 // V a l u e o f qw i s c a l c u l a t e d wrong i n book a s l a m r
comes t o be 0 . 0 8 n o t 0 . 8 9 m
Scilab code Exa 4.9 Rectangular resonator

2 // Example 4 . 9
3 clc ;
4 // Given
114
Figure 4.9: Rectangular resonator
5 c =3 D +8; //m/ s
6 a =2.286; //cm
7 a1 = a /100; //m
8 b =1.106; //cm
9 b1 = b /100; //m
10
11 // For f r 1 =9.3D+9;
12 fr1 =9.3 D +9; // hz
13 lamr1 = c / fr1 ; //m
14 c1 =(2* a1 ) / sqrt ((((2* a1 ) / lamr1 ) ^2) -1) ;
15
16 // For f r 2 =10.2D+9;
17 fr2 =10.2 D +9; // hz
18 lamr2 = c / fr2 ; //m
19 c2 =(2* a1 ) / sqrt ((((2* a1 ) / lamr2 ) ^2) -1) ;
20
21 r = c1 - c2 ;
22 disp ( ’ cm ’ ,r *100 , ’ Range o f p i s t o n movement : ’ ) ;
115
Scilab code Exa 4.10 Cylindrical resonator

2 // Example 4 . 1 0
3 clc ;
4 // Given
5 a =3; //cm
6 a1 = a /100; //m
7 d =10; //cm
8 d1 = d /100; //m
9 df =2.5 D +6;
10 er =2.25;
11 p11 =1.841;
12 c =3 D +8; //m/ s
13
14 // R e s o n a n t f r e q u e n c y
15 fr =( c /2) *( sqrt (( p11 / a1 ) ^2+( %pi / d1 ) ^2) ) ; // hz
16 disp ( ’ Ghz ’ , fr /10^9 , ’ R e s o n a n t f r e q u e n c y : ’ ) ;
17
18 //Q w i t h o u t d i e l e c t r i c
19 q0 = fr / df ;
20 disp ( q0 , ’Q w i r h o u t d i e l e c t r i c c o n s t a n t : ’ ) ;
21
22 // Q w i t h d i e l e c t r i c
23 fr1 = fr / sqrt ( er ) ;
24 qd =1 D +3;
25 q =( q0 * qd ) /( q0 + qd ) ;
26 disp (q , ’Q w i t h d i e l e c t r i c c o n s t a n t : ’ ) ;
116
Figure 4.10: Cylindrical resonator
117
Figure 4.11: Cylindrical resonantor
118
Scilab code Exa 4.11 Cylindrical resonantor

2 // Example 4 . 1 1
3 clc ;
4 // Given
5 f =9.375 D +9; // hz
6 sig =5.8 D +7;
7 eet =377;
8 c =3 D +8; //m/ s
9 w =2* %pi * f ;
10 r =1.5;
11 u =4 D -7* %pi ;
12
13 // R a d i u s
14 a = c /( f *2.62) ; //m
15 disp ( ’ cm ’ ,a *100 , ’ R a d i u s o f r e s o n a n t o r ’ ) ;
16
17 //O
18 rs = sqrt (( w * u ) /(2* sig ) ) ; //ohm
19 x =1.202* eet ;
20 y = rs *(1+(1/ r ) ) ;
21 q=x/y;
22 disp (q , ’Q o f t h e r e s o n a t o r : ’ ) ;
23
24 // Answer f o r Q i s c a l c u l a t e d a s 1 0 8 7 5 i n book but i t
i s 10763.303
Scilab code Exa 4.12 Cylindrical Resonator
119
Figure 4.12: Cylindrical Resonator
120
2 // Example 4 . 1 2
3 clc ;
4 // Given
5 f =5 D +9; // hz
6 sig =5.813 D +7;
7 er =2.25;
8 tandel =4 D -4;
9 c =3 D +8; //m/ s
10 h01 =3.832;
11 u =4 D -7* %pi ;
12
13 // Length o f r e s o n a t o r
14 lamr = c /( f * sqrt ( er ) ) ;
15 d = sqrt ([{(((2*3.832) ^2) +( %pi * %pi ) ) *( lamr * lamr )
}/(2*2* %pi * %pi ) ]) ;
16 disp ( ’ cm ’ ,d *100 , ’ Length o f r e s o n a t o r : ’ ) ;
17
18 //Q o f r e s o n a t o r
19 n =(120* %pi ) / sqrt ( er ) ;
20 Rs = sqrt (( f * u ) / sig ) ;
21 a = d /2;
22 Qw1 = n *[[( h01 / a ) ^2+( %pi / d ) ^2]^(3/2) ];
23 Qw2 =2* Rs *[(( h01 * h01 ) /( a * a * a ) ) +((2* %pi * %pi ) /( d * d * d ) )
];
24 Qw = Qw1 / Qw2 ;
25 Qd =1/ tandel ;
26 Q =( Qw * Qd ) /( Qw + Qd ) ;
27 disp (Q , ’Q o f r e s o n a t o r : ’ ) ;
28
29 // V a l u e o f Qw i s c a l c u l a t e d wrong i n t h e book , i t
s h o u l d be 5 0 0 5 7 . 9 1 i n s t e a d o f 5 3 4 7 3 . 8
30 // Hence t h e v a l u e o f Q a l s o d i f f e r s
121
Figure 4.13: Lossless plane H tee
122
Scilab code Exa 4.13 Lossless plane H tee

2 // Example 4 . 1 3
3 clc ;
4 // Given
5 p =100; //mW
6 // As 2 and 3 a r e matched t e r m i n a l s
7 x =1/2;
8 y =1/ sqrt (2) ;
9 s =[ x -x y ; - x 0 y ; y y 0];
10
11 // Power d e l i v e r e d
12 // P o r t 1
13 p1 = p *(1 - s (1 ,1) ^2) ;
14 disp ( ’mW’ ,p1 , ’ Power a t p o r t 1 : ’ ) ;
15
16 // P o r t 2
17 p2 = p * s (2 ,1) ^2;
19
20 // P o r t 3
21 p3 = p * s (3 ,1) ^2;
Scilab code Exa 4.14 E plane tee

2 // Example 4 . 1 4
3 clc ;
4 // Given
5 p =40; //mW
6 // S i n c e p o r t 3 i s matched
123
Figure 4.14: E plane tee
124
7 x = sqrt (2) ;
8 s =[1 1 x ;1 1 -x ; x -x 0];
9 r1 =40; //ohm
10 r2 =60; //ohm
11 w =50; //ohm
12
13 // R e f l e c t i o n c o e f f i c i e n t s
14 T1 =( w - r1 ) /( w + r1 ) ;
15 T2 =( r2 - w ) /( r2 + w ) ;
16
17 // As power i s f e d i n t o 1 and 2 e q u a l l y
18 pd = p /2;
19
20 // Power d e l i v e r e d
21 // P o r t 1
22 p1 = pd *(1 - T1 ^2) ;
24
25 // P o r t 2
26 p2 = pd *(1 - T2 ^2) ;
Scilab code Exa 4.15 Magic Tee

2 // Example 4 . 1 5
3 clc ;
4 // Given
5 T1 =1/2;
6 T2 =3/5;
7 T3 =0;
8 T4 =4/5;
9 p =500 D -3; //W
125
126
Figure 4.15: Magic Tee

10 // S m a t r i x f o r magic Tee
11 x =1/ sqrt (2) ;
12 s =[0 0 x x ;0 0 x -x ; x x 0 0; x -x 0 0];
13 // U s i n g t h e i n p u t o u t p u t r e l a t i o n
14 // [ b ] = [ s ] ∗ [ a ]
15 b = [0.6565;0.7576;0.5536;0.0892];
16
17 // ( i ) Power t r a n s m i t t e d t h r o u g h p o r t s
18 // P o r t 1
19 p1 =(1/2) * b (1 ,1) ^2*(1 - T1 ^2) ;
20 disp ( ’W’ ,p1 , ’ Power a t p o r t 1 : ’ ) ;
21
22 // P o r t 2
23 p2 =(1/2) *( b (2 ,1) ^2) *(1 -( T2 ^2) ) ;
25
26 // P o r t 4
27 p4 =(1/2) * b (4 ,1) ^2*(1 - T4 ^2) ;
29
30 // ( i i ) Power r e f l e c t e d a t p o r t 3
31 // P o r t 3
32 p3 = p * b (3 ,1) ^2;
34
35 // ( i i i ) Power a b s o r b e d
36 pabs =p -( p1 + p2 + p3 + p4 ) ;
37 disp ( ’W’ , pabs , ’ Power a b s o r b e d : ’ ) ;
38
39 // Answer f o r power a b s o r b e d i s c a l c u l a t e d wrong i n
book
Scilab code Exa 4.18 Directional Coupler
127
Figure 4.16: Directional Coupler

2 // Example 4 . 1 8
3 clc ;
4 // Given
5 C =10; //dB
6 D =30; //dB
7
8 // P a r a m e t e r s
9 bet =10^( - C /20) ;
10 x = bet /(10^( D /20) ) ;
11 a = sqrt (1 -( bet * bet ) ) ;
12 // S c a t t e r i n g m a t r i x
13 // Assuming symmetery
14 s =[0 a x ( bet * %i ) ; a 0 ( bet * %i ) x ; x ( bet * %i ) 0 a ;( bet
* %i ) x a 0];
15 disp (s , ’ S c a t t e r i n g m a t r i x : ’ ) ;
128
Figure 4.17: Directional coupler
129
Scilab code Exa 4.20 Directional coupler

2 // Example 4 . 2 0
3 clc ;
4 // Given
5 vswr =2;
6 D1 =8; //mW
7 D2 =2; //mW
8
9 // R e f l e c t i o n c o e f f i c i e n t a t arm 4
10 T =( vswr -1) /( vswr +1) ;
11 //Powwe d e l i v e r e d t o D1
12 P =( D1 *100) /(1 - T ^2) ;
13 P1 =0.99* P ;
14 // Power r e f l e c t e d a t D1
15 W1 =( P /100) * T * T ;
16 // Power r e f l e c t e d a t l o a d
17 W2 = D2 - W1 ;
18 Tt = sqrt (( W2 *100) /( P1 ) ) ;
19 pt =(1+ Tt ) /(1 - Tt ) ;
20 disp ( pt , ’VSWR: ’ ) ;
21 Pl = P1 *(1 -( Tt * Tt ) ) ;
22 disp ( ’mW’ ,Pl , ’ Power d e l i v e r e d : ’ ) ;
23
24 // Answer f o r P1 s h o u l d be 792 but i t i s g i v e n a s 800
Scilab code Exa 4.21 Isolator Matrix
130
Figure 4.18: Isolator Matrix
131
Figure 4.19: Circulator Matrix
2 // Example 4 . 2 1
3 clc ;
4 // Given
5 I =30; //dB
6 Il =0.4; //dB
7
8 S12 =10^( I / -20) ;
9 S21 =10^( Il / -20) ;
10 s =[0 S12 ; S21 0];
132
Scilab code Exa 4.22 Circulator Matrix

2 // Example 4 . 2 2
3 clc ;
4 // Given
5 I =30; //dB
6 Il =2; //dB
7 p =1.3;
8
9 // E l e l m e n t s
10 T =( p -1) /( p +1) ;
11 S11 = T ;
12 S22 = T ;
13 S33 = T ;
14 S12 =10^( - Il /20) ;
15 S13 =10^( - I /20) ;
16 S21 = S13 ;
17 S32 = S13 ;
18 S23 = S12 ;
19 S31 = S23 ;
20 s =[ S11 S21 S31 ; S12 S22 S32 ; S13 S23 S33 ];

2 // Example 4 . 2 3
3 clc ;
4 // Given
5 f =10 D +9; // Hz
6 u =4 D -7* %pi ;
7 c =3 D +8; //m/ s
133
134
8 a =2.29; //cm
9 a1 = a /100;
10 b =1.02; //cm
11 b1 = b /100;
12
13 //E/H
14 w =2* %pi * f ;
15 EbyH =( w * u ) / sqrt ((( w / c ) ^2) +(( %pi / a1 ) ^2) ) ;
16 lam = c / f ;
17 lamc =2* a1 ;
18 d =(1/4) *( lam / sqrt (1 -(( lam / lamc ) ^2) ) ) ;
19 disp ( ’ cm ’ ,d *100 , ’ P o s i t i o n : ’ ) ;
20
21 // Answer f o r p o s i t o n i s c a l c u l a t e d wrong i n book
Scilab code Exa 4.24 Attenuator matrix

2 // Example 4 . 2 4
3 clc ;
4 // Given
5 // As i t i s p e r f e c t l y matched
6 S12 =1/ sqrt (2) ;
7 S21 = S12 ;
8 s =[0 S12 ; S21 0];
135
Figure 4.21: Attenuator matrix
136
Chapter 5
Microwave Tubes Klystrons
Scilab code Exa 5.1 Two Cavity Klystron

2 // Example 5 . 1
3 clc ;
4 // Given
5 f =10 D +9; // Hz
6 v =9 D +3; //V
7 i =40 D -3; //A
8 l =3; //cm
9 l1 = l /100; //m
10 G =2 D -6; //mho
11 bet =0.92;
12 j1x =0.582;
13 x =1.841;
14 ebym =1.7 D +11; // J
15
16 //Maximum v o l t a g e
17 w =2* %pi * f ;
18 v0x = sqrt (2* ebym ) ;
19 thet =( w * l1 ) /( v0x * sqrt ( v ) ) ;
137
Figure 5.1: Two Cavity Klystron
138
20
21 av =( bet ^2* thet * i * j1x ) /( x * v * G ) ;
22 disp ( ’V ’ ,av , ’ Maximum v o l t a g e : ’ ) ;
23
24 // Power Gain
25 ic =2* i * j1x ;
26 v2 =( bet * ic ) / G ;
27 pout = bet * ic * v2 ;
28 pin =2* i * v ;
29
30 // E f f i c i e n c y
31 eet = pout / pin ;
32 disp ( ’% ’ , eet *100 , ’ Power g a i n : ’ ) ;
33
34 // Answer f o r e f f c i e n c y comes o u t t o be wrong , i t is
c a l c u l t e d w r o n g l y i n book
Scilab code Exa 5.2 Two cavity Klystron

2 // Example 5 . 2
3 clc ;
4 // Given
5 l =2; //cm
6 l1 = l /100; //m
7 f =5 D +9; // Hz
8 i =25 D -3; //A
9 n =21/4;
10 e =1.6 D -19;
11 m =9.1 D -31;
12 thetag =0;
13 bet =1;
14 j1x =0.582;
139
Figure 5.2: Two cavity Klystron
140
15 x =1.841;
16
17 // ( i ) Beam V o l t a g e
18 v0 =( m * l1 * l1 * f * f ) /(2* e * n * n ) ;
19 disp ( ’V ’ ,v0 , ’ Beam v o l t a g e : ’ ) ;
20
21 // ( i i ) I n p u t v o l t a g e
22 v1 = x * v0 /( %pi * bet * n ) ;
23 disp ( ’V ’ ,v1 , ’ I n p u t v o l t a g e : ’ ) ;
24
25 // ( i i i ) Output v o l t a g e
26 v2 =0.25* v0 ;
27 disp ( ’V ’ ,v2 , ’ Output v o l t a g e ’ ) ;
28
29 // ( i v ) Power o u t p u t
30 pmax = i * v0 * j1x ;
31 disp ( ’W’ , pmax , ’ Maximum power o u t p u t : ’ ) ;
32
33 // ( v ) E f f i c i e n c y
34 eet = j1x * bet * v2 / v0 ;
35 disp ( ’% ’ , eet *100 , ’ E f f i c i e n c y : ’ ) ;

2 // Example 5 . 3
3 clc ;
4 // Given
5 r0 =45 D +3; //W
6 j0 =25 D -3; //A
7 V =1500; //V
8 f =5 D +9; // hz
9 d =1; //mm
141
142
10 d1 = d /1000; //m
11 l =3.5; //cm
12 l1 = l /100; //m
13 rsh =32 D +3; // ohms
14 j1x =0.582;
15 x =1.841;
16
17 // ( i ) I n p u t gap v o l t a g e
18 w =2* %pi * f ;
19 v0 =(5.93 D +5* sqrt ( V ) ) ;
20 thetag =( w * d1 ) / v0 ;
21 bet = sin ( thetag /2) /( thetag /2) ;
22 theta0 =( w * l1 ) / v0 ;
23 v1 =(2* V * x ) /( bet * theta0 ) ;
24 disp ( v1 , ’ I n p u t gap v o l t a g e : ’ ) ;
25
26 // ( i i ) V o l t a g e g a i n
27 av =( bet ^2* theta0 * j1x * rsh ) /( r0 * x ) ;
28 disp ( av , ’ V o l t a g e g a i n : ’ ) ;

2 // Example 5 . 4
3 clc ;
4 // Given
5 V =1000; //V
6 r0 =40 D +3; //ohm
7 i0 =25 D -3; //A
8 f =3 D +9; // Hz
9 d =1; //mm
10 d1 = d /1000; //m
11 l =4; //cm
143
144
12 l1 =4/100; //m
13 j1x =0.582;
14 x =1.841;
15 rsh =30 D +3; //ohm
16
17 // ( i ) I n p u t gap v o l t a g e
18 w =2* %pi * f ;
19 v0 =(5.93 D +5* sqrt ( V ) ) ;
20 thetag =( w * d1 ) / v0 ;
22 theta0 =( w * l1 ) / v0 ;
23 vmax =(2* V * x ) /( bet * theta0 ) ;
24 disp ( ’V ’ , vmax , ’ I n p u t gap v o l t a g e : ’ ) ;
25
26 // ( i i ) V o l t a g e g a i n
27 av =( bet * bet * theta0 * j1x * rsh ) /( r0 * x ) ;
29
30 // ( i i i ) E f f i c i e n c y
31 v2 = bet *2* i0 * j1x * rsh ;
32 eet =( bet *2* i0 * j1x * v2 ) /(2* i0 * V ) ;
34
35 // ( i v ) Beam l o a d i n g c o n d u c t a n c e
36 gbl =( i0 /(2* V ) ) *(( bet * bet ) -( bet * cos ( thetag /2) ) ) ;
37 disp ( gbl , ’ Beam l o a d i n g c o n d u c t a n c e : ’ ) ;
38
39 // Ansewr f o r beam l o a d i n g c o n d u c t a n c e i s c a l c u l a t e d
wrong i n book
145
146
2 // Example 5 . 5
3 clc ;
4 // Given
5 f =3 D +9; // hz
6 v =900; //V
7 i =30 D -3; //A
8 d =4; //cm
9 d1 = d /100; //m
10 gap =1; //mm
11 gap1 =1/1000; //m
12 rsh =40 D +3; //ohm
13 x =1.841;
14 j1x =0.582;
15 r =40 D +3; //ohm
16 ebym =1.758 D +11; // J
17
18 // ( i ) E l e c t r o n v e l o c i t y
19 v0 = sqrt (2* ebym * v ) ;
20 disp ( ’m/ s ’ ,v0 , ’ E l e c t r o n v e l o c i t y : ’ ) ;
21
22 // ( i i ) E l e c t r o n t r a n s i t t i m e
23 t = d1 / v0 ;
24 disp ( ’ s ’ ,t , ’ E l e c t r o n t r a n s i t t i m e : ’ ) ;
25
26 // ( i i i ) I n p u t v o l t a g e gap
27 w =2* %pi * f ;
28 theta0 =( w * d1 ) / v0 ;
29 thetag =( w * gap1 ) / v0 ;
31 v2 =(2* v * x ) /( bet * theta0 ) ;
32 disp ( ’V ’ ,v2 , ’ I n p u t v o l t a g e gap : ’ ) ;
33
34 // ( i v ) V o l t a g e g a i n
35 av =( bet ^2* theta0 * j1x * rsh ) /( x * r ) ;
37
38 // V a l u e s o f v and f a r e c h a n g e d i n q u e s t i o n and
answer , h e n c e v a u l e s u s e d i n a n s w e r a r e t a k e n .
147
39 // A l s o s e c o n d p a r t h a s n o t b e e n done i n book

2 // Example 5 . 6
3 clc ;
4 // Given
5 f =8 D +9; // hz
6 i =2.5; //A
7 v =20 D +3; //V
8 bet =1;
9 amp =10* sqrt (2) ; //V
10 rsh =10 D +3; //ohm
11 rsho =30 D +3; //ohm
12 dc =1 D -6; // c /mˆ3
13 rf =0.5;
14 e =1.6 D -19;
15 ee =8.854 D -12;
16 m =9.1 D -31; // kg
17
18 // ( i ) I n d u c e d c u r r e n t
19 w =2* %pi * f ;
20 wq = rf * sqrt (( e * dc ) /( m * ee ) ) ;
21
22 // A m p l i t u d e o f i n d u c e d c u r r e n t
23 ic =( i * w *( bet ^2) * amp ) /(2* v * wq ) ;
24 disp ( ’A ’ ,ic , ’ I n d u c e d c u r r e n t : ’ ) ;
25
26 // I n d u c e d v o l t a g e
27 icrms = ic / sqrt (2) ;
28 v2rms = icrms * rsho ;
29 disp ( ’V ’ , v2rms , ’ I n d u c e d v o l t a g e : ’ ) ;
148
149
30
31 // ( i i ) Power g a i n
32 pg =((( i * w ) ^2) *( bet ^4) * rsh * rsho ) /(4*(( v * wq ) ^2) ) ;
33 pgdb =10* log10 ( pg ) ;
34 disp ( ’ dB ’ , pgdb , ’ Power g a i n : ’ ) ;
35
36 // ( i i i ) E l e c t r o n i c e f f i c i e n c y
37 eeta =(( icrms ^2) * rsho ) /( i * v ) ;
38 disp ( ’% ’ , eeta *100 , ’ E l e c t r o n i c e f f i c i e n c y : ’ ) ;

2 // Example 5 . 7
3 clc ;
4 // Given
5 f =3 D +9; // hz
6 l =4; //cm
7 l1 =4/100; //m
8 d =0.1; //cm
9 d1 = d /100; //m
10 V =900; //V
11 i0 =30 D -3; //A
12 rsh =25 D +3; //ohm
13 x =1.841;
14 j1x =0.582;
15
16 // ( i ) I n p u t v o l t a g e f o r maximum o u t p u t
17 v0 =0.593 D +6* sqrt ( V ) ;
18 w =2* %pi * f ;
19 theta0 = w * l1 / v0 ; // r a d
20 thetag = w * d1 / v0 ; // r a d
150
151
22 v1max =2* V * x /( bet * theta0 ) ; // v
23 disp ( ’V ’ , v1max , ’ I n p u t v o l t a g e f o r maximum o u t p u t : ’ ) ;
24
25 // ( i i ) V o l t a g e g a i n
26 r0 = V / i0 ; //ohm
27 av =(( bet ^2) * theta0 * j1x * rsh ) /( x * r0 ) ; //V
28 disp ( ’V ’ ,av , ’ V o l t a g e g a i n : ’ ) ;
29
30 // ( i i i ) E f f i c i e n c y
31 ic =2* i0 * j1x ; //A
32 v2 = bet * ic * rsh ; //V
33 eet = bet * ic * v2 /(2* i0 * V ) ;
35
36 // ( i v ) Beam l o a d i n g c o n d u c t a n c e
37 gb =( i0 /( V *2) ) *( bet ^2 -( bet * cos ( thetag /2) ) ) ; //ohm
38 disp ( ’ ohm ’ ,gb , ’ Beam l o a d i n g c o n d u c t a n c e : ’ ) ;

2 // Example 5 . 8
3 clc ;
4 // Given
5 f =5 D +9; // hz
6 v0 =10 D +3; //V
7 d =1; //mm
8 d1 = d /1000; //m
9 v1 =100; //V
10
11 // ( i ) Gap t r a n s i t t i m e
12 vv0 =0.593 D +6* sqrt ( v0 ) ; //m/ s e c
13 tau = d1 / vv0 ; // s e c
152
153
14 disp ( ’ s e c ’ ,tau , ’ Gap t r a n s i t t i m e : ’ ) ;
15
16 //Gap t r a n s i t a n g l e
17 w =2* %pi * f ;
18 thetag = w * tau ; // r a d
19 disp ( ’ r a d ’ , thetag , ’ Gap t r a n s i t a n g l e : ’ ) ;
20
21 // ( i i ) Beam c o u p l i n g c o e f f i c i e n t
22 betin = sin ( thetag /2) /( thetag /2) ;
23 disp ( betin , ’ Beam c o u p l i n g c o e f f i c i e n t : ’ ) ;
24
25 // ( i i i ) V e l o c i t y o f e l e c t r o n l e a v i n g b u n c h e r gap
26 vig = vv0 *(1+(( betin * v1 ) /(2* v0 ) ) ) ; //m/ s e c
27 disp ( ’m/ s e c ’ ,vig , ’ V e l o c i t y o f e l e c t r o n l e a v i n g
b u n c h e r gap : ’ ) ;
28
29 // ( i v ) Depth o f m o d u l a t i o n
30 m = betin * v1 / v0 ;
31 disp (m , ’ Depth o f m o d u l a t i o n : ’ ) ;
Scilab code Exa 5.9 Four cavity Klystron

2 // Example 5 . 9
3 clc ;
4 // Given
5 f =10 D +9; // hz
6 v0 =15 D +3; //V
7 i0 =2.5 D -3; //A
8 d =1; //cm
9 d1 = d /100; //m
10 vrms =10; //V
11 bet =1;
154
Figure 5.9: Four cavity Klystron
155
12 p =1 D -8; //C/mˆ3
13 rf =0.6;
14 e =1.6 D -19;
15 m =9.1 D -31;
16 ee =8.854 D -12;
17
18 // ( i ) DC e l e c t r o n beam p h a s e c o b s t a n t
19 vv0 =(0.593 D +6* sqrt ( v0 ) ) ;
20 w =2* %pi * f ;
21 bete = w / vv0 ; // r a d /m
22 disp ( ’ r a d /m ’ , bete , ’DC e l e c t r o n beam p h a s e c o n s t a n t : ’
);
23
24 // ( i i ) Reduced p l a s m a f r e q u e n c y and r e d u c e d p l a s m a
phase constant
25 wq = rf * sqrt ( e * p /( m * ee ) ) ; // r a d /m
26 disp ( ’ r a d /m ’ ,wq , ’ Reduced p l a s m a f r e q u e n c y : ’ ) ;
27 betq = wq / vv0 ; // r a d / s e c
28 disp ( ’ r a d / s e c ’ , betq , ’ Reduced p l a s m a p h a s e c o n s t a n t : ’
);
29
30 // ( i i i ) Gap t r a n s i t t i m e
31 tau = d1 / vv0 ; // s e c
32 vtg = vv0 *(1+( bet * vrms * sin ( w * tau ) /(2* v0 ) ) ) ; //m/ s e c
33 disp ( ’m/ s e c ’ ,vtg , ’ Gap t r a n s i t t i m e : ’ ) ;
Scilab code Exa 5.10 Four cavity Klystron

2 // Example 5 . 1 0
3 clc ;
4 // Given
5 f =4 D +9; // hz
156
Figure 5.10: Four cavity Klystron
157
6 v0 =10 D +3; //V
7 i0 =0.75; //A
8 v1 =2; //V
9 bet =1;
10 rsh =10 D +3; //ohm
11 p =5 D -5; //C/mˆ3
12 r =0.6;
13 rsht =4 D +3; //ohm
14 e =1.6 D -19;
15 m =9.1 D -31;
16 ee =8.854 D -12;
17
18 // ( i ) I n d u c e d c u r r e n t and v o l t a g e i n o u t p u t c a v i t y
19 w1 = sqrt ( e * p /( m * ee ) ) ; // r a d / s e c
20 w =2* %pi * f ;
21 wq =0.5* w1 ; // r a d / s e c
22 rr = w / wq ;
23
24 i4 =[( i0 ^3) *( rr ^3) *( bet ^6) * v1 *( rsh ^2) ]/(8*( v0 ^3) ) ; //
A
25 disp ( ’A ’ ,i4 , ’ I n d u c e d c u r r e n t : ’ ) ;
26 v4 = i4 * rsht ; //V
27 disp ( ’ kV ’ , v4 /1000 , ’ I n d u c e d v o l t a g e : ’ ) ;
28
29 // ( i i ) Power o u t p u t
30 pout =( i4 ^4) * rsht ; //W
31 disp ( ’W’ , pout , ’ Power o u t p u t : ’ ) ;
32
33 // Answer f o r Pout s h o u l d be 1 3 . 4 3 kW but i t i s g i v e n
a s 1 0 . 8 9kW a s v a l u e o f I 4 i s c a l c u l a t e d a s 1 . 2 8 9
but i t comes o u t t o be 1 . 3 5
Scilab code Exa 5.11 Reflex Klystron
158
Figure 5.11: Reflex Klystron
159
2 // Example 5 . 9
3 clc ;
4 // Given
5 f =8 D +9; // hz
6 v0 =500; //V
7 l =1.2; //mm
8 l1 = l /1000; //m
9 rsh =18 D +3; //ohm
10 ebym =1.759 D +11;
11 ee =8.854 D -12;
12
13 // ( i ) R e p e l l e r v o l t a g e
14 n =1+(3/4) ;
15 v11 =( ebym * n * n ) /(8*( l1 ^2) *( f ^2) ) ;
16 vr = sqrt ( v0 / v11 ) - v0 ;
17 disp ( ’V ’ ,vr , ’ R e p e l l e r v o l t a g e : ’ ) ;
18
19 // ( i i ) R e q u i r e d dc c u r r e n t
20 v2 =200; //V
21 j1x =0.582;
22 i = v2 /(2* rsh * j1x ) ; //A
23 disp ( ’mA ’ ,i *1000 , ’ R e q u i r e d dc c u r r e n t : ’ ) ;
24
25 // Answer f o r r e p e l l e r v o l t a g e i s c a l c u l a t e d wrong i n
book

2 // Example 5 . 1 2
3 clc ;
4 // Given
160
161
5 f =9 D +9; // hz
6 v0 =361; //V
7 i0 =30 D -3; //A
8 l =0.1; //cm
9 l1 = l /100; //m
10 x =2.408;
11 j1x =0.582;
12 ebym =1.759 D +11;
13
14 //Maximum power o u t p u t
15 n =1;
16 pout =2* i0 * v0 * x * j1x /(2* %pi *( n +(3/4) ) ) ; //W
17 disp ( ’W’ , pout , ’ Maximum power o u t p u t : ’ ) ;
18
19 // O p e r a t i n g r e p e l l e r v o l t a g e
20 vr =((6.744 D -6* sqrt ( v0 ) * l1 * f ) /( n +(3/4) ) ) - v0 ; // v
21 disp ( ’V ’ ,vr , ’ O p e r a t i n g r e p e l l e r v o l t a g e : ’ ) ;

2 // Example 5 . 1 3
3 clc ;
4 // Given
5 f =9 D +9; // hz
6 v0 =250; //V
7 l =0.5; //cm
8 l1 = l /100; //m
9
10 // Bandwidth
11 n =3;
12 df =( n +(3/4) ) /(6.774 D -6* l1 * sqrt ( v0 ) ) ; // hz
13 disp ( ’ Mhz ’ , df /10^6 , ’ Bandwidth : ’ ) ;
162
163

2 // Example 5 . 1 4
3 clc ;
4 // Given
5 f =10 D +9; // hz
6 v0 =600; //V
7 vr =250; //V
8 ebym =1.759 D +11;
9
10 // R e p e l l e r s p a c e
11 n =1;
12 l = sqrt (( ebym *( n +(3/4) ) ^2*( vr + v0 ) ^2) /(8* f ^2* v0 ) ) ; //m
13 disp ( ’mm’ ,l *1000 , ’ R e p e l l e r s p a c e : ’ ) ;

2 // Example 5 . 1 5
3 clc ;
4 // Given
5 v0 =300; //V
6 i0 =20 D -3; //A
7 v1 =40; //V
8 n =2;
9 x =2.408;
164
165
166
10 j1x =0.52;
11
12 // ( i ) I n p u t power
13 pin = i0 * v0 ; //W
14 disp ( ’W’ ,pin , ’ I n p u t power : ’ ) ;
15
16 // ( i i ) Output power
17 pout =(2* v0 * i0 * x * j1x ) /((2* %pi * n ) -( %pi /2) ) ; //W
18 disp ( ’W’ , pout , ’ Output power : ’ ) ;
19
21 eet = pout / pin ;
23
24 // Answer f o r o u t p u t power i n book i s 0 . 7 which i s
wrong , i t s h o u l d be 1 . 3W
25 // Hence a n s w e r o f e f f i c i e n c y a l s o c h a n g e s

2 // Example 5 . 1 6
3 clc ;
4 // Given
5 f =10 D +9; // hz
6 v0 =600; //V
7 l =0.1; //cm
8 l1 = l /100; //m
9 bet =0.9;
10 ebym =1.759 D +11;
11 n =2;
12 j1x =0.575; // from s t a n d a r d t a b l e
13
167
168
14
15 // ( i ) R e p e l l e r v o l t a g e
16 vr =((6.744 D -6* sqrt ( v0 ) * l1 * f ) /( n -(1/4) ) ) - v0 ; //V
17 disp ( ’V ’ , round ( vr ) , ’ R e p e l l e r v o l t a g e : ’ ) ;
18
19 // ( i i ) Bunching p a r a m e t e r
20 v1 =200; //V
21 x = bet * v1 *2* %pi *( n -(1/4) ) /(2* v0 ) ;
22 disp (x , ’ Bunching p a r a m e t e r : ’ ) ;
23
24 // ( i i i ) R e q u i r e d DC c u r r e n t
25 rsh =20 D +3; //ohm
26 i = v1 /(2* rsh * j1x ) ; //A
27 disp ( ’mA ’ ,i *1000 , ’ R e q u i r e d DC c u r r e n t : ’ ) ;
28
29 // ( i v ) E l e c t r o n i c e f f i c i e n c y
30 eet =2* x * j1x /(2* %pi *( n -(1/4) ) ) ;
31 disp ( ’% ’ , eet *100 , ’ E l e c t r o n i c e f f i c i e n c y : ’ ) ;
Scilab code Exa 5.17 Electron Gun

2 // Example 5 . 1 7
3 clc ;
4 // Given
5 f =10 D +9; // hz
6 v0 =300; //V
7 j0 =0.3; //A/cm
8 i0 =45 D -3; //A
9
10 rb = sqrt ( i0 /( %pi * j0 ) ) ; //mm
11 disp ( ’mm’ , rb *10 , ’ E l e c t r o n beam r a d i u s : ’ ) ;
12 r = rb *(120/100) ; //mm
169
Figure 5.17: Electron Gun
170
13 disp ( ’mm’ ,r *10 , ’ R a d i u s o f c a t h o d e d i s c : ’ ) ;
14 d = sqrt (2.335 D -6*(300) ^(3/2) / j0 ) ; //mm
15 disp ( ’mm’ ,d *10 , ’ Cathode anode s p a c i n g : ’ ) ;
16 // Anode h o l e h a s t o be 15% l a r g e r t h a n c a t h o d e d i s c
17 ra = r *1.15; //mm
18 disp ( ’mm’ , ra *10 , ’ Anode h o l e : ’ ) ;
Scilab code Exa 5.18 Re entrant Coaxial Cavity
1 // Page Number :
2 // Example 5 . 1 8
3 clc ;
4 // Given
5 f =9 D +9; // hz
6 v0 =300; //V
7 vr =125; //V
8 bet =0.9;
9 c =3 D +8; //m/ s
10 w =2* %pi * f ;
11 br =2.18; //mm
12 e0 =8.854 D -12;
13 ebym =1.7 D +11;
14
15 // From s i n ( t h e t a ) / t h e t a t a b l e , t h e t a g i s f o u n d o u t
t o be
16 thetag =0.25* %pi ;
17 d =(2* thetag *0.593 D +6* sqrt ( v0 ) ) / w ;
18 disp ( ’mm’ ,d *1000 , ’ D i s t a n c e : ’ ) ;
19
20 // A x i a l c a v i t y l e n g t h
21 l = c /(10* f ) ; //m
22 disp ( ’mm’ ,l *1000 , ’ A x i a l c a v i t y l e n g t h : ’ );
23
171
172
Figure 5.18: Re entrant Coaxial Cavity

24 // R a t i o o f o u t e r t o i n n e r c o n d u c t o r
25 a =1.5* br ;
26 a1 = a /1000;
27 x = d /( w * e0 * a1 * a1 *60* tan (( w * l ) / c ) ) ;
28 bbya = exp ( x ) ;
29 disp ( bbya , ’ R a t i o o f o u t e r t o i n n e r c o n d u c t o r : ’ ) ;
30
31 // r a d i i o f o u t e r and i n n e r c o n d u c t o r
32 disp ( ’mm’ ,a , ’ R a d i u s o f o u t e r c o n d u c t o r : ’ ) ;
33
34 b =1.52* a ; //mm
35 disp ( ’mm’ ,b , ’ R a d i u s o f i n n e r c o n d u c t o r : ’ ) ;
36
37 // R e p e l l e r s p a c i n g
38 lopt = sqrt ( ebym *(19/4) ^2*( v0 + vr ) ^2/(8* f ^2* v0 ) ) ; //m
39 disp ( ’mm’ , lopt *1000 , ’ R e p e l l e r s p a c i n g : ’ ) ;
40
41 // Answer f o r r a d i i o f o u t e r and i n n e r c o n d u c t o r have
wrong c a l c u l a t i o n s i n book
42 // A l s o r a t i o o f o u t e r t o i n n e r c o n d u c t o r i s a l s o
c a l c u l a t e d wrong
173
Chapter 6
Microwave Travelling Wave

Tubes O type
Scilab code Exa 6.1 TWT

2 // Example 6 . 1
3 clc ;
4 // Given
5 clc ;
6 // Given
7 I0 =30 D -3; //A
8 V0 =3 D +3; //V
9 Z0 =10; //ohm
10 l =0.1624; //m
11 f =10 D +9; // Hz
12
13 // ( i ) Gain p a r a m e t e r
14 C =(( I0 * Z0 ) /(4* V0 ) ) ^(1/3) ;
15 disp (C , ’ Gain p a r a m e t e r : ’ ) ;
16
17 N =( l * f ) /(0.593 D +6* sqrt ( V0 ) ) ;
174
Figure 6.1: TWT
175
18
19 // ( i i ) Power Gain
20 Ap = -9.54+(47.3* C * N ) ;
21 disp ( ’ dB ’ ,Ap , ’ Power g a i n : ’ ) ;
22
23 ve =0.593 D +6* sqrt ( V0 ) ;
24 be =(2* %pi * f ) / ve ;
25
26 // Four p r o p o g a t i o n c o n s t a n t s
27 gam1 =(( - sqrt (3) * be * C ) /2) +( %i * be *(2+ C ) ) /2;
28 gam2 =(( sqrt (3) * be * C ) /2) +( %i * be *(2+ C ) ) /2;
29 gam3 = %i * be *(1 - C ) ;
30 gam4 = - %i * be *(1 -(( C * C * C ) /4) ) ;
31
32 disp ( gam4 , gam3 , gam2 , gam1 , ’ Four p r o p o g a t i o n c o n s t a n t s
: ’ );
33
34 // C a l c u l a t i o n s f o r p r o p o g a t i o n c o n s t a n t s a r e wrong
i n book f o r gam 3 and 4 , h e n c e a n s w e r s d o n t match
Scilab code Exa 6.2 Helix TWT

2 // Example 6 . 2
3 clc ;
4 // Given
5 I0 =20 D -3; //A
6 V0 =4 D +3; //V
7 Z0 =100; //ohm
8 N =30;
9
10 C =(( I0 * Z0 ) /(4* V0 ) ) ^(1/3) ;
11 // Gain
176
Figure 6.2: Helix TWT
177
12 Ap = -9.54+(47.3* C * N ) ;
13 disp ( ’ dB ’ ,Ap , ’ Gain : ’ ) ;
Scilab code Exa 6.3 Helical TWT

2 // Example 6 . 3
3 clc ;
4 // Given
5 c =3 D +8; //m/ s
6 d =2 D -3; //m
7 p =50 D +2; // t u r n s p e r m
8 e =1.6 D -19; // J
9 m =9.1 D -31;
10
11 // A x i a l p h a s e v e l o c i t y
12 vp = c /( %pi * p * d ) ;
13 disp ( ’m/ s ’ ,vp , ’ A x i a l p h a s e v e l o c i t y : ’ ) ;
14
15 // Anode v o l t a g e
16 V0 =( m * vp * vp ) /(2* e ) ;
17 disp ( ’V ’ ,V0 , ’ Anode v o l t a g e : ’ ) ;
Scilab code Exa 6.4 O type TWT

2 // Example 6 . 4
3 clc ;
4 // Given
178
Figure 6.3: Helical TWT
179
Figure 6.4: O type TWT
180
5 a =(4.4* %pi ) /180; // r a d i a n s
6 c =3 D +8 //m/ s
7 f =8 D +9; // Hz
8 al =2; //Np/m
9
10 // Phase v e l o c i t y
11 vp = c * sin ( a ) ;
12
13 // P r o p o g a t i o n c o n s t a n t
14 be =(2* %pi * f ) / vp ;
15
16 gam = al +( %i * be ) ;
17 disp ( gam , ’ P r o p o g a t i o n c o n s t a n t : ’ ) ;
Scilab code Exa 6.5 Cavity coupled

2 // Example 6 . 5
3 clc ;
4 // Given
5 Vc =11 D +3; //V
6 Ir =0.85; //A
7 V0 =31 D +3; //V
8 Pout =50 D +3; //W
9 I =7; //A
10
11 // E l e c t r o n i c e f f i c i e n c y
12 ne = Pout /( V0 * I ) ;
13 disp ( ’% ’ , ne *100 , ’ E l e c t r o n i c e f f i c i e n c y : ’ ) ;
14
15 // O v e r a l l e f f i c i e n c y
16 no = Pout /( Vc *( I - Ir ) ) ;
17 disp ( ’% ’ , no *100 , ’ O v e r a l l e f f i c i e n c y : ’ ) ;
181
Figure 6.5: Cavity coupled
182
Figure 6.6: O Type Backward Wave amplifier
18
19 // Answer f o r e l e c y t r o n i c e f f i c i e n c y s h o u l d be 2 3 . 0 4%
but i t i s g i v e n a s 3 6 . 4 i n book
Scilab code Exa 6.6 O Type Backward Wave amplifier

2 // Example 6 . 6
183
Figure 6.7: Multicavity TWT
3 clc ;
4 // Given
5 I0 =0.95; //A
6 V0 =7 D +3; //V
7 Z0 =20; //ohm
8 N =20;
9
10 C =(( I0 * Z0 ) /(4* V0 ) ) ^(1/3) ;
11 // Gain
12 Ap = -9.54+(47.3* C * N ) ;
13 disp ( ’ dB ’ ,Ap , ’ Gain : ’ ) ;
184
Scilab code Exa 6.7 Multicavity TWT

2 // Example 6 . 7
3 clc ;
4 // Given
5 Vc =12 D +3; //V
6 V0 =30 D +3; //V
7 Pout =60 D +3; //W
8 I =7.5; //A
9
11 ne = Pout /( V0 * I ) ;
13
15 no = Pout /( Vc * I ) ;
16 disp ( ’% ’ , no *100 , ’ O v e r a l l e f f i c i e n c y : ’ ) ;
Scilab code Exa 6.8 Gridded TWT

2 // Example 6 . 8
3 clc ;
4 // Given
5 Vc =20 D +3; //V
6 V0 =32 D +3; //V
7 Pout =75 D +3; //W
8 I =7; //A
9
11 ne = Pout /( V0 * I ) ;
185
Figure 6.8: Gridded TWT
186
Figure 6.9: Helix TWT
13
15 no = Pout /( Vc * I ) ;
16 disp ( ’% ’ , no *100 , ’ O v e r a l l e f f i c i e n c y : ’ );
Scilab code Exa 6.9 Helix TWT
187
2 // Example 6 . 9
3 clc ;
4 // Given
5 I0 =500 D -3; //A
6 V0 =10 D +3; //V
7 Z0 =25; //ohm
8 l =.20; //m
9 f =5.93 D +9; // Hz
10
11 // Gain p a r a m e t e r
12 C =(( I0 * Z0 ) /(4* V0 ) ) ^(1/3) ;
14
15 N =( l * f ) /(0.593 D +6* sqrt ( V0 ) ) ;
16 // Gain
17 Ap = -9.54+(47.3* C * N ) ;
18 disp ( ’ dB ’ ,Ap , ’ Gain o f TWT: ’ ) ;
Scilab code Exa 6.10 Low Power TWT

2 // Example 6 . 1 0
3 clc ;
4 // Given
5 Pout =250; //W
6 n =0.15;
7 V0 =7.5 D +3; //V
8 f =6.15 D +9; // Hz
9 c =3 D +8; //m/ s
10
11 // ( i ) I n p u t Power
12 Pi = Pout / n ;
188
189
Figure 6.10: Low Power TWT

13 disp ( ’W’ ,Pi , ’ I n p u t Power : ’ ) ;
14
15 // ( i i ) Beam c u r r e n t
16 I0 = Pi / V0 ;
17 disp ( ’A ’ ,I0 , ’ Beam c u r r e n t : ’ ) ;
18
19 // ( i i i ) Beam v e l o c i t y
20 vb =0.593 D +6* sqrt ( V0 ) ;
21 disp ( ’m/ s ’ ,vb , ’ Beam v e l o c i t y : ’ ) ;
22
23 // ( i v ) R a d i u s o f h e l i x
24 a =(2* vb ) /(2* %pi * f ) ;
25 disp ( ’m ’ ,a , ’ R a d i u s o f h e l i x : ’ ) ;
26
27 // ( v ) E l e c t r o n beam r a d i u s
28 r =(3* a ) /4;
29 disp ( ’m ’ ,r , ’ E l e c t r o n beam r a d i u s : ’ ) ;
30
31 // ( v i ) P i t c h o f h e l i x
32 p =(2* %pi * a * vb ) / c ;
33 disp ( ’m ’ ,p , ’ P i t c h o f h e l i x : ’ ) ;
34
35 // ( v i i ) C u r r e n t d e n s i t y
36 J0 = I0 /( %pi * r * r ) ;
37 disp ( ’ kA/ msqr ’ , J0 /1000 , ’ C u r r e n t d e n s i t y : ’ ) ;
38
39 // ( v i i i ) M a g n e t i c f i e l d f o r beam c o n f i n e m e n t
40 B =(4*8.3 D -4* sqrt ( I0 /( r * r * sqrt ( V0 ) ) ) ) ;
41 disp ( ’mT ’ , round ( B *1000) , ’ M a g n e t i c f i e l d f o r beam
confi nement : ’ );
190
Figure 6.11: TWT

2 // Example 6 . 1 1
3 clc ;
4 // Given
5 I0 =30 D -3; //A
6 V0 =3 D +3; //V
7 Z0 =10; //ohm
8 l =0.1624; //m
9 f =10 D +9; // Hz
10 C =(( I0 * Z0 ) /(4* V0 ) ) ^(1/3) ;
11 N =( l * f ) /(0.593 D +6* sqrt ( V0 ) ) ;
12
13 // Gain
14 Ap = -9.54+(47.3* C * N ) ;
15 disp ( ’ dB ’ ,Ap , ’ Gain : ’ ) ;
16
17 ve =0.593 D +6* sqrt ( V0 ) ;
18 be =(2* %pi * f ) / ve ;
19
191
21 gam1 =(( - sqrt (3) * be * C ) /2) +( %i * be *(2+ C ) ) /2;
22 gam2 =(( sqrt (3) * be * C ) /2) +( %i * be *(2+ C ) ) /2;
23 gam3 = %i * be *(1 - C ) ;
24 gam4 = - %i * be *(1 -(( C * C * C ) /4) ) ;
25
: ’ );
27
f o r gam 3 and 4 h e n c e a n s w e r s d o n t match

2 // Example 6 . 1 2
3 clc ;
4 // Given
5 I0 =35 D -3; //A
6 V0 =4 D +3; //V
7 Z0 =20; //ohm
8 f =10 D +9; // Hz
9
10 // ( i ) Gain p a r a m e t e r
11 C =(( I0 * Z0 ) /(4* V0 ) ) ^(1/3) ;
13
14 ve =0.593 D +6* sqrt ( V0 ) ;
15 be =(2* %pi * f ) / ve ;
16
18 gam1 =(( - sqrt (3) * be * C ) /2) +( %i * be *(2+ C ) ) /2;
19 gam2 =(( sqrt (3) * be * C ) /2) +( %i * be *(2+ C ) ) /2;
192
Figure 6.12: TWT
193
20 gam3 = %i * be *(1 - C ) ;
21 gam4 = - %i * be *(1 -(( C * C * C ) /4) ) ;
22
: ’ );
24
h e n c e a n s w e r s d o n t match
194
Chapter 7
Cross Field Microwave Tubes

M Type
Scilab code Exa 7.1 X band Magnetron

2 // Example 7 . 1
3 clc ;
4 // Given
5 f =10 D +9; // Hz
6 C =2.5 D -12; //F
7 Gr =2 D -4; //mho
8 Ge =0.025 D -3; //mho
9 Ploss =18.5 D +3; //W
10 V0 =5.5 D +3; //V
11 I0 =4.5; //A
12
13 w =2* %pi * f ;
14
15 // ( i ) Unloaded Q
16 Qun =( w * C ) / Gr ;
17 disp ( Qun , ’ Unloaded q u a l i t y f a c t o r : ’ ) ;
195
Figure 7.1: X band Magnetron
196
18
19 // E x t e r n a l Q
20 Qe =( w * C ) / Ge ;
21 disp ( Qe , ’ E x t e r n a l q u a l i t y f a c t o r : ’ ) ;
22
23 // ( i i ) C i r c u i t e f f c i e n c y
24 n =1/(1+( Qe / Qun ) ) ;
25 disp ( ’% ’ ,n *100 , ’ C i r c u i t e f f c i e n c y : ’ ) ;
26
27 // E l e c t r o n i c e f f c i e n c y
28 ne =1 -( Ploss /( V0 * I0 ) ) ;
29 disp ( ’% ’ , ne *100 , ’ E l e c t r o n i c e f f c i e n c y : ’ ) ;
30
31 // Answer f o r Qe i s g i v e n a s 6 2 8 5 . 6 but i t s h o u l d be
6283.1
Scilab code Exa 7.2 Cylindrical Magnetron

2 // Example 7 . 2
3 clc ;
4 // Given
5 V0 =25 D +3; //V
6 ebym =1.76 D +11;
7 B0 =0.0336; //T
8 a =5 D -2; //m
9 b =10 D -2; //m
10
11 // ( i ) Cut o f f v o l t a g e
12 x =( b /(( b * b ) -( a * a ) ) ) ^2;
13 V =( ebym * B0 * B0 ) /(8* x ) ;
14 disp ( ’KV ’ ,V /1000 , ’ Cut o f f v o l t a g e : ’ ) ;
15
197
Figure 7.2: Cylindrical Magnetron
198
Figure 7.3: Cylindrical Magnetron
16 // ( i i ) Cut o f f m a g n e t i c f i e l d
17 y =((8* V0 * x ) / ebym ) ;
18 B = sqrt ( y ) ;
19 disp ( ’mT ’ ,B *1000 , ’ Cut o f f m a g n e t i c f i e l d : ’ );
Scilab code Exa 7.3 Cylindrical Magnetron

2 // Example 7 . 3
3 clc ;
199
4 // Given
5 Pout =250 D +3; //W
6 V0 =25 D +3; //V
7 I0 =25; //A
8 ebym =1.76 D +11;
9 B0 =0.035; //T
10 a =4 D -2; //m
11 b =8 D -2; //m
12
13
14 // ( i ) E f f i c i e n c y
15 n = Pout /( V0 * I0 ) ;
16 disp ( ’% ’ ,n *100 , ’ E f f i c i e n c y : ’ ) ;
17
18 // ( i i ) C y c l o t r o n f r e q u e n c y
19 f =( ebym * B0 ) /(2* %pi ) ;
20 disp ( ’ Ghz ’ ,f /10^9 , ’ C y c l o t r o n f r e q u e n c y : ’ ) ;
21
22 // ( i i i ) Cut o f f m a g n e t i c f i e l d
23 x =( b /(( b * b ) -( a * a ) ) ) ^2;
24 y =((8* V0 * x ) / ebym ) ;
25 B = sqrt ( y ) ;
26 disp ( ’mT ’ ,B *1000 , ’ Cut o f f m a g n e t i c f i e l d : ’ ) ;
27
28 // ( i v ) Cut o f f v o l t a g e
29 V =( ebym * B0 * B0 ) /(8* x ) ;
30 disp ( ’KV ’ , round ( V /1000) , ’ Cut o f f v o l t a g e : ’ ) ;
31
32 // Answer f o r C y c l o t r o n f r e q u e n c y i s i s g i v e n a s 9 . 8
GHz but i t s h o u l d be 0 . 9 8 GHz a s v a l u e o f B0
=0.035 not 0 . 3 5 as taken i n p a r t 2
Scilab code Exa 7.4 Conventional Magnetron
200
Figure 7.4: Conventional Magnetron
201
2 // Example 7 . 4
3 clc ;
4 // Given
5 Gr =3 D -4; //mho
6 Ge =3 D -5; //mho
7 Ploss =200 D +3; //W
8 V0 =22 D +3; //V
9 I0 =28; //A
10
11 // ( i ) C i r c u i t e f f c i e n c y
12 n =1/(1+( Gr / Ge ) ) ;
14
15 // ( i i ) E l e c t r o n i c e f f c i e n c y
16 ne =1 -( Ploss /( V0 * I0 ) ) ;
Scilab code Exa 7.5 Conventional Magnetron

2 // Example 7 . 5
3 clc ;
4 // Given
5 f =9 D +9; // Hz
6 C =2.5 D -12; //F
7 Gr =2 D -4; //mho
8 Ge =2.5 D -5; //mho
9 Ploss =18.5 D +3; //W
10 V0 =5.5 D +3; //V
11 I0 =4.5; //A
12
13 // ( i ) A n g u l a r r e s o n a n t f r e q u e n c y
202
Figure 7.5: Conventional Magnetron
203
14 w =2* %pi * f ;
15 disp ( ’ r a d / s ’ ,w , ’ A n g u l a r r e s o n a n t f r e q u e n c y : ’ ) ;
16
17 // ( i i ) Unloaded Q
18 Qun = round (( w * C ) / Gr ) ;
19 disp ( Qun , ’ Unloaded q u a l i t y f a c t o r : ’ ) ;
20
21 // ( i i i ) Loaded Q
22 Ql = round (( w * C ) /( Gr + Ge ) ) ;
23 disp ( Ql , ’ Loaded q u a l i t y f a c t o r : ’ ) ;
24
25 // ( i v ) E x t e r n a l Q
26 Qe =( w * C ) / Ge ;
27 disp ( Qe , ’ E x t e r n a l q u a l i t y f a c t o r : ’ ) ;
28
29 // ( v ) C i r c u i t e f f c i e n c y
30 n =1/(1+( Qe / Qun ) ) ;
32
33 // ( v i ) E l e c t r o n i c e f f c i e n c y
34 ne =1 -( Ploss /( V0 * I0 ) ) ;
36
37 // Answer f o r e x t e r n a l Q i s g i v e n a s 5 6 . 5 7 but i t
s h o u l d be 5 6 5 4 . 8
Scilab code Exa 7.6 Carcinotron

2 // Example 7 . 6
3 clc ;
4 // Given
5 f =4 D +9; // Hz
204
Figure 7.6: Carcinotron
205
6 V0 =25 D +3; //V
7 I0 =3; //A
8 B0 =0.3; //T
9 D =0.8;
10 Z0 =50; //ohm
11 ebym =1.76 D +11;
12
13 // ( i ) E l e c t r o n beam p h a s e c o n s t a n t
14 be =(2* %pi * f ) / sqrt (2* ebym * V0 ) ;
15 disp ( ’ r a d / s ’ ,be , ’ E l e c t r o n beam p h a s e c o n s t a n t : ’ ) ;
16
17 // ( i i ) Gain P a r a m e t e r
18 C =(( I0 * Z0 ) /(4* V0 ) ) ^(1/3) ;
19 disp (C , ’ Gain P a r a m e t e r : ’ ) ;
20
21 // ( i i i ) Length f o r o s c i l l a t i o n c o n d i t i o n
22 N =1.25/ D ;
23 l =(2* %pi * N ) / be ;
24 disp ( ’m ’ ,l , ’ Length f o r o s c i l l a t i o n c o n d i t i o n : ’ ) ;
Scilab code Exa 7.7 Frequency Aglile Magnetron

2 // Example 7 . 7
3 clc ;
4 // Given
5 N =20;
6 t =0.2 D -6; // s
7 DC =0.001; // Duty c y c l e
8
9 // ( i ) A g i l e e x c u r s i o n
10 A=N/t;
11 disp ( ’MHz ’ ,A /10^6 , ’ A g i l e e x c u r s i o n : ’ ) ;
206
Figure 7.7: Frequency Aglile Magnetron
207
12
13 // ( i i ) S i g n a l f r e q u e n c y
14 f = DC / t ;
15 disp ( ’ Khz ’ ,f /1000 , ’ S i g n a l f r e q u e n c y : ’ ) ;
16
17 // ( i i i ) A g i l e r a t e
18 R = f /(2* N ) ;
19 disp ( ’ Hz ’ ,R , ’ A g i l e Rate : ’ ) ;
Scilab code Exa 7.8 Cross field amplifier

2 // Example 7 . 8
3 clc ;
4 // Given
5 V0 =1.8 D +3; //V
6 I0 =1.3; //A
7 Pin =70; //W
8 n =0.22;
9
10 // ( i ) Power g e n e r a t e d
11 Pgen = n * I0 * V0 ;
12 disp ( ’W’ , Pgen , ’ Power g e n e r a t e d : ’ ) ;
13
14 // ( i i ) T o t a l RF power g e n e r a t e d
15 Pt = Pin + Pgen ;
16 disp ( ’W’ ,Pt , ’ T o t a l RF power g e n e r a t e d : ’ ) ;
17
18 // ( i i i ) Power g a i n
19 G = Pt / Pin ;
20 Gdb =10* log10 ( G ) ;
21 disp ( ’ dB ’ ,Gdb , ’ Power Gain : ’ ) ;
208
Figure 7.8: Cross field amplifier
209
Figure 7.9: Inverted coaxial Magnetron
Scilab code Exa 7.9 Inverted coaxial Magnetron

2 // Example 7 . 9
3 clc ;
210
4 // Given
5 V0 =10 D +3; //V
6 I0 =2; //A
7 b =4 D -2; //m
8 a =3 D -2; //m
9 B0 =0.01; //Wb/m2
10 ebym =1.759 D +11;
11
12 // Cut o f f v o l t a g e
13 x =1 -(( b * b ) /( a * a ) ) ;
14 V =( ebym *( B0 ^2) *( a ^2) *( x ^2) ) /8;
15 KV = V /1000; // K i l o v o l t s
16 disp ( ’KV ’ ,KV , ’ Cut o f f v o l t a g e : ’ ) ;
17
18 // M a g n e t i c f l u x d e n s i t y
19 y = - sqrt ((8* V0 ) / ebym ) ;
20 B = y /( a * x ) ;
21 disp ( ’T ’ ,B , ’ M a g n e t i c f l u x d e n s i t y : ’ ) ;
Scilab code Exa 7.10 Inverted coaxial Magnetron

2 // Example 7 . 1 0
3 clc ;
4 // Given
5 V0 =10 D +3; //V
6 I0 =2; //A
7 b =4 D -2; //m
8 a =3 D -2; //m
9 B0 =0.01; //Wb/m2
10 ebym =1.759 D +11;
11
211
Figure 7.10: Inverted coaxial Magnetron
212
Figure 7.11: Linear Magnetron
13 x =1 -(( b * b ) /( a * a ) ) ;
14 V =( ebym *( B0 ^2) *( a ^2) *( x ^2) ) /8;
15 disp ( ’KV ’ ,V /1000 , ’ Cut o f f v o l t a g e : ’ ) ;
16
18 y = - sqrt ((8* V0 ) / ebym ) ;
19 B = y /( a * x ) ;
20 disp ( ’T ’ ,B , ’ M a g n e t i c f l u x d e n s i t y : ’ ) ;
Scilab code Exa 7.11 Linear Magnetron
213
2 // Example 7 . 1 1
3 clc ;
4 // Given
5 e =1.6 D -19; // J
6 B0 =0.01; //Wb/m2
7 d =6 D -2; //m
8 V0 =20 D +3; //V
9 ebym =1.759 D +11;
10
11 // ( i ) H u l l c u t o f f v o l t a g e
12 Voc =( B0 * B0 * d * d * ebym ) /2;
13 disp ( ’KV ’ , Voc /1000 , ’ H u l l c u t o f f v o l t a g e : ’ ) ;
14
15 // ( i i ) H u l l m a g n e t i c f i e l d
16 Boc = sqrt ((2* V0 ) / ebym ) / d ;
17 disp ( ’mT ’ , Boc *1000 , ’ H u l l m a g n e t i c f i e l d : ’ ) ;
Scilab code Exa 7.12 Inverted Coaxial Magnetron

2 // Example 7 . 1 2
3 clc ;
4 // Given
5 V0 =10 D +3; //V
6 V01 =5 D +3; //V
7 I0 =2; //A
8 b =3 D -2; //m
9 a =2 D -2; //m
10 B0 =0.01; //Wb/m2
11 ebym =1.759 D +11;
12
214
Figure 7.12: Inverted Coaxial Magnetron
215
14 x =1 -(( b * b ) /( a * a ) ) ;
15 V =( ebym *( B0 ^2) *( a ^2) *( x ^2) ) /8;
16 KV = V /1000; // K i l o v o l t s
17 disp ( ’KV ’ ,KV , ’ Cut o f f v o l t a g e : ’ ) ;
18
20 y = - sqrt ((8* V01 ) / ebym ) ;
21 B = y /( a * x ) ;
22 disp ( ’Wb/m2 ’ ,B , ’ M a g n e t i c f l u x d e n s i t y : ’ ) ;
23
24 // Answer i n book i s wrong f o r M a g n e t i c f l u x d e n s i t y
a s a ∗ a , where a =2 , i s t a k e n a s 5 , which s h o u l d be
4
Scilab code Exa 7.13 Agile coaxial Magnetron

2 // Example 7 . 1 3
3 clc ;
4 // Given
5 N =15;
6 t =0.3 D -6; // s
7 DC =0.0011; // Duty c y c l e
8
9 // ( i ) A g i l e e x c u r s i o n
10 A=N/t;
11 disp ( ’MHz ’ ,A /10^6 , ’ A g i l e e x c u r s i o n : ’ ) ;
12
13 // ( i i ) P u l s e t o p u l s e f r e q u e n c y s e p e r a t i o n
14 fp =1/ t ;
15 disp ( ’ Mhz ’ , fp /10^6 , ’ P u l s e t o p u l s e f r e q u e n c y
s e p e r a t i o n : ’ );
16
216
Figure 7.13: Agile coaxial Magnetron
217
17 // ( i i i ) S i g n a l f r e q u e n c y
18 f = DC / t ;
19 disp ( ’ Khz ’ ,f /1000 , ’ S i g n a l f r e q u e n c y : ’ ) ;
20
21 // ( i v ) A g i l e r a t e
22 Tp = N / f ;
23 R =1/(2* Tp ) ;
24 disp ( ’ p s ’ ,R , ’ A g i l e Rate : ’ ) ;
218
Chapter 8
Microwave Solid State Control

Devices
Scilab code Exa 8.1 Single pole Switch

2 // Example 8 . 1
3 clc ;
4 // Given
5 Rf =0.5; //ohm
6 Rr =1; //ohm
7 Ls =0.3 D -9; //H
8 Cj =0.1 D -12; //F
9 f =3.18 D +9; // Hz
10 Z0 =50; //ohm
11
12 Zf = Rf +( %i * round (2* %pi * f * Ls ) ) ;
13 Zr = Rr +( %i *( round (2* %pi * f * Ls ) -(1/(2* %pi * f * Cj ) ) ) ) ;
14
15 // S e r i e s C o n f i g u r a t i o n
16 disp ( ’ S e r i e s C o n f i g u r a t i o n ’ ) ;
17
219
Figure 8.1: Single pole Switch
220
18 // I n s e r t i o n L o s s
19 x =(2* Z0 ) /((2* Z0 ) + Zf ) ;
20 x1 = sqrt (( real ( x ) ) ^2+( imag ( x ) ) ^2) ;
21 IN = -20* log10 ( x1 ) ;
22 disp ( ’ dB ’ ,IN , ’ I n s e r t i o n L o s s : ’ ) ;
23
24 // I s o l a t i o n L o s s
25 y =(2* Z0 ) /((2* Z0 ) + Zr ) ;
26 y1 = sqrt (( real ( y ) ) ^2+( imag ( y ) ) ^2) ;
27 IS = -20* log10 ( y1 ) ;
28 disp ( ’ dB ’ ,IS , ’ I s o l a t i o n L o s s : ’ ) ;
29
30 // Shunt C o n f i g u r a t i o n
31 disp ( ’ Shunt C o n f i g u r a t i o n ’ ) ;
32
34 a =(2* Zr ) /((2* Zr ) + Z0 ) ;
35 a1 = sqrt (( real ( a ) ) ^2+( imag ( a ) ) ^2) ;
36 INs = -20* log10 ( a1 ) ;
37 disp ( ’ dB ’ ,INs , ’ I n s e r t i o n L o s s : ’ ) ;
38
40 b =(2* Zf ) /((2* Zf ) + Z0 ) ;
41 b1 = sqrt (( real ( b ) ) ^2+( imag ( b ) ) ^2) ;
42 ISs = -20* log10 ( b1 ) ;
43 disp ( ’ dB ’ ,ISs , ’ I s o l a t i o n L o s s : ’ ) ;
44
45 // Answer f o r S e r i e s c o n f i g u r a t i o n i n s e r t i o n l o s s i s
0 . 0 5 8 but i s g i v e n a s 0 . 5 8 db
Scilab code Exa 8.2 Pin diode switches
221
Figure 8.2: Pin diode switches
222
2 // Example 8 . 2
3 clc ;
4 // Given
5 Rf =1; //ohm
6 Rr =4; //ohm
7 Ls =0.3 D -9; //H
8 Cj =0.1 D -12; //F
9 f =3.18 D +9; // Hz
10 Z0 =50; //ohm
11
12 Zf = Rf +( %i * round (2* %pi * f * Ls ) ) ;
13 Zr = Rr +( %i *( round (2* %pi * f * Ls ) -(1/(2* %pi * f * Cj ) ) ) ) ;
14
15 // S e r i e s C o n f i g u r a t i o n
16 disp ( ’ S e r i e s C o n f i g u r a t i o n ’ ) ;
17
19 x =(2* Z0 ) /((2* Z0 ) + Zf ) ;
20 x1 = sqrt (( real ( x ) ) ^2+( imag ( x ) ) ^2) ;
21 IN = -20* log10 ( x1 ) ;
22 disp ( ’ dB ’ ,IN , ’ I n s e r t i o n L o s s : ’ ) ;
23
25 y =(2* Z0 ) /((2* Z0 ) + Zr ) ;
26 y1 = sqrt (( real ( y ) ) ^2+( imag ( y ) ) ^2) ;
27 IS = -20* log10 ( y1 ) ;
28 disp ( ’ dB ’ ,IS , ’ I s o l a t i o n L o s s : ’ ) ;
29
30 // Shunt C o n f i g u r a t i o n
31 disp ( ’ Shunt C o n f i g u r a t i o n ’ ) ;
32
34 a =(2* Zr ) /((2* Zr ) + Z0 ) ;
35 a1 = sqrt (( real ( a ) ) ^2+( imag ( a ) ) ^2) ;
36 INs = -20* log10 ( a1 ) ;
37 disp ( ’ dB ’ ,INs , ’ I n s e r t i o n L o s s : ’ ) ;
38
223
Figure 8.3: Silicon switching diode
40 b =(2* Zf ) /((2* Zf ) + Z0 ) ;
41 b1 = sqrt (( real ( b ) ) ^2+( imag ( b ) ) ^2) ;
42 ISs = -20* log10 ( b1 ) ;
43 disp ( ’ dB ’ ,ISs , ’ I s o l a t i o n L o s s : ’ ) ;
Scilab code Exa 8.3 Silicon switching diode
224
2 // Example 8 . 3
3 clc ;
4 // Given
5 Vbd =1000; //V
6 f =30 D +9; // Hz
7 E =3 D +5; //V/cm
8 Cj =0.3 D -12; //F
9 er =11.8;
10 e0 =8.854 D -12;
11
12 W = Vbd / E ;
13 Wpi = W /100; //mu
14
15 // T o t a l s e r i e s r e s i s t a n c e
16 R =1/(2* %pi * f * Cj ) ;
17 disp ( ’ ohms ’ ,R , ’ T o t a l s e r i e s r e s i s t a n c e : ’ );
18
19 // J u n c t i o n Area
20 A =( Cj * Wpi ) /( e0 * er ) ;
21 disp ( ’ cm2 ’ ,A *10000 , ’ J u n c t i o n Area : ’ ) ;
Scilab code Exa 8.6 Parametric upconverter

2 // Example 8 . 6
3 clc ;
4 // Given
5 MQ =10;
6 M =0.4;
7 r =20;
8 Td =300; //K
9 T =290; //K
10
225
Figure 8.4: Parametric upconverter
226
11 x =( MQ * MQ ) / r ;
12 // Power Gain
13 Ap =( r * x ) /((1+ sqrt (1+ x ) ) ^2) ;
14 Apdb =10* log10 ( Ap ) ;
15 disp ( ’ dB ’ , Apdb , ’ Power g a i n : ’ ) ;
16
17 // N o i s e f i g u r e
18 z =( Td / T ) / sqrt (1+(( MQ * MQ ) / r ) ) ;
19 F =1+ z ;
20 Fdb =10* log10 ( F ) ;
21 disp ( ’ dB ’ ,F , ’ N o s i e f i g u r e : ’ ) ;
22
23 // Bandwidth
24 BW =2* M * sqrt ( r ) ;
25 disp ( BW , ’ Bandwidth : ’ ) ;
Scilab code Exa 8.7 Parametric amplifier

2 // Example 8 . 7
3 clc ;
4 // Given
5 MQ =10;
6 r =10;
7
8 x =( MQ * MQ ) / r ;
9
10 // Gain
11 Ap =( r * x ) /((1+ sqrt (1+ x ) ) ^2) ;
12 Apdb =10* log10 ( Ap ) ;
13 disp ( ’ dB ’ , Apdb , ’ Gain : ’ ) ;
227
Figure 8.5: Parametric amplifier
228
Figure 8.6: Negative resistance parametric amplifier
Scilab code Exa 8.8 Negative resistance parametric amplifier

2 // Example 8 . 8
3 clc ;
4 // Given
229
Figure 8.7: 330 stage CCD
5 Rs =1; //ohm
6 ws =5 D +9; // Hz
7 M =0.25;
8 C0 =2 D -12; //F
9
10 // ( i ) E f f e c t i v e Q
11 Q =1/( Rs * ws * C0 *(1 -( M * M ) ) ) ;
12 disp (Q , ’ E f f e c t i v e Q : ’ ) ;
230
Scilab code Exa 8.9 330 stage CCD

2 // Example 8 . 9
3 clc ;
4 // Given
5 e =0.0001;
6 s =330;
7
8 // Charge t r a n s f e r e f f c i e n c y
9 n =1 - e ;
10
11 // F i n a l c h a r g e p u l s e
12 // x=P/P0
13 x =(1 -( e * s ) ) ;
14 disp (x , ’ F i n a l c h a r g e p u l s e : ’ ) ;
Scilab code Exa 8.10 3 phase CCD

2 // Example 8 . 1 0
3 clc ;
4 // Given
5 Qmax =0.05 D -12; //C
6 f =10 D +6; // Hz
7 V =10; //V
8 n =3;
9
10 // Power d i s s p a t e d p e r b i t
11 P = n * f * V * Qmax ;
231
Figure 8.8: 3 phase CCD
232
Figure 8.9: Surface channel CCD
12 disp ( ’muW ’ ,P *10^6 , ’ Power d i s s p a t e d p e r b i t : ’ ) ;
Scilab code Exa 8.11 Surface channel CCD

2 // Example 8 . 1 1
233
3 clc ;
4 // Given
5 e0 =8.854 D -12;
6 er =3.9;
7 d =0.15 D -6; //m
8 e =1.6 D -19; // J
9 Nmax =2.2 D +16; //m−2
10 A =0.6 D -8; //m
11 P =0.67 D -3; //W
12 n =3;
13
14 // ( i ) J u n c t i o n c a p a c i t a n c e
15 Ci =( e0 * er ) / d ;
16
17 // Gate v o l t a g e
18 V =( Nmax * e ) / Ci ;
19 disp ( ’V ’ ,V , ’ Gate v o l t a g e : ’ ) ;
20
21 // ( i i ) Charge s t o r e d
22 Qmax = Nmax * e * A ;
23
24 // C l o c k f r e q u e n c y
25 f = P /( n * V * Qmax ) ;
26 disp ( ’MHz ’ ,f /10^6 , ’ C l o c k f r e q u e n c y : ’ ) ;
Scilab code Exa 8.12 3 phase CCD

2 // Example 8 . 1 2
3 clc ;
4 // Given
5 Qmax =0.06 D -12; //C
6 f =20 D +6; // Hz
234
Figure 8.10: 3 phase CCD
235
7 V =10; //V
8 n =3;
9
10 // Power d i s s p a t e d p e r b i t
11 P = n * f * V * Qmax ;
12 disp ( ’muW ’ ,P *10^6 , ’ Power d i s s p a t e d p e r b i t : ’ ) ;
Scilab code Exa 8.13 Surface channel CCD

2 // Example 8 . 1 3
3 clc ;
4 // Given
5 e0 =8.854 D -12;
6 er =4;
7 d =0.1 D -6; //m
8 si =0.85;
9 e =1.6 D -19; // J
10 Na =1 D +20;
11
12 Ci =( e0 * er ) / d ;
13 disp ( ’ F/m ’ ,Ci , ’ J u n c t i o n c a p a c i t a n c e : ’ ) ;
14
15 W = sqrt ((2* e0 * er * si ) /( e * Na ) ) ;
16 disp ( ’m ’ ,W , ’ D e p l e t i o n l a y e r w i d t h : ’ ) ;
236
Figure 8.11: Surface channel CCD
237
Chapter 9
Microwave Solid State

Generators and Amplifiers
Scilab code Exa 9.2 Bipolar transistor

2 // Example 9 . 2
3 clc ;
4 // Given
5 fc =5 D +9; // Hz
6 Em =2 D +7; //V/m
7 vs =4 D +3; //ms/ s
8 Xc =1; //ohm
9
10 //Maximum a l l o w a b l e power
11 Pm =(( Em * vs ) ^2) /(((2* %pi * fc ) ^2) * Xc ) ;
12 disp ( ’W’ ,Pm , ’ Maximum a l l o w a b l e power : ’ ) ;
238
Figure 9.1: Bipolar transistor
239
Figure 9.2: Heterojunction transistor
240
Scilab code Exa 9.3 Heterojunction transistor

2 // Example 9 . 3
3 clc ;
4 // Given
5 XeGe =4.0; // eV
6 XeGaAs =4.1; // eV
7 delEgGe =0.78; // eV
8 delEgGaAs =1.42; // eV
9
10 // C o n d u c t i o n band d i f f e r e n t i a l
11 delEc = XeGe - XeGaAs ;
12 disp ( ’ eV ’ , delEc , ’ C o n d u c t i o n band d i f f e r e n t i a l : ’ ) ;
13
14 // V a l e n c e band d i f f e r e n t i a l
15 delEv = delEgGaAs - delEgGe - delEc ;
16 disp ( ’ eV ’ , delEv , ’ V a l e n c e band d i f f e r e n t i a l : ’ ) ;
Scilab code Exa 9.4 GaAs FET

2 // Example 9 . 4
3 clc ;
4 // Given
5 S11 =0.89;
6 S12 =0.02;
7 S21 =3.1;
8 S22 =0.78;
9
10 del =( S11 * S22 ) -( S12 * S21 ) ;
11 K =[1 -( S11 ) ^2 -( S22 ) ^2+( del ) ^2;]/(2* S12 * S21 ) ;
12 if (K <1)
241
Figure 9.3: GaAs FET
242
13 disp ( ’ A m p l i f i e r i s p o t e n t i a l l y u n s t a b l e ’ ) ;
14 else
15 disp ( ’ A m p l i f i e r i s p o t e n t i a l l y s t a b l e ’ ) ;
16 end
Scilab code Exa 9.5 Microwave transistor

2 // Example 9 . 5
3 clc ;
4 // Given
5 S11 =0.40;
6 S12 =0.01;
7 S21 =2.00;
8 S22 =0.35;
9
10 ZL =20; //ohm
11 ZS =30; //ohm
12 Z0 = ZL + ZS ; //ohm
13
14 // R e f l e c t i o n c o e f f i c i e n t s o f s o u r c e and l o a d
15 TL =( ZL - Z0 ) /( ZL + Z0 ) ;
16 TLm = - TL ;
17 TS =( ZS - Z0 ) /( ZS + Z0 ) ;
18 TSm = - TS ;
19
20 // R e f l e c t i o n c o e f f i c i e n t s o f i n p u t and o u t p u t
21 Tin = S11 +(( S12 * S21 * TL ) /(1 -( S22 * TL ) ) ) ;
22 Tout = S22 +(( S12 * S21 * TS ) /(1 -( S22 * TS ) ) ) ;
23
25 x =(1 -( TSm ) ^2) /((1 -( S11 * TSm ) ) ^2) ; // V a l u e o f s h o u l d
be 1 . 1 4 5
243
Figure 9.4: Microwave transistor
244
26 y =( S21 * S21 ) ;
27 z =(1 -( TLm ) ^2) /((1 -( Tout * TLm ) ) ^2) ;
28 GT = x * y * z ;
30
32 z1 =1 -( Tout ) ^2;
33 GA =( x * y ) / z1 ;
35
36 // Power Gain
37 z2 =1 -( Tin ) ^2;
38 GP =( x * y ) / z2 ;
40
41 // A l l t h e end c a l c u l a t i o n s o f f i n d i n g g a i n a r e n o t
a c c u r a t e i n t h e book , h e n c e t h e a n s w e r s d o n t
match
Scilab code Exa 9.6 Transistor Amplifier

2 // Example 9 . 6
3 clc ;
4 // Given
5 S11 =0.60;
6 S12 =0.045;
7 S21 =2.50;
8 S22 =0.50
9 TS =0.5;
10 TL =0.4;
11 Vrms =10; //V
12 Z0 =50; //ohm
245
246
Figure 9.5: Transistor Amplifier

13
14 // ( i ) R e f l e c t i o n c o e f f i c i e n t s o f i n p u t and o u t p u t
15 Tin = S11 +(( S12 * S21 * TL ) /(1 -( S22 * TL ) ) ) ;
16 Tout = S22 +(( S12 * S21 * TS ) /(1 -( S22 * TS ) ) ) ;
17 disp ( Tin , ’ R e f l e c t i o n c o e f f i c i e n t s o f i n p u t : ’ ) ;
18 disp ( Tout , ’ R e f l e c t i o n c o e f f i c i e n t s o f o u t p u t : ’ ) ;
19
20 // ( i i ) G a i n s
22 x =(1 -( TS ) ^2) /((1 -( S11 * TS ) ) ^2) ;
23 y =( S21 * S21 ) ;
24 z =(1 -( TL ) ^2) /((1 -( Tout * TL ) ) ^2) ;
25 GT = x * y * z ;
27
29 z1 =1 -( Tout ) ^2;
30 GA =( x * y ) / z1 ;
32
33 // Power Gain
34 z2 =1 -( Tin ) ^2;
35 GP =( x * y ) / z2 ;
37
38 // C a l c u l a t i o n f o r Tout and G a i n s a r e wrong i n t h e
book , h e n c e t h e a n s w e r s d o n t match
39
40 // ( i i i ) Power a v a i l a b l e
41 Gt =9.4;
42 Pas =( sqrt (2) * Vrms ) ^2/(8* Z0 ) ;
43 Pal = Gt * Pas ;
44 disp ( ’W’ ,Pas , ’ Power a v a i l a b l e at s o u r c e : ’ );
45 disp ( ’W’ ,Pal , ’ Power a v a i l a b l e at l o a d : ’ );
247
248
Figure 9.6: Microwave transistor
Scilab code Exa 9.7 Microwave transistor

2 // Example 9 . 7
3 clc ;
4 // Given
5 S11 =0.90;
6 S12 =0;
7 S21 =2.40;
8 S22 =0.80;
9
10 Gmax =( S21 * S21 ) /((1 -( S11 ) ^2) *(1 -( S22 ) ^2) ) ;
11 Gdb =10* log10 ( Gmax ) ;
12 disp ( Gdb , ’ Maximum g a i n : ’ ) ;
Scilab code Exa 9.8 JEFT

2 // Example 9 . 8
3 clc ;
4 // Given
5 e =1.6 D -19;
6 Nd =1.1 D +23; //m−3
7 a =0.2 D -6; //m
8 er =11.8;
9 e0 =8.854 D -12;
10 mue =800 D -4; //m2/ Vs
11 Z =50 D -6;
12 L =8.5 D -6; //m
13 W0 =1; //V
249
Figure 9.7: JEFT
250
14 Vd =12; //V
15 Vg =1.5; //V
16
17 // ( i ) P i n c h o f f v o l t a g e and p i n c h o f f c u r r e n t
18 Vp =( e * Nd * a * a ) /(2* er * e0 ) ;
19 disp ( ’V ’ ,Vp , ’ P i n c h o f f v o l t a g e : ’ ) ;
20
21 Ip =( mue * e * e * Nd * Nd * Z * a * a ) /( e0 * er * L ) ;
22 disp ( ’A ’ ,Ip , ’ P i n c h o f f c u r r e n t : ’ ) ;
23 // Answer f o r I p i s 5 5 8 0 9 A but i t i s g i v e n a s
0.00558 A
24
25 // ( i i ) D r a i n and maximum d r a i n c u r r e n t
26 // Taking I p =5.58mA a s g i v e n i n book
27 Ip1 =0.00558; //A
28 x =(2/3) *((( Vd + Vg + W0 ) / Vp ) ^(3/2) ) ;
29 y =(2/3) *((( Vg + W0 ) / Vp ) ^(3/2) ) ;
30 Id = Ip1 *[( Vd / Vp ) -x + y ];
31 disp ( ’A ’ ,-Id , ’ D r a i n c u r r e n t : ’ ) ;
32
33 // S a t u r a t i o n C u r r e n t
34 Is = Ip1 *[(1/3) -(( Vg + W0 ) / Vp ) +((2/3) *((( Vg + W0 ) / Vp )
^(3/2) ) ) ];
35 disp ( ’A ’ ,Is , ’ D r a i n s a t u r a t i o n c u r r e n t : ’ ) ;
36
37 // ( i i i ) Cut o f f f r e q u e n c y
38 f =(2* mue * e * Nd * a * a ) /( %pi * er * e0 * L * L ) ;
39 disp ( ’GHz ’ ,f /10^9 , ’ Cutt o f f f r e q e n c y : ’ ) ;
Scilab code Exa 9.9 MESFET

2 // Example 9 . 9
251
Figure 9.8: MESFET
252
3 clc ;
4 // Given
5 e =1.6 D -19;
6 Nd =8 D +23; //m−3
7 a =0.12 D -6; //m
8 er =13.2;
9 e0 =8.854 D -12;
10
11 // P i n c h o f f v o l t a g e
12 Vp =( e * Nd * a * a ) /(2* er * e0 ) ;
13 disp ( ’V ’ ,Vp , ’ P i n c h o f f v o l t a g e : ’ ) ;
Scilab code Exa 9.10 Gunn device

2 // Example 9 . 1 0
3 clc ;
4 // Given
5 vd =2 D +5; //m/ s
6 L =10 D -6; //m
7 Ec =3.2 D +5; //V/m
8
9 // N a t u r a l f r e q u e n c y
10 f = vd / L ;
11 disp ( ’GHz ’ ,f /10^9 , ’ N a t u r a l f r e q u e n c y : ’ ) ;
12
13 // C r i t i c a l v o l t a g e
14 Vc = Ec * L ;
15 disp ( ’V ’ ,Vc , ’ C r i t i c a l v o l t a g e : ’ ) ;
253
Figure 9.9: Gunn device
254
Figure 9.10: Gunn oscillator
Scilab code Exa 9.11 Gunn oscillator

2 // Example 9 . 1 1
3 clc ;
4 // Given
5 n =0.08;
6 A =3 D -8; //m2
7 n0 =1 D +21; //m−3
8 e =1.6 D -19;
9 vd =1.5 D +5; //m/ s
10 M =3.2
11 E =350 D +3; //V
12 L =12 D -6; //m
255
13
14 // Power o u t p u t
15 Pout = n * A * n0 * e * vd * M * L * E ;
16 disp ( ’mW’ , Pout *1000 , ’ Power o u t p u t : ’ ) ;
Scilab code Exa 9.12 Tunnel diode

2 // Example 9 . 1 2
3 clc ;
4 // Given
5 G =15.85;
6 Rn =75; //ohm
7
8 Rl = Rn -( Rn / G ) ;
9 C = Rl +(10* %i ) ;
10 disp ( ’ ohms ’ ,C , ’ C a v i t y i m p e d a n c e : ’ ) ;
Scilab code Exa 9.13 Gunn diode

2 // Example 9 . 1 3
3 clc ;
4 // Given
5 e =1.6 D -19;
6 n1 =1 D +16; //m−3
7 mu1 =8000 D -4; //m2/ Vs
8 nu =1 D +14; //m−3
9 muu =180 D -4; //m2/ Vs
256
Figure 9.11: Tunnel diode
257
Figure 9.12: Gunn diode
258
10
11 // / C o n d u c t i v i t y
12 C = e *(( n1 * mu1 ) +( nu * muu ) ) ;
13 disp ( ’m mho ’ ,C *1000 , ’ C o n d u c t i v i t y : ’ ) ;

2 // Example 9 . 1 4
3 clc ;
4 // Given
5 e0 =8.854 D -12;
259
6 er =13.1;
7 vd =2.5 D +5; //m/ s
8 e =1.6 D -19;
9 mu =0.015; //m2/ Vs
10
11 // C r i t e r i a
12 n0L =( e0 * er * vd ) /( e * mu ) ;
13 disp ( ’mˆ−3 ’ ,n0L , ’ n0L s h o u l d be g r e a t e r t h a n ’ ) ;

2 // Example 9 . 1 5
3 clc ;
4 // Given
260
5 L =10 D -6; //m
6 f =10 D +9; // Hz
7 e =1.6 D -19;
8 n0 =2 D +20; //m3
9 E =3200 D +2; //V/m
10
11 // C u r r e n t d e n s i t y
12 vd = L * f ;
13 J = n0 * e * vd ;
14 disp ( ’A/m s q r ’ ,J , ’ C u r r e n t d e n s i t y : ’ ) ;
15
16 // N e g a t i v e e l e c t r o n m o b i l i t y
17 mu = - vd / E ;
18 disp ( ’ cm s q r / Vs ’ , mu *10000 , ’ N e g a t i v e e l e c t r o n
m o b i l i t y : ’ );
19
20 // Answer f o r N e g a t i v e e l e c t r o n m o b i l i t y i s 3 1 2 5 but
i t i s g i v e n as 3100
Scilab code Exa 9.17 IMPATT diode

2 // Example 9 . 1 7
3 clc ;
4 // Given
5 n =0.15;
6 Vdc =100; //V
7 Idc =200 D -3; //A
8 vd =2 D +5; //m/ s
9 L =6 D -6; //m
10
11 // ( i ) Maximum CW o u t p u t power
12 Pdc = Vdc * Idc ;
261
Figure 9.15: IMPATT diode
262
13 Pout = n * Pdc ;
14 disp ( ’W’ , Pout , ’ Maximum CW power o u t p u t : ’ ) ;
15
16 // ( i i ) R e s o n a n t f r e q u e n c y
17 f = vd /(2* L ) ;
18 disp ( ’GHz ’ ,f /10^9 , ’ R e s o n a n t f r e q u e n c y : ’ ) ;

2 // Example 9 . 1 8
3 clc ;
4 // Given
5 n =0.1;
6 Vdc =100; //V
7 Idc =100 D -3; //A
8 vd =2 D +5; //m/ s
9 L =5 D -6; //m
10 V0 =90; //V
11 k =3;
12
13 // ( i ) Maximum CW o u t p u t power
15 Pout = n * Pdc ;
16 disp ( ’W’ , Pout , ’ Maximum CW power o u t p u t : ’ ) ;
17
18 // ( i i ) R e s o n a n t f r e q u e n c y
19 f = vd /(2* L ) ;
20 disp ( ’ Hz ’ ,f , ’ R e s o n a n t f r e q u e n c y : ’ ) ;
21
22 // ( i i i ) T r a n s i t t i m e
23 T = L / vd ;
24 disp ( ’ s ’ ,T , ’ T r a n s i t t i m e : ’ ) ;
263
264
25
26 // ( i v ) A v a l a n c h e m u l t i p l i c a t i o n f a c t o r
27 M =1/(1 -(( Vdc / V0 ) ^ k ) ) ;
28 disp ( -M , ’ A v a l a n c h e m u l t i p l i c a t i o n f a c t o r : ’ ) ;

2 // Example 9 . 1 9
3 clc ;
4 // Given
5 n =0.1;
6 Vdc =100; //V
265
7 Idc =0.9; //A
8 t =0.01 D -9; // s
9 f =16 D +9; // Hz
10
11 // ( i ) Power o u t p u t
13 Pout = n * Pdc ;
14 disp ( ’W’ , Pout , ’ Power o u t p u t : ’ ) ;
15
16 // ( i i ) Duty c y c l e
17 D =( t /2) +(1/(2* f ) ) ;
18 disp ( ’ s ’ ,D , ’ Duty c y c l e : ’ ) ;

2 // Example 9 . 2 0
3 clc ;
4 // Given
5 Cj =0.5 D -12; //F
6 Lp =0.5 D -9; //H
7 Irf =0.65; //A
8 Rl =2; // ohms
9 Vbd =80; //V
10 Idc =0.08; //A
11
12 // R e s o n a n t f r e q u e n c y
13 f =1/(2* %pi * sqrt ( Cj * Lp ) ) ;
14 disp ( ’ Hz ’ ,f , ’ R e s o n a n t f r e q u e n c y : ’ ) ;
15
17 Pout =( Irf * Irf * Rl ) /2;
18 Pin = Vbd * Idc ;
266
267
19 n =( Pout *100) / Pin ;
20 disp ( ’% ’ ,n , ’ E f f i c i e n c y : ’ ) ;
Scilab code Exa 9.21 TRAPATT diode

2 // Example 9 . 2 1
3 clc ;
4 // Given
5 J =25 D +7; //A/m;
6 Na =2.5 D +21; //m3
7 e =1.6 D -19;
8
9 // A v l a n c e z o n e v e l o c i t y
10 vz = J /( Na * e ) ;
11 disp ( ’m/ s ’ ,vz , ’ A v l a n c h e z o n e v e l o c i t y : ’ ) ;
Scilab code Exa 9.22 BARITT diode

2 // Example 9 . 2 2
3 clc ;
4 // Given
5 e =1.6 D -19;
6 N =4 D +21; //m
7 L =10 D -6; //m
8 e0 =8.854 D -12;
9 er =11;
10
268
Figure 9.19: TRAPATT diode
269
Figure 9.20: BARITT diode
270
Figure 9.21: Laser
11 // Breakdown v o l t a g e
12 Vbd =( e * N * L * L ) /( e0 * er ) ;
13 disp ( ’V ’ , round ( Vbd ) , ’ Breakdown v o l t a g e : ’ ) ;
14
15 // Breakdown e l e c t r i c f i e l d
16 E = Vbd / L ;
17 disp ( ’V/m ’ ,E , ’ Breakdown e l e c t r i c f i e l d : ’ ) ;
Scilab code Exa 9.23 Laser
271
2 // Example 9 . 2 3
3 clc ;
4 // Given
5 lam =8000 D -10; //m
6 a =0.5 D -2; //m
7 D =4 D +8; //m
8
9 // A n g u l a r S p r e a d
10 t =(1.22* lam ) / a ;
11 disp ( ’ r a d ’ ,t , ’ A n g u l a r s p r e a d : ’ ) ;
12
13 // A e r i a l s p r e a d
14 A = %pi *(( D * t ) ^2) ;
15 disp ( ’m s q r ’ ,A , ’ A e r i a l s p r e a d : ’ ) ;
16
17
18 // Answer f o r A i s g i v e n a s 193 m s q r but i t i s 1.915
D+10 m s q r
Scilab code Exa 9.24 Laser

2 // Example 9 . 2 4
3 clc ;
4 // Given
5 E =10; //W
6 T =1 D -9; // s
7 c =3 D +8; //m/ s
8 lam =650 D -9; //m
9
10 // P u l s e Power
11 P=E/T;
272
Figure 9.22: Laser
12 disp ( ’W’ ,P , ’ P u l s e Power : ’ ) ;

13
14 //Q v a l u e
15 Q =( c * T ) / lam ;
16 disp (Q , ’Q v a l u e : ’ ) ;
Scilab code Exa 9.25 Heterojunction laser

2 // Example 9 . 2 5
3 clc ;
4 // Given
5 h =6.626 D -34;
6 c =3 D +8; //m/ s
273
Figure 9.23: Heterojunction laser
7 e =1.6 D -19;
8 Eg =1.85; // eV
9
10 // Wavelenght e m i t t e d
11 lam =( h * c ) /( Eg * e ) ;
12 lamarm = lam *1 D +10;
13 disp ( ’A ’ , round ( lamarm ) , ’ Wavelenght e m i t t e d : ’ ) ;
274
Chapter 10
Striplines and Microstrip Lines
Scilab code Exa 10.1 Copper stripline

2 // Example 1 0 . 1
3 clc ;
4 // Given ,
5
6 z0 =50; //ohm
7 t =0.001; //mm
8 b =0.32; //cm
9 er =2.20;
10 tandel = 0.0005;
11 rs =0.026; //ohm
12 f =10 D +9; // Hz
13 c =3 D +8; //m/ s e c
14
15 p = sqrt ( er ) * z0 ;
16 // As p<120
17 w = b *[((30* %pi ) / p ) -0.441];
18 disp ( ’ cm ’ ,w , ’ Width ’ ) ;
19
275
Figure 10.1: Copper stripline
276
20 // A t t e n u a t i o n
21 k ={(2* %pi * f * sqrt ( er ) ) / c };
22 ad =( k * tandel ) /2;
23
24 // and
25 A =1+{(2* w ) /( b - t ) }+[{( b + t ) /(( b - t ) * %pi ) }* log (((2* b ) -t )
/ t ) ];
26 // Hence
27 ac =(2.7 D -3* rs * er * z0 * A ) /{30* %pi *( b - t ) *1 D -2};
28 // T o t a l a t t e n u t i o n
29 a = ad + ac ;
30
31 // T o t a l a t t e n u t i o n i n db
32 x = exp ( a ) ;
33 alp =20* log10 ( x ) ; // db /m
34
35 // T o t a l a t t e n u t i o n i n db / lambda :
36 lam = c /( sqrt ( er ) * f ) ;
37 lamm = lam *1 D +2;
38 alph = alp / lamm ;
39 disp ( ’ db / lambda ’ , alph , ’ T o t a l a t t e n u t i o n i n db / lambda
: ’ );
40
41
42 // Answer i n book f o r a l p h i s g i v e n a s 0 . 8 5 6 but i t
s h o u l d be 0 . 0 8 5 6 a s v a l u e o f f i s t a k e n a s 10D+10
but i t s h o u l d be 10D+9
Scilab code Exa 10.2 Microstrip line

2 // Example 1 0 . 2
3 clc ;
277
Figure 10.2: Microstrip line
278
4 // Given ,
5 er =9.7;
6 h =0.25; //mm
7 w =0.25; //mm
8 f =5 D +9; // Hz
9 c =3 D +8; //m/ s
10
11 // ( i ) D i e l e c t r i c c o n s t a n t
12 dc =(( er +1) /2) +((( er -1) /2) *(1/ sqrt (1+12* h / w ) ) ) ;
13 disp ( dc , ’ D i e l e c t r i c c o n s t a n t : ’ ) ;
14
15 // ( i i ) Phase c o n s t a n t
16 lam0 = c / f ;
17 pc = sqrt ( dc ) *(2* %pi / lam0 ) ;
18 disp ( ’ r a d /m ’ , pc /100 , ’ Phase c o n s t a n t : ’ ) ;
19
20 // ( i i i ) M i c r o s t r i p w a v e l e n g t h
21 lams = lam0 / sqrt ( dc ) ;
22 disp ( ’ cm ’ , lams *100 , ’ M i c r o s t r i p w a v e l e n g t h : ’ ) ;
23
24 // ( i v ) C a p a c i t a n c e p e r u n i t l e n g t h
25 e0 =8.854 D -12;
26 cap =(2* %pi * e0 ) / log ((8* h / w ) -( w /(4* h ) ) ) ;
27 disp ( ’ F/cm ’ ,cap , ’ C a p a c i t a n c e p e r u n i t l e n g t h : ’ ) ;
28
29 // ( v ) C h a r a c t e r s t i c Impedance
30 ci =(60/ sqrt ( dc ) ) * log ((8* h / w ) +( w /(4* h ) ) ) ;
31 disp ( ’ ohm ’ ,ci , ’ C h a r a c t e r s t i c i m p e d a n c e : ’ ) ;
Scilab code Exa 10.3 Microstrip

2 // Example 1 0 . 3
279
Figure 10.3: Microstrip
280
Figure 10.4: Stripline
3 clc ;
4 // Given ,
5 er =5.23;
6 w =10; // m i l s
7 t =2.8; // m i l s
8 h =7; // m i l s
9
10 dc =(( er +1) /2) +((( er -1) /2) *(1/ sqrt (1+12* h / w ) ) ) ;
12
13 // As w/h>1
14 ci =(120* %pi ) /( sqrt ( dc ) *(( w / h ) +1.393+0.667* log (( w / h )
+1.444) ) ) ;
281
Scilab code Exa 10.4 Stripline

2 // Example 1 0 . 4
3 clc ;
4 // Given ,
5
6 q =2.5;
7 dh =1.58;
8 er =9;
9 f =10;
10 c =3 D +8;
11
12 erff =(( er +1) /2) +((( er -1) /2) *((1+(12/ q ) ) ^( -1/2) ) ) ;
13 vp =( c / sqrt ( erff ) ) * erff ;
14 fe1 = c /( sqrt ( vp ) *2* dh * q ) ;
15 if f < fe1 then
16 disp ( ’ S t r i p s u p p o r t s TEM mode o n l y ’ ) ;
17 else
18 disp ( ’ S t r i p d o e s n o t s u p p o r t TEM mode o n l y ’ ) ;
19 end

2 // Example 1 0 . 5
3 clc ;
4 // Given ,
5
6 er =9.7;
282
283
7 h =0.5; //mm
8 w =0.5; //mm
9 lt =2 D -4;
10 t =0.02; //mm
11 f =5 D +9; // Hz
12 fg =5; //HZ
13 c =3 D +8;
14 rs =8.22 D -3* sqrt ( fg ) ;
15
16 // ( i ) D i e l e c t r i c c o n s t a n t
17 dc =(( er +1) /2) +((( er -1) /2) *(1/ sqrt (1+12* h / w ) ) ) ;
19
20 // ( i i ) C h a r a c t e r s t i c Impedance
21 ci =(60/ sqrt ( dc ) ) * log ((8* h / w ) +( w /(4* h ) ) ) ;
23
24 // ( i i i ) D i e l e c t r i c a t t e n u a t i o n
25 lam0 = c / f ;
26 alphd =( %pi / lam0 ) *( er / sqrt ( dc ) ) *(( dc -1) /( er -1) ) * lt ;
27 disp ( ’ Np/m ’ , alphd , ’ D i e l e c t r i c a t t e n u a t i o n : ’ ) ;
28
29 // C o n d u c t o r a t t e n u a t i o n
30 r1 =[0.94+(0.132*( w / h ) ) -(0.0062*(( w / h ) ^2) ) ]*[(1/ %pi )
+(1/( %pi ^2) ) * log ((4* %pi * w ) / t ) ]*( rs /( w *1 D -3) ) ;
31 r1m = r1 *1 D -2;
32 r2 =( w / h ) /[(( w / h ) +5.8+(0.03*( h / w ) ) ) ]*( rs /( w *1 D -3) ) ;
33 r2m = r2 *1 D -2;
34 alphc =( r1 + r2 ) /(2* ci ) ;
35 disp ( ’ Np/m ’ , alphc , ’ C o n d u c t o r a t t e n u a t i o n : ’ ) ;
36
37 // ( i v ) T o t a l a t t e n u a t i o n
38 A = alphc + alphd ;
39 Adb = A *8.686*1 D -2;
40 disp ( ’ db /cm ’ ,Adb , ’ T o t a l a t t e n u a t i o n : ’ ) ;
284

2 // Example 1 0 . 6
3 clc ;
4 // Given
5
6 sig =5.8 D +7;
7 f =10; //GHz
8 h =0.12 D -2; //m
9
10 q =62.8* h * sqrt ( f * sig ) ;
11 disp ( round ( q ) , ’ c o n d u c t o r Q o f t h e s t r i p l i n e : ’ ) ;
285
Figure 10.7: Parallel stripline
286
Scilab code Exa 10.7 Parallel stripline

2 // Example 1 0 . 7
3 clc ;
4 // Given
5 Er =6;
6 h =4 D -3; //m
7
8 // ( i ) W f o r Z0=50W
9 Z0 =50; //W
10 W =(120* %pi * h ) /( sqrt ( Er ) * Z0 ) ;
11 disp ( ’mm’ ,W *1000 , ’ R e q u i r e d Width : ’ ) ;
12
13 // ( i i ) S t r i p l i n e c a p a c i t a n c e
14 E0 =8.854 D -12;
15 C =( E0 * Er * W ) / h ;
16 disp ( ’ pF/m ’ ,C *10^12 , ’ S t r i p l i n e c a p a c i t a n c e : ’ ) ;
17
18 // ( i i i ) S t r i p l i n e i n d u c t a n c e
19 Mu0 =4* %pi *10 D -7;
20 L =( Mu0 * h ) / W ;
21 disp ( ’muH/m ’ ,L *10^5 , ’ S t r i p l i n e i n d u c t a n c e : ’ ) ;
22
23 // ( i v ) Phase v e l o c i t y
24 c =3 D +8;
25 vp = c / sqrt ( Er ) ;
26 disp ( ’m/ s ’ ,vp , ’ Phase v e l o c i t y ’ ) ;
Scilab code Exa 10.8 Stripline coupler
287
Figure 10.8: Stripline coupler
288
2 // Example 1 0 . 8
3 clc ;
4 // Given
5 cl =3 D +8; //m/ s
6 f =5 D +9; // Hz
7 Er =9;
8 C = -10; // db
9 Z0 =50; //ohm
10 // Length
11 L =( cl / f ) /(4* sqrt ( Er ) ) ;
12 disp ( ’ cm ’ ,L *100 , ’ Length : ’ ) ;
13
14 // C o u p l i n g c o e f f i c i e n t
15 C0 =10^( C /20) ;
16 disp ( C0 , ’ C o u p l i n g c o e f f i c i e n t : ’ ) ;
17
18 // Even and odd mode i m p e d a n c e
19 Z0e =( Z0 * sqrt (1+ C0 ) ) / sqrt (1 - C0 ) ;
20 disp ( ’ ohm ’ ,Z0e , ’ Even mode i m p e d a n c e : ’ ) ;
21
22
23 Z0o =( Z0 * sqrt (1 - C0 ) ) / sqrt (1+ C0 ) ;
24 disp ( ’ ohm ’ ,Z0o , ’ Odd mode i m p e d a n c e : ’ ) ;
Scilab code Exa 10.9 Branch coupler

2 // Example 1 0 . 9
3 clc ;
4 // Given
5 Z0 =50; //ohm
6 C =3; // db
289
Figure 10.9: Branch coupler
290
7
8 // L i n e i m p e d a n c e
9 Z01sqr =(1 -(10^( C / -10) ) ) ;
10 Z01 = sqrt ( Z0 * Z0 * Z01sqr ) ;
11 disp ( ’ ohm ’ ,Z01 , ’ Z01 : ’ ) ;
12
13 Z02 = Z01 /( sqrt (1 -(1/ sqrt (2) ) ^2) ) ;
14 disp ( ’ ohm ’ , round ( Z02 ) , ’ Z02 : ’ ) ;
Scilab code Exa 10.10 Broadside stripline

2 // Example 1 0 . 1 0
3 clc ;
4 // Given
5 W =6; //m
6 s =2.2; //m
7 b =4.8; //m
8 Er =2.2;
9
10 // Even and odd mode i m p e d a n c e
11 Z0e =((120* %pi ) *( b - s ) ) /(2* sqrt ( Er ) * W ) ;
12 disp ( ’ ohm ’ ,Z0e , ’ Even mode i m p e d a n c e : ’ ) ;
13
14
15 Z0o =( Z0e * s ) / b ;
16 disp ( ’ ohm ’ ,Z0o , ’ Odd mode i m p e d a n c e : ’ ) ;
17
18 // Mid band c o u p l i n g
19 x =( Z0e - Z0o ) /( Z0e + Z0o ) ;
20 C = -20* log10 ( x ) ;
21 disp ( ’ db ’ ,C , ’ Mid band c o u p l i n g : ’ ) ;
22
291
Figure 10.10: Broadside stripline
292
Figure 10.11: Paralle stripline
23 // Answer i n book f o r C i s g i v e n a s 5 4 . 2 but i t

s h o u l d be 8 . 6 0
Scilab code Exa 10.11 Paralle stripline

2 // Example 1 0 . 1 1
3 clc ;
4 // Given
5 Er =6;
6 d =3 D -3; //m
7 Z0 =50; //ohm
293
8 E0 =8.854 D -12; //F/m
9 Mu0 =4* %pi *10 D -7; //H/m
10
11 // ( i ) W
12 W =(377* d ) /( sqrt ( Er ) * Z0 ) ;
13 disp ( ’mm’ ,W *1000 , ’ R e q u i r e d Width : ’ ) ;
14
15 // ( i i ) S t r i p l i n e c a p a c i t a n c e
16 C =( E0 * Er * W ) / d ;
17 disp ( ’ pF/m ’ ,C *10^12 , ’ S t r i p l i n e c a p a c i t a n c e : ’ ) ;
18
19 // ( i i i ) S t r i p l i n e i n d u c t a n c e
20 L =( Mu0 * d ) / W ;
21 disp ( ’muH/m ’ ,L *10^6 , ’ S t r i p l i n e i n d u c t a n c e : ’ ) ;
22
23 // ( i v ) Phase v e l o c i t y
24 c =3 D +8;
25 vp = c / sqrt ( Er ) ;
26 disp ( ’m/ s ’ ,vp , ’ Phase v e l o c i t y ’ ) ;
Scilab code Exa 10.12 Shielded stripline

2 // Example 1 0 . 1 2
3 clc ;
4 // Given
5 Er =2.56;
6 w =25; // m i l s
7 t =14; // m i l s
8 d =70; // m i l s
9 E0 =8.854 D -12; //F/m
10
11 // ( i ) K f a c t o r
294
Figure 10.12: Shielded stripline
12 K =1/(1 -( t / d ) ) ;
13 disp (K , ’K f a c t o r : ’ ) ;
14
15 // ( i i ) F r i n g e c a p a c i t a n c e
16 C =[( E0 * Er ) *[2* K * log ( K +1) -(K -1) * log ( K ^2 -1) ]]/ %pi ;
17 disp ( ’ pF/m ’ ,C *10^12 , ’ F r i n g e c a p a c i t a n c e : ’ ) ;
18
19 // ( i i i ) C h a r e c t e r i s t i c Impedance
20 X =1/[(( w * K ) / d ) +( C /( E0 * Er ) ) ];
21 Z0 =(94.15* X ) / sqrt ( Er ) ;
22 disp ( ’ ohm ’ ,Z0 , ’ C h a r e c t e r i s t i c Impedance : ’ ) ;
23
24
25 // Answer i n book f o r Z0 i s g i v e n a s 5 0 . 2 9 but i t
s h o u l d be 5 1 . 7
295
Figure 10.13: Lossless stripline
296
Scilab code Exa 10.13 Lossless stripline

2 // Example 1 0 . 1 3
3 clc ;
4 // Given
5 Z0 =50; //ohm
6 // S i n c r r a t i o o f power i s 2 : 3
7 x1 =5/2;
8 y1 =5/3;
9 // Output Impedance
10 Z1 = x1 * Z0 ;
11 Z2 = y1 * Z0 ;
12 disp ( ’ ohm ’ ,Z1 , ’ Output Impedance 1 : ’ )
13 disp ( ’ ohm ’ ,Z2 , ’ Output Impedance 2 : ’ )
14
15 // I n p u t Impedance
16 Zin =[(( Z2 *2* Z2 ) /3) /(( Z2 +(2* Z2 ) /3) ) ];
17
18 // L o o k i n g i n t o Z1 , Z2 i s | | t o Z0
19 A1 =( Z2 * Z0 ) /( Z2 + Z0 ) ;
20
21 // L o o k i n g i n t o Z , Z2 i s | | t o Z0
22 A2 =( Z1 * Z0 ) /( Z1 + Z0 ) ;
23
24 // R e f l e c t i o n C o e f f c i e n t s
25 R1 =( A1 - Z1 ) /( A1 + Z1 ) ;
26 R2 =( A2 - Z2 ) /( A2 + Z2 ) ;
27
28 disp ( R2 , R1 , ’ R e f l e c t i o n C o e f f c i e n t s : ’ ) ;
297
Chapter 11
Microwave Integrated Circuits
Scilab code Exa 11.1 Costs

2 // Example 1 1 . 1
3 clc ;
4 // Given
5 fabc =10000; // Rs / w a f f e r
6 c =100;
7 y =40/100;
8 coc = fabc /( y * c ) ;
9 // C o s t o f one c h i p
10 disp ( ’ Rs ’ ,coc , ’ C o s t o f one c h i p : ’ ) ;
11
12 // Market C o s t
13 mc =2* coc ;
14 disp ( ’ Rs ’ ,mc , ’ Market c o s t o f one c h i p : ’ ) ;
298
Figure 11.1: Costs
299
Figure 11.2: Yield
300
Scilab code Exa 11.2 Yield

2 // Example 1 1 . 2
3 clc ;
4 // Given
5 c =5000; // Rs
6 S =0.6; //cm
7 // S i d e s
8 x =3; //cm
9 y =2.54; //cm
10 // b r e a k e v e n c o s t
11 bec =250;
12 // hence , c h i p s / w a f f e r s n e e d e d
13 cpw = c / bec ;
14 D=x*y;
15 // For g i v e n Area , a t l e a s t 40 c h i p s a r e r e q u i r e d
16 n =2* cpw ;
17
18 // D i a m e t e r
19 N = D /( sqrt (2) * S ) ;
20 // Lower round o f f
21 NN = floor ( N ) ;
22 // C h i p s p o s s i b l e
23 cp = NN ^2;
24
25 // Y i e l d
26 Y =( n / cp ) *100; // P e r c e n t
27 disp ( ’% ’ ,Y , ’ Y i e l d : ’ ) ;
301
Chapter 12
Microwave Measurements
Scilab code Exa 12.1 Microwave diode

2 // Example 1 2 . 1
3 clc ;
4 // Given
5 Is =0.1*(10^ -6) ; //A
6 Pi =0; //dBm
7 Cs =0.1*(10^ -12) ; //F
8 Ls =2*(10^ -9) ;
9 Cj =0.15*(10^ -12) ; //F
10 Rs =10; //ohm
11 T =293; //K
12 nktbye =25*(10^ -3) ; //V
13
14 // Rj
15 Rj =( nktbye / Is ) ;
16 disp ( ’Kohm ’ , Rj /1000 , ’ Rj : ’ ) ;
17
18 // Bi
19 Bi = nktbye /2;
302
Figure 12.1: Microwave diode
303
Figure 12.2: Detector mismatch
20 Bii = Bi *1000;
21 disp ( ’A/W’ ,Bii , ’ Bi : ’ ) ;
22
23 //Bv
24 Bv = Rj * Bii ;
25 disp ( ’V/W’ ,Bv , ’ Bv : ’ ) ;
Scilab code Exa 12.2 Detector mismatch

2 // Example 1 2 . 2
3 clc ;
4 // Given
5 vswr =4;
304
Figure 12.3: Transmission waveguide
6
7 modT =( vswr -1) /( vswr +1) ;
8 Lm = -10* log10 (1 -( modT * modT ) ) ; //dB
9 disp ( ’ dB ’ ,Lm , ’ Mismatch L o s s : ’ ) ;
10
11 // S e n s i t i v i t y r e d u c e s by a f a c t o r
12 Bvd =(1 -( modT * modT ) ) ;
13 Bvdp = Bvd *100;
14 disp ( ’% ’ , Bvdp , ’ V o l t g e s e n s i t i v i t y r e d u c e s by : ’ ) ;
Scilab code Exa 12.3 Transmission waveguide

2 // Example 1 2 . 3
3 clc ;
4 // Given
5 f =10 D +9; // Hz
6 c =3 D +10; //cm/ s
305
Figure 12.4: VSWR of waveguide
7 a =4; //cm
8 s =0.1; //cm
9 lmb = c / f ; //cm
10 lmbg = lmb /( sqrt (1 -(( lmb /(2* a ) ) ^2) ) ) ;
11 vswr = lmbg /( %pi * s ) ;
12 disp ( vswr , ’VSWR: ’ ) ;
13
14 // Answer i n book f o r lmbg i s g i v e n a s 3 . 4 9 but i t
s h o u l d be 3 . 2 3 and h e n c e t h e a n s w e r w i l l be 1 0 . 3
Scilab code Exa 12.4 VSWR of waveguide

2 // Example 1 2 . 4
3 clc ;
4 // Given
306
Figure 12.5: Directional couplers
5 delx =3.5; //cm

6 s =0.25; //cm
7
8 lmbg =2* delx ;
9 vswr = lmbg /( %pi * s ) ;
10 disp ( vswr , ’VSWR: ’ ) ;
Scilab code Exa 12.5 Directional couplers

2 // Example 1 2 . 5
3 clc ;
4 // Given
5 vswr =2;
6 Pin =4.5 D -3; //W
7
8 modT =( vswr -1) /( vswr +1) ;
307
9 // Power r e f l e c t e d ,
10 Pr =( modT ^2) * Pin ;
11 // As c o u p l e r s a m p l e s o n l y 1 / 1 0 0 0 t h power
12 Prr = Pr *1000;
13 disp ( ’W’ ,Prr , ’ R e f l e c t e d Power : ’ ) ;

2 // Example 1 2 . 6
3 clc ;
4 // Given
5 Z0 =50; //ohm
6 p =2.4;
7 L =0.313;
8 x =2* %pi * L ;
9 y = tan ( x ) ;
308
10
11 Zl =( Z0 *(1+( p * p * %i ) ) ) /( p +( p * %i ) ) ;
12 T =( Zl - Z0 ) /( Zl + Z0 ) ;
13 p = sqrt (( real ( T ) ) ^2+( imag ( T ) ) ^2) ;
14 disp (p , ’ R e f l e c t i o n c o e f f i c i e n t : ’ ) ;
309
2 // Example 1 2 . 7
3 clc ;
4 // Given
5 Zl =25+25* %i ; //ohm
6 Z0 =50; //ohm
7
8 T =( Zl - Z0 ) /( Zl + Z0 ) ;
9 p = sqrt (( real ( T ) ) ^2+( imag ( T ) ) ^2) ;
10 disp (p , ’ R e f l e c t i o n c o e f f i c i e n t : ’ ) ;
11
12 vswrr =(1+ p ) /(1 - p ) ;
13 disp ( vswrr , ’VSWR: ’ ) ;
14
15 // F r a c t i o n o f power d e l i v e r e d
16 Pd =1 -( p ^2) ;
17 Pdp = Pd *100;
18 disp ( ’% ’ ,Pdp , ’ F r a c t i o n o f power d e l i v e r e d : ’ ) ;

2 // Example 1 2 . 8
3 clc ;
4 // Given
5 d =2.4; //cm
6 lmbc =1.8;
7 c =3*10^10; //cm/ s
8
9 lmbg =2* d ;
10 lmb =( lmbg * lmbc ) /( sqrt ( lmbg ^2+ lmbc ^2) ) ;
11 // O p e r a t i n g f r e q u e n c y
12 f = c / lmb ;
13 disp ( ’GHz ’ ,f /10^9 , ’ O p e r a t i n g f r e q u e n c y : ’ ) ;
310
Figure 12.8: Rectangular waveguide
311
Figure 12.9: Three port circulator
Scilab code Exa 12.9 Three port circulator

2 // Example 1 2 . 9
3 clc ;
4 // Given
5 p =1.5;
6 IsL =1; //dB
7 InL =30; //dB
8
9 S21 =10^( - IsL /20) ;
10
11 // Assuming t g r e e p o r t s t o be i d e n t i c a l
12 S32 = S21 ;
312
Figure 12.10: Air filled cavity
13 S13 = S21 ;
14
15 // I s o l a t i o n s a r e a l s o t h e same
16 S31 =10^( - InL /20) ;
17 S23 = S31 ;
18 S12 = S31 ;
19
20 // R e f e l c t i o n c o e f f i c i e n t s a r e a l s o t h e same
21 T =( p -1) /( p +1) ;
22 S11 = T ;
23 S22 = T ;
24 S33 = T ;
25
26 S =[ S11 S12 S13 ; S21 S22 S23 ; S31 S32 S33 ];
27 disp (S , ’ M a t r i x i s : ’ ) ;
313
Scilab code Exa 12.10 Air filled cavity

2 // Example 1 2 . 1 0
3 clc ;
4 // Given
5 R1 =10.6; //GHz
6 R2 =8.30; //GHz
7 Q0 =8200;
8 Q0d =890;
9
10 Er =( R1 / R2 ) ^2;
11 disp ( Er , ’ D i e l e c t r i c c o n s t a n t ’ ) ;
12
13 Qd =( Q0 - Q0d ) /( Q0 * Q0d ) ;
14 disp ( Qd , ’ L o s s t a n g e n t o f d i e l e c t r i c ’ ) ;

2 // Example 1 2 . 1 1
3 clc ;
4 // Given
5 l0 =0.15; //cm
6 lmbg =2*2.24; //cm
7 le =1.14; //cm
8 a =2.286; //cm
9 d =2;
10
11 B0 =(2* %pi ) / lmbg ;
314
12 x = tan ( B0 * l0 ) /( B0 * l0 ) ;
13 // A l s o
14 x1 =( l0 * x ) / le ;
15 // C o r r e c t v a l u e s e e m s t o be
16 Bele =2.786;
17 e1 =(((( a / %pi ) ^2) *( Bele / le ) ^2) +1) ;
18 e2 =(((2* a ) / lmbg ) ^2) +1;
19 Er = e1 / e2 ;
20 disp ( Er , ’ Er : ’ ) ;
21
22
23 // Answer i n book f o r Er i s g i v e n a s 2 . 0 6 2 but i t
s h o u l d be 2 . 0 3 8
315

Microwave Engineering - G. S. Raghuvanshi

Uploaded by

Copyright:

Available Formats

Microwave Engineering - G. S. Raghuvanshi

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Microwave Engineering - G. S. Raghuvanshi

Uploaded by

Copyright:

Available Formats

What is the numbering policy used in this document?

What is the numbering policy used in this document?

What are the contents covered in this document?

What are the contents covered in this document?

Scilab Textbook Companion for

1 Funded by a grant from the National Mission on Education through ICT,

Title: Microwave Engineering

Publisher: Cengage Learning, New Delhi

Exa Example (Solved example)

Eqn Equation (Particular equation of the above book)

AP Appendix to Example(Scilab Code that is an Appednix to a particular

List of Scilab Codes 4

3 Microwave Network Analysis 83

4 Microwave Resonators and Waveguide Components 92

5 Microwave Tubes Klystrons 129

6 Microwave Travelling Wave Tubes O type 166

7 Cross Field Microwave Tubes M Type 187

8 Microwave Solid State Control Devices 211

9 Microwave Solid State Generators and Amplifiers 230

10 Striplines and Microstrip Lines 267

11 Microwave Integrated Circuits 290

12 Microwave Measurements 294

Exa 1.2 Lossless line . . . . . . . . . . . . . . . . . . . . . . . . 6

1.1 Lossless line . . . . . . . . . . . . . . . . . . . . . . . . . . . 7

2.1 Dominant mode . . . . . . . . . . . . . . . . . . . . . . . . . 36

3.1 Scattering matrix . . . . . . . . . . . . . . . . . . . . . . . . 84

4.1 Rectangular cavity resonator . . . . . . . . . . . . . . . . . . 93

5.1 Two Cavity Klystron . . . . . . . . . . . . . . . . . . . . . . 130

6.1 TWT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 167

7.1 X band Magnetron . . . . . . . . . . . . . . . . . . . . . . . 188

8.1 Single pole Switch . . . . . . . . . . . . . . . . . . . . . . . . 212

9.1 Bipolar transistor . . . . . . . . . . . . . . . . . . . . . . . . 231

10.1 Copper stripline . . . . . . . . . . . . . . . . . . . . . . . . . 268

11.1 Costs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 291

12.1 Microwave diode . . . . . . . . . . . . . . . . . . . . . . . . 295

Scilab code Exa 1.2 Lossless line

Figure 1.1: Lossless line

Scilab code Exa 1.3 Microwave line

Scilab code Exa 1.4 Quater wave transformer

Scilab code Exa 1.5 Parallel resonant circuit

Scilab code Exa 1.6 Lossless line

Scilab code Exa 1.7 Lossy cable

Scilab code Exa 1.8 Transmission Line

Scilab code Exa 1.10 Transmission Line

Scilab code Exa 1.11 Microwave line

Scilab code Exa 1.12 Active Device

Scilab code Exa 1.13 Transmission line

Scilab code Exa 1.15 Microwave line

Scilab code Exa 1.16 EM Plane

Scilab code Exa 1.17 Polyethylene

Scilab code Exa 1.18 Electromagnetic wave

Scilab code Exa 1.19 Range in sea water

Scilab code Exa 2.1 Dominant mode

Scilab code Exa 2.2 Rectangular Waveguide

Scilab code Exa 2.3 Rectangular Waveguide

Scilab code Exa 2.4 Wider dimension

Scilab code Exa 2.7 Rectangular waveguide

Scilab code Exa 2.8 Lossless Rectangular Waveguide

Scilab code Exa 2.9 Dimensions