Microwave Engineering M Kulkarni - Compress PDF

Scilab Textbook Companion for
Microwave and Radar Engineering

by M. Kulkarni1
Created by
Chandawar Saichander
ECE
Electronics Engineering
SASTRA UNIVERSITY
College Teacher
N. Raju
Cross-Checked by
K. V. P. Pradeep
June 13, 2014
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 and Radar Engineering
Author: M. Kulkarni
Publisher: Umesh Publications, New Delhi
Edition: 3
Year: 2008
ISBN: 81-88114-00-6
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
3 Transmission Lines 7
4 Microwave Transmission Lines 18
5 Cavity Resonators 45
6 Microwave Components 50
7 Microwave Measurements 60
8 Microwave Tubes and Circuits 64
9 Solid State Microwave Devices 83
10 Microwave Communication Systems 95
11 Radars 107
3
List of Scilab Codes
Exa 3.1 Terminating Impedance . . . . . . . . . . . . . . . . . 7

Exa 3.2 power . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
Exa 3.3 phase velocity . . . . . . . . . . . . . . . . . . . . . . 9
Exa 3.4 power . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
Exa 3.5 VSWR and Reflection coefficient . . . . . . . . . . . . 11
Exa 3.6 point of attachment . . . . . . . . . . . . . . . . . . . 12
Exa 3.7 Terminating Impedance . . . . . . . . . . . . . . . . . 13
Exa 3.8 VSWR and Impedance . . . . . . . . . . . . . . . . . 14
Exa 3.9 Reflection loss . . . . . . . . . . . . . . . . . . . . . . 15
Exa 3.10 Characteristic Impedance . . . . . . . . . . . . . . . . 16
Exa 4.2 attenuation and phase constants . . . . . . . . . . . . 19
Exa 4.3 breakdown power . . . . . . . . . . . . . . . . . . . . . 20
Exa 4.6 ratio of areas . . . . . . . . . . . . . . . . . . . . . . . 23
Exa 4.7 breadth of the guide . . . . . . . . . . . . . . . . . . . 25
Exa 4.8 guide wavelength . . . . . . . . . . . . . . . . . . . . . 26
Exa 4.10 area . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
Exa 4.11 possible modes . . . . . . . . . . . . . . . . . . . . . . 30
Exa 4.13 frequency . . . . . . . . . . . . . . . . . . . . . . . . . 32
Exa 4.15 possible modes . . . . . . . . . . . . . . . . . . . . . . 34
Exa 4.18 proof . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
4
Exa 4.19 amount of attenuation . . . . . . . . . . . . . . . . . . 38
Exa 4.20 max power handling capacity . . . . . . . . . . . . . . 39
Exa 4.21 maximum power . . . . . . . . . . . . . . . . . . . . . 40
Exa 4.22 peak value of electric field . . . . . . . . . . . . . . . . 41
Exa 5.1 minimum distance . . . . . . . . . . . . . . . . . . . . 45
Exa 5.2 lowest resonant frequency . . . . . . . . . . . . . . . . 46
Exa 5.3 resonant frequency . . . . . . . . . . . . . . . . . . . . 47
Exa 5.4 resonant frequency . . . . . . . . . . . . . . . . . . . . 48
Exa 6.2 distance to be shifted . . . . . . . . . . . . . . . . . . 50
Exa 6.3 Scattering parameters . . . . . . . . . . . . . . . . . . 51
Exa 6.4 powers in the remaining ports . . . . . . . . . . . . . . 52
Exa 6.5 power . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
Exa 6.6 Reflected power . . . . . . . . . . . . . . . . . . . . . 53
Exa 6.7 Scattering matrix . . . . . . . . . . . . . . . . . . . . 54
Exa 6.9 Scattering matrix . . . . . . . . . . . . . . . . . . . . 55
Exa 6.10 output powers . . . . . . . . . . . . . . . . . . . . . . 56
Exa 6.11 coupling isolation directivity . . . . . . . . . . . . . . 57
Exa 6.12 VSWR . . . . . . . . . . . . . . . . . . . . . . . . . . 57
Exa 6.13 phase shift . . . . . . . . . . . . . . . . . . . . . . . . 58
Exa 7.1 measured distance . . . . . . . . . . . . . . . . . . . . 60
Exa 7.2 VSWR and Reflected power . . . . . . . . . . . . . . . 61
Exa 7.3 VSWR . . . . . . . . . . . . . . . . . . . . . . . . . . 62
Exa 7.4 Reflected power . . . . . . . . . . . . . . . . . . . . . 62
Exa 8.1 electron velocity and etc . . . . . . . . . . . . . . . . . 64
Exa 8.2 power and etc . . . . . . . . . . . . . . . . . . . . . . 65
Exa 8.3 efficiency and etc . . . . . . . . . . . . . . . . . . . . . 66
Exa 8.6 electronic efficiency and etc . . . . . . . . . . . . . . . 70
Exa 8.8 cyclotron frequency and etc . . . . . . . . . . . . . . . 73
Exa 8.9 phase velocity and anode voltage . . . . . . . . . . . . 74
Exa 8.11 dc electron velocity and etc . . . . . . . . . . . . . . . 76
Exa 8.12 gap transit angle . . . . . . . . . . . . . . . . . . . . 77
5
Exa 8.14 cyclotron angular frequency and etc . . . . . . . . . . 80
Exa 8.16 repeller voltage and etc . . . . . . . . . . . . . . . . . 81
Exa 9.1 frequency . . . . . . . . . . . . . . . . . . . . . . . . . 83
Exa 9.2 threshold electric field . . . . . . . . . . . . . . . . . . 84
Exa 9.3 power gain . . . . . . . . . . . . . . . . . . . . . . . . 84
Exa 9.4 breakdown voltage and etc . . . . . . . . . . . . . . . 85
Exa 9.5 Avalanche zone velocity . . . . . . . . . . . . . . . . . 86
Exa 9.6 power gain . . . . . . . . . . . . . . . . . . . . . . . . 87
Exa 9.7 minimum voltage . . . . . . . . . . . . . . . . . . . . . 88
Exa 9.8 rational frequency . . . . . . . . . . . . . . . . . . . . 88
Exa 9.10 drift time . . . . . . . . . . . . . . . . . . . . . . . . . 90
Exa 9.11 breakdown voltage and etc . . . . . . . . . . . . . . . 91
Exa 9.12 max power gain and etc . . . . . . . . . . . . . . . . . 92
Exa 9.13 gain and etc . . . . . . . . . . . . . . . . . . . . . . . 93
Exa 10.1 radio horizon . . . . . . . . . . . . . . . . . . . . . . . 95
Exa 10.2 value of factor . . . . . . . . . . . . . . . . . . . . . . 96
Exa 10.3 power . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
Exa 10.4 power . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
Exa 10.5 antenna beam angle . . . . . . . . . . . . . . . . . . . 99
Exa 10.6 round trip time . . . . . . . . . . . . . . . . . . . . . . 100
Exa 10.7 figure of merit . . . . . . . . . . . . . . . . . . . . . . 101
Exa 10.8 CNR . . . . . . . . . . . . . . . . . . . . . . . . . . . 102
Exa 10.9 system noise temperature . . . . . . . . . . . . . . . . 102
Exa 10.10 HPBW . . . . . . . . . . . . . . . . . . . . . . . . . . 103
Exa 10.11 gain . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104
Exa 10.12 gain . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
Exa 10.13 power gain . . . . . . . . . . . . . . . . . . . . . . . . 106
Exa 11.1 max range . . . . . . . . . . . . . . . . . . . . . . . . 107
Exa 11.2 max range . . . . . . . . . . . . . . . . . . . . . . . . 108
Exa 11.3 RCS . . . . . . . . . . . . . . . . . . . . . . . . . . . . 109
Exa 11.4 Duty cycle and etc . . . . . . . . . . . . . . . . . . . . 110
Exa 11.5 max range . . . . . . . . . . . . . . . . . . . . . . . . 111
Exa 11.6 factor . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
6
Chapter 3
Transmission Lines
Scilab code Exa 3.1 Terminating Impedance
1 // c h a p t e r −3 p a g e 47 e x a m p l e 3 . 1
2 //
==================================================================
3 clc ;
4 clear ;
5
6 Z0 =100; // C h a r a c t e r i s t i c Impedance i n ohms
7 S =5; // V o l t a g e S t a n d i n g Wave R a t i o (VSWR)
8
9 //CALCULATION
10 Zm = Z0 * S ; // T e r m a i n a t i n g i m p e d a n c e a t a max o f t h e
v o l t a g e s t a n d i n g wave
11 Zl = Zm ; // L o a d i n g Impedance
12
13 //OUTPUT
14 mprintf ( ’ T e r m i n a t i n g i m p e d a n c e a t a maximum o f t h e
v o l t a g e s t a n d i n g wave i s Z l= %3 . 0 f ohms ’ , Zl ) ;
15
16 //====================END OF PROGRAM
========================================
7
Scilab code Exa 3.2 power
2 //
==================================================================
3 clc ;
4 clear ;
5
6 R =8; // R e s i s t a n c e o f a t r a n s m i s s i o n l i n e i n ohm/km
7 L =0.002; // I n d u c t a n c e o f a t r a n s m i s s i o n l i n e i n h e n r y
/km
8 C =0.002*(10^( -6) ) ; // C a p a c i t a n c e o f a t r a n s m i s s i o n
l i n e in Farads
9 G =0.07*(10^( -6) ) ; // C o n d u c t a n c e o f a t r a n s m i s s i o n
l i n e i n s i e m e n s /km
10 f =2000; // F r e q u e n c y i n Hz
11 w =2*( %pi ) * f ; // A n g u l a r F r e q u e n c y i n r a d / s e c
12 Vs =2; // I n p u t V o l t a g e i n v o l t s
13 l =500; // L e n g th o f T r a n s m i s s i o n l i n e i n km
14
15 //CALCULATIONS
16 Z0 = sqrt (( R +( w * L *( %i ) ) ) /( G +( w * C *( %i ) ) ) ) ; //
C h a r a c t e r i s t i c Impedance
17 x = real ( Z0 ) ;
18 y = imag ( Z0 ) ;
19 disp ( ’ C h a r a c t e r i s t i c Impedance i n ohms i s ’ ) ;
20 disp ( Z0 ) ;
21 g = sqrt (( R +( w * L *( %i ) ) ) *( G +( w * C *( %i ) ) ) ) ; // P r o p a g a t i o n
Constant
22 a = real ( g ) ; // A t t e n u a t i o n C o n s t a n t i n NP/km
23 b = imag ( g ) ; // Phase C o n s t a n t i n r a d /km
24 Is = Vs / Z0 ;
25 I0 = Is * exp ( -( g * l ) ) ; // Load c u r r e n t
8
26 m = sqrt (( real ( I0 ) ) ^2+( imag ( I0 ) ^2) ) ;
27 P =( m ^2) * x ; // Power d e l i v e r e d t o t h e l o a d i n w a t t s
28
29 //OUTPUT
30 mprintf ( ’ \ n A t t e n u a t i o n C o n s t a n t i s a=%1 . 6 f NP/km \
nPhase C o n s t a n t i s b=%1 . 6 f r a d /km \ nPower
d e l i v e r e d t o t h e l o a d i s P=%1 . 6 f w a t t s ’ ,a ,b , P ) ;
31
32 //===============END OF PROGRAM
================================
Scilab code Exa 3.3 phase velocity
2 //
==================================================================
3 clc ;
4 clear ;
5
6 w =4*( %pi ) ; // A n g u l a r F r e q u e n c y i n r a d / s e c
7 b =0.02543; // Phase C o n s t a n t i n r a d /km
8
9 //CALCULATION
10 Vp = w / b ; // Phase V e l o c i t y i n km/ s e c
11
12 //OUTPUT
13 mprintf ( ’ Phase V e l o c i t y i s Vp=%3 . 2 f km/ s e c ’ , Vp ) ;
14
15 //=========END OF PROGRAM=========================
16
17 //NOTE:CHECK THE CALCULATION PART GIVEN IN THE
TEXTBOOK
18 //GIVEN ANSWER 4 9 4 . 2 2 KM/SEC
19 //GETTING ANSWER 4 9 4 . 1 6 KM/SEC
9
2 //
==================================================================
3 clc ;
4 clear ;
5
6 f =37.5*10^6; // F r e q u e n c y i n Hz
7 c =3*10^8; // V e l o c i t y o f L i g h t i n m/ s e c
8 l1 =10; // L e n g th o f l i n e i n met
9 Vg =200; // G e n e r a t o r V o l t a g e i n v o l t s ( rms )
10 Zint =200; // I n t e r n a l R e s i s t a n c e o f G e n e r a t o r i n ohms
12 Zl =100; // Load i m p e d a n c e i n ohms
13
14 //CALCULATIONS
15 w = c / f ; //Wave L e n g th i n met
16 b =2*( %pi ) / w ;
17 l1 =(5/4) * w ; // For L o s s l e s s L i n e
18 Zi = Z0 *(( Zl +( Z0 *( %i ) * tan ( b * l1 ) ) ) /( Z0 +( Zl *( %i ) * tan ( b *
l1 ) ) ) ) ; // I n p u t Impedance a t G e n e r a t o r end
19 Vs = Vg *( Zi /( Zi + Z0 ) ) ; // V o l t a g e i n l i n e i n v o l t s
20 Is = Vg /( Zi + Z0 ) ; // C u r r e n t i n L i n e drawn from G e n e r a t o r
i n amps
21 Ps = Vs * Is ; // Power drawn i n l i n e
22 Pl = Ps ; // For L o s s l e s s L i n e s Power d e l i v e r e d t o l o a d
i s e q u a l t o t h e Power drawn i n l i n e
23 Il = sqrt (( Pl / Zl ) ) ; // C u r r e n t f l o w i n g i n t h e l o a d
24 m = real ( Il ) ; // Magnitude o f C u r r e n t f l o w i n g i n t h e
load
25 p = imag ( Il ) ; // Phase o f C u r r e n t f l o w i n g i n t h e l o a d
26
10
27 //CALCULATIONS
28 mprintf ( ’ \ n C u r r e n t drawn from G e n e r a t o r i s I s=%1 . 3 f
amps \ nMagnitude o f C u r r e n t f l o w i n g i n t h e l o a d
i s m=%1 . 3 f \ nPhase o f C u r r e n t f l o w i n g i n t h e l o a d
i s p=%2 . 2 f deg \ nPower d e l i v e r e d t o l o a d i s Pl=
%2 . 2 f w a t t s ’ ,Is ,m ,p , Pl ) ;
29
30 //=========================END OF PROGRAM
==============================================================
Scilab code Exa 3.5 VSWR and Reflection coefficient
2 //
==================================================================
3 clc ;
4 clear ;
5
7 f =300*10^6; // F r e q u e n c y i n Hz
8 Zl =50+(50*( %i ) ) ; // T e r m i n a t i n g l o a d i m p e d a n c e i n ohms
9 w =((3*10^8) / f ) ; //Wave L e n g th
10
11 //CALCULATIONS
12 p =(( Zl - Z0 ) /( Zl + Z0 ) ) ; // R e f l e c t i o n C o e f f i c i e n t ( Complex
Form )
13 y = real ( p ) ;
14 z = imag ( p ) ;
15 x = sqrt ( y ^2+ z ^2) ; // R e f l e c t i o n C o e f f i c i e n t V a l u e
16 s =((1+ x ) /(1 - x ) ) ; // V o l t a g e S t a n d i n g Wave R a t i o (VSWR)
17
18 //OUTPUT
19 mprintf ( ’ \ n R e f l e c t i o n C o e f f i c i e n t i s x=%1 . 4 f \
11
n V o l t a g e S t a n d i n g Wave R a t i o (VSWR) i s s=%1 . 2 f ’ ,x ,
s);
20
21 //===================END OF PROGRAM
=====================================
Scilab code Exa 3.6 point of attachment
2 //
==================================================================
3 clc ;
4 clear ;
5
6 Zl =100; // Pure Load r e s i s t a n c e o f a d i p o l e antenna in
ohms
7 Z0 =600; // C h a r a c t e r i s t i c Impedance o f a w i r e f e e d e r
i n ohms
8 f =100*10^6; // F r e q u e n c y i n Hz
10
11 //CALCULATIONS
13 l =(( w /(2*( %pi ) ) ) * atan ( sqrt ( Zl / Z0 ) ) ) ; // The p o s i t i o n
o f t h e Stub i n met
14 x = atand ( sqrt (( Zl * Z0 ) ) /( Zl - Z0 ) ) ;
15 y =180+ x ; // I n D e g r e e s
16 l1 =(( w /(2*( %pi ) ) ) * y ) ; // L e n g th o f S h o r t C i r c u i t e d
Stub i n met
17 l0 = l1 *(( %pi ) /180) ;
18
19 //OUTPUTS
20 mprintf ( ’ \ nThe P o i n t o f Attachment i s l=%1 . 3 f met \
nLength o f SC Stub i s l 1=%1 . 2 f met ’ ,l , l0 ) ;
12
21
22 //=========================END OF PROGRAM
==============================================================
Scilab code Exa 3.7 Terminating Impedance
2 //
==================================================================
3 clc ;
4 clear ;
5
7 S =3.2; // V o l t a g e S t a n d i n g Wave R a t i o (VSWR)
8
9 // I t i s p o s s i b l e t o m e a s u r e t h e l o a d i m p e d a n c e if
t h e l i n e i s assumed l o s s l e s s , by m e a s u r i n g t h e
VSWR, w a v e l e n g t h and t h e d i s t a n c e from t h e l o a d t o
t h e n e a r e s t v o l t a g e minimum
10 //CALCULATIONS
11 w =1; // Assume Wavelength i n met
12 Xmin =0.23* w ; // D i s t a n c e from t h e l o a d t o t h e n e a r e s t
v o l t a g e minimum i n met
13 b =(2*( %pi ) ) / w ;
14 Zl = Z0 *((1 -( S *( %i ) * tan ( b * Xmin ) ) ) /( S -(( %i ) * tan ( b * Xmin )
) ) ) ; // Load i m p e d a n c e i n ohms
15 disp ( ’ Load i m p e d a n c e i n ohms i s ’ ) ;
16 disp ( Zl ) ;
17
18
19 //=========================END OF PROGRAM
===================================================
13
20
21 // Note : Check t h e a n s w e r g i v e n i n Text book o n c e . I
t h i n k i t i s wrong i n t e x t book . .
Scilab code Exa 3.8 VSWR and Impedance
2 //
==================================================================
3 clc ;
4 clear ;
5
8 f =300*10^3; // F r e q u e n c y i n Hz
9 Pl =0.05; // Load Power i n w a t t s
11
12 //CALCULATIONS
14 p =(( Zl - Z0 ) /( Zl + Z0 ) ) ; // R e f l e c t i o n C o e f f i c i e n t
15 S =((1+ p ) /(1 - p ) ) ; // V o l t a g e S t a n d i n g Wave R a t i o (VSWR)
16
17 // S i n c e Zl >Z0 , f i r s t Vmax i s l o c a t e d a t t h e l o a d and
f i r s t Vmin i s l o c a t e d a t Wavelength /4
18 x1max =0; // P o s i t i o n o f f i r s t Vmax ( l o c a t e d a t t h e
l o a d ) from l o a d i n met
19 x1min = w /4; // P o s i t i o n o f f i r s t Vmin from l o a d i n met
20 Vmax = sqrt ( Pl * Zl ) ; // V a l u e o f maximum v o l t a g e i n v o l t s
21 Vmin = Vmax / S ; // V a l u e o f minimum v o l t a g e i n v o l t s
22 Zmax = Z0 * S ; // Impedance a t Vmax i n ohms
23 Zmin = Z0 / S ; // Impedance a t Vmin i n ohms
24
25 //OUTPUTS
14
26 mprintf ( ’ \ n V o l t a g e S t a n d i n g Wave R a t i o (VSWR) i s S=%1
. 0 f \ n P o s i t i o n o f f i r s t Vmax from l o a d i s x1max=
%d met ( l o c a t e d a t t h e l o a d ) \ n P o s i t i o n o f f i r s t
Vmin from l o a d i s x1min=%3 . 0 f met \ nValue o f
maximum v o l t a g e i s Vmax=%1 . 2 f v o l t s \ nValue o f
minimum v o l t a g e i s Vmin=%1 . 2 f v o l t s \ nImpedance
a t Vmax i s Zmax=%3 . 0 f ohms \ nImpedance a t Vmin i s
Zmin=%2 . 0 f ohms ’ ,S , x1max , x1min , Vmax , Vmin , Zmax ,
Zmin ) ;
27
28 //=========================END OF PROGRAM
==============================================================
Scilab code Exa 3.9 Reflection loss
2 //
==================================================================
3 clc ;
4 clear ;
5
7 Zs =50; // G e n e r a t o r i m p e d a n c e i n ohms
8 l =200; // L e n g th o f t r a n s m i s s i o n l i n e i n met
10
11 //CALCULATIONS
12 p =(( Zl - Z0 ) /( Zl + Z0 ) ) ; // R e f l e c t i o n C o e f f i c i e n t
13 x = abs ( p ) ;
14 Lr =10* log10 (1/(1 - x ^2) ) ; // R e f l e c t i o n l o s s i n dB
15 La =0; // S i n c e t h e l i n e i s l o s s l e s s , a t t e n u a t i o n l o s s
i s z e r o dB
16 Lt = La + Lr ; // T r a n s m i s s i o n l o s s i n dB
15
17 Lrt =10* log10 ( x ) ; // R e tu rn l o s s i n dB
18
19 //OUTPUT
20 mprintf ( ’ \ n R e f l e c t i o n l o s s i s Lr=%1 . 3 f dB \
n T r a n s m i s s i o n l o s s i s Lt=%1 . 3 f dB \ nReturn l o s s
i s L r t=%2 . 3 f dB ’ ,Lr , Lt , Lrt ) ;
21
22 //=========================END OF PROGRAM
==============================================================
Scilab code Exa 3.10 Characteristic Impedance
1 // c h a p t e r −3 p a g e 52 e x a m p l e 3 . 1 0
2 //
==================================================================
3 clc ;
4 clear ;
5
6 f =1000; // F r e q u e n c y i n Hz
7 l =10000; // L e n g th o f open w i r e t r a n s m i s s i o n l i n e in
met
8 z1 =2930; // Magnitude o f a s h o r t c i r c u i t i m p e d a n c e i n
ohms
9 p1 =26; // Phase o f a s h o r t c i r c u i t i m p e d a n c e i n deg
10 z2 =260; // Magnitude o f a open c i r c u i t i m p e d a n c e i n
ohms
11 p2 = -32; // Phase o f a open c i r c u i t i m p e d a n c e i n deg
12 //CALCULATIONS
13 Zsc =(( z1 * cosd ( p1 ) ) +(( %i ) *( z1 * sind ( p1 ) ) ) ) ;
14 Zoc =(( z2 * cosd ( p2 ) ) +(( %i ) *( z2 * sind ( p2 ) ) ) ) ;
15 Z0 = sqrt ( Zsc * Zoc ) ; // C h a r a c t e r i s t i c Impedance i n ohms
16 disp ( ’ C h a r a c t e r i s t i c Impedance i n ohms i s ’ ) ;
17 [ ro , theta ]= polar ( Z0 )
16
18 disp ( ro ) ;
19 disp ( theta *180/ %pi ) ;
20 g =((1/ l ) *( atanh ( sqrt ( Zsc / Zoc ) ) ) ) ; // P r o p a g a t i o n
Constant
21 disp ( g )
22 b = imag ( g ) ; // Phase C o n s t a n t
23 w =2* f *( %pi ) ; // A n g u l a r F r e q u e n c y i n r a d / s e c
24 Vp = w / b ; // Phase V e l o c i t y i n m/ s e c
25 disp ( Vp )
26 //OUTPUT
27 mprintf ( ’ \ nPhase V e l o c i t y i s Vp=%5 . 2 f m/ s e c ’ , Vp ) ;
28
29 //=========================END OF PROGRAM
==============================================================
30
31
32 // Note : Check t h e c a l c u l a t i o n o n c e
17
Chapter 4
Microwave Transmission Lines
2 //
==================================================================
3 clc ;
4 clear ;
5
6 d =0.0049; // D i a m e t e r o f i n n e r c o n d u c t o r i n met
7 D =0.0110; // I n n e r D i a m e t e r o f o u t e r c o n d u c t o r i n met
8 er =2.3; // P o l y e t h y l e n e d i e l e c t r i c
10
11 //CALCULATIONS
12 x = log ( D / d ) ;
13 L =(2*10^( -1) * x ) ; // I n d u c t a n c e p e r u n i t l e n g t h s i n
microH /m
14 C =(55.56*( er / x ) ) ; // The C a p a c i t a n c e p e r u n i t l e n g t h s
i n p i c o F /m
15 R0 =( x *(60/ sqrt ( er ) ) ) ; // The C h a r a c t e r i s t i c Impedance
i n ohms
16 V =( c / sqrt ( er ) ) /(10^3) ; // The V e l o c i t y o f p r o p a g a t i o n
18
i n Km/ s
17
18 //OUTPUT
19 mprintf ( ’ \ n I n d u c t a n c e p e r u n i t l e n g t h s i s L=%1 . 5 f
microH /m \ nThe C a p a c i t a n c e p e r u n i t l e n g t h s i s C=
%2 . 2 f p i c o F /m \ nThe C h a r a c t e r i s t i c Impedance i s
R0=%2 . 2 f ohms \ nThe V e l o c i t y o f p r o p a g a t i o n i s V=
%6 . 2 f Km/ s ’ ,L ,C , R0 , V ) ;
20
21 //=========================END OF PROGRAM
===================================================
Scilab code Exa 4.2 attenuation and phase constants
2 //
==================================================================
3 clc ;
4 clear ;
5
6 R =0.05; // R e s i s t a n c e i n ohm/m
7 L =0.16173*10^( -6) ; // I n d u c t a n c e p e r u n i t lengths in H
/m
8 C =0.15802*10^( -6) ; // The C a p a c i t a n c e p e r u n i t l e n g t h s
i n F/m
9 V =197814.14; // The V e l o c i t y o f p r o p a g a t i o n i n Km/ s
10 l =50; // L e n g th o f C o a x i a l L i n e i n met
11 Pin =480; // I n p u t Power t o t h e System i n w a t t s
12 f =3*10^9; // F r e q u e n c y i n Hz
13 c =3*10^5; // V e l o c i t y o f L i g h t i n Km/ s e c
14 e0 =8.854*10^( -12) ; // P e r m i t t i v i t y i n f r e e s p a c e i n F/
m
15
19
16 //CALCULATIONS
17 Z0 = sqrt ( L / C ) ;
18 A =( R /(2* Z0 ) ) ; // A t t e n u a t i o n C o n s t a n t i n NP/m
19 w =(2*( %pi ) * f ) ; // A n g u l a r F r e q u e n c y i n r a d / s e c
20 B =( w * sqrt ( L * C ) ) ; // Phase C o n s t a n t i n r a d /m
21 Vp =(1/ sqrt ( L * C ) ) /(10^3) ; // Phase V e l o c i t y i n Km/ s
22 er =((( c / V ) ^2) / e0 ) ; // R e l a t i v e P e r m i t t i v i t y
23 Pl =(2* Pin * l ) ; // Power L o s s i n w a t t s
24
25 //OUTPUT
26 mprintf ( ’ \ n A t t e n u a t i o n C o n s t a n t i s A=%1 . 4 f NP/m \
nPhase C o n s t a n t i s B=%4 . 3 f r a d /m \ nPhase V e l o c i t y
i s Vp=%4 . 3 f Km/ s \ n R e l a t i v e P e r m i t t i v i t y i s e r=
%12 . 2 f \ nPower L o s s i s Pl=%5 . 0 f w a t t s ’ ,A ,B , Vp , er ,
Pl ) ;
27
28 //=========================END OF PROGRAM
===========================================
Scilab code Exa 4.3 breakdown power
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For an a i r f i l l e d c o a x i a l c a b l e
7 f =9.375*10^9; // o p e r a t i n g f r e q u e n c y i n Hz
8 c =3*10^10; // V e l o c i t y o f L i g h t i n cm/ s e c
9 disp ( ’ Assuming a r a t i o o f ( b / a ) =2.3 and ( b+a ) <(w/ p i )
t o e x c l u d e h i g h e r o r d e r modes and a dominant
mode p r o p a g a t i n g ’ ) ;
10 a =0.36432; // l e n g t h o f c o a x i a l c a b l e i n cm
20
11 x =2.3; // r a t i o o f b / a
12
13 //CALCULATION
14 w0 =( c / f ) ; // f r e e s p a c e w a v e l e n g t h i n cm
15 Pbd =(3600*( a ^2) * log ( x ) ) ; // Breakdown power o f a
c o a x i a l c a b l e i n kW
16
17 //OUTPUT
18 mprintf ( ’ \ nBreakdown power o f a c o a x i a l c a b l e i s Pbd
=%3 . 0 f kW ’ , Pbd ) ;
19
20 //=========================END OF PROGRAM
===============================
2 //
==================================================================
3 clc ;
4 clear ;
5
6 b =0.3175; // D i s t a n c e b e t w e e n g r o u n d p l a n e s of strip
l i n e i n cm
7 d =0.0539; // D i a m e t e r o f c i r c u l a r c o n d u c t o r i n cm
8 er =2.32; // D i e l e c t r i c C o n s t a n t
10
11 //CALCULATION
12 Z0 =((60/ sqrt ( er ) ) * log ((4* b ) /( d *( %pi ) ) ) ) ; //
C h a r a c t e r i s t i c Impedance i n ohms
13 V =( c / sqrt ( er ) ) /(10^3) ; // The V e l o c i t y o f p r o p a g a t i o n
i n Km/ s
14
21
15 //OUTPUT
16 mprintf ( ’ \ n C h a r a c t e r i s t i c Impedance i s Z0=%2 . 2 f ohms
\ nThe V e l o c i t y o f p r o p a g a t i o n i s V=%5 . 2 f Km/ s ’ ,
Z0 , V ) ;
17
18 //=========================END OF PROGRAM
===================================================
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a m i c r o s t r i p t r a n s m i s s i o n l i n e
7 er =9.7; // r e l a t i v e d i e l e c t r i c c o n s t a n t o f an a l u m i n a
substrate
8 x1 =0.5; //w/ h r a t i o i n f i r s t t r a n s m i s s i o n l i n e
9 x2 =5; //w/ h r a t i o i n s e c o n d t r a n s m i s s i o n l i n e
11
12 //CALCULATION
13 disp ( ’ For c a s e 1 : w/ h =0.5 ’ ) ;
14 disp ( ’ S i n c e x1 =0.5 <1 , f o r t h i s we u s e h i g h i m p e d a n c e
a n a l y s i s ’ );
15 Eeff1 =((( er +1) /2) +(( er -1) /2) *(1/(( sqrt (1+(12/ x1 ) ) )
+(0.04*(1 - x1 ) ^2) ) ) ) ; // E f f e c t i v e d i e l e c t r i c
constant
16 Zo1 =((60/ sqrt ( Eeff1 ) ) * log ((8/ x1 ) +( x1 /4) ) ) ; //
C h a r a c t e r i s t i c i m p e d a n c e i n ohms
17 V1 =( c / sqrt ( Eeff1 ) ) /10^8; // V e l o c i t y o f p r o p a g a t i o n i n
22
1 0 ˆ 8 m/ s e c
18 mprintf ( ’ \ n E f f e c t i v e d i e l e c t r i c c o n s t a n t i s E e f f 1=%1
. 2 f \ n C h a r a c t e r i s t i c i m p e d a n c e i s Zo1=%2 . 2 f ohms
\ n V e l o c i t y o f p r o p a g a t i o n i s V1=%1 . 1 f ∗ 1 0 ˆ 8 m/ s e c
’ , Eeff1 , Zo1 , V1 ) ;
19
20 disp ( ’ For c a s e 2 : w/ h=5 ’ ) ;
21 disp ( ’ h e r e x2 >1 ’ ) ;
22 Eeff2 =((( er +1) /2) +(( er -1) /2) *(1/( sqrt (1+(12/ x2 ) ) ) ) ) ;
// E f f e c t i v e d i e l e c t r i c c o n s t a n t
23 Zo2 =((120*( %pi ) / sqrt ( Eeff2 ) ) *(1/( x2 +1.393+(0.667* log
(1.444+ x2 ) ) ) ) ) ; // C h a r a c t e r i s t i c i m p e d a n c e i n ohms
24 V2 =( c / sqrt ( Eeff2 ) ) /10^8; // V e l o c i t y o f p r o p a g a t i o n i n
1 0 ˆ 8 m/ s e c
25 mprintf ( ’ \ n E f f e c t i v e d i e l e c t r i c c o n s t a n t i s E e f f 2=%1
. 2 f \ n C h a r a c t e r i s t i c i m p e d a n c e i s Zo2=%2 . 2 f ohms
\ n V e l o c i t y o f p r o p a g a t i o n i s V2=%1 . 2 f ∗ 1 0 ˆ 8 m/ s e c
’ , Eeff2 , Zo2 , V2 ) ;
26
27 //=========================END OF PROGRAM
===============================
Scilab code Exa 4.6 ratio of areas
2 //
==================================================================
3 clc ;
4 clear ;
5
6 //To c a l c u l a t e the r a t i o o f c i r c u l a r waveguide cross
−s e c t i o n a l area to the r e c t a n g u l a r waveguide
cross section
7 disp ( ’ Assuming t h a t b o t h t h e s e w a v e g u i d e s have
23
s i m i l a r or equal c u t o f f f r e q u e n c i e s / wavelengths ’ )
;
8
9 disp ( ’ Case1 : When TE wave i s p r o p a g a t e d ’ ) ;
10 disp ( ’ For s t a n d a r d r e c t a n g u l a r w a v e g u i d e s a=2b and
For TE11 dominant mode i n c i r c u l a r w a v e g u i d e wc1
=(2( p i ) r ) / 1 . 8 4 1 ’ ) ;
11 disp ( ’ where r i s t h e r a d i u s o f t h e c i r c u l a r
w a v e g u i d e and wc1 i s t h e c u t o f f w a v e l e n g t h f o r
c i r c u l a r waveguide ’ );
12 disp ( ’ I t i s g i v e n wc1=wc2 where wc2 i s t h e c u t o f f
wavelength f o r r e c t a n g u l a r waveguide ’ );
13 disp ( ’ For TE10 ( dominant mode ) o f p r o p a g a t i o n i n
r e c t a n g u l a r w a v e g u i d e wc2=2a ’ ) ;
14 disp ( ’ S i n c e wc2 =(2 ab ) / ( s q r t ( ( mb) ˆ2+( nb ) ˆ 2 ) ) a s m=1; n
=0 f o r TE10 wc2=2ab / b=2a ’ ) ;
15 disp ( ’ By e q u a t i n g wc1=wc2 , we g e t a = 1 . 7 0 6 4 5 r ’ ) ;
16 disp ( ’ For a s t a n d a r d w a v e g u i d e a=2b t h e r e f o r e , b=a /2
’ );
17 disp ( ’Now t h e a r e a o f r e c t a n g u l a r w a v e g u i d e=a ∗b=a ∗ a
/ 2 = 1 . 7 0 6 4 5 r ∗ 1 . 7 0 6 4 5 r / 2 = 1 . 4 5 6 r ˆ2 ’ ) ;
18 disp ( ’ Area o f r e c t a n g u l a r w a v e g u i d e = 1 . 4 5 6 r ˆ2 , Area
o f c i r c u l a r w a v e g u i d e =( p i ) ∗ r ˆ2 ’ ) ;
19 disp ( ’ R a t i o o f a r e a o f c i r c u l a r t o a r e a o f
r e c t a n g u l a r w a v e g u i d e =( Area o f c i r c u l a r w a v e g u i d e
/ Area o f r e c t a n g u l a r w a v e g u i d e ) =( p i ∗ r ˆ 2 ) / ( 1 . 4 5 6 r
ˆ2) =2.1576873=2.2 ’ );
20 disp ( ’ T h i s c l e a r l y shows t h a t t h e s p a c e o c c u p i e d by
a r e c t a n g u l a r w a v e g u i d e s y s t e m i s l e s s compared
t o t h a t f o r a c i r c u l a r w a v e g u i d e s y s t e m . Hence
c i r c u l a r w a v e g u i d e s a r e n o t p r e f e r r e d i n some
a p p l i c a t i o n s ’ );
21
22 disp ( ’ Case2 : When TM wave i s p r o p a g a t e d ’ ) ;
23 disp ( ’ For TM01 mode wc1 =(2∗ p i ∗ r ) / (Pnm) min =(2∗ p i ∗ r ) /
Pnm=(2∗ p i ∗ r ) / 2 . 4 0 5 where r i s t h e r a d i u s o f
c i r c u l a r w a v e g u i d e wc1 = 2 . 6 1 5 5 r ’ ) ;
24 disp ( ’Now i f wc2 i s t h e w a v e l e n g t h f o r TM11 wave
24
propagating in a standard r e c t a n g u l a r waveguide
wc2=wc1 but wc2 =(2 ab ) / s q r t ( aˆ2+b ˆ 2 ) ’ ) ;
25 disp ( ’ For s t a n d a r d w a v e g u i d e s , we know a=2b , wc2 =(2∗2
b∗b ) / s q r t ( 4 bˆ2+b ˆ 2 ) =(4b ˆ 2 ) / s q r t ( 5 b ˆ 2 ) =4b / s q r t ( 5 ) ’
);
26 disp ( ’ By e q u a t i n g wc1=wc2 , we g e t 2 . 6 1 5 5 r =4b / s q r t ( 5 )
=>b = 1 . 4 6 2 1 r ’ ) ;
27 disp ( ’ Area o f r e c t a n g u l a r w a v e g u i d e=b∗b=b ˆ2 but b
= 1 . 4 6 2 1 r , s o Area o f r e c t a n g u l a r w a v e g u i d e
= ( 1 . 4 6 2 1 r ) ˆ 2 = 2 . 1 3 2 r ˆ2 and Area o f c i r c u l a r
w a v e g u i d e= p i ∗ r ˆ2 ’ ) ;
28 disp ( ’ R a t i o o f a r e a o f c i r c u l a r t o a r e a o f
r e c t a n g u l a r w a v e g u i d e =( Area o f c i r c u l a r w a v e g u i d e
/ Area o f r e c t a n g u l a r w a v e g u i d e ) =( p i ∗ r ˆ 2 ) / ( 2 . 1 3 2 r
ˆ 2 ) =1.5 ’ ) ;
29
30 //=========================END OF PROGRAM
===============================
Scilab code Exa 4.7 breadth of the guide
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a r e c t a n g u l a r w a v e g u i d e
7 disp ( ’ For a r e c t a n g u l a r w a v e g u i d e t h e dominant mode
i s t h e TE10 mode . TE10 mode can p r o p a g a t e a t a
l o we r f r e q u e n c y ’ );
8 f =9*10^9; // f r e q u e n c y i n Hz
9 wg =4; // g u i d e w a v e l e n g t h i n cm
25
11 disp ( ’ For TE10 mode wc=2a ’ ) ;
12
13 //CALCULATION
15 wc =( w0 / sqrt (1 -( w0 / wg ) ^2) ) ; // C u t o f f w a v e l e n g t h f o r
TE10 mode i n cm
16 disp ( ’ F r e e s p a c e w a v e l e n g t h w0 i n cm i s ’ ) ;
17 disp ( w0 ) ;
18 disp ( ’ C u t o f f w a v e l e n g t h wc i n cm i s ’ ) ;
19 disp ( wc ) ;
20 disp ( ’ S i n c e wc>w0 , t h e wave p r o p a g a t e s ’ ) ;
21 a =( wc /2) ; // l e n g t h o f t h e g u i d e i n cm
22 b =( wc /4) ; // b r e a d t h o f t h e g u i d e i n cm
23
24 //OUTPUT
25 mprintf ( ’ \ n l e n g t h o f t h e g u i d e i s a=%1 . 0 f cm \
n b r e a d t h o f t h e g u i d e i s b=%1 . 1 f cm ’ ,a , b ) ;
26
27 //=========================END OF PROGRAM
===============================
Scilab code Exa 4.8 guide wavelength
2 //
==================================================================
3 clc ;
4 clear ;
5
6 a =10; // b r e a d t h o f a r e c t a n g u l a r w a v e g u i d e i n cm
7 f =2.5*10^9; // F r e q u e n c y i n Hz i n TE10 mode
9
10 //CALCULATION
26
11 wc =2* a ; // C u t o f f w a v e l e n g t h f o r TE10 mode i n cm
12 w0 =( c / f ) ; // F r e e s p a c e w a v e l e n g t h i n cm
13 x = sqrt (1 -( w0 / wc ) ^2) ;
14 wg =( w0 / x ) ; // Guide w a v e l e n g t h i n cm
15 Vp =( c / x ) /10^5; // Phase V e l o c i t y i n Km/ s e c
16 Vg =(( c ^2) / Vp ) /10^10; // Group V e l o c i t y i n Km/ s e c
17
18 //OUTPUT
19 mprintf ( ’ \ n C u t o f f w a v e l e n g t h f o r TE10 mode i s wc=%2
. 0 f cm \ nGuide w a v e l e n g t h i s wg=%2 . 0 f cm \ nPhase
V e l o c i t y i s Vp=%7 . 2 f Km/ s e c \ nGroup V e l o c i t y i s
Vg=%6 . 2 f Km/ s e c ’ ,wc , wg , Vp , Vg ) ;
20
21 //=========================END OF PROGRAM
===============================
2 //
==================================================================
3 clc ;
4 clear ;
5
6 f =8.6*10^9; // f r e q u e n c y i n Hz
8 a =2.5; // L e n g th o f a Waveguide i n cm
9 b =1; // Width o f a Waveguide i n cm
10
11 //CALCULATION
12 disp ( ’ The c o n d i t i o n f o r t h e wave t o p r o p a g a t e a l o n g
a g u i d e i s t h a t wc>w0 . ’ ) ;
13 w0 = c / f ; // f r e e s p a c e w a v e l e n g t h i n cm
27
15 disp ( w0 ) ;
16 disp ( ’ For TE waves , wc=(2 ab / s q r t ( ( mb) ˆ2+( na ) ˆ 2 ) ) ’ ) ;
17 disp ( ’ For TE01 waves ’ ) ;
18 m1 =0;
19 n1 =1;
20 wc1 =((2* a * b ) /( sqrt (( m1 * b ) ^2+( n1 * a ) ^2) ) ) ; // C u t o f f
w a v e l e n g t h f o r TE01 mode i n cm
21 disp ( ’ C u t o f f w a v e l e n g t h f o r TE01 mode i n cm i s ’ ) ;
22 disp ( wc1 ) ;
23 disp ( ’ S i n c e wc f o r TE01=2cm i s n o t g r e a t e r t h a n w0
TE01 , w i l l n o t p r o p a g a t e f o r TE01 mode . ’ ) ;
25 m2 =1;
26 n2 =0;
29 disp ( wc2 ) ;
30 disp ( ’ S i n c e wc TE10 > w0 TE10 i s a p o s s i b l e mode . ’ ) ;
31 fc =( c / wc2 ) /10^9; // C u t o f f f r e q u e n c y i n GHz
32 disp ( ’ For TE11 and TM11 waves ’ ) ;
33 m3 =1;
34 n3 =1;
36 disp ( ’ C u t o f f w a v e l e n g t h f o r TE11 and TM11 modes i n
cm i s ’ ) ;
37 disp ( wc3 ) ;
38 disp ( ’ As wc f o r TE11 and TM11 i s < w0 b o t h TE11 and
TM11 do n o t p r o p a g a t e a s h i g h e r modes . ’ ) ;
39 wg =( w0 / sqrt (1 -( w0 / wc2 ) ^2) ) ; // Guide w a v e l e n g t h i n cm
40 disp ( ’ From t h e a b o v e a n a l y s i s we c o n c l u d e t h a t o n l y
TE10 mode i s p o s s i b l e ’ ) ;
41
42 //OUTPUT
43 mprintf ( ’ \ n C u t o f f f r e q u e n c y i s f c=%1 . 0 f GHz \ nGuide
w a v e l e n g t h i s wg=%1 . 3 f cm ’ ,fc , wg ) ;
44
28
45 //=========================END OF PROGRAM
===============================
Scilab code Exa 4.10 area
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For an a i r f i l l e d c i r c u l a r Waveguide i n t h e
dominant mode
8 disp ( ’ For an a i r f i l l e d c i r c u l a r Waveguide TE11 i s
t h e dominant mode i e p r o p a g a t e d ’ ) ;
9 wc =10; // c u t o f f wave l e n g t h i n cm
10
11 //CALCULATION
12 r =((1.841* wc ) /(2*( %pi ) ) ) ; // r a d i u s of circular
Waveguide i n cm
13 A =( %pi ) * r ^2; // C r o s s s e c t i o n a l a r e a o f t h e g u i d e i n
s q . cms
14 fc =( c / wc ) /10^9; // C u t o f f f r e q u e n c y f o r TE11 mode i n
GHz
15 disp ( ’ C u t o f f f r e q u e n c y f o r TE11 mode i n GHz i s ’ ) ;
16 disp ( fc ) ;
17 disp ( ’ F r e q u n c y a b o v e 3GHz can be p r o p a g a t e d t h r o u g h
the waveguide ’ );
18
19 //OUTPUT
20 mprintf ( ’ \ n C r o s s s e c t i o n a l a r e a o f t h e g u i d e i s A=%2
. 2 f s q . cms ’ ,A ) ;
21
29
22 //=========================END OF PROGRAM
===============================
Scilab code Exa 4.11 possible modes
2 //
==================================================================
3 clc ;
4 clear ;
5
7 f =5*10^9; // f r e q u e n c y i n Hz
9 a =4; // Le n g th o f R e c t a n g u l a r Waveguide i n cm
10 b =3; // Width o f R e c t a n g u l a r Waveguide i n cm
11
12 //CALCULATION
13 disp ( ’ The c o n d i t i o n f o r t h e wave t o p r o p a g a t e a l o n g
a g u i d e i s t h a t wc>w0 . ’ ) ;
14 w0 = c / f ; // f r e e s p a c e w a v e l e n g t h i n cm
16 disp ( w0 ) ;
17 disp ( ’ For TE waves , wc=(2 ab / s q r t ( ( mb) ˆ2+( na ) ˆ 2 ) ) ’ ) ;
19 m1 =0;
20 n1 =1;
23 disp ( wc1 ) ;
24 disp ( ’ S i n c e wc f o r TE01=6cm i s n o t g r e a t e r t h a n w0
TE01 , w i l l n o t p r o p a g a t e f o r TE01 mode . ’ ) ;
30
26 m2 =1;
27 n2 =0;
30 disp ( wc2 ) ;
31 disp ( ’ S i n c e wc TE10 > w0 TE10 i s a p o s s i b l e mode . ’ ) ;
33 m3 =1;
34 n3 =1;
37 disp ( wc3 ) ;
38 disp ( ’ As wc TE11 < w0 TE11 d o e s n o t p r o p a g a t e . ’ ) ;
39
40 //=========================END OF PROGRAM
===============================
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For an a i r f i l l e d c i r c u l a r Waveguide i n t h e
dominant mode
7 D =4; // I n n e r d i a m e t e r o f an a i r f i l l e d c i r c u l a r
Waveguide i n cm
9
10 //CALCULATION
31
11 disp ( ’ The dominant mode i n t h e c i r c u l a r w a v e g u i d e
would be l i k e TE11 , wc i s maximum ’ ) ;
12 r = D /2; // r a d i u s i n cm
13 wc =((2*( %pi ) * r ) /1.841) ; // C u t o f f w a v e l e n g t h i n cms
14 fc =( c / wc ) /10^9; // C u t o f f f r e q u e n c y i n GHz
15 mprintf ( ’ \ n C u t o f f w a v e l e n g t h i s wc=%1 . 4 f cms \
n C u t o f f f r e q u e n c y i s f c=%1 . 3 f GHz ’ ,wc , fc ) ;
16 disp ( ’ S i n c e cut − o f f f r e q u e n c y i s 4 . 3 9 5 GHz ,
f r e q u e n c i e s h i g h e r t h a n f c w i l l be p r o p a g a t e d .
Assume a s i g n a l o f f r e q u e n c y o f 5 GHz i s b e i n g
propagated ’ );
17 f =5*10^9; // f r e q u e n c y o f s i g n a l i n Hz
19 wg =( w0 / sqrt (1 -( w0 / wc ) ^2) ) ; // Guide w a v e l e n g t h i n cm
20 mprintf ( ’ \nWave l e n g t h i n t h e g u i d e i s wg=%2 . 2 f cm ’ ,
wg ) ;
21
22 //=========================END OF PROGRAM
===============================
Scilab code Exa 4.13 frequency
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a r e c t a n g u l a r w a v e g u i d e i n TE10 mode
10 x =4.55; // d i s t a n c e b e t w e e n maximum and minimum i n cm
11
32
12 //CALCULATIONS
13 wc =2* a ; // C u t o f f w a v e l e n g t h f o r a TE10 mode i n cms
14 wg =4* x ; // Guide Wavelength i n cm
15 w0 =( wg / sqrt (1+( wg / wc ) ^2) ) ; // // F r e e s p a c e w a v e l e n g t h
i n cm
16 f =( c / w0 ) /10^9; // F r e q u e n c y o f t h e wave i n GHz
17
18 //OUTPUT
19 mprintf ( ’ \ n F r e q u e n c y o f t h e wave i s f=%1 . 3 f GHz ’ ,f ) ;
20
21 //=========================END OF PROGRAM
===============================
2 //
==================================================================
3 clc ;
4 clear ;
5
7 b =2.5; // L e n g th o f R e c t a n g u l a r Waveguide i n cm
8 a =5; // b r e a d t h o f R e c t a n g u l a r Waveguide i n cm
10 w0 =4.5; // F r e e s p a c e w a v e l e n g t h i n cm
11
12 //CALCULATION
13 disp ( ’ For a TE10 mode which i s t h e dominant mode ’ ) ;
14 wc =2* a ; // C u t o f f w a v e l e n g t h i n cm
16 Vp =( c / sqrt (1 -( w0 / wc ) ^2) ) /10^10; // Phase V e l o c i t y i n
1 0 ˆ 1 0 cm/ s e c
17 B =((2*( %pi ) * sqrt ( wc ^2 - w0 ^2) ) /( w0 * wc ) ) ; // Phase
33
constant in radians
18
19 //OUTPUT
20 mprintf ( ’ \ nGuide w a v e l e n g t h i s wg=%1 . 5 f cm \ nPhase
c o n s t a n t i s B=%1 . 3 f r a d i a n s \ nPhase V e l o c i t y i s
Vp=%1 . 2 f ∗ 1 0 ˆ 1 0 cm/ s e c ’ ,wg ,B , Vp ) ;
21
22 //=========================END OF PROGRAM
===========================================
23
24 // Note : Check t h e a n s w e r s o n c e
25 // C o r r e c t a n s w e r s a r e
26 // Guide w a v e l e n g t h i s wg = 5 . 0 3 9 0 3 cm
27 // Phase c o n s t a n t i s B= 1 . 2 4 7 r a d i a n s
28 // Phase V e l o c i t y i s Vp=3.36 ∗ 1 0 ˆ 1 0 cm/ s e c
Scilab code Exa 4.15 possible modes
2 //
==================================================================
3 clc ;
4 clear ;
5
6 wcTE10 =16; // C r i t i c a l w a v e l e n g t h o f TE10 mode i n cm
7 wcTM11 =7.16; // C r i t i c a l w a v e l e n g t h o f TM11 mode i n cm
8 wcTM21 =5.6; // C r i t i c a l w a v e l e n g t h o f TM21 mode i n cm
9 disp ( ’ For any wave t o be p r o p a g a t e d , t h e c o n d i t i o n
t o be met i s wc>wo ’ ) ;
10 wo1 =10; // F r e e s p a c e w a v e l e n g t h i n cm
11 wo2 =5; // F r e e s p a c e w a v e l e n g t h i n cm
12 disp ( ’ C r i t i c a l w a v e l e n g t h o f TE10 mode i n cm i s ’ ) ;
13 disp ( wcTE10 ) ;
14 disp ( ’ C r i t i c a l w a v e l e n g t h o f TM11 mode i n cm i s ’ ) ;
34
15 disp ( wcTM11 ) ;
16 disp ( ’ C r i t i c a l w a v e l e n g t h o f TM21 mode i n cm i s ’ ) ;
17 disp ( wcTM21 ) ;
18 disp ( ’ For wo1=10cm , The mode t h a t p r o p a g a t e s o n l y
TE10 . B e c a u s e wcTE10>wo1 and a l l o t h e r modes t h a t
i s TM11 TM21 d o n o t p r o p a g a t e ’ ) ;
19 disp ( ’ For wo2=5cm ’ ) ;
20 disp ( ’ wcTE10>wo2 , s o TE10 mode p r o p a g a t e s ’ ) ;
21 disp ( ’ wcTM11>wo2 , s o TE11 mode p r o p a g a t e s ’ ) ;
22 disp ( ’ wcTE21>wo2 , s o TE21 mode p r o p a g a t e s ’ ) ;
2 //
==================================================================
3 clc ;
4 clear ;
5
6 n =120*( %pi ) ; // I n t r i n s i c Impedance
9 f =10^10; // F r e q u e n c y i n Hz
11
12 //CALCULATION
13 wc =((2* a * b ) / sqrt ( a ^2+ b ^2) ) ; // C u t o f f w a v e l e n g t h i n
TM11 mode i n cms
14 w0 =( c / f ) ; // F r e e s p a c e w a v e l e n g t h i n cms
15 ZTM =( n * sqrt (1 -( w0 / wc ) ^2) ) ; // C h a r a c t e r i s t i c Wave
Impedance i n ohms
16
17 //OUTPUT
18 mprintf ( ’ \ n C h a r a c t e r i s t i c Wave Impedance i s ZTM=%2 . 3
35
f ohms ’ , ZTM ) ;
19
20
21 //=========================END OF PROGRAM
=================================
22
23 // Note : Check t h e g i v e n a n s w e r o n c e i t i s wrong
24 // c u r r e c t a n s w e r i s 1 6 3 . 2 4 2 ohms
2 //
==================================================================
3 clc ;
4 clear ;
5
7 f =6*10^9; // F r e q u e n c y i n Hz
8
9 //CALCULATION
10 fc =(0.8* f ) ; // Given C u t o f f f r e q u e n c y f o r TE11 mode i n
Hz
11 wc =( c / fc ) ; // C u t o f f w a v e l e n g t h i n cms
12 D =((1.841* wc ) /( %pi ) ) ; // D i a m e t e r o f w a v e g u i d e i n cm
13 w0 =( c / f ) ; // F r e e s p a c e w a v e l e n g t h i n cm
15
16 //OUTPUT
17 mprintf ( ’ \ n D i a m e t e r o f t h e w a v e g u i d e i s D=%1 . 4 f cm \
nGuide w a v e l e n g t h i s wg=%1 . 3 f cm ’ ,D , wg ) ;
18
19 //=========================END OF PROGRAM
===============================
36
Scilab code Exa 4.18 proof
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a TE10 mode
7 a =1.5; // L e n g th o f an a i r f i l l e d s q u a r e Waveguide i n
m
8 b =1; // b r e a d t h o f an a i r f i l l e d s q u a r e Waveguide i n
cm
10 f =6*10^9; // I m p r e s s e d F r e q u e n c y i n Hz
11 er =4; // d i e l e c t r i c c o n s t a n t
12
13 //CALCULATION
14 wc =2* a ; // C u t o f f w a v e l e n g t h i n cm
15 fc =( c / wc ) /10^9; // C u t o f f f r e q u e n c y i n GHz
16 disp ( ’ C u t o f f f r e q u e n c y i n GHz i s ’ ) ;
17 disp ( fc ) ;
18 disp ( ’ The i m p r e s s e d f r e q u e n c y o f 6 GHz i s l e s s than
t h e C u t o f f f r e q u e n c y and h e n c e t h e s i g n a l w i l l
not pass through the gu i d e ’ );
19 w =( c / f ) ; // Wavelength i n cm
20 disp ( ’ A l t e r n a t i v e l y , t h e w a v e l e n g t h o f t h e i m p r e s s e d
s i g n a l i n cm i s ’ ) ;
21 disp ( w ) ;
22 wair = w ;
23 disp ( ’ which i s l o n g e r t h a n t h e c u t o f f w a v e l e n g t h ( 3
cm ) and h e n c e no p r o p a g a t i o n o f t h e wave ’ ) ;
24 w1 = wair / sqrt ( er ) ; // Wavelength i n cm
37
25 disp ( ’ I f t h e w a v e g u i d e i s l o a d e d w i t h d i e l e c t r i c o f
e r =4 , t h e n t h e w a v e l e n g t h i n cm i s ’ ) ;
26 disp ( w1 ) ;
27 disp ( ’ which i s l e s s t h a n w a i r ’ ) ;
28 disp ( ’Now t h e s i g n a l w i t h 6 GHz f r e q u e n c y w i l l p a s s
through the d i e l e c t r i c loaded waveguide ’ );
29
30 //=========================END OF PROGRAM
===============================
Scilab code Exa 4.19 amount of attenuation
2 //
==================================================================
3 clc ;
4 clear ;
5
6 a =0.015; // L e n g th o f h o l l o w R e c t a n g u l a r Waveguide i n
m
7 b =1; // b r e a d t h o f h o l l o w R e c t a n g u l a r Waveguide i n cm
8 f =6*10^9; // F r e q u e n c y i n Hz i n TE10 mode
10 m =1; // V a l u e o f m i n TE10 mode
11 n =0; // V a l u e o f n i n TE10 mode
12 u =4*( %pi ) *10^( -7) ; // P e r m e a b i l i t y i n f r e e s p a c e i n
Henry
13 e =8.854*10^( -12) ; // P e r m i t t i v i t y i n f r e e s p a c e i n F/m
14
15 //CALCULATION
16 wc =2* a ; // C u t o f f w a v e l e n g t h f o r TE10 mode i n m
17 fc = c / wc ; // C u t o f f f r e q u e n c y i n Hz
18 w =2*( %pi ) * f ; // A n g u l a r f r e q u e n c y i n r a d / s e c
19
38
20 // So 6GHz s i g n a l w i l l n o t p a s s t h r o u g h w a v e g u i d e but
w i l l get attenuated
21 A =( sqrt (( m *( %pi ) / a ) ^2+( n *( %pi ) / b ) ^2 -( w ^2* u * e ) ) ) ; //
A t t e n u a t i o n i n NP/m
22 AdB = A *(20/ log (10) ) ; // A t t e n u a t i o n i n dB/m
23
24 //OUTPUT
25 mprintf ( ’ \Amount o f A t t e n u a t i o n i s A=%3 . 1 f NP/m \
n A t t e n u a t i o n i s AdB=%4 . 2 f dB/m ’ ,A , AdB ) ;
26
27 //=========================END OF PROGRAM
===============================
Scilab code Exa 4.20 max power handling capacity
2 //
==================================================================
3 clc ;
4 clear ;
5
10 Emax =3000; //Max p o t e n t i a l g r a d i e n t i n V/cm
11
12 //CALCULATION
14 disp ( ’ F r e e s p a c e Wavelength i n cm i s ’ ) ;
15 disp ( w0 ) ;
16 wc =2* a ; // C u t o f f w a v e l e n g t h i n TE10 mode i n cms
17 wg =( w0 / sqrt (1 -( w0 / wc ) ^2) ) ; // Guide w a v e l e n g t h i n cms
18 disp ( ’ Guide Wavelength i n cm i s ’ ) ;
39
19 disp ( wg ) ;
20 P =((6.63*10^( -4) ) *( Emax ^2) * a * b *( w0 / wg ) ) /1000; // Power
h a n d l i n g c a p a b i l i t y o f t h e w a v e g u i d e i n kW
21
22 //OUTPUT
23 mprintf ( ’ \ nPower h a n d l i n g c a p a b i l i t y o f t h e
w a v e g u i d e i s P=%2 . 3 f kW ’ ,P ) ;
24
25
26 //=========================END OF PROGRAM
=================================
Scilab code Exa 4.21 maximum power
2 //
==================================================================
3 clc ;
4 clear ;
5
6 d =5; // I n t e r n a l D i a m e t e r o f c i r c u l a r w a v e g u i d e i n cm
9 Emax =300; //Max f i e l d s t r e n g t h i n V/cm
10
11 //CALCULATION
13 wc =(( d *( %pi ) ) /1.841) ; // C u t o f f w a v e l e n g t h i n TE11
mode i n cms
15 Pmax =(0.498*( Emax ^2) *( d ^2) *( w0 / wg ) ) /1000; //Maximum
power i n kWatts
16
17 //OUTPUT
40
18 mprintf ( ’ \nMaximum power i s Pmax=%4 . 2 f kWatts ’ , Pmax )
;
19
20
21 //=========================END OF PROGRAM
===============================
Scilab code Exa 4.22 peak value of electric field
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For an a i r f i l l e d s q u a r e w a v e g u i d e
7 a =0.01; // L e n g th o f an a i r f i l l e d s q u a r e Waveguide i n
m
8 b =0.01; // b r e a d t h o f an a i r f i l l e d s q u a r e Waveguide
in m
11 Pmax =746; //Max power =1 h o r s e p o w e r in W
12 n =120*( %pi ) ; // Impedance o f f r e e s p a c e i n ohms
13
14 //CALCULATION
15 w0 =( c / f ) ; // F r e e s p a c e w a v e l e n g t h i n m
16 wc =2* a ; // C u t o f f w a v e l e n g t h i n m
17 ZTE =( n / sqrt (1 -( w0 / wc ) ^2) ) ; // Impedance i n ohms
18 Emax =( sqrt (( Pmax *4* ZTE ) /( a * b ) ) ) /1000; // The Peak
value of E l e c t r i c f i e l d occuring in the guide in
kV/m
19 // From P= ( 1 / 2 ) ∗ I n t e g r a t i o n ( Re (E∗H) ) da
20 // and Pmax= ( 1 / ( 4 ∗ZTE) ) ∗Emaxˆ2∗ a ∗b
41
21
22 //OUTPUT
23 mprintf ( ’ \ nThe Peak v a l u e o f E l e c t r i c f i e l d o c c u r i n g
i n t h e g u i d e i s Emax=%3 . 2 f kV/m ’ , Emax ) ;
24
25 //=========================END OF PROGRAM
===============================
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For an a i r f i l l e d r e c t a n g u l a r w a v e g u i d e
7 a =0.023; // L e n g th o f an a i r f i l l e d Rectangular
Waveguide i n m
8 b =0.01; // b r e a d t h o f an a i r f i l l e d R e c t a n g u l a r
Waveguide i n m
10 f =9.375*10^9; // F r e q u e n c y i n Hz i n TE11 mode
11 w0 =0.01; // F r e e s p a c e w a v e l e n g t h i n m
12 wc =0.02; // C u t o f f w a v e l e n g t h i n m
13 Pmax =746; //Max power =1 h o r s e p o w e r in W
14
15 //CALCULATION
16 wo =( c / f ) ; // F r e e s p a c e w a v e l e n g t h i n cm
17 Pbd =(597* a * b * sqrt (1 -( wo /(2* a ) ) ^2) ) ; // The breakdown
power f o r t h e dominant mode i e TE11 i n W
18 wg =( w0 / sqrt (1 -( w0 / wc ) ^2) ) ; // Guide w a v e l e n g t h i n m
19 Emax =( sqrt (( Pmax * wg ) /(6.63*10^( -4) * w0 ) ) ) /1000; //Max
e l e c t r i c f i e l d i n kV/m
42
20
21 //OUTPUT
22 mprintf ( ’ \ nThe breakdown power f o r t h e dominant mode
i e TE11 i s Pbd=%1 . 5 f W \nMax e l e c t r i c f i e l d i s
Emax=%1 . 4 f kV/m ’ ,Pbd , Emax ) ;
23
24 //=========================END OF PROGRAM
===========================================
25
26
29 // The breakdown power f o r t h e dominant mode i e TE11
i s Pbd = 0 . 0 9 8 6 4 W
30 //Max e l e c t r i c f i e l d i s Emax = 1 . 1 3 9 8 kV/m
2 //
==================================================================
3 clc ;
4 clear ;
5
6 a =2.5; // R a d i u s o f c i r c u l a r w a v e g u i d e i n cm
7 d =5; // I n t e r n a l D i a m e t e r o f c i r c u l a r w a v e g u i d e i n cm
10
11 //CALCULATION
13 wc =(( d *( %pi ) ) /1.841) ; // C u t o f f w a v e l e n g t h i n TE11
mode i n cms
14 fc =( c / wc ) ; // C u t o f f f r e q u e n c y i n Hz
43
15 Pbd =(1790*( a ^2) * sqrt (1 -( fc / f ) ^2) ) /1000; // Breakdown
Power i n TE11 mode i n kW
16
17 //OUTPUT
18 mprintf ( ’ \ nBreakdown Power i n TE11 mode i s Pbd=%5 . 3 f
kW ’ , Pbd ) ;
19
20 //=========================END OF PROGRAM
===============================
44
Chapter 5
Cavity Resonators
Scilab code Exa 5.1 minimum distance
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a c i r c u l a r w a v e g u i d e
7 a =3; // r a d i u s i n cm
8 f0 =10*10^9; // r e s o n a n t f r e q u e n c y o f a c i r c u l a r
r e s o n a t o r i n Hz
9 disp ( ’ Given t h e mode o f o p e r a t o r i s TM011 s o h e r e n
=0 ,m=1 , p=1 ’ ) ;
10 c =3*10^10; // V e l o c i t y o f l i g h t i n cm/ s e c
11 m =1;
12 n =0;
13 p =1;
14 Pnm =2.405; // dominant mode v a l u e [ TM01 ]
15
16 //CALCULATION
17 d =(( p *( %pi ) ) /( sqrt ((2*( %pi ) * f0 / c ) ^2 -( Pnm / a ) ^2) ) ) ; //
45
The minimum d i s t a n c e b e t w e e n t h e two end p l a t e s
i n cms
18
19 //OUTPUT
20 mprintf ( ’ \ nThe minimum d i s t a n c e b e t w e e n t h e two end
p l a t e s o f a c i r c u l a r w a v e g u i d e i s d=%1 . 2 f cms ’ ,d )
;
21
22 //=========================END OF PROGRAM
===============================
Scilab code Exa 5.2 lowest resonant frequency
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a r e c t a n g u l a r c a v i t y r e s o n a t o r
7 a =2; // b r e a d t h i n cm
8 b =1; // h e i g h t i n cm
9 l =3; // l e n g t h o f r e c t a n g u l a r w a v e g u i d e i n cm
10 disp ( ’ L owest r e s o n a n t f r e q u e n c y i s o b t a i n e d f o r t h e
dominant mode TE10 [ f=c /w where w i n c r e a s e s a s f
d e c r e a s e s . I n dominant mode wc i s maximum ] ’ ) ;
11 disp ( ’ So t h e dominant mode i s TE101 s o h e r e m=1 , n =0 ,
p=1 ’ ) ;
13 m =1;
14 n =0;
15 p =1;
16
17 //CALCULATION
46
18 f0 =(( c /2) * sqrt (( m / a ) ^2+( n / b ) ^2+( p / l ) ^2) ) /10^9; // The
resonant frequency of a rectangular cavity
r e s o n a t o r i n GHz
19
20 //OUTPUT
21 mprintf ( ’ \ nThe r e s o n a n t f r e q u e n c y o f a r e c t a n g u l a r
c a v i t y r e s o n a t o r i s f 0=%1 . 0 f GHz ’ , f0 ) ;
22
23 //=========================END OF PROGRAM
===============================
Scilab code Exa 5.3 resonant frequency
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a c i r c u l a r r e s o n a t o r
7 D =12.5; // d i a m e t e r i n cm
8 l =5; // l e n g t h o f c i r c u l a r w a v e g u i d e i n cm
9 disp ( ’ Given t h e mode o f o p e r a t o r i s TM012 s o h e r e n
=0 ,m=1 , p=2 ’ ) ;
11 m =1;
12 n =0;
13 p =2;
14 Pnm =2.405; // dominant mode v a l u e [ TM01 ]
15
16 //CALCULATION
17 a = D /2; // r a d i u s i n cm
18 f0 =(( c /(2*( %pi ) ) ) * sqrt (( Pnm / a ) ^2+(( p *( %pi ) ) / l ) ^2) )
/10^9; // The r e s o n a n t f r e q u e n c y o f a c i r c u l a r
47
r e s o n a t o r i n GHz
19
20 //OUTPUT
21 mprintf ( ’ \ nThe r e s o n a n t f r e q u e n c y o f a c i r c u l a r
r e s o n a t o r i s f 0=%1 . 2 f GHz ’ , f0 ) ;
22
23 //=========================END OF PROGRAM
===============================
Scilab code Exa 5.4 resonant frequency
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a c i r c u l a r r e s o n a t o r
7 a =3; // r a d i u s i n cm
8 b =2; // d i m e n s i o n i n cm
9 l =4; // l e n g t h o f c i r c u l a r w a v e g u i d e i n cm
10 disp ( ’ Given t h e mode o f o p e r a t o r i s TE101 s o h e r e m
=1 , n =0 , p=1 ’ ) ;
12 m =1;
13 n =0;
14 p =1;
15
16 //CALCULATION
17 f0 =(( c /2) * sqrt (( m / a ) ^2+( n / b ) ^2+( p / l ) ^2) ) /10^9; // The
r e s o n a n t f r e q u e n c y o f a c i r c u l a r r e s o n a t o r i n GHz
18
19 //OUTPUT
20 mprintf ( ’ \ nThe r e s o n a n t f r e q u e n c y o f a c i r c u l a r
48
r e s o n a t o r i s f 0=%1 . 2 f GHz ’ , f0 ) ;
21
22 //=========================END OF PROGRAM
===============================
49
Chapter 6
Microwave Components
Scilab code Exa 6.2 distance to be shifted
1 // Chapter −6 , Example 6 . 2 , Page 234

