Lab - Solution Fundamentals of Electromagnetics With Engineering Applications by Stuart M. Wentworth (Solution
Lab - Solution Fundamentals of Electromagnetics With Engineering Applications by Stuart M. Wentworth (Solution
1. General Properties
P8.1: In free space, a wave propagating radially away from an antenna at the origin has
Is
Hs cos2 a ,
r
where the driving current phasor Is = Ioej. Determine (a) Es, (b) P(r,) and (c) Rrad.
� I �
E s h a P �H s ho ar �� s cos 2 a � ,
�r �
hI
(a )E s o s cos 2 a
r
1 1 ho I o e j
� I o e j �
P Re � �E s �H *
s�
� Re � cos 2
a � cos 2 a �
2 2 � r r �
2
1 �I �
(b) P ( r , , ) ho �o �cos 4 a r
2 �r �
Now to find Rrad:
1
Prad � P (r , , )gdS I o2 Rrad ,
2
p 2p
1 cos 4
1
Prad ho I o2 � 2 a r gr 2 sin d d a r ho I o2 � cos 4 sin d � d
2 r 2 0 0
pho I o2 2 1
( cos5 ) pho I o2 I o2 Rrad
p
Prad
5 0 5 2
Solving:
2
p ( 120p ) I o2
Rrad 5 96p 2 W
1 2
I
2 o
(c ) Rrad 950W
P8.2: What is the pattern solid angle and the directivity for an isotropic antenna? How
about for a semi-isotropic antenna, that radiates equally in all directions above = p/2,
but is zero otherwise?
p 2p
WP �
�Pn ( , ) d W �
(1)sin d d 4p
�
0 0
4p
W P 4p , and Dmax 1.
WP
For a semi-isotropic antenna:
8-2
p 2 2p
WP �(1)sin d �
d 2p , and Dmax 2.
0 0
P8.3: Sketch an appropriate cross-section of the radiation pattern and determine the
beamwidth, pattern solid angle and directivity for the following normalized radiation
intensities:
(a ) Pn ( , ) cos for 0 p 2, 0 otherwise.
(b) Pn ( , ) cos 2 for 0 p 2, 0 otherwise.
(c ) Pn ( , ) cos3 for 0 p 2, 0 otherwise.
% MLP0803
% Polar Plots for P8.03
%
% The polar plot function doesn't allow multiple
% plots. So we have to let the angle theta loop
% around several times, changing the rho function
% each time.
%
for i=1:100
theta(i)=-pi/2+i*pi/100;
rho(i)=cos(theta(i));
end
for j=101:200
theta(j)=(-pi/2)+j*pi/100;
rho(j)=0;
end
for i=201:300
theta(i)=-pi/2+i*pi/100;
rho(i)=(cos(theta(i)))^2;
end
for j=301:400
theta(j)=-pi/2+j*pi/100;
rho(j)=0;
end
for i=401:500
theta(i)=-pi/2+i*pi/100;
rho(i)=(cos(theta(i)))^3;
end Fig. P8.3
for j=501:600
8-3
theta(j)=-pi/2+j*pi/100;
rho(j)=0;
end
polar(theta,rho)
(a) Beamwidth: solving for cos = 0.5, we find = 60°, and beamwidth =2 = 120°.
WP �
�Pn ( , ) d W �
�cos sin d d 2p �cos sin d
0
u cos , du sin d , so
p 2
W P 2p �
udu p u 2 p cos 2 p sr.
0
4p
Dmax 4
WP
(b) Beamwidth: solving for cos2 = 0.5, we find = 45°, and beamwidth =2 = 90°.
WP �
�cos 2 sin d d 2p �cos 2
sin d
0
u cos , du sin d , so
p 2
2p 3 2p 2p
W P 2p �
u 2 du u cos3 sr
3 3 0 3
4p
Dmax 6
WP
(c) Beamwidth: solving for cos3 = 0.5, we find = 37.5°, and beamwidth =2 = 75°.
P8.4: Sketch an appropriate cross-section of the radiation pattern and determine the
beamwidth, pattern solid angle and directivity for the following normalized radiation
intensities:
(a ) Pn ( , ) sin
(b) Pn ( , ) sin 2
(c ) Pn ( , ) sin 3
1 �1 �
(a) Pn sin , sin 1 � � 30o.
2 �2 �
Now, since the beam is maximum at = 90°, the beamwidth is from 30° to 150°, or BW =
120°.
p 2p
WP �
�Pn d W �
�sin sin d d �
sin d �
d 2
0 0
p
1 4p
so W P 2p �( 1 cos 2 ) d p sr , Dmax 1.27
2
0
2 WP
1 �1� o
(b) Pn sin 2 , sin 1 �
� �
� 45 , BW 135o 45o 90 o
2 �2�
p
WP �
�Pn d W �
� ( 1 cos2 ) sin d
sin 2 sin d d 2p �
0
p p
W P 2p �
sin d 2p �
cos 2 sin d
0 0
p 2p p 4p 8p
2p ( cos ) 0 + cos3 4p sr
3 0 3 3
4p
so Dmax 1.5
WP
(c) Pn
1
2 ( )
sin 3 , sin 1 ( 1 2 ) 3 52.5o, BW ( 180o 52.5o ) 52.5o 75o
1
p p
1 p p� 1 �
sin 4 d 2p �( 1 cos 2 ) d � 1 2 cos 2 + ( 1 + cos 4 ) �
2
W P 2p � � d
0 0
4 2 0� 2 �
Plots are generated using MLP0804:
% MLP0804
%
% Generate polar plots
%
clc
clear
for i=1:100
theta(i)=i*pi/50;
Fig. P8.4
8-5
rho(i)=abs(sin(theta(i)));
end
for i=101:200
theta(i)=i*pi/50;
rho(i)=(sin(theta(i)))^2;
end
for i=201:300
theta(i)=i*pi/50;
rho(i)=abs(sin(theta(i))^3);
end
polar(theta,rho)
P8.5: (JustAsk): You are given the following normalized radiation intensity:
Pn ( , ) sin 2 sin 3 for 0 p ,
0 otherwise.
Find the beamwidth, pattern solid angle, and directivity.
1
The beam is pointing in the ay direction, and we have BW
2
( BW + BW ) .
To find BW, we fix = p/2 and set sin2 equal to ½. Then,
� 1�
sin 1 � � 45 , so BW ( 180 45 ) 45 90 .
o o o o o
� 2�
� �
To find BW, we fix = p/2, and set sin3 = ½, giving us
(
sin 1 ( 1 2 )
1
3
) 52.5 ,o
so BW ( 180o 52.5o ) 52.5o 75o.
1
Finally, BW
2
( 90o + 75o) 82.5o.
The pattern solid angle is
WP �� Pn d W � � ( sin 2 sin 3 ) sin d d ,
p p
sin d �
WP � sin 3 d , (note limits on )
3
0 0
p
1 p 1 2
so � cos 2 x sin xdx cos3 x (1 1) .
0
3 0 3 3
So we have
p
2 4
y� sin 3 xdx 2 ,
0
3 3
p p
�4 ��4 � 16
and W P � sin 3 d � sin 3 d � �� � 1.78sr.
0 0 �3 ��3 � 9
4p 4p
Dmax 7.1
W P 1.78
1
BW
2
( BW + BW ) ,
BW: Fix = p, sin2 = 1/2, = 45°, BW = (180° – 45°) – 45° =90°.
BW: Fix = p/2, sin(/2)=1/2, BW = (360° – 60°) – 60° = 240°
BW = (90°+240°)/2=165°
By inspection, the direction of maximum radiation is at = p and = p/2. (i.e. the –ax
direction).
2p
p 3
WP � �sin sin sin d d �
2
sin d �
sin d
2 0
2 0
Do each integral separately:
2p 2p
� sin d 2cos 2( 1 1) 4
0
2 20
p p p p
� ( 1 cos2 ) sin dx �
sin 3 d � sin d �
cos 2 sin d
0 0 0 0
p 1 p 1 2 4
cos 0 + cos3 (1 1) + (1 1) 2
3 0 3 3 3
�4 � 16 4p 3p
So W P (4) � � , and Dmax 2.4
�3 � 3 WP 4
BOS AOS
y
( o Io e j y ) ax j o Ioe j y ax .
BOS
HOS j I oe j y a x .
o
EOS hoa y HOS hoa y j I oe j y a x jho I oe j y a z .
1
P Re EOS H*OS , so
2
1 1
P Re jho I o e j y a z + j I o e + j y a x ho ( I o ) a y .
2
2 2
Note that this problem is not very realistic, but is good to teach the mechanics of going
from AOS to P.
P8.9: (JustAsk): Suppose a Hertzian dipole antenna is 1.0 cm long and is excited by a 10.
mA amplitude current source at 100. MHz. What is the maximum power density radiated
by this antenna at a 1.0 km distance? What is the antenna’s radiation resistance?
c 3 x108 m s
clf, l 3m.
f 100 x106 1 s
ho 2 I o2 l 2 120p ( 2p ) ( 0.010 ) ( 0.010 )
2 2 2
pW
Pmax 0.052 2
32p r 2 2
32 32
p 2
1000 2
m
2 2
�l � �0.01 �
Rrad 80p 2 � � 80p 2 � � 8.8mW
�l � � 3 �
8-9
P8.10: A 1.0 cm long, 1.0 mm diameter copper wire is used as a Hertzian dipole radiator
at 1.0 GHz. (a) Find Rrad. (b) Estimate Rdiss by considering the skin effect resistance of
the wire. (c) Find efficiency, e. (d) Find the maximum power gain Gmax.
c 3 x108 m s
l 0.3m
f 1x109 1 s
2 2
�l � �0.01 �
Rrad 80p � � 80p 2 �
2
� 0.877W
�l � � .3 �
From Example 8.2 we have d Cu 1GHz 2.09 x10 m
6
Rrad 0.877
e 0.97
Rrad + Rdiss 0.877 + 0.026
Gmax eDmax 0.97 ( 1.5 ) 1.46
P8.11: Evaluate the curl of Aos (equation (8.59)) to find Hos. Now apply a far-field
approximation to verify (8.60).
o I s S
2 (
Eqn. (8.59): A OS 1 + j r ) e j r sin a
4p r
Since AOS only has an A component, we have
1 � 1 �
� A OS
r sin �
( sin A ) a r
r �r
( rA ) a
A A jA j r
We’ll let A 2 ( 1 + j r ) e
j r
sin 2 e j r sin + e sin ,
r r r
I S
where A o S .
4p
We break up the derivative into two parts:
� � A j r jA j r ��sin � A j r jA j r �
2
1.
�
( sin A) � 2
�r
e +
r
e �
� �
� 2 e
�r
+
r
e �2sin cos
�
or
� 2A j 2 A j r
�
( sin A ) 2 e j r sin cos +
r r
e sin cos
1 � �2 A j 2 A �
r sin �
( sin A ) a r � 3 e j r cos + 2 e j r cos �
�r r �
ar
� ��A �
2.
�
r
( rA ) � e j r sin + jA e j r sin �
�r �r
,
�
8-10
P8.12: Neglecting resistive losses in the wire, how much current must drive a loop
antenna of radius 2.0 cm at 60 MHz to radiate 1.0 W of power? Repeat for a 20 turn
loop.
8-11
c 3 x108
At 60 MHz, l 5m . So at 2 cm radius, we have the small loop situation.
f 60 x106
2
4h p 3 I 2 �S �
For a small loop: Prad o o � 2 �,
3 �l �
Here,
S p a 2 p ( 0.02m ) 1.257 x103 m 2 .
2
ho I s e j r
EOS ho a r �HOS sin cos a
4p r
2
1 1 � I o � 2
P (r , , ) Re � EOS �H OS � *
� 2 ho �4p r �sin cos a r
2
2 � � �
Note also that P ( r , , ) P P
max n ( , )a r , where here
2
1 � I �
Pmax ho � o �, and Pn ( , ) sin 2 cos 2 .
2 �4p r �
Then, W P �
�Pn ( , )d W.
p
2 4
sin 3 d 2
Referring to P8.5, � , and
0
3 3
2p
1
�
cos d �
2
� d + � cos 2 d �
� p .
0
2�
4p
So, W P �� Pn ( , )d W sr
3
� � I o �
2 1
2
�
�4p �
Prad r Pmax W P r � ho � ��
2
� � .
�2 �4p r �� �3 �
� �
Using ho 120p and 2p l , we find
8-12
20p 2 I o2
Prad .
l2
1 2 2P 40p 2
Finally, Prad I o Rrad , so Rrad rad ,
2 I o2 l2
and since for this problem, l c f 3m, Rrad ( 40p 9 ) 44W
2
P8.14: Suppose in the far-field for a particular antenna at the origin, the electric field is
e j r
Eos ho I o sin a .
pr
What is the radiation resistance of this antenna?
1 2
We’ll use: Prad I o Rrad r 2 Pmax W P , so we must find Pmax and Wp.
2
1 e j r e j r
H S a r �ho Io sin a I o sin a
ho pr pr
1 1 1
P Re � EOS �H*OS �
� � ho I o2 sin 2 a r ,
( pr)
2
2 2
so
1 1
Pmax ho I o2 , and Pn sin 2 .
(pr)
2
2
p 2p p
WP �
� sin d �
Pn d W � 3
( 1 cos2 ) sin d
d 2p �
0 0 0
p p
� � � p 1 p � 8p
W P 2p �� sin d � cos 2 sin d � 2p � cos 0 + cos 2 � sr
�
0 0 � � 3 0
� 3
2 2 1 1 8p 8 (120p )8
Rrad 2 r 2 Pmax W P 2 r 2 ho I o2 ho 320W
(pr) 3 3p 3p
2
Io Io 2
P8.15: Derive the expressions for radiated power (equation (8.64)) and radiation
resistance (equation (8.65)) for a small loop antenna.
1 2
We’ll use: Prad I o Rrad r 2 Pmax W P
2
2 o 2 I o2 S 2 2
From (8.63) we have Pmax
32hop 2 r 2
and
p 2p
8p
WP ��sin 2
d W � sin 3
d � d sr (see integral solution of P8.14)
0 0
3
8-13
Now,
� 2 2 I 2 S 2 2 �
�8p �
Prad r 2 Pmax W P r 2 � o o 2 2 � � �
� 32hop r � �3 �
o
Using the conversions: o o , 2p l , and ho o
we arrive at:
2
4 3 2 �S � 1 2
Prad hop I o � 2 � I o Rrad
3 �l � 2
Solving for Rrad,
2
4 �S �
Rrad 320p � 2 �W
�l �
3. Dipole Antennas
P8.16: MATLAB: Develop a routine to calculate the beamwidth for a dipole antenna of
arbitrary length between 0.1land 1l.
% MLP0816
%
% Determine beamwidth for an arbirary length
% dipole antenna. Equation (8.74) is used.
%
% 2/10/03 Wentworth
%
% Variables
% L dipole length in wavelengths
% th(i) theta angle (degrees)
% num,den calculation variables
% F(i) function F from (8.74)
% Fmax maximum F(i)
% thmax angle where Fmax occurs
% diff(i) calculation variable used to find BW
% diffmin calculation variable used to find BW
% thhalf theta at the half-power point
%
clc
clear
% dipole length
L=1; %dipole length in wavelengths
% initialize settings
diffmin=.1;
8-14
% perform calculations
bL=2*pi*L;
% Find Fmax
for i=1:1:180
th(i)=i*pi/360;
num=cos((bL/2)*cos(th(i)))-cos(bL/2);
den=sin(th(i));
F(i)=(num/den)^2;
end
Fmax=max(F);
% Find thmax
for i=1:1:180
if F(i)>=Fmax
Fmax=F(i);
thmax=180*th(i)/pi;
end
end
for i=1:1:180
Pn(i)=F(i)/Fmax;
diff(i)=abs(Pn(i)-0.5);
if diff(i)<diffmin
diffmin=diff(i);
thhalf=180*th(i)/pi;
end
end
BW=2*abs(thhalf-thmax)
P8.17: How long is a 1.5l long dipole antenna at 1.0 GHz? Suppose this antenna is
constructed using AWG#20 (0.406 mm radius) copper wire. Determine Rdiss, e, and Gmax.
8-15
c 3 x108
l 0.3m, L 1.5l 0.45m
f 1x109
1l
Rdiss
s S
From example 8.2, the skin depth for this wire at 1 GHz is 2.09x10 -6m. Then, the cross-
sectional surface over which we consider the current to be conducted is:
S 2p rd Cu 5.33x10 9 m 2
Then:
1 0.45m
Rdiss 7 1
1.456W
5.8 x10 5.33x109 m 2
Wm
30
Now we need radiation resistance, Rrad F ( ) max W P , and we use Matlab 0804 to find
p
WP = 8.08 (and Dmax = 1.55), and Fmax = 1.366. Therefore, Rrad = 105W.
The efficiency is
Rrad
e 0.986
Rrad + Rdiss
Finally, Gmax = e Dmax = 1.53.
One approach is to carefully plot the pattern and then estimate the beamwidth (see Figure
P8.18). A more exact method is as follows.
�p �
cos 2 � cos �
Here we have
Pn �2 � 1 , or 2 cos 2 �p cos � sin 2
� �
sin
2
2 �2 �
We can define a function
�p �
F ( ) 2 cos 2 � cos � sin 2 , and then
�2 �
dF ( ) �p � � p � �p �
4cos � cos � sin � cos � sin � 2sin cos .
�
d �2 � �2 � �2 �
Rearranging, we have
dF �p � �p �
F '( ) 2p sin cos � cos � sin � cos � 2sin cos .
d �2 � �2 �
Now we can apply the Newton Raphson routine to converge onto a solution:
F (i )
i +1 i .
F '(i )
Newton-Raphson routine:
% MLP0818
%
8-16
for i=1:N
F=2*cos((pi/2)*cos(thr(i)))^2-sin(thr(i))^2;
dF=2*pi*sin(thr(i))*cos((pi/2)*cos(thr(i)))*sin((pi/2)*sin(t
hr(i)))-2*sin(thr(i))*cos(thr(i));
thr(i+1)=thr(i)-F/dF;
end
th=180*thr(N+1)/pi;
BW=2*(90-th)
So BW = 78°
Fig. P8.18
P8.19: A 2.45 GHz l/2 dipole
antenna is driven by a 2.0 A
amplitude current source. Find the maximum power density at a distance of 1.0 km.
15 I o2 15(2) 2 W
Pmax 19 2
pr 2
p (1000) 2
m
P8.20: Given a z-polarized half-wave dipole antenna at the origin, and a driving current
i(t) = 10cos(2px109t) A, find the instantaneous electric and magnetic fields at a point 2.0
km distant and angle =60.
� �p �
�
j r � cos � cos ��
jI e 2
H os o � � �
� a
2p r � sin �
�
� �
�
� �p ��
cos � cos 60o �
j (20.944)(2000) � �
j (10) e 2 A
� � � o
�a 6.5 x104 e j156 a
2p 2000 � sin 60
o
� m
�
� �
�
o V
Eos ho a r �H os 0.245e j156 a
m
Converting to instantaneous form:
V
E 0.245cos ( t 156o ) a
m
mA
H 0.65cos ( t 156o ) a
m
These equations can also be written in terms of sin as:
V
E 0.245sin ( t + 114o ) a
m
mA
H 0.65sin ( t + 114o ) a
m
P8.21: MATLAB: Modify MATLAB 8.4 to calculate directivity and radiation resistance
for an arbitrary length dipole antenna. Evaluate these properties for a 0.75l dipole
antenna.
% M-File: MLP0821
% Modify ML0804 to calculate D and radiation resistance.
% All that is needed is a line to calculate Rrad
% using Equation (8.79).
%
clc %clears the command window
clear %clears variables
% Initialize variables
L=.75;
bL2=pi*L;
N=90;
% Perform calculations
i=1:1:N;
dth=pi/N;
8-18
th(i)=i*pi/N;
num(i)=cos(bL2.*cos(th(i)))-cos(bL2);
den(i)=sin(th(i));
F(i)=((num(i)).^2)./den(i);
Fmax=max(F);
Pn=F./Fmax;
omegaP=2*pi*dth*sum(Pn)
Dmax=4*pi/omegaP
Fmax
% Calculate Rrad
Rrad=(30/pi)*Fmax*omegaP
Dmax =
1.8821
Fmax =
2.9142
Rrad =
185.8086
>>
c 3 x108 m s 3m
(a) l 1.2m, L = 2.5l
6
f 250 x10 / s 1.2m l
Now we use this information in MLP0821. Plugging in L = 2.5, we have
omegaP =
7.4529
D=
1.6861
Fmax =
8-19
1.6969
Rrad =
120.7662
>>
Following the same approach for the other two frequencies, we arrive at the following
table of results:
f(MHz) L(l) Dmax Rrad(W)
250 2.5 1.69 121
500 5.0 2.37 342
750 7.5 2.23 154
At 50 MHz,l = 6 m so the antenna is a l/2 dipole and Zant = 73.2+j42.5 W. We can find
the VSWR with the Smith chart. Or we can calculate it as follows:
Z Z o 73.2 + j 42.5 50 o
G L ant 0.3715e j 42.3
Z ant + Z o 73.2 + j 42.5 50
1 + GL
VSWR 2.18
1 GL
The stub matching solution uses the approach of chapter 6. We first locate the
normalized load (z = 1.46+j0.85) at point a, convert it to a normalized admittance at point
b, and move along the constant gamma circle to point c (distance traveled is 0.215l)
where the admittance is y = 1+j1.8. Then we move from a short in the admittance chart
to the point 0-j1.8, moving a distance 0.142l.
8-20
(b)
Fig. P8.23
P8.24: MATLAB: Use MATLAB 8.2 to generate plots like those of Figure 8.19 for a
dipole antenna of length 3l.
Fig. P8.24
8-21
P8.25: A 0.485l dipole antenna is constructed for operation at 4.0 GHz. (a) How long is
the antenna? (b) What impedance is required of a quarter-wave transformer to match this
antenna to a 50 W impedance line?
c 3x108 �0.075m �
l 0.075m, L ( 0.485l ) � � 0.0364m, L 3.6cm
� l �
9
f 4 x10
From section 8.3 describing the half-wave dipole, we know that a 0.485l dipole has Zant =
73W. A quarter wavelength matching section will have an impedance:
Z (73)(50) 60W
P8.26: MATLAB: Modify MATLAB 8.3 to run the movie from 0.1l up to 4l.
% M-File: MLP0826
%
% Modifies ML0803 to extend movie for L
% up to 4 wavelengths
%
clc %clears the command window
clear %clears variables
% Initialize variables
N=360;
th=1:1:N;
thr=th*pi./180;
T=num2str(L(n));
S=strvcat('Length',T,'wavelengths');
text(10,10,S)
hold on
num=cos(pi*L(n)*cos(thr))-cos(pi*L(n));
den=sin(thr);
F=(num./den).^2;
polar(thr,F)
hold off
M(:,1)=getframe;
end
The figure shows a snapshot of the movie when the angle reaches 4l.
Fig. P8.26
P8.27: MATLAB: Using MATLAB 8.4, generate data of the pattern solid angle versus
number of increments N to see the function convergence. Consider a 1.25l dipole. Try
N = 2, 4, 8, 16, 32, 64, 128.
10
8
Pattern Solid Angle (sr)
0
0 5 10 15 20 25 30 35
N iterations
Fig. P8.27
The function converges above N = 8, so data from the N = 64 and 128 runs was omitted.
4. Monopole Antennas
P8.28: Consider a 1.0 nC charge at (0.0, 0.0, 5.0m) above a conductive sheet occupying
the x-y plane at z = 0. Use image theory to find the electric field intensity at the point
(0.0, 5.0m, 5.0m).
36p m ( )
5 m (2 x
)
V
0.36a x .
m
from Q2:
1x109 C ( 5a x + 10a z )
E2
( )( )
9 3
4p 10 F 125m
36p m Fig. P8.28
V
0.032a x 0.064a z .
m
The total field is then EP = 0.33ax – 0.064az V/m.
8-24
In problem P8.18, we found the bandwidth for a l /2 dipole antenna was 78°. For the
l/4 monopole, the bandwidth will be half that of the l /2 dipole, or 39°.
P8.30: MATLAB: Devise a routine to give a polar plot of the normalized power radiated
for an arbitrary length monopole antenna. Use your program to generate the polar plot
for a half-wave monopole.
% MLP0830
% Modify ML0802 to plot the normalized power radiated
% for an arbitrary length dipole.
% We remove the current distribution plot, and plot for
% theta from -90° to +90°.
%
clc %clears the command window
clear %clears variables
clf %clear figure
% Initialize variables
Lmono=0.5;
L=2*Lmono;
bL2=pi*L;
N=180;
% Generate Plots
polar(0,1)
hold on
polar(thr,Pn)
T=num2str(Lmono);
S=strvcat('Monopole length',T,'wavelengths');
text(1.0,.8,S)
8-25
Fig. P8.30
P8.31: Determine the pattern solid angle, directivity and radiation resistance for a half-
wave monopole antenna.
P8.32: How long is a 0.75lmonopole antenna at 1.0 GHz? Suppose this antenna is
constructed using AWG#20 (0.406 mm radius) copper wire. Determine Rdiss, e and Gmax.
Compare your results with the 1.5l dipole antenna of problem P8.17.
c 3 x108 m s 3l
l 9
0.30m, L 0.225m 22.5cm
f 1x10 / s 4
8-26
So
Rrad 52.7
e 0.986 (same as a 1.5l dipole)
Rrad + Rdiss 52.7 + 0.73
and
Gmax = eDmax = 3.04 (twice that of a 1.5l dipole)
P8.33: What is the VSWR looking into a quarter-wave monopole antenna if the feed line
has a 50 W impedance? Design an open-ended shunt stub matching network to match
this antenna to the line.
Fig. P8.33
P8.34: Given a 1 GHz quarter-wave monopole antenna at the origin, excited by a 1.0 A
amplitude current, find the amplitudes for the electric and magnetic field intensities at a
point 1.0 km distant at an angle = 80.
�p �
j r cos � cos �
jI o e �2 �a and E jh H a
H os os
2p r sin
os o
�p �
2p ( 1x109 / s ) cos � cos80o �
�2 � 0.978
20.94 rad m , r 1000m, o
c 3 x108 m s sin 80
H os
( 1A ) 0.978 156 A , E 120p H 59 mV
2p 1000m
os os
m m
5. Antenna Arrays
P8.35: (JustAsk): Find and plot the far-field radiation pattern at = p/2 for a two element
dipole antenna array given the following:
1. the dipoles are driven in-phase
2. each dipole is 1l in length oriented in the z-direction
3. the pair of dipoles are 1l apart on the x axis.
Also find the maximum time-averaged power density, in W/m 2, 1.0 km away from the
array if each antenna is driven by a 1.0 A amplitude current source at 1.0 GHz.
% M-File: MLP0835a
% application of ML0804
8-28
%
clc %clears the command window
clear %clears variables
% Initialize variables
L=1.0;
bL2=pi*L;
N=90;
% Perform calculations
i=1:1:N;
dth=pi/N;
th(i)=i*pi/N;
num(i)=cos(bL2.*cos(th(i)))-cos(bL2);
den(i)=sin(th(i));
F(i)=((num(i)).^2)./den(i);
Fmax=max(F);
Pn=F./Fmax;
omegaP=2*pi*dth*sum(Pn)
D=4*pi/omegaP
Fmax
omegaP = 5.2121
D = 2.4110
Fmax = 4
>>
%MLP0835b
clc
clear
phi=.5:.5:360;
phir=phi*pi./180;
Psi=2*pi*cos(phir);
Pn=(cos(Psi./2)).^2;
polar(phir,Pn)
Fig. P8.35
P8.36: Repeat problem P8.35 if the dipoles are 180 out of phase.
Fig. P8.36
P8.37: Repeat P8.35 for the case where the dipoles are 90 out of phase, 1.5l in length,
and separated by l/2.
% M-File: MLP0837a
% modify ML0804
clc %clears the command window
clear %clears variables
% Initialize variables
L=1.5;
bL2=pi*L;
N=90;
% Perform calculations
i=1:1:N;
dth=pi/N;
th(i)=i*pi/N;
num(i)=cos(bL2.*cos(th(i)))-cos(bL2);
den(i)=sin(th(i));
F(i)=((num(i)).^2)./den(i);
Fmax=max(F);
Pn=F./Fmax;
omegaP=2*pi*dth*sum(Pn)
D=4*pi/omegaP
8-31
Fmax
Executing the program we find Wp= 8.087, D =1.554 and Fmax= 1.366.
Now:
15I o2 15I o2
Pmax F ( ) ( 1.366 ) Funit
p r2 max
p r2
Then, for the array we have
�Y � p
Farray 4 cos 2 � � , where Y d cos + p cos +
�2 � 2
� p � 15I o �p p�
2
\P � r , , � (1.366)4 cos 2 � cos + �
� 2 � pr
2
�2 4�
or
15 ( 1)
2
W
Pmax (1.366)4 26 2 .
p ( 1000 )
2
m
Then we plot:
�p p�
Pn cos 2 � cos + �
,
�2 4�
Fig. P8.37
P8.38: Two z-polarized l/2 dipole antennas are spaced l/4 apart, centered at the origin on
the x-axis. (a) If the dipole located at x = -l/8 is driven by Is1 = Ioej0º, what phase shift
would you employ on the other dipole (Is2 = Ioej) to get maximum power at a far-field
point on the +x axis? (b) If the dipole antennas are each driven by 1.0 A amplitude
currents at 500 MHz, with the phase shift from part (a), find the time-averaged power
density vector at 2.0 km on the x-axis.
(a)
�Y �
Farray 4 cos2 � �, where Y d cos + .
�2 �
2p l p
On the x-axis, = 0 so Y + +
l 4 2
2�Y�
We want maximum Farray, or Farray = 4, when cos � � 1.
�2 �
8-32
Y p
This occurs when + 0, p , 2p ...,
2 4 2
satisfied when = -p/2 or 3p/2. So we employ a -90° or a +270° phase shift.
2
� �p ��
2 �cos � cos ��
15 I o
(b) P (r , ) unit � �2 ��a .
p r � sin
2 r
�
� �
� �
2
15I o
At = p/2, we then have Funit .
p r2
The radiated power vector is then
15 ( 1)
2
�15I o2 �
P(r , ) array � 2 � 4a r 4a r 4.78Wa r
�p r � p ( 2000 )
2
P8.39: Two small loop antennas, each oriented in the x-y plane, are centered at ±l/2 on
the x-axis. They each have a 1.0 cm radius and are driven in-phase by a 10. mA current
source at 500. MHz. Find and plot the radiation pattern at = p/2 and determine the
maximum time-averaged power density at a distance 100. m from the array.
�Y � 2p
Farray 4 cos 2 � �, Y d cos + l cos + 0o 2p cos
�2 � l
so Farray 4 cos ( p cos ) , 4
2
Farray
max
so
�( 2p ) ( 500 x106 ) ( 4p x107 ) ( 0.01) ( 2p 0.6 ) p ( 0.01) 2 �
2
1
Pmax 1 loop � �
32 ( 120p ) � p ( 100 ) �
� �
Pmax 1 loop 14.2 pW
m2
( )
Pmax Pmax 1 loop 4 57 pW 2
m
2
A plot of cos (pcos) gives the same result at P8.35.
8-33
P8.40: Given a pair of dipole antennas separated by l/4 and driven in-phase, determine,
for = p/2, (a) the values for at the nulls in the radiation pattern, and (b) the values of
where the radiated power is maximum.
2p l
(b) Y d cos + cos + 0o p cos
l 4
�Y � �p �
Pn cos � � cos � cos � ,
�2 � �2 �
and maximum Pn occurs at = p/2, 3p /2.
(a)To find the location for the minimum Pn, we must take a derivative of the Pn function:
dPn p �p �
sin sin � cos � 0 at 0, p (for Pn min )
d 2 �2 �
P8.41: MATLAB: Create a movie to plot the radiated power pattern in the x-y plane for
the pair of dipoles in Example 8.7 as the separation distance varies from l/10 to 4l.
% MLP0841
% Show in a movie how the radiated power pattern
% in the x-y plane varies with separation distance
% between a pair of in-phase dipoles.
%
% Wentworth, 2/11/03
%
clc %clears the command window
clear %clears variables
% Initialize variables
N=360;
phi=1:1:N;
phir=phi*pi./180;
P8.42: Plot the normalized radiation pattern at = p/2 for 3 dipole antenna elements
spaced l/2 apart with progressive phase steps of 90º.
%MLP0842
clc
clear
8-35
phi=.5:.5:360;
phir=phi*pi./180;
Psi=pi*cos(phir)+pi/2;
N=3;
num=(sin(N.*Psi./2)).^2;
den=(sin(Psi./2)).^2;
Pn=(num./den)./(N^2);
polar(phir,Pn)
Fig. P8.42
P8.43: A particular broadside antenna array consists of 10 l/2 dipole antenna elements
spaced l/2 apart with all currents driven at the same phase. Plot the radiation pattern, and
find the maximum broadside power density (i.e. at = = p/2) at a distance of 10. km if
the antenna is driven by 10.A current sources at 2.45 GHz.
Pn 2
1 sin
2 NY
( )
2 , Y d cos + p cos
N sin Y 2
2 ( )
So we plot
1 sin ( 5p cos )
2
Pn
100 �p �
sin 2 � cos �
�2 �
%MLP0843
clc
clear
phi=.5:.5:360;
phir=phi*pi./180;
Psi=pi*cos(phir);
N=10;
num=(sin(N.*Psi./2)).^2;
den=(sin(Psi./2)).^2;
Pn=(num./den)./(N^2);
Fig. 8.43
8-36
polar(phir,Pn)
P8.44: A particular endfire antenna array consists of 10 l/2 dipole antenna elements
spaced l/2 apart with a progressive phase shift of 90º to each antenna. Plot the radiation
pattern, and find the maximum endfire power density (i.e. at = p/2 and = 0º) at a
distance of 10. km if the antenna is driven by 10A current sources at 2.45 GHz.
c 3x108
At f = 2.45 GHz we have l 0.122m
f 2.45 x109
Also, for half-wave dipoles the maximum directivity is 1.64.
Applying the Friis Equation:
2 2
Prec �l � 2 �0.122 � 9 W
Dmax1 Dmax 2 � � ( 1.64 ) ��4p ( 50 ) �
� 102 x10
Prad �4p R � � � W
or -70dB.
Now, with Prec = -35 dBm, we have Prad = -35 dBm + 70 dB = 35 dBm, or 3.2 W.
Maximum radiation occurs at = 90°, so we choose, for instance, a point on the y-axis.
Then, we know for a small loop antenna the polarization will be in a plane containing the
loop. For maximum power transfer the polarizations must match. Therefore the loop can
be parallel to the y-z plane or to the x-z plane, but not to the x-y plane. Figure P8.46
shows a suggested orientation.
Fig. P8.46
Now, to find the power transfer ratio we apply the Friis equation. We first find
c 3x108
l 0.375m
f 0.800 x109
and we also know for a half-wavelength dipole Dmax = 1.64. Since the loop diameter is
small compared to the wavelength, we have for the small loop Dmax = 1.5. So
2 2
Prec �l � � 0.375 � 9 W
Dmax1 Dmax 2 � � ( 1.64 ) ( 1.5 ) �
�4p ( 100 ) �
� 219 x10
Prad �4p R � � � W
So the power transfer ratio is Prec/Prad = -67 dB.
P8.47: A pair of z-polarized dipole antenna with lengths indicated is shown in Figure
8.50. If the 3.0 m dipole is driven by a 50. MHz source, calculate the power transfer
ratio.
2
Prec �l �
D ( , ) D ( , ) � �
Prad �4p R �
where l c f 6m and R = 4000m.
P ( , ) P ( , ) 4p
D ( , ) n n Pn ( , ) Dmax Pn ( , )
Pn ( , ) avg W p 4p W p
For the transmitter,
8-38
�p �
cos 2 � cos �
Pnt ( 60o, ) �2 � 0.667
sin
2
P8.48: Consider a pair of half-wave dipole antennas operating at 1.0 GHz and separated
by 100. m on the y-axis. Initially, both antennas are aligned in the z-direction for
maximum power transfer. Now to test the effect of polarization, the antenna at the origin
is allowed to rotate an angle in the x-z plane as shown in Figure 8.51. Plot the power
transfer ratio versus from =0º (maximum transfer case) to = 90º.
2
Pout �l �
e p D1 ( , ) D2 max � �
Pin �4p R �
�p �
cos 2 � cos �
D1 ( , ) Dmax Pn ( , ) , Pn ( , ) �2 �
sin
2
Inspecting the problem’s geometry, we see that the angle in our Pn equation is
referenced to the axis of the dipole. For simplification, we’ll let be the angle the dipole
makes with the z-axis and then let be the angle from the dipole axis to a line drawn to
the second antenna. Figure P8.48a shows this situation. Then we have
�p �p ��
cos 2 � cos � � �
Pn ( , ) �2 �2 ��
�p �
sin 2 � �
�2 �
8-39
Fig. P8.48a
Fig. P8.48b
% MLP0848
clc
clear
D=1.64;
lambda=0.3;
R=100;
B=(lambda/(4*pi*R))^2;
alpha=1:1:88;
alphar=alpha*pi./180;
A=cos((pi/2)*cos((pi/2)-alphar));
Pn=A.^2./(sin((pi/2)-alphar)).^2;
D1=D.*Pn;
ep=(cos(alphar)).^2;
Prat=ep.*D1.*D.*B;
PdB=10*log10(Prat);
plot(alpha,PdB)
xlabel('angle(degrees)')
ylabel('(Pout/Pin) dB')
grid on
8-40
P8.49: Design an open-ended shunt stub matching network to match a half-wave dipole
transmitting antenna to a source with 50 W impedance. Now suppose this antenna
network is to be used as a receiver. Use a Smith Chart to determine the impedance
looking into the matching network from the antenna.
The impedance matching network is solved similar to approach used in P8.23, only now
the shunt stub is open-ended instead of shorted. The solution is shown in FigP8.49a.
(a) (b)
Now we need to find the input impedance looking into the matching network from the
antenna, as indicated in Fig P8.49c.
First, looking into the open-ended stub of length 0.392l we see normalized admittance ya
= -j0.8. Adding this to the load admittance (yload = 1+j0) we have point b: yb = 1-j0.8. We
move from this point (at 0.344l on the WTG scale) a distance 0.215l towards the
generator to the point 0.0559l. Then we move to the impedance chart and see zin = 1.5 –
j0.8. De-normalizing, we find: Zin = 75 – j40 W. This is close the the theoretical Zin =
Zant* = 73 – j42 W.
P8.50: Referring to Figure 8.52, suppose a source voltage with amplitude 12.V and
source resistance 50 W drives a half-wave dipole transmitting antenna at 500 MHz. An
identical receiving antenna, 100. meters away and aligned for maximum power transfer,
is coupled to a 50 W load resistance. Clearly neither antenna is impedance-matched to
the transmit and receive circuitry. Calculate the voltage amplitude across this load
resistor.
P8.51: Design open-ended shunt stub matching networks for both the transmitter and
receiver of problem P8.50. Now recalculate the voltage amplitude across the load
resistor.
7. Radar
P8.52: Manipulate (8.125) using (8.113) to arrive at (8.126).
P8.53: (JustAsk): Suppose a 2 GHz radar antenna of effective area 6.0 m 2 transmits 100
kW. If a target with a 12 m 2 radar cross section is 100 km away, (a) what is the round-trip
travel time for the radar pulse? (b) What is the received power? (c) What is the
maximum detectable range if the radar system has a minimum detectable power of 2.0
pW?
ss = 12m2
R=100 km
R
(a) t 2 0.67 ms
c
ss 12m 2
( ) 2 (
6m ) 1.53 pW
2
P
(b) rec1 P A 2
100 kW
4p R l 4p ( 100km ) ( 0.15m )
rad 1 4 2 e 4
R = 93 km.
P8.54: A half-wave dipole antenna is used in a radar system to determine range to a target
that has a 1.0 m2 radar cross section. Consider that 1.0 kW is available to drive the
antenna at 300 MHz. What power is received if the target is (a) 100 m distant? (b) 1.0
km distant?
( 1.64 )
2 2
Prad 1 ( 4p ) R
3 4
( 4p ) R
3 4
(b)
Prec1 ss 6.64 x1011 m2
( )
2
A 2
100 m 27 x1018
Prad 1 4p R l 4p ( 384 x10 m ) ( 0.030m )
4 2 e 4 2
6
Prec1 27 x1012
Prad 1 1MW
27 x1018 27 x1018