2 //
==================================================================
3 // I n p u t p a r a m e t e r s
4 // [ s ] = [ 0 , ( 0 . 3 + ( %i ) ∗ ( 0 . 4 ) ) ; ( 0 . 3 + ( %i ) ∗ ( 0 . 4 ) ) , 0 ] ; / /
s c a t t e r i n g m a t r i x o f a two p o r t
5 // C a l c u l a t i o n s
6 // t o f i n d l s u c h t h a t S12 and S21 w i l l be r e a l when
p o r t 1 i s s h i f t e d lm t o t h e l e f t
7 // l e t p o r t 1 be s h i f t e d by p h i 1 d e g r e e t o t h e l e f t
and p o r t 2 p o s i t i o n be r e m a i n e d unchanged i . e . ,
p h i 2=d e l t a
8 // Then [ p h i ] = [ e ˆ−( j ∗ p h i 1 ) , 0 ; 0 , 1 ]
9 // [ S ’ ] = [ p h i ] ∗ [ s ] ∗ [ p h i ]
10 // f o r S12 and S21 t o be r e a l
11 phi1 =53.13; // i n d e g r e e s
12 phi1 = phi1 *( %pi /180) ; // p h i i n r a d i a n s
13 b =34.3; // me a s u r e d i n r a d /m
14 l =( phi1 ) / b ; // d i s t a n c e o f s h i f t i n m
15 // Output
50
16 mprintf ( ” d i s t a n c e t h a t t h e p o s i t i o n o f p a r t 1 s h o u l d
be s h i f t e d t o t h e l e f t s o t h a t S21 and S12 w i l l
be r e a l numbers i s %1 . 4 f m” ,l )
17 //=================================END OF PROGRAM
==============================
Scilab code Exa 6.3 Scattering parameters

2 //
==================================================================
3 clc ;
5 D =30; // d i r e c t i v i t y i n dB
6 VSWR =1; //VSWR a t e a c h p o r t u n d e r matched c o n d i t i o n s
7 C =10; // c o u p l i n g f a c t o r
9 S41 = sqrt (0.1) ;
10 S14 = S41 ; // u n d e r matched and l o s s l e s s c o n d i t i o n s
11 S31 = sqrt ((( S41 ) ^2) /(10) ^( D /10) ) ;
12 S13 = S31 ;
13 S11 =( VSWR -1) /( VSWR +1) ;
14 S22 = S11 ;
15 S33 = S22 ;
16 S44 = S33 ;
17 // l e t i n p u t power i s g i v e n a t p o r t 1
18 // p1=p2+P3+p4
19 S21 = sqrt (1 -( S41 ) ^2 -( S31 ) ^2) ;
20 S12 = S21 ;
21 S34 = sqrt ((0.5) *(1+( S12 ) ^2 -0.1 -0.0001) ) ;
22 S43 = S34
23 S23 = sqrt (1 -10^ -4 -( S34 ) ^2)
24 S32 = S23 ;
25 S24 = sqrt (1 -0.1 -( S34 ) ^2)
51
26 S42 = S24 ;
27 [ S ]=[ S11 , S12 , S13 , S14 ; S21 , S22 , S23 , S24 ; S31 , S32 , S33 , S34
; S41 , S42 , S43 , S44 ];
28 // Output
29 mprintf ( ” The s c a t t e r i n g m a t r i x i s ” ) ;
30 disp ([ S ])
31 //=================================END OF PROGRAM
==============================
Scilab code Exa 6.4 powers in the remaining ports

2 //
==================================================================
3 clc ;
5 a1 =32*10^ -3; // power i n w a t t s
6 a2 =0;
7 a3 =0;
9 [ S ]=[0.5 , -0.5 ,0.707; -0.5 ,0.5 ,0.707;0.707 ,0.707 ,0]; //
S−m a t r i x f o r H−p l a n e t e e
10 // [ B] = [ b1 , b2 , b3 ]
11 [ B ]=[ S ].*[ a1 ,0 ,0;0 ,0 ,0;0 ,0 ,0];
12 b1 =(0.5) ^2* a1 ; // power a t p o r t 1
13 b2 =( -0.5) ^2* a1 ; // power a t p o r t 2
14 b3 =(0.707) ^2* a1 ; // power a t p o r t 3
15 // Output
16 mprintf ( ” Thus b1 , b2 , b3 a r e %g W, %g W, %g W
r e s p e c t i v e l y ” ,b1 , b2 , b3 ) ;
17 //=================================END OF PROGRAM
==============================
52

2 //
==================================================================
3 clc ;
5 [ S ]=[0.5 , -0.5 ,0.707; -0.5 ,0.5 ,0.707;0.707 ,0.707 ,0];
6 R1 =60; // l o a d a t p o r t 1 i n ohms
7 R2 =75; // l o a d a t p o r t 2 i n ohms
8 R3 =50; // c h a r a c t e r i s t i c i m p e d a n c e i n ohms
9 P3 =20*10^ -3; // power a t p o r t 3 i n Watts
10 // c a l c u l a t i o n s
11 p1 =( R1 - R3 ) /( R1 + R3 ) ;
12 p2 =( R2 - R3 ) /( R2 + R3 ) ;
13 P1 =0.5* P3 *(1 -( p1 ) ^2) ; // power d e l i v e r e d t o t h e p o r t 1
i n Watts
14 P2 =0.5* P3 *(1 -( p2 ) ^2) ; // power d e l i v e r e d t o t h e p o r t 2
i n Watts
15 // Output
16 mprintf ( ” Thus power d e l i v e r e d t o t h e p o r t 1 and p o r t 2
a r e %g W, %g W r e s p e c t i v e l y ” ,P1 , P2 ) ;
17 //=================================END OF PROGRAM
==============================
Scilab code Exa 6.6 Reflected power

2 //
==================================================================
53
3 clc ;
5 p1 =0.5; // r e f l e c t i o n c o e f f i c i e n t a t p o r t 1
7 p3 =1; // r e f l e c t i o n c o e f f i c i e n t a t p o r t 3
9 // [ S
]=[0 ,0 ,0.707 ,0.707;0 ,0 ,0.5 , −0.707;0.707 ,0.707 ,0 ,0; −0.707 ,0.707 ,0 ,
S m a t r i x o f magic Tee
10 // s o l v i n g f o r b1 , b2 , b3 , b4 we g e t i t a s
11 // c a l c u l a t i o n s
12 b1 =0.6566;
13 b2 =0.7576;
14 b3 =0.6536;
15 b4 =0.0893;
16 P1 =( b1 ) ^2; // power a t p o r t 1 i n w a t t s
17 disp ( P1 ) ;
19 disp ( P2 ) ;
21 disp ( P3 ) ;
23 disp ( P4 ) ;
24 //=================================END OF PROGRAM
==============================
Scilab code Exa 6.7 Scattering matrix

2 //
==================================================================
3 clc ;
5 ins =0.5; // i n s e r t i o n l o s s i n db
54
6 iso =30; // i s o l a t i o n l o s s i n db
8 S21 =10^ -( ins /20) ; // i n s e r t i o n l o s s =0.5=−20∗ l o g [ S21 ]
9 S12 =10^ -( iso /20) ; // i s o l a t i o n l o s s =30=−20∗ l o g [ s 1 2 ]
10 S11 =0;
11 S22 =0;
12 [ S ]=[ S11 , S12 ; S21 , S22 ];
13 disp ( S ) ;
14 //=================================END OF PROGRAM
==============================
Scilab code Exa 6.9 Scattering matrix

2 //
==================================================================
3 clc ;
5 ins =0.5; // i n s e r t i o n l o s s i n db
6 iso =20; // i s o l a t i o n l o s s i n db
7 S =2; //VSWR
9 S21 =10^ -( ins /20) ; // i n s e r t i o n l o s s =0.5=−20∗ l o g [ S21 ]
10 S13 = S21 ;
11 S32 = S13 ;
12 S12 =10^ -( iso /20) ; // i s o l a t i o n l o s s =30=−20∗ l o g [ s 1 2 ]
13 S23 = S12 ;
14 S31 = S23 ;
15 p =( S -1) /( S +1) ;
16 S11 = p ;
17 S22 = p ;
18 S33 = p ;
19 [ S ]=[ S11 , S12 , S13 ; S21 , S22 , S23 ; S31 , S32 , S33 ];
20 disp ( S ) ;
55
21 // f o r a p e r f e c t l y matched , non−r e c i p r o c a l , l o s s l e s s 3−
p o r t c i r c u l a t o r , [ S ] i s g i v e n by
22 // [ S ] = [ 0 , 0 , S13 ; S21 , 0 , 0 ; , 0 , S32 , 0 ]
23 // i . e . , S13=S21=S32=1
24 // [ S ] = [ 0 , 0 , 1 ; 1 , 0 , 0 ; 0 , 1 , 0 ]
25 //=================================END OF PROGRAM
==============================
Scilab code Exa 6.10 output powers
1 // Chapter −6 , Example 6 . 1 0 , Page 242

2 //
==================================================================
3
4 clc ;
5 clear ;
6 close ;
7 In_loss =0.5; // i n s e r t i o n l o s s ( i n dB )
8 C =20; // c o u p l i n g c o e f f i c i e n t i n dB
9 D =35; // d i r e c t i v i t y i n dB
10 Pi_Pf =10^( C /10) ;
11 Pi =90; // i n Watts
12 Pf = Pi / Pi_Pf ;
13 Pf_Pb =10^( D /10) ;
14 Pb = Pf / Pf_Pb ;
15 P_rec =( Pi - Pf - Pb ) ; // Power r e c e i v e d ( i n
Watts )
16 P_rec_dB =10* log ( Pi / P_rec ) / log (10) ;
17 P_rec_eff = P_rec_dB - In_loss ; // E f f e c t i v
e power r e c e i v e d ( i n dB )
18 disp ( Pf , ’ Output power t h r o u g h c o u p l e d p o r t ( i n
Watts )= ’ ) ;
19 disp ( Pb , ’ Output power t h r o u g h i s o l a t e d p o r t ( i
n Watts )= ’ ) ;
56
20 disp ( P_rec_dB , ’ Power r e c e i v e d ( i n dB )=
’ );
21 disp ( P_rec_eff , ’ E f f e c t i v e power r e c e
i v e d ( i n dB )= ’ ) ;
22
23
24 //=================================END OF PROGRAM
==============================
Scilab code Exa 6.11 coupling isolation directivity
1 // Chapter −6 , Example 6 . 1 1 , Page 242

2 //
==================================================================
3 clc ;
5 S13 =0.1*( cos (90* %pi /180) +( %i ) * sin (90* %pi /180) ) ; //
c o n v e r s i o n from p o l a r t o r e c t a n g u l a r
6 S13 = abs ( S13 ) ;
7 C = -20* log10 ( S13 ) ; // c o u p l i n g c o e f f i c i e n t i n dB
8 S14 =0.05*( cos (90* %pi /180) +( %i ) * sin (90* %pi /180) ) ; //
c o n v e r s i o n from p o l a r t o r e c t a n g u l a r
9 S14 = abs ( S14 ) ;
10 D =20* log10 ( S13 / S14 ) ; // d i r e c t i v i t y i n dB
11 I = -20* log10 ( S14 ) ; // i s o l a t i o n i n dB
12 mprintf ( ” Thus c o u p l i n g , d i r e c t i v i t y and i s o l a t i o n a r e
%1 . 0 f dB , %1 . 2 f dB and %2 . 2 f dB r e s p e t i v e l y ” ,C ,D
,I);
13 //=================================END OF PROGRAM
==============================
Scilab code Exa 6.12 VSWR
57
2 //
==================================================================
3 clc ;
4 clear ;
5
6 x =3.5; // d i s t a n c e b e t w e e n two minimas i n cm
7 y =0.25; // d i s t a n c e b e t w e e n t w i c e minimum power p o i n t s
i n cm
8
9 //CALCULATION
10 wg =2* x ; // g u i d e d w a v e l e n g t h i n cm
11 S =( wg /( y *( %pi ) ) ) ; // V o l t a g e S t a n d i n g Wave R a t i o (VSWR)
12
13 //OUTPUT
14 mprintf ( ’ \ n V o l t a g e S t a n d i n g Wave R a t i o (VSWR) i s S=%1
. 4 f ’ ,S ) ;
15
16 //=========================END OF PROGRAM
===============================
Scilab code Exa 6.13 phase shift
2 //
==================================================================
3 clc ;
4 clear ;
5
6 wg =7.2; // g u i d e w a v e l e n g t h i n cm
7 x =10.5; // P o s i t i o n o f r e f e r e n c e n u l l w i t h o u t t h e
w a v e g u i d e component i n cm
8 y =9.3; // P o s i t i o n o f r e f e r e n c e n u l l with the
58
w a v e g u i d e component i n cm
9
10 //CALCULATION
11 z =x - y ; // Path d i f f e r e n c e
i n t r o d u c e d due t o t h e
component i n cm
12 p =(2*( %pi ) *( z / wg ) ) ; // Phase d i f f e r e n c e i n t r o d u c e d i n
rad
13 Pd =( p *180) /( %pi ) ; // Phase s h i f t i n t r o d u c e d i n deg
14
15 //OUTPUT
16 mprintf ( ’ \ nPhase s h i f t i n t r o d u c e d i s Pd=%2 . 0 f deg ’ ,
Pd ) ;
17
18 //=========================END OF PROGRAM
===============================
59
Chapter 7
Microwave Measurements
Scilab code Exa 7.1 measured distance
2 //
==================================================================
3 clc ;
4 clear ;
5
6 a =4; // Le n g th o f Waveguide i n cm
7 b =2.5; // b r e a d t h Waveguide i n cm
8 f =10^10; // F r e q u e n c y i n Hz
9 x =0.1; // d i s t a n c e b e t w e e n t w i c e minimum power p o i n t s
i n cm
11
12 //CALCULATION
13 wc =2* a ; // C u t o f f w a v e l e n g t h i n TE10 mode i n cms
16 S =( wg /( x *( %pi ) ) ) ; // V o l t a g e S t a n d i n g Wave R a t i o (VSWR)
f o r d o u b l e minimum method
17
60
18 //OUTPUT
19 mprintf ( ’ \ nFor d o u b l e minimum method , V o l t a g e
S t a n d i n g Wave R a t i o (VSWR) i s S=%2 . 1 f ’ ,S ) ;
20
21 //=========================END OF PROGRAM
===============================
Scilab code Exa 7.2 VSWR and Reflected power
2 //
==================================================================
3 clc ;
4 clear ;
5
6 x =3; //O/P i n c i d e n t power from first directional
c o u p l e r i n mW
7 y =0.1; //O/P r e f l e c t e d power from s e c o n d d i r e c t i o n a l
8
9 //CALCULATION
10 Pi = x *100; // I n c i d e n t Power i n mW
11 Pr = y *100; // R e f l e c t e d Power i n mW
12 p = sqrt ( Pr / Pi ) ; // R e f l e c t i o n C o e f f i c i e n t
14
15 //OUTPUT
16 mprintf ( ’ \ n V o l t a g e S t a n d i n g Wave R a t i o (VSWR) i n t h e
main w a v e g u i d e i s S=%1 . 2 f \ n R e f l e c t e d Power i s Pr
=%2 . 0 f mW’ ,S , Pr ) ;
17
18 //=========================END OF PROGRAM
===============================
61
Scilab code Exa 7.3 VSWR
2 //
==================================================================
3 clc ;
4 clear ;
5
6 Pi =2.5; // I n c i d e n t Power from one d i r e c t i o n a l coupler
i n mW
7 Pr =0.15; // R e f l e c t e d Power from o t h e r d i r e c t i o n a l
8
9 //CALCULATION
10 p = sqrt ( Pr / Pi ) ; // R e f l e c t i o n C o e f f i c i e n t
12
13 //OUTPUT
14 mprintf ( ’ \ n V o l t a g e S t a n d i n g Wave R a t i o (VSWR) i n t h e
w a v e g u i d e i s S=%1 . 2 f ’ ,S ) ;
15
16 //=========================END OF PROGRAM
===============================
Scilab code Exa 7.4 Reflected power
2 //
==================================================================
3 clc ;
62
4 clear ;
5
6 S =2; // V o l t a g e S t a n d i n g Wave R a t i o (VSWR)
7 C =30; // C o u p l i n g Power o f a D i r e c t i o n a l C o u p l e r i n dB
8 Pf =4.5; // C o u p l e r I n c i d e n t S a m p l i n g Power i n mW
9
10 //CALCULATION
11 p =(( S -1) /( S +1) ) ; // R e f l e c t i o n C o e f f i c i e n t
12 Pi = Pf *10^( C /10) ; // I n c i d e n t Power i n mW [ From C=10 l o g
( Pi / Pf ) ]
13 Pr =( Pi *( p ^2) ) /10^3; // R e f l e c t e d Power i n W [ From p=
s q r t ( Pr / Pi ) ]
14
15 //OUTPUT
16 mprintf ( ’ \ nValue o f R e f l e c t e d Power i s Pr=%1 . 2 f W’ ,
Pr ) ;
17
18 //=========================END OF PROGRAM
===============================
63
Chapter 8
Microwave Tubes and Circuits
Scilab code Exa 8.1 electron velocity and etc
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a f o u r c a v i t y K l y s t r o n
7 V0 =14500; //Beam v o l t a g e i n V
8 I =1.4; //Beam c u r r e n t i n A
9 f =10^10; // O p e r a t i o n f r e q u e n c y i n Hz
10 p0 =10^( -6) ; // dc e l e c t r o n c h a r g e d e n s i t y i n C/mˆ3
11 p =10^( -8) ; //RF c h a r g e d e n s i t y i n C/mˆ3
12 V =10^5; // V e l o c i t y p e r t u r b a t i o n s i n m/ s e c
13 e0 =8.854*10^( -12) ; // P e r m i t t i v i t y o f f r e e s p a c e i n F/
m
14 R =0.4;
15
16 //CALCULATION
17 v0 =(0.593*10^6* sqrt ( V0 ) ) /10^8; // The dc e l e c t r o n
v e l o c i t y i n 1 0 ˆ 8 m/ s e c
64
18 w =2*( %pi ) * f ; // a n g u l a r f r e q u e n c y i n r a d / s e c
19 v = v0 *10^8;
20 c =( w / v ) ; // The dc Phase C o n s t a n t
21 wp =( sqrt (1.759*10^11*( p0 / e0 ) ) ) /10^8; // The Plasma
Frequency i n 10ˆ8 rad / s e c
22 wp1 = wp *10^8;
23 wq =( R * wp1 ) /10^8; // The Reduced Plasma F r e q u e n c y i n
10ˆ8 rad / s e c
24 J0 = p0 * v ; // The dc beam c u r r e n t d e n s i t y i n A/sqm
25 J =( p * v ) +( p0 * V ) ; // The i n s t a n t a n e o u s beam c u r r e n t
d e n s i t y i n A/sqm
26
27 //OUTPUT
28 mprintf ( ’ \ nThe dc e l e c t r o n v e l o c i t y i s v0=%2 . 3 f
∗ 1 0 ˆ 8 m/ s e c \ nThe dc Phase C o n s t a n t i s c=%1 . 2 f
r a d / s e c \ nThe Plasma F r e q u e n c y i s wp=%1 . 2 f ∗ 1 0 ˆ 8
r a d / s e c \ nThe Reduced Plasma F r e q u e n c y i s wq=%1 . 3
f ∗ 1 0 ˆ 8 r a d / s e c \ nThe dc beam c u r r e n t d e n s i t y i s
J0=%2 . 1 f A/sqm \ nThe i n s t a n t a n e o u s beam c u r r e n t
d e n s i t y i s J=%1 . 3 f A/sqm ’ ,v0 ,c , wp , wq , J0 , J ) ;
29
30 //=========================END OF PROGRAM
===============================
Scilab code Exa 8.2 power and etc
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a 2 c a v i t y k l y s t r o n a m p l i f i e r
7 Av =15; // V o l t a g e g a i n i n dB
65
8 Pin =0.005; // I /P power i n W
9 Rin =30000; // Rsh o f i / p c a v i t y i n ohms
10 R0 =40000; // Rsh o f o / p c a v i t y i n ohms
11 Rl =40000; // l o a d i m p e d a n c e i n ohms
12 R =20000; // P a r a l l e l r e s i s t a n c e o f R0 and Rl ( R0// Rl )
i n ohms
13
14 //CALCULATION
15 Vin = sqrt ( Pin * Rin ) ; // The i n p u t rms v o l t a g e in V [
From Pin=Vin ˆ2/ Rin ]
16 V0 = Vin *10^( Av /20) ; // The o u t p u t rms v o l t a g e i n V [
From Av=20 l o g ( V0/ Vin ) ]
17 P0 =( V0 ^2) / R ; // The Power d e l i v e r e d t o t h e l o a d i n W
18
19 //OUTPUT
20 mprintf ( ’ \ nThe i n p u t rms v o l t a g e i s Vin=%2 . 2 f V \
nThe o u t p u t rms v o l t a g e i s V0=%2 . 2 f V \ nThe Power
d e l i v e r e d t o t h e l o a d i s P0=%1 . 4 f W’ ,Vin , V0 , P0 ) ;
21
22 //=========================END OF PROGRAM
===============================
Scilab code Exa 8.3 efficiency and etc
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a r e f l e x k l y s t r o n
7 n =2; // peak mode v a l u e
8 V0 =300; // beam v o l t a g e i n V
9 I0 =0.02; // beam c u r r e n t i n A
66
10 Vs =40; // s i g n a l v o l t a g e i n V
11 J1 =1.25; // b e s s e l c o e f f i c i e n t f o r n=2
12
13 //CALCULATION
14 Pdc = V0 * I0 ; // The i n p u t power i n w a t t s
15 Pac =((2* Pdc * J1 ) /((2* n *( %pi ) ) -(( %pi ) /2) ) ) ; // The
o u t p u t power i n w a t t s
16 n =( Pac / Pdc ) *100; // E f f i c i e n c y i n p e r c e n t a g e
17
18 //OUTPUT
19 mprintf ( ’ \ nThe i n p u t power i s Pdc=%1 . 0 f w a t t s \ nThe
o u t p u t power i s Pac=%1 . 2 f w a t t s \ n E f f i c i e n c y i s n
=%2 . 1 f p e r c e n t a g e ’ ,Pdc , Pac , n ) ;
20
21 //=========================END OF PROGRAM
===============================
2 //
==================================================================
3 clc ;
4 clear ;
5
8 I0 =0.03; //Beam c u r r e n t i n A
9 f =8*10^9; // f r e q u e n c y i n Hz
10 d =0.001; // gap s p a c i n g i n e i t h e r c a v i t y i n m
11 L =0.04; // s p a c i n g b e t w e e n c e n t e r s o f c a v i t i e s i n m
12 Rsh =49000; // E f f e c t i v e s h u n t i m p e d a n c e i n ohms
13 J1 =0.582; // v a l u e o f J1 (X)
14 X =1.841; // b u n c h i n g p a r a m e t e r
67
15
16 //CALCULATION
17 v0 =(0.593*10^6* sqrt ( V0 ) ) /10^6; // v e l o c i t y of electron
i n 1 0 ˆ 6 m/ s e c
18 w =2*( %pi ) * f ; // a n g u l a r f r e q u e n c y i n r a d
19 v = v0 *10^6;
20 T0 =( L / v ) /10^( -8) ; // dc t r a n s i t t i m e o f e l e c t r o n s in
10ˆ( −8) s e c
21 a = w * T0 *10^( -8) ; // t r a n s i t a n g l e i n r a d
22 tg = w * d / v ; // a v e r a g e gap t r a n s i t a n g l e i n r a d
23 Tg = tg *(180/( %pi ) ) ;
24 Bi =( sind ( Tg /2) ) /( tg /2) ; // beam c o u p l i n g c o e f f i c i e n t
25 Bo = Bi ; // o u t p u t c a v i t y c o u p l i n g c o e f f i c i e n t
26 V1max =((3.68* V0 ) /( Bi * a ) ) ; // I n p u t v o l t a g e f o r Maximum
output v o l t a g e in V
27 R0 = V0 / I0 ; // i m p e d a n c e i n ohms
28 Av =( Bo ^2* a * Rsh * J1 ) /( R0 * X ) ; // V o l t a g e g a i n
29 AvdB =20* log10 ( Av ) ; // V o l t a g e g a i n i n dB
30
31 //OUTPUT
32 mprintf ( ’ \ n V e l o c i t y o f e l e c t r o n i s v0=%2 . 2 f ∗ 1 0 ˆ 6 m/
s e c \ nThe dc t r a n s i t t i m e o f e l e c t r o n s i s T0=%1 . 3
f ∗10ˆ( −8) s e c \ n I n p u t v o l t a g e f o r Maximum o u t p u t
v o l t a g e i s V1max=%2 . 3 f V \ n V o l t a g e g a i n i s Av=%2
. 2 f \ nThe V o l t a g e g a i n i n dB i s AvdB=%2 . 2 f dB ’ ,v0
, T0 , V1max , Av , AvdB ) ;
33
34 //=========================END OF PROGRAM
===============================
35
36 // Note : Check t h e c a l c u l a t i o n g i v e n i n t e x t book f o r
v o l t a g e g a i n Rsh=49 kohms
37 // but , t a k e n a s 40 kohms
38 // c o r r e c t a n s w e r s a r e V o l t a g e g a i n i s Av= 2 8 . 5 2
39 // The V o l t a g e g a i n i n dB i s AvdB= 2 9 . 1 0 dB
68
2 //
==================================================================
3 clc ;
4 clear ;
5
9 f =8*10^9; // f r e q u e n c y i n Hz
13 J1 =0.582; // v a l u e o f J1 (X)
15
16 //CALCULATION
18 v0 =0.593*10^6* sqrt ( V0 ) ; // v e l o c i t y o f e l e c t r o n i n m/
sec
19 Vomax =((3.68* V0 * v0 ) /( w * L ) ) ; //max o u t p u t power i n V
20 tg =( w * d ) / v0 ; // avg gap t r a n s i t a n g l e i n r a d
21 Tg = tg *(180/( %pi ) ) ;
22 Bi =( sind ( Tg /2) ) /( tg /2) ; // beam c o u p l i n g c o e f f i c i e n t
23 Bo = Bi ; // o u t p u t c a v i t y c o u p l i n g c o e f f i c i e n t
24 Vimax = Vomax / Bi ; // The i n p u t m i c r o w a v e v o l t a g e i n
o r d e r t o g e n e r a t e maximum o u t p u t v o l t a g e i n V
25 t0 = w * L / v0 ; // t r a n s i t a n g l e i n r a d
26 R0 = V0 / I0 ; // i m p e d a n c e i n ohms
27 Av =(( Bo ^2* J1 * t0 * Rsh ) /( R0 * X ) ) ; // V o l t a g e g a i n
28 I2 =2* I0 * J1 ;
69
29 V2 = Bo * I2 * Rsh ;
30 disp ( ’ n e g l e c t i n g beam l o a d i n g ’ ) ;
31 Eff =0.58*( V2 / V0 ) *100; // E f f i c i e n c y i n %
32 G0 =1/ R0 ;
33 GB =( G0 /2) *( Bo *( Bo - cos ( Tg /2) ) ) ; //Beam l o a d i n g
c o n d u c t a n c e i n mhos
34 RB =(1/ GB ) /1000; //Beam l o a d i n g r e s i s t a n c e i n Kohms
35 disp ( ’ Beam l o a d i n g r e s i s t a n c e i n Kohms i s ’ ) ;
36 disp ( RB ) ;
37 disp ( ’ The v a l u e 73 kohms i s v e r y much c o m p a r a b l e t o
Rsh and c a n n o t be n e g l e c t e d b e c a u s e Tg i s q u i t e
high ’ );
38
39 //OUTPUT
40 mprintf ( ’ \ nThe i n p u t m i c r o w a v e v o l t a g e i n o r d e r t o
g e n e r a t e maximum o u t p u t v o l t a g e i s Vimax=%2 . 2 f V
\ nThe v o l t a g e g a i n i s Av=%2 . 2 f p e r c e n t a g e \nBeam
l o a d i n g c o n d u c t a n c e i s GB=%1 . 1 0 f mhos ’ , Vimax , Av ,
GB ) ;
41
42 //=========================END OF PROGRAM
===============================
Scilab code Exa 8.6 electronic efficiency and etc
2 //
==================================================================
3 clc ;
4 clear ;
5
70
9 Rsh =20000; // Shunt r e s i s t a n c e i n ohms
10 L =0.001; // d i s t a n c e i n m
11 f =8*10^(9) ; // // O p e r a t i o n f r e q u e n c y i n Hz
12 V1 =200; // m i c r o w a v e gap v o l t a g e i n V
13 x =1.759*10^11; // e /m v a l u e i n C/ kg
14 J1 =0.582;
15
16 //CALCULATION
17 disp ( ’ Assume t h e gap t r a n s i t t i m e and beam l o a d i n g
a r e n e g l e c t e d ’ );
19 VR =( V0 +(( sqrt (8* V0 / x ) * w * L ) /((2*( %pi ) * n ) -(( %pi ) /2) ) ) )
; // R e p e l l e r v o l t a g e i n V
20 disp ( ’ Assuming o u t p u t c o u p l i n g c o e f f i c i e n t Bo=1 ’ ) ;
21 I0 =( V1 /(2* J1 * Rsh ) ) /10^( -3) ; //Beam c u r r e n t n e c e s s a r y
t o o b t a i n an m i c r o w a v e gap v o l t a f e o f 200V i n mA
22 v0 =0.593*10^6* sqrt ( V0 ) ; // v e l o c i t y o f e l e c t r o n i n m/
sec
23 t0 =(( w *2* L * v0 ) /( x *( VR + V0 ) ) ) ; // t r a n s i t a n g l e i n r a d
24 Bi =1; // beam c o u p l i n g c o e f f i c i e n t [ assume ]
25 X =(( Bi * V1 * t0 ) /(2* V0 ) ) ;
26 disp ( ’ S i n c e X= 1 . 5 1 , from graph , J1 (X) =0.84 ’ ) ;
27 XJ1 =0.84;
28 Eff =((2*( XJ1 ) ) /((2* n *( %pi ) ) -(( %pi ) /2) ) ) *100 //
E f fi c i e n cy in %
29
30 //OUTPUT
31 mprintf ( ’ \ n R e p e l l e r v o l t a g e i s VR=%3 . 2 f V \ nThe dc
n e c e s s a r y t o g i v e an m i c r o w a v e gap v o l t a f e o f 200
V i s I 0=%1 . 2 f mA \ n E l e c t r o n i c E f f i c i e n c y i s E f f=
%2 . 2 f p e r c e n t a g e ’ ,VR , I0 , Eff ) ;
32
33 //=========================END OF PROGRAM
===============================
34
35 // Note : Check t h e a n s w e r f o r VR o n c e
36 // C o r r e c t a n s w e r i s R e p e l l e r v o l t a g e i s VR= 1 1 8 9 . 3 6 V
71
2 //
==================================================================
3 clc ;
4 clear ;
5
7 n =1; // mode v a l u e
8 Pi =0.04; // t h e dc i n p u t power i n W
9 x =0.278; // r a t i o o f V1 o v e r V0
10
11 //CALCULATION
12 X = x *3*( %pi ) /4;
13 J1 =0.3205; // b e s s e l c o e f f i c i e n t v a l u e [ J I (X ’ ) ]
14 ef =((2* X * J1 ) /((2*( %pi ) * n ) -(( %pi ) /2) ) ) *100; //
E f f i c i e n c y of the r e f l e x k l y s t r o n in %
15 Pout =(( ef /100) * Pi ) /10^( -3) ; // T o t a l power o u t p u t i n
mW
16 p =20; // p e r c e n t a g e o f t h e power d e l i v e r e d by t h e
e l e c t r o n beam d i s s i p a t e d i n t h e c a v i t y w a l l s
17 Pd = Pout *(100 - p ) /100; // Power d e l i v e r e d t o l o a d i n mW
18
19 //OUTPUT
20 mprintf ( ’ \ n E f f i c i e n c y o f t h e r e f l e x k l y s t r o n i s e f=
%1 . 2 f p e r c e n t a g e \ n T o t a l power o u t p u t i s Pout=%1 . 3
f mW \ n I f t h e 20 p e r c e n t a g e o f t h e power
d e l i v e r e d by t h e e l e c t r o n beam i s d i s s i p a t e d i n
t h e c a v i t y w a l l s t h e n t h e Power d e l i v e r e d t o l o a d
i s Pd=%1 . 2 f mW’ ,ef , Pout , Pd ) ;
21
22 //=========================END OF PROGRAM
72
===============================
Scilab code Exa 8.8 cyclotron frequency and etc
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a c i r c u l a r magnetron
7 a =0.15; // i n n e r r a d i u s i n m
8 b =0.45; // o u t e r r a d i u s i n m
9 B =1.2*10^( -3) ; // m a g n e t i c f l u x d e n s i t y i n Wb/sqm
10 x =1.759*10^11; // V a l u e o f e /m i n C/ kg
11 V =6000; // beam v o l t a g e i n V
12
13 //CALCULATION
14 V0 =(( x /8) *( B ^2) *( b ^2) *(1 -( a / b ) ^2) ^2) /1000; // H u l l cut
− o f f v o l t a g e i n kV
15 Bc =(( sqrt (8*( V / x ) ) ) /( b *(1 -( a / b ) ^2) ) ) *1000; // Cut− o f f
m a g n e t i c f l u x d e n s i t y i n mWb/sqm
16 fc =(( x * B ) /(2*( %pi ) ) ) /10^9; // C y c l o t r o n f r e q u e n c y i n
GHz
17
18 //OUTPUT
19 mprintf ( ’ \ n H u l l cut − o f f v o l t a g e i s V0=%2 . 3 f kV\ nCut−
o f f m a g n e t i c f l u x d e n s i t y i s Bc=%1 . 6 f mWb/sqm \
n C y c l o t r o n f r e q u e n c y i s f c=%1 . 4 f GHz ’ ,V0 , Bc , fc ) ;
20
21 //=========================END OF PROGRAM
===============================
22
23
73
24 // Check t h e a n s w e r s o n c e
26 // H u l l cut − o f f v o l t a g e i s V0 = 5 . 0 6 6 kV
27 // Cut− o f f m a g n e t i c f l u x d e n s i t y i s Bc = 1 . 3 0 5 9 5 3 mWb/
sqm
28 // C y c l o t r o n f r e q u e n c y i s f c = 0 . 0 3 3 6 GHz
Scilab code Exa 8.9 phase velocity and anode voltage
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a h e l i c a l TWT
7 c =3*10^8; // V e l o c i t y o f l i g h t i n m/ s e c
8 d =0.002; // d i a m e t e r i n m
9 x =5000; // no . o f t u r n s p e r m
10 m =9.1*10^( -31) ; // mass o f an e l e c t r o n i n kg
11 e =1.6*10^( -19) ; // c h a r g e o f an e l e c t r o n i n C
12
13 //CALCULATION
14 y =( %pi ) * d ; // c i r c u m f e r e n c e i n m
15 p =1/ x ; // p i t c h i n m
16 Vp =( c * p ) / y ; // A x i a l p h a s e v e l o c i t y i n m/ s e c
17 V0 =(( m * Vp ^2) /(2* e ) ) ; // The Anode v o l t a g e a t which t h e
TWT can be o p e r a t e d f o r u s e f u l g a i n i n V
18
19 //OUTPUT
20 mprintf ( ’ \ n A x i a l p h a s e v e l o c i t y i s Vp=%6 . 2 f m/ s e c \
nThe Anode v o l t a g e a t which t h e TWT can be
o p e r a t e d f o r u s e f u l g a i n i s V0=%2 . 2 f V ’ ,Vp , V0 ) ;
21
74
22 //=========================END OF PROGRAM
===============================
2 //
==================================================================
3 clc ;
4 clear ;
5
9 f =8*10^9; // f r e q u e n c y i n Hz
13 y =0.582; // v a l u e o f J1 (X)
14 X =1.841;
15
16 //CALCULATION
17 v0 =(0.593* sqrt ( V0 ) *10^6) /10^7; // The e l e c t r o n
v e l o c i t y i n 1 0 ˆ 7 m/ s e c
18 v = v0 *10^7;
19 t0 =( d / v ) /10^( -10) ; // T r a n s i t t i m e i n 10ˆ( −10) s e c
20 t = t0 *10^( -10) ;
21 a =2*( %pi ) * f * t ; //Gap t r a n s i t a n g l e i n r a d
22 Bi =( sin ( a /2) ) /( a /2) ; //Beam c o u p l i n g c o e f f i c i e n t
23 Bo = Bi ;
24 to =(2*( %pi ) * f * L ) / v ; // dc t r a n s i t a n g l e i n r a d
25 disp ( ’ For maximum o u t o u t v o l t a g e , V2 J1 (X) = 0 . 5 8 2 ,X
=1.841 ’ );
26 V1 =((2* V0 * X ) /( Bo * to ) ) // The i n p u t v o l t a g e f o r maximum
75
output v o l t a g e in V
27 Ro =( V0 / I0 ) ;
28 Av =(( Bo ^2* to * y * Rsh ) /( Ro * X ) ) ; // V o l t a g e g a i n
30
31 //OUTPUT
32 mprintf ( ’ \ nThe e l e c t r o n v e l o c i t y i s v0=%1 . 1 f ∗ 1 0 ˆ 7 m
/ s e c \ nThe dc e l e c t r o n T r a n s i t t i m e i s t 0=%1 . 2 f
∗10ˆ( −10) s e c \ nThe i n p u t v o l t a g e f o r maximum
o u t p u t v o l t a g e i s V1=%2 . 2 f V \ n V o l t a g e g a i n i s Av
=%2 . 2 f \ n V o l t a g e g a i n i n dB i s AvdB=%2 . 2 f dB ’ ,v0 ,
t0 , V1 , Av , AvdB ) ;
33
34 //=========================END OF PROGRAM
===============================
Scilab code Exa 8.11 dc electron velocity and etc
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a f o u r c a v i t y K l y s t r o n
8 I =2; //Beam c u r r e n t i n A
9 f =9*10^9; // O p e r a t i o n f r e q u e n c y i n Hz
10 p0 =10^( -6) ; // dc e l e c t r o n c h a r g e d e n s i t y i n C/mˆ3
11 p =10^( -8) ; //RF c h a r g e d e n s i t y i n C/mˆ3
12 V =10^5; // V e l o c i t y p e r t u r b a t i o n s i n m/ s e c
13 e0 =8.854*10^( -12) ; // P e r m i t t i v i t y o f f r e e s p a c e i n F/
m
14 R =0.5;
76
15
16 //CALCULATION
17 v0 =(0.593*10^6* sqrt ( V0 ) ) /1000; // The dc e l e c t r o n
v e l o c i t y i n Km/ s e c
18 w =2*( %pi ) * f ; // a n g u l a r f r e q u e n c y i n r a d / s e c
19 v = v0 *1000;
20 c =( w / v ) ; // The dc Phase C o n s t a n t
21 wp =( sqrt (1.759*10^11*( p0 / e0 ) ) ) /10^8; // The Plasma
Frequency i n 10ˆ8 rad / s e c
22 wp1 = wp *10^8;
23 wq =( R * wp1 ) /10^8; // The Reduced Plasma F r e q u e n c y i n
10ˆ8 rad / s e c
24 J0 = p0 * v ; // The dc beam c u r r e n t d e n s i t y i n A/sqm
25 J =( p * v ) -( p0 * V ) ; // The i n s t a n t a n e o u s beam c u r r e n t
d e n s i t y i n A/sqm
26
27 //OUTPUT
28 mprintf ( ’ \ nThe dc e l e c t r o n v e l o c i t y i s v0=%4 . 2 f Km/
s e c \ nThe dc Phase C o n s t a n t i s c=%3 . 2 f r a d / s e c \
nThe Plasma F r e q u e n c y i s wp=%1 . 2 f ∗ 1 0 ˆ 8 r a d / s e c \
nThe Reduced Plasma F r e q u e n c y i s wq=%1 . 3 f ∗ 1 0 ˆ 8
r a d / s e c \ nThe dc beam c u r r e n t d e n s i t y i s J0=%2 . 2 f
A/sqm \ nThe i n s t a n t a n e o u s beam c u r r e n t d e n s i t y
i s J=%1 . 4 f A/sqm ’ ,v0 ,c , wp , wq , J0 , J ) ;
29
30 //=========================END OF PROGRAM
===============================
Scilab code Exa 8.12 gap transit angle
2 //
==================================================================
3 clc ;
77
4 clear ;
5
7 f =5*10^9; // F r e q u e n c y o f o p e r a t i o n i n hz
8 V0 =1000; // anode v o l t a g e i n V
9 d =0.002; // c a v i t y gap i n m
10 Vr = -500; // r e p e l l e r v o l t a g e i n V
11
12 //CALCULATION
13 N =7/4; // mode v a l u e
14 VR = abs ( Vr ) ;
15 L =((( VR + V0 ) * N ) /(6.74*10^( -6) * f * sqrt ( V0 ) ) ) /10^( -3) ; //
Optimum l e n g t h o f t h e d r i f t r e g i o n i n mm
16 u =5.93*10^5* sqrt ( V0 ) ; // i n m/ s e c
18 Tg =( w * d ) / u ; //Gap t r a n s i t a n g l e i n r a d
19
20 //OUTPUT
21 mprintf ( ’ \nOptimum l e n g t h o f t h e d r i f t r e g i o n i s L=
%1 . 3 f mm \nGap t r a n s i t a n g l e i s Tg=%1 . 3 f r a d ’ ,L ,
Tg ) ;
22
23 //=========================END OF PROGRAM
===============================
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a 2 c a v i t y klystron amplifier
78
9 f =10*10^9; // f r e q u e n c y i n Hz
13 J1 =0.582; // v a l u e o f J1 (X)
15
16 //CALCULATION
17 v0 =0.593*10^6* sqrt ( V0 ) ; // v e l o c i t y o f r e f e r e n c e
e l e c t r o n i n m/ s e c
19 a = w * L / v0 ; // t r a n s i t a n g l e w i t h o u t RF v o l t a g e i n r a d
20 tg = a * d / L ; // a v e r a g e gap t r a n s i t a n g l e i n r a d
21 Bi =( sin ( tg /2) ) /( tg /2) ; // beam c o u p l i n g c o e f f i c i e n t
22 V1max =((2* X * V0 ) /( Bi * a ) ) ; // I n p u t RF v o l t a g e f o r
Maximum o u t p u t v o l t a g e i n V
23 B0 = Bi ; // o u t p u t c a v i t y c o u p l i n g c o e f f i c i e n t
24 V2 =2* B0 * I0 * J1 * Rsh ; // i n V
25 Av = V2 / V1max ; // V o l t a g e g a i n
27 n =0.58*( V2 / V0 ) *100; //Maximum e f f i c i e n c y i n %
28
29 //OUTPUT
30 mprintf ( ’ \ n I n p u t RF v o l t a g e f o r Maximum o u t p u t
v o l t a g e i s V1max=%2 . 2 f V \ nThe V o l t a g e g a i n i s
AvdB=%2 . 2 f dB \nMaximum e f f i c i e n c y i s I 0=%2 . 2 f
p e r c e n t a g e ’ , V1max , AvdB , n ) ;
31
32 //=========================END OF PROGRAM
===============================
33
35 // There a r e s l i g h t c h a n g e s i n v a l u e s
36 // I n p u t RF v o l t a g e f o r Maximum o u t p u t v o l t a g e i s
V1max = 5 5 . 2 8 V
37 // The V o l t a g e g a i n i s AvdB= 2 4 . 3 5 dB
79
38 //Maximum e f f i c i e n c y i s I0 =44.11 p e r c e n t a g e
Scilab code Exa 8.14 cyclotron angular frequency and etc
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For aa X−band c y l i n d r i c a l magnetron
7 a =0.04; // i n n e r r a d i u s i n m
8 b =0.08; // o u t e r r a d i u s i n m
9 B =0.01; // m a g n e t i c f l u x d e n s i t y i n Wb/sqm
10 x =1.759*10^11; // V a l u e o f e /m i n C/ kg
11 V =30000; // beam v o l t a g e i n V
12
13 //CALCULATION
14 w =( x * B ) /10^9; // C y c l o t r o n a n g u l a r f r e q u e n c y i n 1 0 ˆ 9
rad / s e c
15 VHC =(( x /8) *( B ^2) *( b ^2) *(1 -( a / b ) ^2) ^2) /1000; // H u l l
cut − o f f v o l t a g e i n kV
16 Bc =(( sqrt (8*( V / x ) ) ) /( b *(1 -( a / b ) ^2) ) ) *1000; // Cut− o f f
m a g n e t i c f l u x d e n s i t y i n mWb/sqm
17
18 //OUTPUT
19 mprintf ( ’ \ n C y c l o t r o n a n g u l a r f r e q u e n c y i s w=%1 . 3 f
∗ 1 0 ˆ 9 r a d / s e c \ n H u l l cut − o f f v o l t a g e i s VHC=%1 . 4 f
kV \ nCut− o f f m a g n e t i c f l u x d e n s i t y i s Bc=%2 . 3 f
mWb/sqm ’ ,w , VHC , Bc ) ;
20
21 //=========================END OF PROGRAM
===============================
80
2 //
==================================================================
3 clc ;
4 clear ;
5
9 I0 =0.022; // beam c u r r e n t i n A
10 Vs =30; // s i g n a l v o l t a g e i n V
11 J1 =1.25; // b e s s e l c o e f f i c i e n t f o r n=2
12
13 //CALCULATION
14 Pdc = V0 * I0 ; // The i n p u t power i n w a t t s
15 Pac =((2* Pdc * J1 ) /((2* n *( %pi ) ) -(( %pi ) /2) ) ) ; // The
o u t p u t power i n w a t t s
16 n =( Pac / Pdc ) *100; // E f f i c i e n c y i n p e r c e n t a g e
17
18 //OUTPUT
19 mprintf ( ’ \ nThe i n p u t power i s Pdc=%1 . 2 f w a t t s \ nThe
o u t p u t power i s Pac=%1 . 1 f w a t t s \ n E f f i c i e n c y i s n
=%2 . 2 f p e r c e n t a g e ’ ,Pdc , Pac , n ) ;
20
21 //=========================END OF PROGRAM
===============================
Scilab code Exa 8.16 repeller voltage and etc
81
2 //
==================================================================
3 clc ;
4 clear ;
5
9 Rsh =20000; // Shunt r e s i s t a n c e i n ohms
10 L =0.001; // d i s t a n c e i n m
11 J1 =0.582; // b e s s e l c o e f f i c i e n t v a l u e [ J I (X ’ ) ]
12 f =8*10^(9) ; // // O p e r a t i o n f r e q u e n c y i n Hz
13 V1 =200; //RF gap v o l t a g e i n V
14 x =1.759*10^11; // e /m v a l u e i n C/ kg
15
16 //CALCULATION
17 disp ( ’ Assume t h e gap t r a n s i t t i m e and beam l o a d i n g
a r e n e g l e c t e d ’ );
19 VR =( V0 +(( sqrt (8* V0 / x ) * w * L ) /((2*( %pi ) * n ) -(( %pi ) /2) ) ) )
; // R e p e l l e r v o l t a g e i n V
20 disp ( ’ Assuming o u t p u t c o u p l i n g c o e f f i c i e n t Bo=1 ’ ) ;
21 I0 =( V1 /(2* J1 * Rsh ) ) /10^( -3) ; //Beam c u r r e n t n e c e s s a r y
t o o b t a i n an RF gap v o l t a f e o f 200V i n mA
22
23 //OUTPUT
24 mprintf ( ’ \ nThe R e p e l l e r v o l t a g e i s VR=%3 . 2 f V \nBeam
c u r r e n t n e c e s s a r y t o o b t a i n an RF gap v o l t a f e o f
200V i s I 0=%1 . 2 f mA ’ ,VR , I0 ) ;
25
26 //=========================END OF PROGRAM
===============================
82
Chapter 9
Solid State Microwave Devices
Scilab code Exa 9.1 frequency
2 //
==================================================================
3 clc ;
4 clear ;
5
6 L =2*10^( -6) ; // D r i f t L e n g th o f a IMPATT d i o d e i n m
7 Vd =(10^7) *(10^( -2) ) ; // D r i f t V e l o c i t y f o r S i i n m/ s e c
8
9 //CALCULATION
10 f =( Vd /(2* L ) ) /10^9; // O p e r a t i n g F r e q u e n c y i n GHz
11
12 //OUTPUT
13 mprintf ( ’ \ n O p e r a t i n g F r e q u e n c y o f t h e IMPATT d i o d e
i s f=%2 . 0 f GHz ’ ,f ) ;
14
15 //=========================END OF PROGRAM
===============================
83
Scilab code Exa 9.2 threshold electric field
2 //
==================================================================
3 clc ;
4 clear ;
5
6 L =75*10^( -6) ; // D e v i c e L e n g th i n m
7 V =25; // V o l t a g e P u l s e A m p l i f i e d i n V
8 f =10*10^9; // O p e r a t i n g F r e q u e n c y i n Hz
9
10 //CALCULATION
11 Eth =( V / L ) /10^5; // T h r e s h o l d E l e c t r i c F i e l d i n kV/cm
12
13 //OUTPUT
14 mprintf ( ’ \ n T h r e s h o l d E l e c t r i c F i e l d i s Eth=%1 . 2 f kV/
cm ’ , Eth ) ;
15
16 //=========================END OF PROGRAM
===============================
Scilab code Exa 9.3 power gain
2 //
==================================================================
3 clc ;
4 clear ;
5
84
6 fs =2*10^9; // S i g n a l F r e q u e n c y i n Hz
7 fp =12*10^9 //Pump F r e q u e n c y i n Hz
8 Ri =16; //O/P r e s i s t a n c e o f s i g n a l g e n e r a t o r i n ohms
9 Rs =1000; //On t y p e s r e s i s t a n c e o f s i g n a l g e n e r a t o r i n
ohms
10
11 //CALCULATION
12 P =10* log10 (( fp - fs ) / fs ) ; // Power g a i n i n dB
13 Pusb =10* log10 (( fp + fs ) / fs ) ; // Power g a i n a s USB
c o n v e r t e r i n dB
14
15 //OUTPUT
16 mprintf ( ’ \ nPower g a i n i s P=%1 . 2 f dB \ nPower g a i n a s
USB c o n v e r t e r i s Pusb=%1 . 2 f dB ’ ,P , Pusb ) ;
17
18 //=========================END OF PROGRAM
===============================
19
20
21 // Note : Answer g i v e n i n t e x t b o o k i s wrong Check i t
once . .
22 // C o r r e c t a n s w e r s a r e Power g a i n i s P=6.99 dB
23 // Power g a i n a s USB c o n v e r t e r i s
Pusb =8.45 dB
Scilab code Exa 9.4 breakdown voltage and etc
2 //
==================================================================
3 clc ;
4 clear ;
5
6 es =12.5; // R e l a t i v e D i e l e c t r i c constant
85
7 e0 =8.854*10^( -12) ; // P e r m i t t i v i t y i n F r e e S p a c e i n F/
m
8 N =3.2*10^22; // Donor C o n c e n t r a t i o n p e r mˆ3
9 L =8*10^( -6) ; // Le n g th o f S i BARITT d i o d e i n m
10 q =1.6*10^( -19) ; // Charge o f an E l e c t r o n i n C
11
12 //CALCULATION
13 Vc =(( q * N * L ^2) /(2* es * e0 ) ) /10^3; // C r i t i c a l V o l t a g e i n
kV
14 Vbd =2* Vc ; // Breakdown V o l t a g e i n kV
15 Ebd =( Vbd / L ) /100; // Breakdown E l e c t r i c F i e l d i n kV/cm
16
17 //OUTPUT
18 mprintf ( ’ \ n C r i t i c a l V o l t a g e i s Vc=%1 . 2 f kV \
nBreakdown V o l t a g e i s Vbd=%1 . 2 f kV \ nBreakdown
E l e c t r i c F i e l d i s Ebd=%6 . 2 f kV/cm ’ ,Vc , Vbd , Ebd ) ;
19
20 //=========================END OF PROGRAM
===============================
Scilab code Exa 9.5 Avalanche zone velocity
2 //
==================================================================
3 clc ;
4 clear ;
5
6 J =33000; // C u r r e n t d e n s i t y i n A/ sqcm
7 Na =2.5*10^16; // Doping C o n c e n t a t i o n i n TRAPATT d i o d e
p e r c u b i c cm
9
10 //CALCULATION
86
11 Vz =( J /( q * Na ) ) /10^5; // A v a l a n c h e Zone V e l o c i t y i n Km/
sec
12
13 //OUTPUT
14 mprintf ( ’ \ n A v a l a n c h e Zone V e l o c i t y i s Vz=%2 . 1 f Km/
s e c ’ , Vz ) ;
15
16 //=========================END OF PROGRAM
===============================
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For an IMPATT d i o d e power a m p l i f i e r
7 Rd =25; // N e g a t i v e R e s i s t a n c e i n ohms
8 Rl =50; // Load R e s i s t a n c e i n ohms
9
10 //CALCULATION
11 x = abs ( Rd ) ;
12 G =(( - x - Rl ) /( - x + Rl ) ) ^2; // Power g a i n o f an IMPATT
diode
13
14 //OUTPUT
15 mprintf ( ’ \ nPower g a i n o f an IMPATT d i o d e i s G=%1 . 0 f ’
,G ) ;
16
17 //=========================END OF PROGRAM
===============================
87
Scilab code Exa 9.7 minimum voltage
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a Gunn Dio d e
7 L =5*10^( -4) ; // D r i f t L e n g th i n cm
8 Vg =3300; // V o l t a g e g r a d i e n t i n V/cm [ Vg>3.3 kV/cm ]
9
10 //CALCULATION
11 Vmin = Vg * L ; // Minimum V o l t a g e n e e d e d t o i n i t i a t e Gunn
e f f e c t in volts
12
13 //OUTPUT
14 mprintf ( ’ \nMinimum V o l t a g e n e e d e d t o i n i t i a t e Gunn
e f f e c t i s Vmin=%1 . 2 f v o l t s ’ , Vmin ) ;
15
16 //=========================END OF PROGRAM
===============================
Scilab code Exa 9.8 rational frequency
2 //
==================================================================
3 clc ;
4 clear ;
88
5
6 // For a Gunn Dio d e
7 L =20*10^( -4) ; // A c t i v e L e n g th i n cm
8 Vd =2*10^7; // D r i f t V e l o c i t y o f E l e c t r o n s i n cm/ s e c
9 Ec =3.3*10^3; // C r i t i c l F i e l d f o r GaAs i n V/cm
10
11 //CALCULATION
12 fn =( Vd / L ) /10^9; // N a t u r a l ( R a t i o n a l ) F r e q u e n c y i n GHz
13 Vc = L * Ec ; // C r i t i c a l V o l t a g e o f t h e d i o d e i n v o l t s
14
15 //OUTPUT
16 mprintf ( ’ \ n N a t u r a l ( R a t i o n a l ) F r e q u e n c y i s f n=%2 . 0 f
GHz \ n C r i t i c a l V o l t a g e o f t h e d i o d e i s Vc=%1 . 1 f
v o l t s ’ ,fn , Vc ) ;
17
18 //=========================END OF PROGRAM
===============================
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For an IMPATT d i o d e
7 Lp =0.5*10^( -9) ; // I n d u c t a n c e i n Henry
8 Cj =0.5*10^( -12) ; // C a p a c i t a n c e i n Farad
9 Ip =0.8; //RF peak c u r r e n t i n A
10 Rl =2; // Load R e s i s t a n c e i n ohms
11 Vbd =100; // Breakdown V o l t a g e i n V
12 Ib =0.1; // dc B i a s c u r r e n t i n A
13
89
14 //CALCULATION
15 f =(1/(2*( %pi ) * sqrt ( Lp * Cj ) ) ) /10^9; // R e s o n a n t
F r e q u e n c y i n GHz
16 n =(( Rl * Ip ^2) /(2* Vbd * Ib ) ) *100; // E f f i c i e n c y i n
Percentage
17
18 //OUTPUT
19 mprintf ( ’ \ n R e s o n a n t F r e q u e n c y i s f=%2 . 0 f GHz \
n E f f i c i e n c y i s n=%1 . 1 f p e r c e n t a g e ’ ,f , n ) ;
20
21 //=========================END OF PROGRAM
===============================
Scilab code Exa 9.10 drift time
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For an IMPATT d i o d e
7 L =2*10^( -6) ; // D r i f t L e n g th i n m
8 Vd =10^5; // C a r r i e r D r i f t V e l o c i t y ( Assume / C o n s i d e r )
i n m/ s e c
9
10 //CALCULATION
11 t =( L / Vd ) /10^( -9) ; // D r i f t Time o f t h e C a r r i e r i n nano
s e c [ From f =(1/2 t ) =(Vd/2L ) ]
12 t1 = t *10^( -9) ;
13 f =(1/(2* t1 ) ) /10^9; // O p e r a t i n g F r e q u e n c y o f t h e d i o d e
i n GHz
14
15 //OUTPUT
90
16 mprintf ( ’ \ n D r i f t Time o f t h e C a r r i e r i s t=%1 . 2 f nano
s e c \ n O p e r a t i n g F r e q u e n c y o f t h e d i o d e i s f=%2 . 0
f GHz ’ ,t , f ) ;
17
18 //=========================END OF PROGRAM
===============================
Scilab code Exa 9.11 breakdown voltage and etc
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For an M−S i −M B a s i t t d i o d e
7 er =11.8; // R e l a t i v e d i e l e c t r i c c o n s t a n t o f S i
8 e0 =8.854*10^( -12) ; // P e r m i t t i v i t y o f f r e e s p a c e i n F/m
9 N =3*10^(21) ; // Donor C o n c e n t r a t i o n p e r mˆ3
10 L =6.2*10^( -6) ; // S i L e n g th i n m
12
13 //CALCULATION
14 Vbd =(( q * N * L ^2) /( er * e0 ) ) ; // Breakdown V o l t a g e i n V
15 Ebd =( Vbd / L ) /10^3; // Breakdown E l e c t r i c F i e l d i n kV/m
16
17 //OUTPUT
18 mprintf ( ’ \ nBreakdown V o l t a g e i s Vbd=%3 . 1 f V \
nBreakdown E l e c t r i c F i e l d i s Ebd=%5 . 0 f kV/m ’ ,Vbd ,
Ebd ) ;
19
20 //=========================END OF PROGRAM
===============================
91
Scilab code Exa 9.12 max power gain and etc
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For an u p c o n v e r t e r p a r a m e t r i c a m p l i f i e r
7 rQ =8; // f i g u r e o f m e r i t f o r a d i o d e n o n l i n e a r
capacitor
8 r =0.2;
9 y =8; // r a t i o
of output frequency over s i g n a l
frequency ( f0 / fs )
10 Td =300; // Dio d e T e m p e r a t u r e i n K
11 T0 =300; // Ambient T e m p e r a t u r e i n K
12
13 //CALCULATION
14 X =(( rQ ) ^2) / y ;
15 G =(( y * X ) /(1+ sqrt (1+ X ) ) ^2) ; //Max power g a i n
16 GdB =10* log10 ( G ) ; //Maximum Power Gain i n dB
17 F =(1+((2* Td / T0 ) *((1/ rQ ) +(1/ rQ ) ^2) ) ) ; // N o i s e F i g u r e
18 FdB =10* log10 ( F ) ; // N o i s e F i g u r e i n dB
19 BW =2* r * sqrt ( y ) ; // Bandwidth
20
21 //OUTPUT
22 mprintf ( ’ \nMaximum Power Gain i s GdB=%1 . 2 f dB\ n N o i s e
F i g u r e i s FdB=%1 . 2 f dB \ nBandWidth i s BW=%1 . 2 f ’ ,
GdB , FdB , BW ) ;
23
24 //=========================END OF PROGRAM
===============================
92
Scilab code Exa 9.13 gain and etc
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a n e g a t i v e r e s i s t a n c e p a r a m e t r i c a m p l i f i e r
7 fs =2*10^9; // S i g n a l F r e q u e n c y i n Hz
8 fp =12*10^9; //pump F r e q u e n c y i n Hz
9 fi =10*10^9; // i d l e r F r e q u e n c y i n Hz
10 fd =5*10^9; // F r e q u e n c y i n Hz
11 Ri =1000; // o / p r e s i s t a n c e o f i d l e r g e n e r a t o r i n ohms
12 Rg =1000; // o / p r e s i s t a n c e o f s i g n a l g e n e r a t o r i n ohms
13 RTs =1000; // t o t a l s e r i e s r e s i s t a n c e a t f s i n ohms
14 RTi =1000; // t o t a l s e r i e s r e s i s t a n c e a t f i i n ohms
15 r =0.35;
16 rQ =10; // f i g u r e o f m e r i t
17 Td =300; // Avg Diod e T e m p e r a t u r e i n K
18 T0 =300; // Ambient T e m p e r a t u r e i n K
19 C =0.01*10^( -12) ; // C a p a c i t a n c e i n F
20
21 //CALCULATION
22 ws =2*( %pi ) * fs ;
23 wi =2*( %pi ) * fi ;
24 R =(( r ^2) /( ws * wi * RTi * C ^2) ) ; // E q u i v a l e n t n o i s e
r e s i s t a n c e i n ohms
25 a =( R / RTs ) ;
26 G =((4* fi * a * Rg * Ri ) /( fs * RTs * RTi *(1 - a ) ^2) ) ; // Gain
27 GdB =10* log10 ( G ) ; // Gain i n dB
28 F =(1+((2* Td / T0 ) *((1/ rQ ) +(1/ rQ ) ^2) ) ) ; // N o i s e F i g u r e
29 FdB =10* log10 ( F ) ; // N o i s e F i g u r e i n dB
93
30 BW =(( r /2) *( sqrt ( fd /( fs * G ) ) ) ) ;
31
32 //OUTPUT
33 mprintf ( ’ \ n E q u i v a l e n t n o i s e r e s i s t a n c e
i s R=%4 . 1 f
ohms \ nGain i s GdB=%2 . 1 f dB\ n N o i s e F i g u r e i s FdB=
%1 . 2 f dB \ nBandWidth i s BW=%1 . 3 f ’ ,R , GdB , FdB , BW ) ;
34
35 //=========================END OF PROGRAM
===============================
36
37 // Note : Check t h e Bandwidth a n s w e r i s o n c e I t s h o u l d
be 0 . 0 2 7
94
Chapter 10
Microwave Communication
Systems
Scilab code Exa 10.1 radio horizon
1 // c h a p t e r −10 p a g e 486 e x a m p l e 1 0 . 1
2 //
==================================================================
3 clc ;
4 clear ;
5
6 ht =144; //TV t r a n s m i t t e r a n t e n n a h e i g h t i n m
7 hr =25; //TV r e c e i v e r a n t e n n a h e i g h t i n m
8 // Radio h o r i z o n i s a b o u t 4/3 a s f a r a s t h e o p t i c a l
horizon
9
10 //CALCULATION
11 dr =4* sqrt ( hr ) ; // d i s t a n c e i n km
12 dt =4* sqrt ( ht ) ; // Radio H o r i z o n i n km
13 d = dt + dr ; // The Maximum d i s t a n c e o f P r o p a g a t i o n o f t h e
TV s i g n a l i n km
14
15 //OUTPUT
95
16 mprintf ( ’ \ nThe Maximum d i s t a n c e o f P r o p a g a t i o n o f
t h e TV s i g n a l i s d=%2 . 0 f km \ nRadio H o r i z o n i s d t
=%2 . 0 f km ’ ,d , dt ) ;
17
18 //=========================END OF PROGRAM
===============================
Scilab code Exa 10.2 value of factor
2 //
==================================================================
3 clc ;
4 clear ;
5
6 r =6370*10^3; // r a d i u s o f t h e e a r t h i n m
7 x = -0.05*10^( -6) ; // t h e g r a d i e n t o f r e f r a c t i v e index
o f a i r n e a r t h e g r o u n d p e r m [ du / dh ]
8
9 //CALCULATION
10 k =1/(1+( r * x ) ) ; // The v a l u e o f t h e f a c t o r by which t h e
h o r i z o n d i s t a n c e o f a t r a n s m i t t e r w i l l be
modified
11
12 //OUTPUT
13 mprintf ( ’ \ nThe v a l u e o f t h e f a c t o r by which t h e
h o r i z o n d i s t a n c e o f a t r a n s m i t t e r w i l l be
m o d i f i e d i s k=%1 . 4 f ’ ,k ) ;
14
15 //=========================END OF PROGRAM
===================================
96
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a m i c r o w a v e LOS l i n k
7 f =2*10^9; // f r e q u e n c y o f o p e r a t i o n i n Hz
9 r =50000; // r e p e a t e r s p a c i n g i n m
10 PrdBm = -20; // r e q u i r e d c a r r i e r power a t t h e r e c e i v e r i
/ p t o a v o i d d e t e r i o r a t i o n due t o f a d i n g and n o i s e
i n dBm
11 GtdB =34; // a n t e n n a g a i n o f t r a n s m i t t e r i n dB
12 GrdB =34; // a n t e n n a g a i n o f r e c e i v e r i n dB
13 LdB =10; // c o u p l i n g and w a v e g u i d e l o s s i n t r a n s m i t t e r
i n dB
14
15 //CALULATION
16 w = c / f ; // w a v e l e n g t h i n m
17 x =( w ^2) /(4*( %pi ) ) ;
18 y =(4*( %pi ) * r ^2) ;
19 PtdBm = PrdBm +(10* log10 ( y ) ) - GtdB -(10* log10 ( x ) ) + LdB -
GrdB ; // The r e q u i r e d C a r r i e r T r a n s m i t t e r power i n
dBm
20
21 //OUTPUT
22 mprintf ( ’ \ nThe r e q u i r e d C a r r i e r T r a n s m i t t e r power i s
PtdBm=%2 . 1 f dBm ’ , PtdBm ) ;
23
24 //=========================END OF PROGRAM
===================================
97
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a g e o s t a t i o n a r y c o m m u n i c a t i o n s a t e l l i t e
7 f =6*10^(9) ; // u p l i n k f r e q u e n c y i n Hz
8 Pt =1000; // T r a n s m i t t e r power i n W
9 x =36000*10^3; // v e r t i c a l d i s t a n c e b e t w e e n s u r f a c e o f
e a r t h and s a t e l l i t e i n m
10 a =5; // a n t e n n a e l e v a t i o n a n g l e i n deg
11 GtdB =60; // a n t e n n a g a i n o f t r a n s m i t t e r i n dB
12 GrdB =0; // a n t e n n a g a i n o f r e c e i v e r i n dB
14
15 //CALCULATION
16 Gt =10^( GtdB /10) ; // a n t e n n a g a i n o f t r a n s m i t t e r
17 Gr =10^( GrdB /10) ; // a n t e n n a g a i n o f r e c e i v e r
19 Ar =( w ^2) *( Gr /(4*( %pi ) ) ) ; // a r e a i n sqm
20 r = x /( sind ( a ) ) ; // d i s t a n c e b e t w e e n t r a n s m i t t e r and
r e c e i v e r i n m [ From S i n e f o r m u l a and d i a g r a m ]
21 Pr =(( Pt * Gt * Ar ) /(4*( %pi ) * r ^2) ) /10^( -12) ; // The
r e c e i v e d power a t t h e i n p u t o f t h e s a t e l l i t e
r e c e i v e r in pico watts
22
23 //OUTPUT
24 mprintf ( ’ \ nThe r e c e i v e d power a t t h e i n p u t o f t h e
s a t e l l i t e r e c e i v e r i s Pr=%1 . 2 f p i c o w a t t s (pW) ’ , Pr
);
98
25
26 //=========================END OF PROGRAM
===============================
Scilab code Exa 10.5 antenna beam angle
2 //
==================================================================
3 clc ;
4 clear ;
5
6 x1 =35855; // D i s t a n c e b e t w e e n g e o s t a t i o n a r y o r b i t to
s u r f a c e o f e a r t h i n km
7 x2 =6371; // D i s t a n c e b e t w e e n c e n t e r o f e a r t h t o
8
9 //CALCULATION
10 x = x1 + x2 ; // d i s t a n c e o f s a t e l l i t e from c e n t e r o f e a r t h
i n km
11 y = x2 *( %pi ) ; // C i r c u m f e r e n c e o f h a l f c i r c l e a r c i n km
12 b = y / x ; //Beam a n g l e i n r a d
13 Bdeg =( b *180) /( %pi ) ; //Beam a n g l e i n deg
14
15 //OUTPUT
16 mprintf ( ’ \ nAntenna Beam a n g l e r e q u i r e d by a
s a t e l l i t e antenna to provide f u l l g l o b a l coverage
from a g e o s t a t i o n a r y o r b i t i s Bdeg=%2 . 2 f deg ’ ,
Bdeg ) ;
17
18 //=========================END OF PROGRAM
===============================
99
Scilab code Exa 10.6 round trip time
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a s a t e l l i t e c o m m u n i c a t i o n s y s t e m
7 h =35855; // D i s t a n c e b e t w e e n g e o s t a t i o n a r y o r b i t to
8 r =6371; // D i s t a n c e b e t w e e n c e n t e r of earth to s u r f a c e
o f e a r t h i n km
9 a =5; // e a r t h s t a t i o n e l e v a t i o n a n g l e wrt t h e
g e o s t a t i o n a r y s a t e l l i t e i n deg
10 b =5; // a n g l e i n deg
11 c =3*10^5; // V e l o c i t y o f l i g h t i n km/ s e c
12 b1 =90; // a n g l e f o r v e r t i c a l t r a n s m i s s i o n i n deg
13 a1 =0;
14
15 //CALCULATION
16 d =( sqrt (( r + h ) ^2 -( r * cosd ( a ) ) ^2) ) - sind ( b ) ; // d i s t a n c e
i n km
17 T =2*( d / c ) ; // The round t r i p t i m e b e t w e e n t h e e a r t h
s t a t i o n and t h e s a t e l l i t e i n s e c
18 d1 =( sqrt (( r + h ) ^2 -( r * cosd ( a ) ) ^2) ) - sind ( b ) ; // d i s t a n c e
i n km
19 Tv =(2/ c ) *( d1 - r ) ; // The round t r i p t i m e f o r v e r t i c a l
t r a n s m i s s i o n b e t w e e n t h e e a r t h s t a t i o n and t h e
s a t e l l i t e in sec
20
21 //OUTPUT
22 mprintf ( ’ \ nThe round t r i p t i m e b e t w e e n t h e e a r t h
100
s t a t i o n and t h e s a t e l l i t e i s T=%1 . 3 f s e c \ nThe
round t r i p t i m e f o r v e r t i c a l t r a n s m i s s i o n b e t w e e n
t h e e a r t h s t a t i o n and t h e s a t e l l i t e i s Tv=%1 . 3 f
s e c ’ ,T , Tv ) ;
23
24 //=========================END OF PROGRAM
===============================
Scilab code Exa 10.7 figure of merit
2 //
==================================================================
3 clc ;
4 clear ;
5
6 Tant =25; // e f f e c t i v e n o i s e t e m p e r a t u r e i n K
7 Tr =75; // r e c e i v e r n o i s e t e m p e r a t u r e i n K
8 GdB =45; // I s o t r o p i c power g a i n o f t h e a n t e n n a i n dB
9
10 //CALCULATION
11 T = Tant + Tr ; // The t o t a l n o i s e i n K
12 TdB =10* log10 ( T ) ; // The t o t a l n o i s e i n dB
13 MdB = GdB - TdB ; // F i g u r e o f m e r i t o f e a r t h s t a t i o n i n dB
14
15 //OUTPUT
16 mprintf ( ’ \ n F i g u r e o f m e r i t o f e a r t h s t a t i o n i s MdB=
%2 . 0 f dB ’ , MdB ) ;
17
18 //=========================END OF PROGRAM
===============================
101
Scilab code Exa 10.8 CNR
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a S a t e l l i t e c o m m u n i c a t i o n l i n k
7 EIRPdB =55.5; // S a t e l l i t e ESM i n dBW
8 MdB =35; //G/T r a t i o o f e a r t h s t a t i o n i n dB
9 LfsdB =245.3 // F r e e s p a c e l o s s i n dB
10
11 //CALCULATION
12 CNRdB = EIRPdB + MdB - LfsdB +228.6; // C a r r i e r t o N o i s e
R a t i o a t t h e e a r t h s t a t i o n r e c e i v e r i n dB
13
14 //OUTPUT
15 mprintf ( ’ \ n C a r r i e r t o N o i s e R a t i o a t t h e e a r t h
s t a t i o n r e c e i v e r i s CNRdB=%2 . 1 f dB ’ , CNRdB ) ;
16
17 //=========================END OF PROGRAM
===============================
Scilab code Exa 10.9 system noise temperature
2 //
==================================================================
3 clc ;
4 clear ;
5
6 D =30; // D i a m e t e r o f a d i s h a n t e n n a w i t h circular
102
aperture in m
7 f =4*10^9; // down l i n k f r e q u e n c y i n Hz
8 MdB =20; //G/T r a t i o o f e a r t h s t a t i o n i n dB
10
11 //CALCULATION
12 A =(( %pi ) /4) * D ^2; // a r e a i n sqm
14 G =(4*( %pi ) * A ) / w ^2; // Gain
15 GdB =10* log10 ( G ) ; // Gain i n dB
16 TsdB = GdB - MdB ; // The System N o i s e T e m p e r a t u r e i n dB
17
18 //OUTPUT
19 mprintf ( ’ \ nThe System N o i s e T e m p e r a t u r e i s TsdB=%2 . 2
f dB ’ , TsdB ) ;
20
21 //=========================END OF PROGRAM
===============================
Scilab code Exa 10.10 HPBW
1 // c h a p t e r −10 p a g e 489 e x a m p l e 1 0 . 1 0
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a p a r a b o l i c a n t e n n a
7 Gp =1500; // Power g a i n
8 w =0.1; // w a v e l e n g t h i n m
9
10 //CALCULATION
11 D = sqrt ( Gp ) *( w /( %pi ) ) ; // D i a m e t e r o f t h e c i r c u l a r
mouth o f a p a r a b o l i c a n t e n n a i n m
103
12 HPBW =58*( w / D ) ; // H a l f Power BeamWidth o f t h e a n t e n n a
i n deg
13
14 //OUTPUT
15 mprintf ( ’ \ n D i a m e t e r o f t h e c i r c u l a r mouth o f a
p a r a b o l i c a n t e n n a i s D=%1 . 4 f m \ n H a l f Power
BeamWidth o f t h e a n t e n n a i s HPBW=%1 . 3 f deg ’ ,D ,
HPBW ) ;
16
17 //=========================END OF PROGRAM
===============================
Scilab code Exa 10.11 gain
2 //
==================================================================
3 clc ;
4 clear ;
5
6 D =1; // Assume d i a m e t e r o f t h e p a r a b o l i c r e f l e c t o r s in
the o r i g i n a l system in m
7 w =1; // Assume w a v e l e n g t h i n m
8
9 //CALCULATION
10 D1 =2* D ; // d i a m e t e r o f t h e p a r a b o l i c r e f l e c t o r s in the
modified system in m
11 G =6*( D / w ) ^2; // g a i n i n o r i g i n a l s y s t e m
12 G1 =6*( D1 / w ) ^2; // g a i n i n m o d i f i e d s y s t e m
13 GdB =10* log10 ( G1 / G ) ; // O v e r a l l g a i n t h a t can be
e x p e c t e d i n dB
14 GdBo =2* GdB ; // O v e r a l l g a i n o f t h e s y s t e m ( c o m b i n i n g
t h e two a n t e n n a s one a t t h e Tx and o t h e r a t t h e
Rx ) i n dB
104
15
16 //OUTPUT
17 mprintf ( ’ \ n O v e r a l l g a i n t h a t can be e x p e c t e d i s GdB=
%1 . 0 f dB \ n O v e r a l l g a i n o f t h e s y s t e m ( c o m b i n i n g
t h e two a n t e n n a s one a t t h e Tx and o t h e r a t t h e
Rx ) i s GdBo=%1 . 0 f dB ’ ,GdB , GdBo ) ;
18
19 //=========================END OF PROGRAM
===============================
20
21 // Note : Check t h e a n s w e r o n c e . . i t s h o u l d be GdB=10
l o g ( 4 ) =6 dB and GdBo=12dB
Scilab code Exa 10.12 gain
2 //
==================================================================
3 clc ;
4 clear ;
5
6 D =3; // d i m e n s i o n o f a p a r a b o l o i d i n m
7 f =3*10^9; // f r e q u e n c y ( S band ) i n Hz
9
10 //CALCULATION
11 w = c / f ; // wave l e n g t h i n m
12 BWFN =140*( w / D ) ; // BeamWidth b e t w e e n F i r s t N u l l s i n
deg
13 BWHP =70*( w / D ) ; // BeamWidth b e t w e e n H a l f P o w e r p o i n t s
i n deg
14 G =6*( D / w ) ^2; // Gain o f t h e a n t e n n a
15
16 //OUTPUT
105
17 mprintf ( ’ \nBeamWidth b e t w e e n F i r s t N u l l s i s BWFN=%1
. 2 f deg \nBeamWidth b e t w e e n H a l f P o w e r p o i n t s i s
BWHP=%1 . 2 f deg \ nGain o f t h e Antenna i s G=%4 . 0 f ’
, BWFN , BWHP , G ) ;
18
19 //=========================END OF PROGRAM
===============================
2 //
==================================================================
3 clc ;
4 clear ;
5
6 l =1; // ( Assume )−d i m e n s i o n ( w a v e l e n g t h ) i n cm
7
8 //CALCULATION
9 x =5* l ; // g i v e n s q u a r e a p e r t u r e o f an optimum h o r n
a n t e n n a a s a s i d e d i m e n s i o n i n cm
10 A = x * x ; // Area i n s q . cm
11 Gp =4.5*( A / l ^2) ; // Power g a i n o f an optimum h o r n
antenna
12
13 //OUTPUT
14 mprintf ( ’ \ nPower g a i n o f an optimum h o r n a n t e n n a i s
Gp=%3 . 1 f ’ , Gp ) ;
15
16 //=========================END OF PROGRAM
===============================
106
Chapter 11
Radars
Scilab code Exa 11.1 max range
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a r a d a r s y s t e m
7 Pt =600000; // peak p u l s e power i n W
8 Smin =10^( -13) ; // minimum d e t e c t a b l e s i g n a l i n W
9 Ae =5; // c r o s s s e c t i o n a l a r e a o f t h e r a d a r a n t e n n a i n
sq m
10 w =0.03; // w a v e l e n g t h i n m
11 s =20; // r a d a r c r o s s s e c t i o n a l a r e a i n s q m
12
13 //CALCULATION
14 Rmax =((( Pt * s * Ae ^2) /(4*( %pi ) * Smin * w ^2) ) ^(1/4) ) /1000;
//Maximum r a n g e o f a r a d a r s y s t e m i n km
15 RMax = Rmax /1.853; // I n n a u t i c a l m i l e s ; 1 nm= 1 . 8 5 3 km
16
17 //OUTPUT
107
18 mprintf ( ’ \nMaximum r a n g e o f a r a d a r s y s t e m i s Rmax=
%3 . 3 f km ’ , Rmax ) ;
19 disp ( ’ I n n a u t i c a l m i l e s ; 1 nm= 1 . 8 5 3 km ’ ) ;
20 mprintf ( ’ \nMaximum r a n g e o f a r a d a r s y s t e m i n
n a u t i c a l m i l e s i s RMax=%3 . 0 f nm ’ , RMax ) ;
21
22 //=========================END OF PROGRAM
===============================
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a r a d a r s y s t e m
7 Pt =250000; // peak t r a n s m i t t e d power i n W
8 G =2500; // power g a i n o f t h e a n t e n n a
9 Smin =10^( -14) ; // minimum d e t e c t a b l e s i g n a l i n W
10 Ae =10; // c r o s s s e c t i o n a l a r e a o f t h e r a d a r a n t e n n a i n
sq m
11 f =10*10^9; // f r e q u e n c y o f r a d a r i n Hz
14
15 //CALCULATION
17 Rmax =((( Pt * G * Ae * s ) /( Smin *(4*( %pi ) ) ^2) ) ^(1/4) ) /1000;
//Maximum r a n g e o f a r a d a r s y s t e m i n km
18
19 //OUTPUT
20 mprintf ( ’ \nMaximum r a n g e o f a r a d a r s y s t e m i s Rmax=
108
%3 . 2 f km ’ , Rmax ) ;
21
22 //=========================END OF PROGRAM
===============================
Scilab code Exa 11.3 RCS
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a m a r i n e r a d a r s y s t e m
7 Pt =250000; // peak t r a n s m i t t e d power i n W
8 G =4000; // power g a i n o f t h e a n t e n n a
9 R =50000; //maximum r a n g e o f r a d a r i n m
10 Pr =10^( -11) ; // minimum d e t e c t a b l e s i g n a l i n W
11 f =10*10^9; // f r e q u e n c y o f r a d a r i n H
13
14 //CALCULATION
16 Ae =(( G * w ^2) /(4*( %pi ) ) ) ; // c r o s s s e c t i o n a l a r e a o f t h e
radar antenna in sq m
17 s =(( Pr *(4*( %pi ) * R ^2) ^2) /( Pt * G * Ae ) ) ; // The c r o s s
s e c t i o n o f t h e t a r g e t t h e r a d a r can s i g h t i n s q m
18
19 //OUTPUT
20 mprintf ( ’ \ nThe c r o s s s e c t i o n o f t h e t a r g e t t h e r a d a r
can s i g h t i s s=%2 . 2 f s q m ’ ,s ) ;
21
22 //=========================END OF PROGRAM
===============================
109
Scilab code Exa 11.4 Duty cycle and etc
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a g u i d e d m i s s i l e t r a c k i n g r a d a r
7 Pt =400000; // t r a n s m i t t e d power i n W
8 prf =1500; // p u l s e r e p i t i t i o n f r e q u e n c y i n pps ( p u l s e
per sec )
9 tw =0.8*10^( -6) ; // p u l s e w i d t h i n s e c
11
12 //CALCULATION
13 Runamb =( c /(2* prf ) ) /1000; // Unambiguous r a n g e i n km
14 dc = tw /(1/ prf ) ; // Duty c y c l e
15 Pav = Pt * dc ; // A v e r a g e power i n W
16 n1 =1;
17 BW1 =( n1 / tw ) /10^6; // S u i t a b l e BW i n MHz f o r n=1
18 n2 =1.4;
19 BW2 =( n2 / tw ) /10^6; // S u i t a b l e BW i n MHz f o r n =1.4
20
21 //OUTPUT
22 mprintf ( ’ \ nUnambiguous r a n g e i s Runamb=%3 . 0 f km \
nDuty c y c l e i s dc=%1 . 4 f \ n A v e r a g e power i s Pav=%3
. 0 f W’ , Runamb , dc , Pav ) ;
23 disp ( ’ For e f f i c i e n c y n =1 , s u i t a b l e bandwidth i n MHz
i s ’ );
24 disp ( BW1 ) ;
25 disp ( ’ For e f f i c i e n c y n = 1 . 4 , s u i t a b l e bandwidth i n MHz
i s ’ );
110
26 disp ( BW2 ) ;
27
28 //=========================END OF PROGRAM
===============================
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a m i l i t a r y r a d a r
7 Pt =2500000; // power o u t p u t i n W
8 f =5*10^9; // f r e q u e n c y o f r a d a r i n H
10 D =5; // a n t e n n a d i a m e t e r i n m
11 B =1.6*10^6; // r e c e i v e r bandwidth i n Hz
13 NF =12; // n o i s e f i g u r e i n dB
14
15 //CALCULATION
17 F =10^( NF /10) ; // n o i s e f i g u r e
18 Rmax =(48*(( Pt * s * D ^4) /( B *( F -1) * w ^2) ) ^(1/4) ) ; //Maximum
d e t e c t i o n r a n g e i n km
19
20 //OUTPUT
21 mprintf ( ’ \nMaximum d e t e c t i o n r a n g e o f a r a d a r i s
Rmax=%3 . 0 f km ’ , Rmax ) ;
22
23 //=========================END OF PROGRAM
===============================
111
Scilab code Exa 11.6 factor
2 //
==================================================================
3 clc ;
4 clear ;
5
6 // For a c i v i l i a n r a d a r
7 Rmax =30; //maximum r a n g e i n kms
8 x =50;
9 y =2;
10 disp ( ’ Maximum r a n g e w i t h an e q u i v a l e n t echoing area
o f 50 t i m e s i n kms i s ’ ) ;
11 R = Rmax * x ^(1/4) ;
12 disp ( R ) ;
13 disp ( ’ Range would be i n c r e a s e d i f Tx power i s
d o u b l e d by a f a c t o r o f ’ ) ;
14 f = y ^(1/4) ;
15 disp ( f ) ;
16
17 //=========================END OF PROGRAM
===============================
112

Microwave Engineering M Kulkarni - Compress PDF

Uploaded by

Copyright:

Available Formats

Microwave Engineering M Kulkarni - Compress PDF

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 M Kulkarni - Compress PDF

Uploaded by

Copyright:

Available Formats

Scilab Textbook Companion for

Microwave and Radar Engineering

June 13, 2014

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

Title: Microwave and Radar Engineering

Publisher: Umesh Publications, 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

4 Microwave Transmission Lines 18

8 Microwave Tubes and Circuits 64

9 Solid State Microwave Devices 83

10 Microwave Communication Systems 95

Exa 3.1 Terminating Impedance . . . . . . . . . . . . . . . . . 7

Scilab code Exa 3.1 Terminating Impedance

Scilab code Exa 3.3 phase velocity

Scilab code Exa 3.5 VSWR and Reflection coefficient

Scilab code Exa 3.6 point of attachment

Scilab code Exa 3.7 Terminating Impedance

Scilab code Exa 3.8 VSWR and Impedance

Scilab code Exa 3.9 Reflection loss

Scilab code Exa 3.10 Characteristic Impedance

Microwave Transmission Lines

Scilab code Exa 4.1 Characteristic Impedance

Scilab code Exa 4.2 attenuation and phase constants

Scilab code Exa 4.3 breakdown power

Scilab code Exa 4.4 Characteristic Impedance

Scilab code Exa 4.5 Characteristic Impedance

Scilab code Exa 4.6 ratio of areas

Scilab code Exa 4.7 breadth of the guide

Scilab code Exa 4.8 guide wavelength

Scilab code Exa 4.9 guide wavelength

Scilab code Exa 4.10 area

Scilab code Exa 4.11 possible modes

Scilab code Exa 4.12 guide wavelength

Scilab code Exa 4.13 frequency

Scilab code Exa 4.14 guide wavelength

Scilab code Exa 4.15 possible modes

Scilab code Exa 4.16 Characteristic Impedance

Scilab code Exa 4.17 guide wavelength

Scilab code Exa 4.19 amount of attenuation

Scilab code Exa 4.20 max power handling capacity

Scilab code Exa 4.21 maximum power

Scilab code Exa 4.22 peak value of electric field

Scilab code Exa 4.23 breakdown power

Scilab code Exa 4.24 breakdown power

Scilab code Exa 5.1 minimum distance

Scilab code Exa 5.2 lowest resonant frequency

Scilab code Exa 5.3 resonant frequency

Scilab code Exa 5.4 resonant frequency

Scilab code Exa 6.2 distance to be shifted

1 // Chapter −6 , Example 6 . 2 , Page 234

Scilab code Exa 6.3 Scattering parameters

1 // Chapter −6 , Example 6 . 3 , Page 236

Scilab code Exa 6.4 powers in the remaining ports

1 // Chapter −6 , Example 6 . 4 , Page 238

1 // Chapter −6 , Example 6 . 5 , Page 239

Scilab code Exa 6.6 Reflected power