ECE351 Notes PDF
ECE351 Notes PDF
1
TABLE OF CONTENTS
1.3.1. Example 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.3.2. Example 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.5.1. Example 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
1.5.2. Example 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
1.5.3. Example 3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2
1.7.1. High frequency analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
1.9.1. Example 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
1.9.2. Example 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
1.11.1. Example 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
1.11.2. Example 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
1.12.1. Decibels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
1.12.2. Phasors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
1.12.3. Example 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
1.12.4. Example 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
1.13.1. Example 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
1.14.1. Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
3
1.14.6. Load impedance example 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
2.1.4. Gradient . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
4
2.4. Electric Flux Density . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
2.5.1. Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
2.6. Divergence . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
2.8.1. Work . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
2.8.4. Potential . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
5
CHAPTER 3. PROPERTIES OF DIELECTRICS AND CONDUCTORS . . . . . . . . . . . . 77
3.1. Current . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
3.1.1. Definition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
3.5. Capacitance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
3.5.1. Definition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
3.5.3. Capacitance between two parallel plates with two dielectrics example 86
4.1.1. Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
6
4.2. Magnetic Flux Density . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
4.4. Curl . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 94
7
5.4. Maxwell’s Equations for Time-Varying Fields . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
8
CHAPTER 1. TRANSMISSION LINES
Transmission lines (TL) are used to deliver power from a source to a load. This can be done
by using coaxial lines or conductors made out of wire. When distances are large enough between
the source and the load, we use TL theory to find the voltage (V) and current (I) along the TL-line.
These TLs are usually categorized as electrically small or electrically large structures. This then
results in the following two methods to analyze the propagation of a wave along a TL:
Lumped Elements - if the time delay between the source and load is negligible (L << λ/10).
e
Distributed Elements - if the time delay between the source and load is not negligible (L>λ/10).
Lossless Transmission Line
V0
S1
{ {
I
I+
_
V + = V0
S2
a)
R
S1 S2
I+ L1 L2 L3
V0 R
C1 C2 C3
_
I
b)
Figure 1. a) Propagation of a wave along a TL; the equivalent circuit of the lossless TL.
9
1.1.1. Transmission Line Equations
Next, consider the per-unit equivalent (lumped element model) circuit of a short TL shown
in Fig. 2. R represents the conductor loss, L represents the line inductance, C represents the line
capacitance and G represents the loss in the dielectric between the conductors. Now we want to
derive expressions for V (z) and I(z) on the TL shown in Fig. 2 in terms of R, L, C and G.
∆z
Figure 2. Per-unit equivalent circuit of a TL.
1 1 ∂I 1 ∂ 1
V = R∆zI + L∆z + L∆z (I + ∆I) + R∆z(I + ∆I) + V + ∆V. (1.1)
2 2 ∂t 2 ∂t 2
⇒
( )
V 1 1 ∂I 1 ∂I ∂∆I 1 V ∆V
= RI + L + L + + R(I + ∆I) + + . (1.2)
∆z 2 2 ∂t 2 ∂t ∂t 2 ∆z ∆z
⇒
∆V ∂I 1 1 ∂∆I
= −RI − L − R∆I − L . (1.3)
∆z ∂t 2 2 ∂t
( )
∂V ∂I
= − RI + L . (1.4)
∂z ∂t
10
Next, using KCL on the circuit shown in Fig. 2:
I = IG + IC + I + ∆I
( ) ( )
∆V ∂ ∆V
= G∆z V + + C∆z V + + I + ∆I. (1.5)
2 ∂t 2
⇒
( ) ( )
I ∆V ∂ ∆V I ∆I
=G V + +C V + + + . (1.6)
∆z 2 ∂t 2 ∆z ∆z
⇒
( ) ( )
∆I ∆V ∂ ∆V
= −G V + −C V + . (1.7)
∆z 2 ∂t 2
( )
∂I ∂V
= − GV + C . (1.8)
∂z ∂t
Equations 1.4 and 1.8 are referred to as the telegraphist’s equations. Their solutions lead to the
wave equations on the TL. Next, differentiating (1.4) w.r.t. z and (1.8) w.r.t. t we get:
∂ 2V ∂I ∂2I
= −R − L (1.9)
∂z 2 ∂z ∂t∂z
and
∂ 2I ∂V ∂ 2V
= −G −C 2 . (1.10)
∂z∂t ∂t ∂t
∂ 2V ∂ 2V ∂V
2
= LC 2
+ (LG + RC) + RGV. (1.11)
∂z ∂t ∂t
11
Similar substitutions result in the following expression:
∂ 2I ∂ 2I ∂I
2
= LC 2
+ (LG + RC) + RGI. (1.12)
∂z ∂t ∂t
Equations (1.11) and (1.12) represent the general wave equations for the TL in Fig. 2.
to the following:
∂2V ∂ 2V
= LC (1.13)
∂z 2 ∂t2
and
∂ 2I ∂ 2I
= LC . (1.14)
∂z 2 ∂t2
Solving the second order partial differential equations in (1.13) and (1.14) results in the following
assumed solutions:
where ν represents the wave velocity, V + represents the forward traveling wave and V − represents
the backward traveling wave or reflected wave. The t − z/ν represents the forward traveling wave.
As t increases, z must also increase to sustain f (0) (the wave front). To solve for the wave velocity
∂2V 2
we consider the solution of V (z, t) to ∂z 2
= LC ∂∂tV2 . Without loss of generality we consider only f1 .
⇒
∂V (z, t) ∂f1 ∂f1 ∂(t − z/ν) −1 ′
= = = f (1.16)
∂z ∂z ∂(t − z/ν) ∂z v 1
′
where f1 is the partial derivative w.r.t. the argument. Similarly we have
∂ 2 f1 1 ′′
2
= 2 f1 (1.18)
∂z ν
12
and
∂ 2 f1 ′′
2
= f1 . (1.19)
∂t
Notice that (1.18) and (1.19) are simply derivatives of f1 w.r.t. z and t, respectively. Substituting
(1.18) and (1.19) into (1.13) results in the following expression:
∂2V 1 ′′ ′′ ∂ 2V
= f = LCf = LC (1.20)
∂z 2 ν2 1 1
∂t2
or
1 ′′ ′′
2
f1 = LCf1 . (1.21)
ν
′′
Canceling f1 results in the following expression for the wave velocity ν:
1
ν=√ . (1.22)
LC
[ ]
1 1
I(z, t) = f1 (t − z/ν) − f2 (t + z/ν) = I + + I − = V (z, t) (1.23)
Lν Z0
where √
L
Z0 = . (1.24)
C
Z0 is referred to as the characteristic impedance of the TL. Notice the expressions for the wave
velocity and characteristic impedance are written entirely in terms of L and C. This indicates
that ν and Z0 are dependent only on the dimensions of the physical structure and not time and
frequency.
13
⇒
⇒
V + = Vf (z, t) = |V0 | cos[ωt + βz + ϕ] (1.26)
and
V − = Vb (z, t) = |V0 | cos[ωt − βz + ϕ] (1.27)
where
ω
β= . (1.28)
νp
β is called the phase constant of the TL. If t = 0 (i.e., fix time and look at the spatial variation)
then we have Vf (z, 0) = |V0 | cos[βz]. Vf (z, 0) represents a periodic function that repeats w.r.t. a
value of z. Denote this value as λ and call it the wavelength of the wave. This then gives βλ = 2π.
⇒
2π νp
λ= = . (1.29)
β f
Now, if z = 0 (fix position and look at time variation) then we have V (0, t) = |V0 | cos[ωt]. This is
illustrated in Fig. 3. Notice that the sinusoid repeats every 2π.
⇒
1
T = (1.30)
f
14
V(0,t)
+
V0
ωt
2π
⇒
[ ] [ ]
±jx 1 jx −jx
cos(x) = Re e = e +e (1.31)
2
and
[ ] [ ]
±jx 1 jx −jx
sin(x) = ±Im e =± e −e (1.32)
2j
[ ]
−jϕ
Note that V0 was used to represent the complex voltage magnitude 1
|V |
2 0
jϕ
e +e in (1.33).
and
Vs (z) = V0 e±jβz (1.35)
where Vc (z, t) is the complex instantaneous voltage and Vs (z) is the phasor voltage. Again, from
15
the general wave equation (1.11) we have
∂ 2V ∂ 2V ∂V
2
= LC 2
+ (LG + RC) + RGV. (1.36)
∂z ∂t ∂t
⇒
d2 Vs
2
= −ω 2 LCVs + jω(LG + RC)Vs + RGVs . (1.37)
dz
⇒
d2 Vs
= (R + jωL)(G + jωC)Vs = γ 2 Vs (1.38)
dz 2
√
where γ = (R + jωL)(G + jωC) = α + jβ is the propagation constant along the TL. Now
dVs
= −(R + jωL)Is (1.39)
dz
and
dIs
= −(G + jωC)Vs . (1.40)
dz
[ ]
⇒ −γV0+ e−γz + γV0− e+γz + −γz
= −Z I0 e + γz
+ I0 e where Z0 = R + jωL. Equating the coefficients
of e−γz and eγz gives −γV0+ = −ZI0+ and γV0− = −ZI0− . Therefore, the characteristic impedance
16
Z0 is
V0+
Z0 =
I0+
V−
= − 0−
I0
Z
=
γ
Z
= √
ZY
√
Z
= (1.41)
Y
where Y = G + jωC.
⇒ √
R + jωL
Z0 = = |Z0 |ejθ . (1.42)
G + jωC
1.3.1. Example 1
Consider an 80 cm long lossless TL with a source connected to one end operating at 600 MHz.
The lumped element values of the TL are L = 0.25µH/m and C = 100pF/m. Find Z0 , β and νp .
Solution: Using (1.42) we have
√ √
R + jωL 0 + jω.25 × 10−6
Z0 = = = 50Ω.
G + jωC 0 + jω100 × 10−12
√ √ √
Also from (1.38) we have γ = (R + jωL)(G + jωC) = (0 + jωL)(0 + jωC) = 0+jβ = jω LC.
√
⇒ β = 2π(600 × 106 ) LC = 18.85rad/m. Finally, we have νp = ω/β = 2 × 108 m/s.
1.3.2. Example 2
A transmission line constructed of two parallel wires in air has a conductance of G = 0. The
two parallel wires are made of good conductors, therefore it is assumed that R = 0. If Z0 = 50Ω,
β = 20 rad/m and f = 700 MHz, find the per-unit inductance and capacitance of the TL.
Solution:
17
√ √ √
Since R = G = 0, β = ω LC = 20 = 2π700M Hz LC and Z0 = CL = 50. Then the ratio
of β and Z0 results in the following:
β
= ωC.
Z0
β
C =
ωZ0
20
=
2π ∗ 700M Hz ∗ 50
= 90.9pF/m.
L = Z02 C
= 227nH/m.
and
Vr (z) = V0r e+αz e+jβz . (1.44)
⇒
VL (0) = Vi (0) + Vr (0) = V0i + V0r (1.45)
18
and
1
IL (0) = I0i + I0r = [V0i − V0r ]. (1.46)
Z0
V0r ZL − Z0
Γ= = = |Γ|ejϕΓ . (1.47)
V0i ZL + Z0
Now using Γ we can write the voltage at the load in terms of the reflection coefficient and the
incident wave in the following manner:
Z0
+
Vi(z)
VL ZL=RL+jXL
Vr(z)
-
z=0
Figure 4. Voltage reflection at the load of a TL.
VL 2ZL
τ= =1+Γ= = |τ |ejϕt . (1.49)
V0i Z0 + ZL
( ∗
)
1 ∗ 1 −αz −jβz V0 −αz +jβz 1 |V0 |2 −2αz
< P >= Re(Vs Is ) = Re V0 e e e e = e cos θ. (1.50)
2 2 |Z0 |ejθ 2 |Z0 |
Note that θ in (1.50) refers to the angle on the characteristic impedance Z0 . Then for z = L we
19
have the following expression for the incident power:
1 |V0 |2 −2αL
< Pi >= e cos θ. (1.51)
2 |Z0 |
Then, to find the reflected power at the load substitute in the reflected wave in (1.48):
( ∗
)
1 −αL −jβL (ΓV0 ) −αL jβL
< Pr > = Re ΓV0 e e e e
2 |Z0 |ejθ
1 |Γ|2 |V0 |2 −2αL
= e cos θ. (1.52)
2 |Z0 |
This then leads to the following expression for the reflected power
< Pr >
= ΓΓ∗ = |Γ|2 . (1.53)
< Pi >
< Pt >
= 1 − |Γ|2 . (1.54)
< Pi >
1.5.1. Example 1
a) ZL = Z0 = 50Ω ⇒ Γ = 0 and τ = 1.
b) ZL = 0, Z0 = 50Ω ⇒ Γ = −1 and τ = 0.
c) ZL → ∞ (open), Z0 = 50Ω ⇒ Γ = 1 and τ = 2.
Notice from the previous derivations and examples that we have −1 ≤ |Γ| ≤ 1 (RL ≥ 0) and
0 ≤ |τ | ≤ 2 (RL ≥ 0).
1.5.2. Example 2
A TL with Z0 = 100Ω is loaded with a series connected 50 Ω resistor and a 10 pF capacitor.
Find the reflection coefficient at the load at 100 MHz.
Solution:
20
The load impedance at 100 MHz is ZL = 50 − j159Ω. Using (1.47) gives
ZL − Z0 50 − j159 − 100
Γ= = = 0.762∠ − 60.78.
ZL + Z0 50 − j159 + 100
1.5.3. Example 3
√
ZL − Z0 jXL − Z0 −(Z0 − jXL ) − Z02 + XL2 e−jθ
Γ= = = = √ 2 = −e−j2θ .
ZL + Z0 jXL + Z0 Z0 + jXL 2 jθ
Z0 + XL e
A cross-section of a coaxial TL is shown in Fig. 5. The radius of the inner conductor is a, the
radius of the inner wall on the outer conductor is b and the radius of the outer wall on the outer
conductor is c. It can be shown that the per unit values of the coaxial TL are:
( )
2πε′ F
C= , (1.55)
ln( ab ) m
( )
2πσ S
G= , (1.56)
ln( ab ) m
( ) ( )
µ b H
Lext = ln , (1.57)
2π a m
( ) ( )
1 1 1 Ω
R= + , (1.58)
2πδσc a b m
and
21
√ √ ( )
Lext 1 µ b
Z0 = = ln (Ω) (1.59)
C 2π ε′ a
where σ is the conductivity of the TL, ε′ is the permittivity of the TL and µ is the permeability of
the TL.
outer
conductor
a (σc)
dielectric
b (σ,ε,µ)
inner
conductor c
(σc)
Figure 5. Cross-section of a coaxial TL.
Also for the LF analysis of the coaxial TL we have the following expressions:
( )
2πε′ F
C= , (1.60)
ln( ab ) m
( )
2πσ S
G= , (1.61)
ln( ab ) m
( )( )
1 1 1 Ω
R= + 2 , (1.62)
πσc a 2 c −b 2 m
and
[ ( ) ( ( ))]( )
µ b 1 1 4c4 c H
Lext = ln + + b − 3c + 2
2 2
ln . (1.63)
2π a 2 2
4 4(c + b ) c −b 2 b m
22
a and the center of each wire is separated by a distance d. It is assumed that the two wires are
emersed in a material with properties (σ, µ, ε′ ). It can be shown that the per unit values of the
two-wire TL are:
( )
πε′ F
C= , (1.64)
ln( ad ) m
( )
πσ S
G= , (1.65)
cosh−1 ( 2a
d
) m
( ) ( )
µ d H
L = ln , (1.66)
π a m
( )
1 Ω
R= (1.67)
πaδσc m
and
√
L
Z0 = (Ω). (1.68)
C
dielectric
σc σc
(σ,ε,µ)
a a
( )
πε′ F
C= , (1.69)
cosh−1 ( 2a
d
) m
23
( )
πσ S
G= , (1.70)
cosh−1 ( 2a
d
) m
[ ( )] ( )
µ 1 −1 d H
L= + cosh , (1.71)
π 4 2a m
and
( )
2 Ω
R= 2 . (1.72)
πa σc m
where ϕ is the angle on the reflection coefficient. Converting (1.73) to the time domain we get:
[ ]
jωt
VsT (z, t) = Re VsT (z)e
The first term in the last expression in (1.74) describes the traveling wave and the second term
in the last expression describes the standing wave along the TL. An example of a standing wave
is illustrated in Fig. 7. The oscillations of the standing wave in Fig. 7 are described by the
trigonometric terms in (1.74). We can also derive expressions for the position of the voltage
maximums and minimums along the TL. The position of the voltage maximum is denoted as zmax
24
|VsT|
λ/2
(1+|Γ|)V0
(1−|Γ|)V0
z
−1 −1 −1 0
(φ+5π) (φ+3π) (φ+π)
−1 2β Γ 2β Γ 2β Γ
−1 −1 −φΓ
(φ+6π) (φ+4π) (φ+2π)
2β Γ
2β Γ
2β Γ
2β
Figure 7. Standing wave along a TL. Note that ϕ in this figure is the phase of the reflection coefficient
Γ = |Γ|ejϕΓ .
and the position of the voltage minimum is denoted as zmin . Thus it can be shown that
1
zmin = − [ϕΓ + (2m + 1)π] (1.75)
2β
and
1
zmax = − [ϕΓ + 2mπ] (1.76)
2β
where m = 0, 1, 2.... Next, we define the VSWR (denoted as s) in the following manner:
A few examples of a standing wave are illustrated in Fig. 8. To understand the maximum
and minimum values of the voltage along the TL we start with the following expression:
( )
−jβz
VsT (z) = V0 e + |Γ|e j(βz+ϕΓ )
. (1.78)
25
We have a voltage minimum when βz = 0 and ϕΓ = π. This then reduces (1.78) to VsT (zmin ) =
V0 (1 − |Γ|). Similarly, we have a voltage maximum when βz = 0 and ϕ = 0. This then reduces
(1.78) to VsT (zmax ) = V0 (1 + |Γ|). Then, if the TL is matched (i.e., Z0 = ZL ) then Γ = 0 and
V0 (1 + |Γ|) = V0 (1 − |Γ|) = V0 . This results in the constant voltage along the TL shown in Fig.
8 a). If ZL = 0 (i.e., short circuit), then |Γ| = −1. This then gives a maximum and minimum
voltage of (1 + 1)V0 = 2V0 and (1 − 1)V0 = 0, respectively. This is illustrated in Fig. 8 b). Finally,
if ZL = ∞ (i.e., open circuit) then |Γ| = 1. This then gives a maximum and minimum voltage of
(1 + 1)V0 = 2V0 and (1 − 1)V0 = 0, respectively. This is illustrated in Fig. 8 c).
We can also see that the standing waves in Fig. 8 b) and c) have a period of λ/2. For example,
to calculate the position of the first minimums and maximums we use the expressions in (1.75) and
(1.76), respectively. First, if ϕΓ = π, then (1.75) and (1.76) reduce to −λ/2 and −λ/4, respectively.
These computations are illustrated in Fig. 8 b). Second, if ϕΓ = 0, then (1.75) and (1.76) reduce
|V(z)|
Matched line +
|V0|
z
0
−λ −3λ −λ −λ
4 2 4
|V(z)|
a)
Short circuit +
2|V0|
z
0
−λ −3λ −λ −λ
4 2 4
b) |V(z)|
Open circuit +
2|V0|
z
0
−λ −3λ −λ −λ
4 2 4
c)
Figure 8. a) Standing wave along the TL for ZL = Z0 ; b) standing wave along the TL for ZL = 0; c)
standing wave along the TL for ZL = ∞.
26
1.9. Examples: Standing Wave Examples
1.9.1. Example 1
A 50 Ω transmission line is terminated with a load of ZL = 100 + j50 Ω. Find the voltage
ZL − Z0 100 + j50 − 50
Γ= = = .45∠26.6.
ZL + Z0 100 + j50 + 50
1.9.2. Example 2
A 140 Ω lossless transmission line is terminated with a load impedance of ZL = 280 + j182 Ω.
If λ = 72 cm, find (a) Γ, (b) s and (c) the first zmin and zmax .
Solution: From (1.47) we have
Now we need to analyze the TL while including the numerous forward and backward reflected
waves. To do this consider the TL in Fig. 9. From before, we have VsT (z) = V0+ e−jβz +V0− ejβz which
represents the total voltage along the TL. We also have IsT (z) = I0+ e−jβz + I0− ejβz which represents
27
the total current along the TL. Using these two expressions, we define the wave impedance as
[ ]
ZL cos βl + jZ0 sin βl
Zin = ZW (−l) = Z0 . (1.80)
Z0 cos βl + jZL sin βl
2π mλ 2π
Next, if βl = λ 2
= mπ where m = 0, 1, ... then Zin (l = mλ/2) = ZL . Also, if βl = λ
(2m + 1) λ4 =
(2m + 1)π/2 where m = 0, 1, ... then Zin (l = λ/4) = Z02 /ZL . The last expression is used for the
Zg
{
+ +
VS Vin Z0 VL ZL
- -
Zin z = -l z=0
Figure 9. General lossless TL in steady state.
Zs=40 Ω
+ +
VS Vin Z0=60+j40 Ω VL ZL=20+j50 Ω
- -
Zin z = -l = 2m z=0
Figure 10. TL for example 1.
28
Z0 =60+j40 Ω and is 2m long. If Vs = 10∠0, Zs = 40 Ω and ZL = 20 + j50 Ω find:
a) Zin
b) Iin
[ ]
ZL cos βl + jZ0 sin βl
Zin = ZW (−l) = Z0
Z0 cos βl + jZL sin βl
[ ]
(20 + j50) cos(2) + j(60 + j40) sin(2)
= (60 + j40)
(60 + j40) cos(2) + j(20 + j50) sin(2)
= 57.28 − j2.11.
b)
Vs
I(−l) =
Zin + Zs
10∠0
=
40 + 57.28 − j2.11
= 102.7∠1.24mA.
1.11.2. Example 2
Now consider the TL with Z0 = 300Ω. The load is two 300 Ω resistors and one capacitor with
Zc = −j300Ω all connected in parallel. Calculate Zin , s, ΓL and PL for l = 2m, νp = 2.5 × 108 m/s,
f = 100M Hz, Zg = 300Ω and Vs =60V.
and
1 + |Γ|
s= = 2.616.
1 − |Γ|
29
βl = 5.02rad = 287.6◦ . Solving for Zin gives
[ ]
ZL cos 287.6 + jZ0 sin 287.6
Zin = Z0
Z0 cos 287.6 + jZL sin 287.6
= 760.1 − j127.6Ω.
Vs
Iin =
Zg + Zin
60
=
300 + Zin
= 56.1∠6.86mA.
1 2
Pin = I R
2
1
= (56.1mA)2 760
2
= 1.199W.
⇒ PL ≈ 1.2W .
In this section we review the use of decibels. The number of decibels is the logarithmic ratio
of the powers of interest or:
( )
P2
P ower = 10 log . (1.81)
P1
The unit of the ratio in (1.81) is dB and referred to as decibels. Equation (1.81) can be used to
express the power gain of a circuit. For example, if P1 = 10W and P2 = 5W in Fig. 11 then the
power gain = 10 log(P2 /P1 ) = −3.01dB.
30
I1 I2
+ +
V1 Z1 Z2 V2
- -
1.12.2. Phasors
Now we switch our focus to phasors. In this course we use the following notation:
[ ]
V (z, t) = Re Ve e jωt
= |Vm |∠θv = Ve = phasor voltage (1.82)
where V (z, t) represents the voltage in the time domain and Ve represents the voltage in the frequency
[ ]
domain. For example, if Ve = 3∠30 = 3e then V (z, t) = Re 3e e
j30 jωt j30
= 3 cos(ωt + 30).
Next, if Ve1 and Ie1 are rms quantities,
[ ]
∗
P1 e
= Re V1 I1e
[ ]
∗
= Re (Ie1 Z1 )(Ie1 )
[ ]
= Re |Ie1 | (R1 + jX1 )
2
= |Ie1 |2 R1
|Ve1 |2
= .
R1
( ) ( )
P2 V22 R1
10 log = 10 log
P1 V 2 R2
( 12 ) ( )
V2 R1
= 10 log 2
+ 10 log
V R
( 1) ( 2)
V2 R1
= 20 log + 10 log .
V1 R2
31
Therefore,
( ) ( )
P2 V2
10 log = 20 log (1.83)
P1 V1
( )
P
P (dBm) = 10 log . (1.84)
1mW
The values written in terms of dBm simply mean that the power value of interest is referenced to
1 mW. ⇒ P = (10P (dBm)/10 )1mW .
1.12.3. Example 1
Next, as another example, express 20 dB across a 50 Ω resistor in terms of dBV, dBµV, dBA
and dBmA.
( )
70.71V
20 log = 37dBV
1V
( )
70.71V
20 log = 157dBµV
1µV
( )
1.414A
20 log = 3dBA
1A
( )
1.414A
20 log = 63dBmA.
1mA
32
1.13. Scattering Matrices or S-parameters
The scattering matrices are used to describe the response of an N-port network to various
incident and reflected voltages on the network. For illustration, consider the N-port network in Fig.
12. The incident voltage (current) and reflected voltage (current) on the N th port is denoted as VN+
(IN+ ) and VN− (IN− ), respectively.
-
V3, I3- V+3,I3+ - -
V 4, I4
V+4, I4+
- -
V2, I2 V+5,I5+
N-port
network - -
V 5 , I5
V+2,I2+
- -
V1, I1 - -
VN,IN
V+1, I1+ V+N, IN+
Figure 12. An N-port network for illustration of the S-parameters.
V1− S11 S12 S13 . . . S1N V1+
V2− S21 V2+
..
V3− .. V3+
= S31 . . )
.. .. ..
. . .
− +
VN SN 1 ... SN N VN
33
and then “zero-out” the incident voltages that are not of interest. This then in general gives
Vi−
Sij = .
Vj+ +
Vk =0 f or k̸=j
The incident voltages can be “zeroed-out” by terminating the ports not of interest with a matched
load. This then reduces the reflection coefficient Γ to zero and minimizes the reflected or scattered
waves at that port.
1.13.1. Example 1
Consider the two port circuit (N = 2) in Fig. 13(a). Since we have a two port network we
have the following scattering matrices:
V1− S11 S12 V1+
= .
V2− S21 S22 V2+
Solving for the reflected voltage at port 1 gives V1− = V1+ S11 + V2+ S12 . To compute S11 we match
V1− Zin − Zo
S11 = + = Γ = .
V1 Zin + Zo
Therefore, if a 50 Ω resistor is connected across port 2 (as shown in Fig. 13(b)), the input impedance
of port 1 can be shown to be 50 Ω also. Thus, S11 = 0.
Next, to compute S21 , the following equation is considered: V2− = V1+ S21 + V2+ S22 . For this
computation, the incident voltage V2 + on port 2 must be reduced to zero by a matched 50 Ω load.
This circuit is shown in Fig. 13(c). Then, if port 1 is driven with V1+ , the voltage on the 50 Ω load
can be computed as V2− = 0.67V1+ . This then gives
0.67V1+
S21 = = 0.67
V1+
34
or in decibels
S21 = 20log0.67 = −3.47dB.
a) b)
10Ω 10Ω
+ +
+ -
Port 1 V1 120Ω V2 50Ω
- -
c)
Figure 13. (a) two port attenuator circuit; (b) port 2 terminated with a match load to compute the input
impedance at port 1 and (c) port 2 terminated with a match load to compute S21 .
All impedance values on the Smith chart are normalized and have the following notation:
ZL RL + jXL
zL = r + jx = = . (1.85)
Z0 Z0
⇒
zL − 1
Γ= (1.86)
zL + 1
or
1+Γ
zL = . (1.87)
1−Γ
35
Next, using rectangular coordinates, we can write Γ = Γr + jΓi . ⇒
1 + Γr + jΓi
zL = r + jx = . (1.88)
1 − Γr − jΓi
and
2Γi
x= . (1.90)
(1 + Γr )2 + Γ2i
Equations (1.89) and (1.90) are the real and imaginary parts of (1.88), respectively. Rearranging
( )2 ( )2 (
gives the following family of circle expressions: Γr − 1+r + Γi = 1+r
r 2 1
and (Γr − 1) + Γi −
2
)2 ( )2
1
x
= x1 . The family of circles are shown in Figs. 14 a) and b). For example, if x = ∞, then
Γ = 1 + j0 and the radius of the circle is zero. If x = +1, then the circle is centered at 1 and 1.
⇒ Γ = 1 + j1 and the radius is 1. Also, if x = 2, then the circle is centered at 1 and 1/2. ⇒
Γ = 1 + j1/2 and the radius of the circle is 1/2. Combining both circles we get the smith chart
Γi Γi
r=0
r = 0.5
x = 0.5 x=2
r=1 r= x=1 x=
r=2 x=0
Γr Γr
x = -2
x = -0.5
x = -1
|Γ| = 1
a) b)
Figure 14. a) r-circles in the Γr and Γi plane and b) x-circles in the Γr and Γi plane.
36
The Smith Chart
0.12 0.13
0.11 0.14
0.38 0.37 0.15
0.1 0.39 0.36
90
0.4 100 80 0.35 0.1
9
0.0 6
45
50
1 110 40 70 0.3
1.0
0.4 4
0.9
1.2
0.1
55
8
0.8
0.0 35
7
1.4
2 0.3
0.7
0.4 0 60 3
12
0.6 60
/Y o)
1.6
0.1
0.0
7 ( +jB 30 8
CE 0.3
3 AN
1.8
0.4 E PT 0.2 2
65
0 SC 50
13 SU
2.0
VE
0.5
0.1
6
TI
0.0
CI 25
9
0.3
4
PA
0.4
1
CA
70
R
, O 0.4
0
o)
40
14
5
0.4
0.2
0.0
Z
5
20
0.3
0.4
jX
(+
T
3.0
75
EN
0.6
N
PO
0.2
4
0.0
0.2
6
OM
0.3
1
30
15
0.4
9
EC
>
0.8 15
R—
80
4.0
NC
TO
TA
0.22
1.0
AC
ERA
0.28
0.47
5.0
RE
1.0
GEN
160
0.2
20
85
IV E
10
UCT
ARD
0.8
0.23
IN D
0.27
0.48
S TOW
ANG
90
ANG
0.6
L E OF
N GTH
10
L E OF
170
0.1
0.4
T R A N SM
0.24
0.0 —> W A V E L E
0.26
0.49
R E FL E C T I O N C O E FFI C I E N T I N D E G
20
0.2
I SSI O N C O E F F I C I E N T I N
50
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
1.2
1.4
1.6
1.8
2.0
3.0
4.0
5.0
10
20
50
0.25
0.25
± 180
0.0
50
RESISTANCE COMPONENT (R/Zo), OR CONDUCTANCE COMPONENT (G/Yo)
A D <—
0.2
20
0.24
ARD LO
0.26
0.49
-170
0.1 0.4
DEGR
10
S TOW
REES
EES
0.6
-90
0.23
)
0.27
0.48
/ Yo
TH
EN G
( -jB
0.8 -10
-160
CE
V EL
-85
-20
0.2
AN
1.0
5.0
0.22
WA
0.28
0.47
T
EP
1.0
<—
SC U
-15 -80
4.0
ES
0.8 -15
IV
0.2
4
0
-30
T
0.0
0.3
UC
0.2
6
1
0.4
9
D
IN
0.6
-75
3.0
O
),
5
Zo
0.2
-20
0.0
X/
5
0.3
0.4
0.4
40
( -j
-4
-1
T 0.4
EN
-70
N
PO
0.1
6
0.0
OM
9
-25
0.3
4
EC
0.4
0.5
1
2.0
30 NC -5
TA 0
-65
-1 C 0.1
7 EA
1.8
0.2
0.0 R
IV E
8
0.6
0.3
3 IT -30
0.4 PA C 2
1.6
CA
-60
0 -60 0.1
8 -12
0.7
0.0 7
1.4
-35 0.3
0.8
2
0.4 3
1.2
-55
0.9
1.0
9 -70 0.1
0.0 -110 -4
0 6
0
-5
0.3
-4
1
0.4 0.1 -100 -80 0.15 4
-90
0.11 0.14 0.35
0.4 0.12 0.13
0.39 0.36
0.38 0.37
T O W A R D L O A D —> <— T O W A R D G E N E R A T O R
SW d [dB F, P or I
RT
)
N. RFL CO
S A
, P ST
LO . C EF
SS OE F, E
RF
S R
I
or
1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.05 0.01 0 0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.5 3 4 5 10 ∞
,E
FF
1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 1 0.99 0.95 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0
TR
OE
.C
SM
CENTER
AN
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2
TR
ORIGIN
37
Z0 = 50Ω. Find ΓL .
Solution: First, normalize the load impedance in the following manner:
ZL 25 + j50
zl = = = .5 + j1.
Z0 50
38
The Smith Chart
0.12 0.13
Angle of reflection
0.11 0.14
0.1 0.39
0.38 0.37
0.36 0.15 coefficient
90
0.4 100 80 0.35 0.1
9
0.0 6
45
50
1 110 40 70 0.3
1.0
0.4 4
0.9
1.2
0.1
55
8
0.8
0.0 35
7
1.4
2 0.3
0.7
0.4 0 60 3
12
0.6 60
)
/ Yo
1.6
0.1
0.0
7 ( +jB 30 8
CE 0.3
3 AN
1.8
0.4 E PT 0.2 2
65
0 SC 50
13 SU
2.0
VE
0.5
0.1
6
TI
0 .0
CI 25
9
0.3
4
PA
0.4
1
CA
70
R
, O 0.4
0
o)
40
14
5
0.4
0.2
0.0
Z
5
20
0.3
0.4
jX
A
(+
3.0
T
75
EN
0.6
N
PO
0.2
4
0.0
0.2
6
0.3
OM
1
30
15
0.4
9
EC
15
>
0.8
R—
80
4.0
NC
TO
TA
1.0
0.22
AC
ERA
0.28
0.47
5.0
RE
1.0
GEN
0.2
20
85
IV E
1 6
10
UCT
ARD
0.8
0.23
IN D
0.27
0.48
S TOW
ANG
90
0.6
ANG
L E OF
N GTH
10
L E OF
170
0.1
0.4
T R A N SM
0.24
0.0 —> W A V E L E
0.26
0.49
R E FL E C T I O N C O E FFI C I E N T I N D E G
20
0.2
I SSI O N C O E F F I C I E N T I N
50
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
1.2
1.4
1.6
1.8
2.0
3.0
4.0
5.0
10
20
50
0.25
0.25
180
0.0
Origin
±
50
RESISTANCE COMPONENT (R/Zo), OR CONDUCTANCE COMPONENT (G/Yo)
L O A D <—
0.2
20
0.24
0.26
0.49
D
70
0.1 0.4
R
DE
-1
A
10
TOW
GREE
REES
0.6
-90
0.23
S
0.27
0.48
S
/ Yo
T H
NG
( -jB
0.8 -10
E
0
CE
L
-85
-20
6
VE
0.2
-1
1.0
5.0
0.22
A
TA
0.28
7
W
.4
EP
1.0
<—
0
SC
SU
-15 -80
4.0
0.8 -15
VE
0.2
4
0
-30
CT
0.0
0.3
0.2
6
1
0.4
DU
9
IN
0.6
-75
3.0
O ,
o)
5
0.2
-20
0.0
Z
X/
5
0.3
0.4
0.4
( -j
-4
4
-1
T 0.4
EN
-70
ON
0.1
6
P
0.0
OM
9
-25
0.3
4
EC
0.4
-65 .5
1
2.0
0
C -5
30 AN 0
-1 T
7 AC 0.2 0.1
1.8
0.0 RE 8
VE
0.6
ITI 0.3
3
0.4 AC -30 2
1.6
CAP
-60
0 -60 0.1
8 -12
0.7
0.0 7
1.4
-35 0.3
0.8
2
0.4 3
1.2
-55
0.9
1.0
9 -70 0.1
0.0 -110 -4
0 6
0
-5
0.3
-4
1
0.4 0.1 -100 -80 0.15 4
-90
0.11 0.14 0.35
0.4 0.12 0.13
0.39 0.36
0.38 0.37
T O W A R D L O A D —> <— T O W A R D G E N E R A T O R
SW d [dB F, P or I
RT
)
N. RFL CO
S A
, P ST
LO . C EF
SS OE F, E
RF
S R
I
or
1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.05 0.01 0 0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.5 3 4 5 10 ∞
,E
FF
1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 1 0.99 0.95 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0
TR
OE
.C
SM
CENTER
AN
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2
TR
ORIGIN
39
1.14.3. Input impedance example
We can also determine the input impedance of the TL. If the length of the TL is l = 60 cm,
λ = 2 m, ZL = 25 + j50 and Z0 = 50Ω, find zin .
ZL 25 + j50
zl = = = .5 + j1.
Z0 50
Next, we need to determine the length of the TL in terms of the wavelength of the source. ⇒
l/λ = 0.6/2 = 0.3. ⇒ l = 0.3λ. This means that we have to move 0.3λ down the TL from the load
to the source. This then places us at the port of the TL. To do this on the Smith chart, we need
to extend the line through point A to intersect the circle labeled “WAVELENGTHS TOWARD
GENERATOR”. This is shown in Fig. 17. This is then our reference value. This intersection
happens at l1 = 0.1345λ. Next, we need move along the circle labeled “WAVELENGTHS TOWARD
GENERATOR” a distance of 0.3λ (which is the length of the TL). Adding the reference value to
the length of the TL gives 0.3λ + 0.1345λ = 0.4345λ. This added value tells us what value we
need to move to along the “WAVELENGTHS TOWARD GENERATOR” circle to be at the port
of the TL. Next, rotate the line going through A around the Smith chart in Fig. 17 in the clockwise
direction. We are now at point B. Notice that a line extending from the origin through B intersects
the “WAVELENGTHS TOWARD GENERATOR” circle at 0.4345λ. Reading the normalized input
impedance value at point B gives zin = 0.28 − j0.4. ⇒ Zin = 50(zin ) = 14 − j20Ω. Note that the
analytical solution using (1.80) gives 13.7-j20.2.
40
The Smith Chart
0.1345 λ
0.12 0.13
0.11 0.14
0.38 0.37 0.15
0.1 0.39 0.36
90
0.4 100 80 0.35 0.1
9
0.0 6
45
50
1 110 40 70 0.3
1.0
0.4 4
0.9
1.2
0.1
55
8
0.8
0.0 35
7
1.4
2 0.3
0.7
0.4 0 60 3
12
0.6 60
)
/ Yo
1.6
0.1
0.0
7 ( +jB 30 8
CE 0.3
3 AN
1.8
0.4 EP
T 0.2 2
65
0 C 50
13 S US
2.0
VE
0.5
0.1
6
TI
0.0
CI 25
9
0.3
4
A
0.4
AP
1
C
70
R
, O 0.4
0
o)
40
14
5
0.4
0.2
A
0.0
Z
5
20
0.3
0.4
jX
(+
3.0
T
75
EN
0.6
N
PO
0.2
4
0.0
0.2
6
0.3
1
30
15
0.4
9
EC
15
>
0.8
R—
80
4.0
NC
TO
TA
1.0
0.22
AC
ERA
0.28
0.47
5.0
RE
1.0
GEN
160
0.2
20
85
IV E
10
UCT
ARD
0.8
0.23
IN D
0.27
0.48
TOW
ANG
90
ANG
0.6
GTH S
L E OF
10
L E OF
170
0.1
0.4
N
T R A N SM
0.24
0.0 —> W A V E L E
0.26
0.49
R E FL E C T I O N C O E FFI C I E N T I N D E G
20
0.2
I SSI O N C O E F F I C I E N T I N
50
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
1.2
1.4
1.6
1.8
2.0
3.0
4.0
5.0
10
20
50
0.25
0.25
± 180
0.0
0.2
20
0.24
0.2
0.49
6
ARD
-170
0.1 0.4
DEGR
10
S TOW
REES
EES
0.6
-90
0.23
)
0.27
0.48
/ Yo
GTH
( -jB
0.8
LEN
-10
-160
E
-85
-20
NC
VE
0.2
1.0
5.0
0.22
A
TA
—W
0.28
0.47
EP
1.0
B
SC
<
U
-15 -80
4.0
ES
0.8 -15
V
0.2
4
0
TI
-30
0.0
0.3
C
0.2
6
1
0.4
DU
9
IN
0.6
-75
3.0
O
),
5
Zo
0.2
-20
0.0
X/
5
0.3
0.4
0.4
40
( -j
-4
-1
T 0.4
EN
-70
ON
0.1
6
P
0.0
M
9
CO -25
0.3
4
0.4
0.5
CE
1
2.0
30 N -5
TA 0
-65
-1
7 AC 0.2 0.1
1.8
0.0 RE
VE
8
0.3 λ rotation
0.6
I 0.3
3
0.4 C IT A -30 2
1.6
CAP
-60
-60
0.4345 λ
0 0.1
8 -12
0.7
0.0 7
1.4
-35 0.3
0.8
2
0.4 3
1.2
-55
0.9
1.0
9 -70 0.1
0.0 -110 0
-4
6
0
-5
0.3
-4
1
0.4 0.1 -100 -80 0.15 4
-90
0.11 0.14 0.35
0.4 0.12 0.13
0.39 0.36
0.38 0.37
T O W A R D L O A D —> <— T O W A R D G E N E R A T O R
SW d [dB F, P or I
RT
)
N. RFL CO
S A
, P ST
LO . C EF
SS OE F, E
RF
S R
I
or
1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.05 0.01 0 0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.5 3 4 5 10 ∞
,E
FF
1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 1 0.99 0.95 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0
TR
OE
.C
SM
CENTER
AN
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2
TR
ORIGIN
41
1.14.4. Voltage minimum and maximum example
We can also determine the voltage minimum and maximum values along the TL using the
Smith chart. Again, let ZL = 25 + j50 and Z0 = 50Ω and determine zmax or zmin closest to the
load. Note that a zmax occurs when r > 1 and zmin occurs when r < 1.
Solution: For our example we have
ZL 25 + j50
zl = = = .5 + j1.
Z0 50
Again, we plot point A (i.e., zl ) on the Smith chart in Fig. 18. Using the line drawn from the origin
through point A as a reference, we move from the load to the source along the “WAVELENGTHS
TOWARD GENERATOR” circle. When the line has rotated to point B we have encountered
our first voltage extremum. Since r = 4.3 > 1 at point B, then we have a voltage maximum.
This maximum occurs at 0.25λ − 0.1345λ = 0.1155λ from the load. This implies that the voltage
minimum occurs λ/4 further down the TL.
42
The Smith Chart
0.1345 λ
0.12 0.13
0.11 0.14
0.38 0.37 0.15
0.1 0.39 0.36
90
0.4 100 80 0.35 0.1
9
0.0 6
45
50
1 110 40 70 0.3
1.0
0.4 4
0.9
1.2
0.1
55
8
0.8
0.0 35
7
1.4
2 0.3
0.7
0.4 0 60 3
12
0.6 60
)
/ Yo
1.6
0.1
0.0
7
E ( +jB 30 8
3 NC 0.3
TA
1.8
0.4 EP 0.2 2
65
0 C 50
13 S US
2.0
VE
0.5
6
0.1
TI
0.0
CI 25
9
4
0.3
PA
0.4
1
CA
70
R
,O 0.4
0
o)
40
14
5
0.4
0.2
A
0 .0
Z
5
20
0.3
0.4
jX
(+
3.0
T
75
EN
0.6
N
PO
4
0.2
0.0
0
6
0.2
M
0.3
1
30
15
0.4
9
EC
15
>
0.8
R—
80
4.0
NC
TO
TA
1.0
0.22
AC
ERA
0.28
0.47
5.0
RE
1.0
GEN
160
0.2
20
85
IV E
10
UCT
ARD
0.8
0.23
IN D
0.48
0.27
S T OW
ANG
90
ANG
0.6
L E OF
N GT H
10
L E OF
170
0.1
0.4
T R A N SM
0.24
0.0 —> W A V E L E
0.49
0.26
R E FL E C T IO N O E FFI C I E N T I N D E G
20
0.2
I SSI O N C O E F F I C I E N T I N
50
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
1.2
1.4
1.6
1.8
2.0
3.0
4.0
5.0
10
20
50
0.25
0.25
± 180
0.0
C
RESISTANCE COMPONENT (R/Zo), OR CONDUCTANCE COMPONENT (G/Yo)
Origin x=0 B 50
L O A D <—
0.2
20
0.24
0.49
0.26
ARD
-170
0.1 0.4
DEGR
10
S T OW
REES
EES
0.6
-90
0.23
o)
0.48
0.27
GTH
/Y
( -jB
0.8
LEN
-10
-160
E
-85
-20
NC
VE
0.2
1.0
5.0
0.22
A
A
—W
0.28
0.47
PT
1.0
E
SC
<
U
-15 -80
4.0
ES
0.8 -15
IV
4
0.2
0
-30
T
0.0
0.3
UC
6
0.2
1
0.4
9
D
IN
0.6
-75
3.0
,O )
5
Zo
0.2
-20
0.0
X/
5
0.3
0.4
0.4
( -j
-4
4
-1
T 0.4
EN
-70
N
PO
6
0.1
0.0
M
9
O -25
4
0.3
EC
0.4
0.5
1
2.0
30 NC A
-5
0
-65
-1 CT 0.1
7 A 0.2
1.8
0.0 RE 8
IV E
0.6
0.3
3
0.4
IT -30
PA C 2
1.6
CA
-60
0 -60 0.1
8 -12
0.7
0.0 7
1.4
2 -35 0.3
0.8
0.4 3
1.2
-55
0.9
1.0
9 -70 0.1
0.0 -110 0
-4
6
0
-5
0.3
-4
1
0.4 0.1 -100 -80 0.15 4
-90
0.11 0.14 0.35
0.4 0.12 0.13
0.39 0.36
0.38 0.37
T O W A R D L O A D —> <— T O W A R D G E N E R A T O R
SW d [dB F, P or I
)
N. RFL CO
S A
, P ST
LO . C EF
SS OE F, E
RF
S R
I
or
1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.05 0.01 0 0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.5 3 4 5 10 ∞
,E
FF
1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 1 0.99 0.95 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0
TR
OE
.C
SM
CENTER
AN
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2
TR
ORIGIN
43
1.14.5. Stub matching example
Next we want to match a load to a 50 Ω TL by placing a short-circuit stub a distance of d
from the load. This is shown in Fig. 19. Assume the length of the short circuit stub is d1 and
has a characteristic impedance of Z0s = 50 Ω (same as the TL). We need to determine d and d1
in Fig. 19. To do this we work with admittances. For this example let zL = 2.1 + j.8. Next, plot
zL on the Smith chart in Fig 20 (point A). Then draw a line from A through the origin to the
intersection of the circle about the origin with a radius of point A. This intersection with the circle
is shown as point B in Fig. 20. The intersection is the admittance of zL and has the normalized
value yL = .42 − j.16. Also note that we use .47 λ as our reference value on the “WAVELENGTHS
TOWARDS GENERATOR” scale. We will use this value to measure how far we need to move
towards the generator to observe an input admittance value of 1 ± jb. Next, we rotate towards the
generator from point B around the circle with a radius of point A to point C on the Smith chart in
Fig. 20. Point C has an input admittance of 1 + j.95. We can also rotate around to point D which
has an input admittance of 1 − j.95. The points C and D are important because both points have
a real value of unity and a reactance that can be canceled with a reactive load.
d1
Z0 ZL
Zin = Z0+j0
(with stub) short-circuit stub
Figure 19. Short-circuit stub matching example.
Next, the distance from point B to point C along the TL is .16λ + (.5 − .47)λ = .19λ and the
distance from point B to point D along the TL is .34λ + (.5 − .47)λ = .37λ. Choosing the shortest
44
distance gives d=.19 λ. At this point (i.e., point C) the input admittance is 1 + j.95. Therefore, we
need to add a stub with an admittance of −j.95 at this point. Since the stub in Fig. 19 is a short
circuit then zs = 0 ⇒ ys = ∞ which is denoted as point E on the Smith chart in Fig. 20. Next, we
need to determine how long the stub needs to be to have an input admittance of −j.95. Plotting
−j.95 on the Smith chart results in the point denoted as F in Fig. 20. Next, we move from point
E to point F towards the generator and this distance is the length of the stub. This gives a value
of d1 = .379λ − .25λ = .129λ. Another way to think about d1 is to view d1 as the distance needed
to travel from the short-circuit load down the TL to observe an input admittance of −j.95. This
distance traveled is then the length of the short-circuit stub. Now adding the admittances at the
junction of the TL and short-circuit stub in Fig. 19 gives yin = 1 + j.95 − j.95 = 1. ⇒ zin = 1. ⇒
Zin = 50 ∗ 1 = 50Ω = Z0 .
45
The Smith Chart
0.12 0.13
0.11 0.14
0.38 0.37
0.1 0.39
90
0.36
80 0.35
0.15
0.16 λ
9 0.4 100 0.1
0.0 6
45
50
1 110 40 70 0.3
1.0
0.4 4
0.9
1.2
0.1
55
8
0.8
0.0 35
7
1.4
0.3
.42
0.7
0 60 3
0 12
0.6 60
)
/ Yo
1.6
0.1
0.0
7 ( +jB 30 8
CE 0.3
3 AN
1.8
0.4 EP
T 0.2 2
65
0 C 50
13 S US
2.0
VE
0.5
6
0.1
TI
0.0
CI 25
9
4
0.3
PA
0.4
1
CA
70
R
,O 0.4
0
o)
40
14
5
0.4
0.2
0.0
Z
5
20
0.3
0.4
jX
(+
3.0
T
75
EN
0.6
N
PO
4
0.2
0.0
0
6
0.2
M
0.3
1
30
15
0.4
9
EC
15
>
0.8
R—
80
4.0
NC
TO
TA
1.0
0.22
AC
ERA
0.28
0.47
C 5.0
RE
1.0
GEN
160
0.2
20
85
IV E
10
UCT
ARD
0.8
0.23
IN D
0.48
0.27
S T OW
ANG
90
ANG
0.6
L E OF
N GT H
10
L E OF
170
0.1
0.4
T R A N SM
0.24
0.0 —> W A V E L E
A
0.49
0.26
R E FL E C T IO N C O E FFI C I E N T
20
0.2
I SSI O N C O E F F I C I E N T I N
50
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
1.2
1.4
1.6
1.8
2.0
3.0
4.0
5.0
10
20
50
0.25
0.25
± 180
0.0
0.2
20
0.2
0.49
0.2
I
4
N
6
B
ARD
-170
0.4
D
0.1
DEGR
E
E
GREE
10
S T OW
EES
0.6
-90
S
0.23
s = 2.5 circle
)
0.48
0.27
/ Yo
GTH
jB
0.8
D
LEN
E (-
-10
-160
-85
-20
NC
VE
0.2
1.0
5.0
0.22
WA
0.28
0.47
PT
1.0
<—
E
SC U
-15 -80
4.0
ES
0.8 -15
0.47 λ
IV
4
0.2
0
-30
T
0.0
0.3
UC
6
0.2
1
0.4
9
ND
0.6
I
-75
3.0
,O
o)
5
0.2
-20
0.0
Z
X/
5
0.3
0.4
0.4
40
( -j
-4
-1
T 0.4
EN
-70
N
PO
6
0.1
0.0
M
9
CO -25
4
0.3
0.4
0.5
CE
1
2.0
30 N -5
TA 0
-65
-1
7 AC 0.2 0.1
1.8
0.0 RE 8
VE
0.6
I 0.3
3
0.4 C IT A -30 2
1.6
CAP
-60
0 -60 0.1
8 -12
0.7
0.0 7
1.4
2 -35 0.3
0.8
0.4
F
3
1.2
-55
0.9
1.0
9 -70 0.1
0.0 -110 0
-4
6
0
-5
0.3
-4
1
0.4 0.1 -100 -80 0.15 4
-90
0.4 0.11
0.12 0.13
0.14
0.36
0.35
0.34 λ
0.39
0.38 0.37
0.379 λ
RADIALLY SCALED PARAMETERS
OE (CO dB CO ]
F
. C K SS [ SS [dB
EF
SM PEA LO LO EN.
T O W A R D L O A D —> <— T O W A R D G E N E R A T O R
SW d [dB F, P or I
RT R FL. C
. T
)
TR S.W RF S.W AT
N. FL O
, P ST
LO . C EF
SS OE F, E
R
I
or
1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.05 0.01 0 0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.5 3 4 5 10 ∞
,E
FF
1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 1 0.99 0.95 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0
OE
.C
SM
CENTER
AN
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2
TR
ORIGIN
46
1.14.6. Load impedance example 1
Let s = 3 on a 50 Ω TL. We have zmin = −5cm from the load and the distance between the
minimum values is 20cm. Find ZL .
Solution: We know the value between the minimum values is λ/2 = 20cm. ⇒ λ = 40cm. This
then gives zmin = −5cm/40cm = −0.125λ. Point A on the Smith chart in Fig. 21 corresponds to
the s = 3 circle and point B corresponds to the zmin value (i.e., r<1). Now, move towards the load
0.125λ. ⇒ λ′ = (0+0.125)λ = 0.125λ. Point C in Fig. 21 is zL = .6−j.8. ⇒ ZL = 50zL = 30−j40.
47
The Smith Chart
0.12 0.13
0.11 0.14
0.38 0.37 0.15
0.1 0.39 0.36
90
0.4 100 80 0.35 0.1
9
0.0 6
45
50
1 110 40 70 0.3
1.0
0.4 4
0.9
1.2
0.1
55
.08
0.8
0 7
35
1.4
2 0.3
0.7
0.4 0 60 3
12
0.6 60
)
/ Yo
1.6
0.1
0.0
7
E ( +jB 30 8
3 NC 0.3
TA
1.8
0.4 P 0.2 2
CE
65
0 50
13 S US
2.0
VE
0.5
6
0.1
TI
0.0
CI 25
9
4
0.3
PA
0.4
1
CA
70
R
,O 0.4
0
o)
40
14
5
0.4
0.2
0.0
Z
5
20
0.3
0.4
jX
(+
3.0
T
75
EN
0.6
N
PO
4
0.2
0.0
0
6
0.2
OM
0.3
1
30
15
0.4
9
EC
15
>
0.8
R—
80
4.0
NC
TO
TA
1.0
0.22
AC
ERA
0.28
0.47
5.0
RE
1.0
GEN
160
0.2
20
85
IV E
10
UCT
ARD
0.8
0.23
IN D
0.48
0.27
S T OW
ANG
90
ANG
0.6
L E OF
N GT H
10
L E OF
170
0.1
0.4
T R A N SM
0.24
0.0 —> W A V E L E
0.49
0.26
R E FL E C T IO N O E FFI C I E N T I N D E G
20
0.2
I SSI O N C O E F F I C I E N T I N
B 50
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
1.2
1.4
1.6
1.8
2.0
3.0
4.0
5.0
10
20
50
0.25
0.25
± 180
0.0
C
RESISTANCE COMPONENT (R/Zo), OR CONDUCTANCE COMPONENT (G/Yo)
Origin A 50
A D <—
0.2
20
0.24
ARD LO
0.49
0.26
-170
0.1 0.4
DEG
E 10
T OW
REES
R EES
0.6
-90
0.23
TH S
)
0.48
0.27
s = 3.0 circle
/ Yo
NG
jB
0.8
E (-
-10
ELE
-160
-85
-20
C
0.2
V
AN
1.0
5.0
0.22
WA
0.28
0.47
T P
1.0
<—
CE S
SU
-15 -80
4.0
0.8 -15
VE
C
4
0.2
0
-30
CT
0.0
0.3
.46
0.2
1
DU
0
9
IN
0.6
-75
OR
3.0
,
o)
5
0.2
-20
0.0
/Z
5
0.3
0.4
X
0.4
40
( -j
-4
-1
T 0.4
EN
-70
N
PO
6
0.1
0.0
M
9
CO -25
4
0.3
0.4
0.5
CE
1
2.0
30 N -5
TA 0
-65
-1
7 AC 0.2 0.1
1.8
0.0 RE 8
VE
0.6
I 0.3
3
0.4 C IT A -30 2
1.6
CAP
-60
0 -60 0.1
8 -12
0.7
0.0 7
1.4
2 -35 0.3
0.8
0.4 3
1.2
-55
0.9
1.0
9 -70 0.1
0.0 -110 -4
0 6
0
-5
0.3
-4
1
0.4 0.1 -100 -80 0.15 4
-90
0.11 0.14 0.35
0.4 0.12 0.13
0.39 0.36
0.38 0.37
0.125λ
RADIALLY SCALED PARAMETERS
OE (CO dB CO ]
F
. C K SS [ SS [dB
EF
SM PEA LO LO EN.
T O W A R D L O A D —> <— T O W A R D G E N E R A T O R
SW d [dB F, P or I
RT
. T
)
TR S.W RF S.W AT
N. RFL CO
, P ST
LO . C EF
SS OE F, E
RF
AN . L.
I
or
1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.05 0.01 0 0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.5 3 4 5 10 ∞
,E
FF
1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 1 0.99 0.95 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0
OE
.C
M
CENTER
A NS
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2
TR
ORIGIN
48
1.14.7. Load impedance example 2
A lossless 100Ω TL of length 3λ/8 is terminated with an unknown impedance. If Zin = −j2.5
then calculate
a) ZL
b) What length of open and short circuit stubs could be used to replace ZL ?
Solution:
a) To calculate the unknown load we first need to normalize the input impedance. ⇒ zin =
Zin /Z0 = −j.025. This point is plotted as point A in Fig. 22. Next, we need to express the length
of the TL in terms of a wavelength. ⇒ L = 3λ/8 = .375λ. This means we need to move from the
source on the TL to the load that is .375λ down the TL. The scale we use for this rotation is labeled
“WAVELENGTHS TOWARDS LOAD”. Rotating .375λ from point A at .004λ puts us at point B
on the Smith chart in Fig. 22. Point B is at .379λ. The load impedance at point B is zL = +j.95.
⇒ ZL = +j95Ω.
b) Next, to find the length of the the open circuit stub that has an input impedance of ZL
we start by plotting point C (zL = ∞) on the Smith chart in Fig. 22. Next, moving towards the
generator on the “WAVELENGTHS TOWARDS GENERATOR” scale to point B we get an input
impedance for the stub of zL = +j.95. The length of the stub is then do = .25λ + .121λ = .371λ.
Next, we plot the short circuit load (zl = 0) on the stub at point D on the Smith chart in Fig. 22.
Then moving towards the generator to point B gives a stub length of ds = .121λ.
What we have shown in part b) is that an open circuit stub of .371λ or a short circuit stub of
.121λ could be used as an equivalent load to the TL as the lumped load ZL determined in part a).
49
Finding the
length of the The Smith Chart
short circuit
stub
0.12 0.13
0.14
0.1
0.11
0.38 0.37
0.36 0.15
Distance
0.39
90
0.0
9 0.4 100 80 0.35 0.1
6 to the load
45
50
1 110 40 70 0.3
1.0
0.4 4
0.9
1.2
0.1
55
8
B
0.8
0.0 35
7
1.4
2 0.3
0.7
0.4 0 60 3
12
0.6 60
/Y o)
1.6
0.1
0.0
7
E ( +jB 30 8
3 NC 0.3
TA
1.8
0.4 P 0.2 2
CE
65
0 S 50
13 SU
2.0
VE
0.5
6
0.1
TI
0.0
CI 25
9
4
0.3
PA
0.4
1
CA
70
R
,O 0.4
0
o)
40
14
5
0.4
0.2
0.0
Z
5
20
0.3
0.4
jX
(+
3.0
T
75
EN
0.6
N
PO
4
0.2
0.0
0
6
0.2
OM
0.3
1
30
15
0.4
9
EC
15
>
0.8
R—
80
4.0
NC
TO
TA
1.0
0.22
AC
ERA
0.28
0.47
5.0
RE
1.0
GEN
160
0.2
20
85
IV E
10
UCT
ARD
0.8
0.23
IN D
0.48
0.27
S T OW
ANG
90
ANG
0.6
L E OF
N GT H
10
L E OF
0
0.1
17
0.4
T R A N SM
0.24
0.0 —> W A V E L E
0.49
0.26
R E FL E C T IO N O E FFI C I E N T I N D E G
20
0.2
I SSI O N C O E F F I C I E N T I N
50
D
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
1.2
1.4
1.6
1.8
2.0
3.0
4.0
5.0
10
20
50
0.25
0.25
± 180
0.0
C
RESISTANCE COMPONENT (R/Zo), OR CONDUCTANCE COMPONENT (G/Yo)
Origin 50
A D <—
A 0.2
20
0.2
ARD LO
0.49
0.2
4
6
-170
0.1 0.4
D
C
EGRE
10
W
REES
T O
0.6
-90
0.23
ES
TH S
)
0.48
0.27
/ Yo
G
( -jB
0.8
N
-10
ELE
-160
CE
-85
-20
0.2
V
1.0
5.0
0.22
A
TA
W
0.28
0.47
1.0
<—
CES U
-15 -80
4.0
ES
0.8 -15
IV
4
0.2
0
-30
T
0.0
0.3
UC
6
0.2
1
0.4
9
ND
0.6
I
-75
3.0
,O
o)
5
0.2
-20
0.0
Z
X/
5
0.3
0.4
0.4
( -j
-4
4
-1
T 0.4
EN
-70
N
PO
6
0.1
0.0
M
9
O -25
4
0.3
EC
0.4
0.5
1
2.0
30 NC A
-5
0
-65
-1 CT 0.1
7 A 0.2
1.8
0.0 RE 8
IV E
0.6
0.3
3
0.4
IT -30
PA C 2
1.6
CA
-60
0 -60 0.1
8 -12
0.7
0.0 7
1.4
2 -35 0.3
0.8
0.4 3
1.2
-55
0.9
1.0
9 -70 0.1
0.0 -110 -4
0 6
0
Finding the
-5
0.3
-4
T O W A R D L O A D —> <— T O W A R D G E N E R A T O R
SW d dB F, P or I
RT
. T
)
TR S.W RF S.W AT
N. RFL CO
, P ST
LO . C EF
SS OE F, E
RF
AN . L.
I
or
1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.05 0.01 0 0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.5 3 4 5 10 ∞
,E
FF
1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 1 0.99 0.95 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0
OE
.C
SM
CENTER
AN
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2
TR
ORIGIN
50
1.15. Transient Analysis
1.15.1. Analytical solution
In this section we study the transient characteristics of a TL. To do this we consider the TL in
Fig. 23 a). The problem has a forward traveling wave with amplitude V + . The wave is propagating
at a wave velocity of ν (group velocity). At time t = 0, the DC battery is connected to the TL
through the switch. At this point a wave front starts propagating down the TL. The closing of the
switch can be modeled as the step function in Fig 23 b). If the TL is lossless, then the current
can be written as IL = V0 /RL . Now, if RL ̸= Z0 , then the reflected wave adds with the forward
traveling wave giving a total wave along the TL as V1+ + V1− were V1+ is the first forward traveling
wave and V1− is the first reflected wave. ⇒
RL − Z0 V−
ΓL = = 1+ (1.91)
RL + Z0 V1
at the load. Since, Rg = 0, we have a reflected wave at the source. This then gives
Zg − Z0 Z0 V+
Γg = =− = −1 = 2− . (1.92)
Zg + Z0 Z0 V1
Now V2+ propagates to the load. ⇒ V2− = ΓL V2+ . This returns to the battery and V3+ = −V2− .
If Rg ̸= 0 then
V0 Z0
V1+ = . (1.94)
Rg + Z0
51
t=0 +
V
+ + +
I
V0
Z0 VL RL = Z0
- -
V0
t
t’ = l/ν
b)
Figure 23. a) A general transmission line with a transient source and b) the step function used to model
the switch
the length of the TL and on the ordinate (y-axis) we measure time. As the wave is launched down
the TL at time zero and at position zero. This point is shown in the bottom left-hand corner of
Fig. 24 a). The first forward wave propagating down the TL reaches the end of the TL at time
l/ν. Then the wave reflects back to the source and arrives at the source at time 2l/ν...and so on.
This behavior results in the back and forth illustration in Fig. 24 a). If we consider the voltage at
z = 3l/4 along the TL then we have the voltage response w.r.t. time shown in Fig. 24 b). We can
see that the voltage slowly converges to V0 RL /(Rg + RL ).
Now say the line was initially charged. This then results in the problem described in Fig. 25.
We know the voltage to the left of the dotted line is less than V0 . ⇒ V1+ has a sign change. Also,
52
t
4l/ν
V2- = Γ
Γ V1+
2
g L
13l/4ν
3l/ν
+ 11l/4ν
+ ΓgΓLV1
V2 =
2l/ν
V1- = Γ +
LV
1
5l/4ν
l/ν
+ 3l/4ν
V1
z
0 3l/4 l
V3/4 a)
V0RL
Rg+RL + - + -
+ - + V1+V1+V2+V2
V1+V1+V2
+ -
V1+V1
+
V1
t
3l/4ν 5l/4ν 11l/4ν 13l/4ν 19l/4ν 21l/4ν
b)
Figure 24. a) Voltage reflections on the transmission line and b) line voltage on the transmission line over
time
t=0 +
V1
V0
+ V0+V+1 +
I1
Rg VR Z0 (open)
-
IR
{
53
IR = −I1+ = −V1+ /Z0 . ⇒ VR = V0 + V1+ = IR Rg = −I1+ Rg = −V1+ Rg /Z0 . ⇒
−V0 Z0
V1+ = . (1.95)
Z0 + Rg
54
CHAPTER 2. ELECTROSTATIC FIELDS
In this section we introduce the notation used for the remainder of the notes. Vectors are
denoted as Ā. In rectangular coordinates, Ā = Ax âx + Ay ây + Az âz . From this we define the
√
magnitude of Ā as |Ā| = A2x + A2y + A2z . The unit vector of Ā is defined as
Ā Ax Ay Az
âA = = âx + ây + âz .
|Ā| |Ā| |Ā| |Ā|
z z
A az
y
ay y
ax
x x
Figure 26. Vector definitions.
Ā · B̄ = |Ā||B̄| cos(θAB )
= Ax Bx + Ay By + Az Bz
55
y
A
B
θΑΒ
x
Figure 27. Dot product.
AXB
B
y
A
x
Figure 28. Cross product.
2.1.4. Gradient
Next, we discuss the gradient of a scalar function. If T (x, y, z) is our scalar function. Then
the gradient of T is defined as:
∂T ∂T ∂T
∇T = âx + ây + âz
∂x ∂y ∂z
56
where the gradient operation is referred to as “del” and the operation is in cartesian coordinates.
After the operation of del on a scalar function T we have a vector whose magnitude is equal to the
maximum rate of change of the physical quantity (i.e., volt) per unit distance and whose direction
is along the direction of maximum increase. Now dot ∇T with a unit vector âl , ∇T · âl , then we
have a directional derivative.
Aρ cos ϕ sin ϕ 0 Ax = z
A = − sin ϕ cos ϕ 0 A = −x = −ρ cos ϕ .
ϕ y
Az 0 0 1 Az = y = ρ sin ϕ
⇒
57
determined by Colonel Charles Coulomb. The result was the following:
Q1 Q2
F =k (2.1)
R2
where Q1 and Q2 are the positive or negative charges, R is the separation distance between the
charges, F is the force between the charges and k is a constant. An illustration of the problem is
shown in Fig. 29. We also have k = 1/(4πε0 ) where ε0 = 8.854 × 10−12 F/m which is called the
free-space permittivity. Using (2.1) gives
Q1 Q2
F = . (2.2)
4πε0 R2
The units for Q are Coulombs (C), for R are meters (m) and F are Newtons (N). In vector form
we have
Q1 Q2 Q1 Q2
F̄1 = −F̄2 = 2
â21 = − 2
â12 . (2.3)
4πε0 R12 4πε0 R12
Q1 a12 R12 = r2 - r1 Q2
r1 r2
origin
Figure 29. The force between two charges as described by Coulomb’s equation.
Now, fix Q1 and move a test charge Qt in the space around Q1 . This force per-unit charge is
F̄t Q1
= â .
2 1t
(2.4)
Qt 4πε0 R1t
58
Equation (2.4) describes a vector field. This field is called the electric field intensity. Therefore,
F̄t Q1
Ē = = â .
2 1t
(2.5)
Qt 4πε0 R1t
The units of (2.5) is Volts/meter or V/m. Equation (2.5) describes the electric field from a single
point charge at a point anywhere in space at vector R1t . This is described in Fig. 30. This then
gives
Q r̄ − r̄′
Ē(r̄) = . (2.6)
4πε0 |r̄ − r̄′ |2 |r̄ − r̄′ |
P(x‘,y‘,z‘) R= r - r‘ P(x,y,z)
Q E
r‘ r
origin
Figure 30. Electric field due to a single point charge.
∫
ρv (r̄′ )dV ′ r̄ − r̄′
Ē(r̄) = (2.8)
V 4πε0 |r̄ − r̄′ |2 |r̄ − r̄′ |
59
Solution: Using a volume integral we have the following:
∫ π ∫ 0.1 ∫ 4
Q = ρ2 z 2 sin(0.6ϕ)ρdzdρdϕ
0 0 2
∫ π ∫ 0.1 4
3 1 3
= ρ sin(0.6ϕ) z dρdϕ
0 0 3 2
∫ π 0.1
56 1 4
= sin(0.6ϕ)dϕ ρ
3 0 4 0
π
56 100 × 10 − cos(0.6ϕ)
−6
=
3 4 0.6 0
= −777.7 × 10−6 (cos(0.6π) − 1)
= 1.0181mC.
Next, consider the infinite line of electric charge on the z-axis in Fig. 31. An incremental unit
of the electric field in the region around the line can be written as
ρl dz ′ (r̄ − r̄′ )
dĒ = (2.9)
4πε0 |r̄ − r̄′ |3
where r̄ = yây = ρâρ , r̄′ = z ′ âz and r̄ − r̄′ = ρâρ − z ′ âz . This then gives
ρl dz ′ (ρâρ − z ′ âz )
dĒ = . (2.10)
4πε0 (ρ2 + (z ′ )2 )3/2
⇒
ρl ρdz ′
dEρ = . (2.11)
4πε0 (ρ2 + (z ′ )2 )3/2
⇒
∫ ∞
ρl ρdz ′
Eρ = . (2.12)
−∞ 4πε0 (ρ2 + (z ′ )2 )3/2
⇒
ρl
Eρ = . (2.13)
2πε0 ρ
60
z
ρl
dQ = ρldz’
(0,0,z’) aR
R = r - r’
r’
θ
P dEρ
y
r
dEz
dE
x
the electric field anywhere above the sheet of charge can be written as
ρs
Ē = âN . (2.14)
2ε0
Notice that (2.14) has no dependence on ρ. Think of this problem as an infinite uniform light with
no attenuation.
z
ρs
x
Figure 32. Electric field due to an infinite sheet of charge in the y-z plane.
61
2.3.6. Streamlines and sketches of fields
Stream lines are used to help us visualize the field lines from a source. To do this consider
ρl
Ē = âρ .
2πε0 ρ
The field lines from a single point charge are sketched in Fig. 33.
a) b)
c) d)
Figure 33. a) A stream line sketch that is not good; b) and c) two good stream line sketches that show
the intensity of the field and d) the traditional stream line sketch indicating the fields in the region around
the charge.
62
charge we have the following electric flux density:
Q
D̄ = ε0 Ē = âr . (2.15)
4πr2
⇒
∫
ρv dV
D̄ = âr . (2.16)
V 4πR2
Insulating
or dielectric
Metal r=b -Q
material
conducting
spheres +Q
r=a
Electric
flux
lines
Figure 34. Electric field density between two conducting spheres with a dielectric insulator.
charge enclosed by that surface. To understand this law, consider the image in Fig. 35. The total
I
Ψ= D̄s · dS̄. (2.17)
S
Equation (2.17) is adding up all of the flux passing through the Gaussian surface S. Then Gauss’s
law states Ψ = Qenc or
I
D̄s · dS̄ = Qenc . (2.18)
S
63
∆S
Ds,normal Gaussian
surface S
θ ∆S
Q
P
Figure 35. Definition of Gauss’s law with a Gaussian surface.
To demonstrate Gauss’s law, consider a point charge Q in Fig. 36. We know the field from
that point charge is Ē = Q
â
4πε0 r2 r
and D̄ = ε0 Ē. ⇒ D̄ = Q
â
4πr2 r
and at the surface, D̄s = Q
â .
4πa2 r
We
know that dS̄ = r2 sin θdθdϕâr = a2 sin θdθdϕâr . This then gives D̄s · dS̄ = Q
â
4πa2 r
· a2 sin θdθdϕâr =
Q
4π
sin θdθdϕ. This then gives
∫ ∫ ∫ π ∫ 2π
2π π
Q 2π
Q Q
sin θdθdϕ = (− cos θ) dϕ = dϕ = Q.
0 0 4π 0 4π 0 0 2π
Ds
θ dS
a
r=
Q y
x
Figure 36. Example of Gauss’s law.
64
1) D̄s is either normal or tangential to the closed surface everywhere, so that D̄s · dS̄ becomes
either Ds dS or zero.
2) On the portion of the closed surface for which D̄s · dS̄ is not zero, Ds = constant.
If D̄s · dS̄ = Ds dS then, Ds dS is a scalar and Ds can be moved outside of the integration.
2.5.4. Point charge example
Again, consider the point charge Q at the origin in Fig. 36. For this example, we choose the
closed surface S to be a sphere. This then gives D̄s = Ds as a constant over the sphere because
the only component of the electric flux density is the r component. Therefore, if the radius of the
sphere is r then
I I I
Q= D̄s · dS̄ = Dr dS = Dr dS = 4πr2 Dr .
S S S
Q
D̄s = âr . (2.21)
4πr2
Next, consider the line charge ρl extending to infinity in the ±z-direction in Fig. 37. On the
surface of the cylinder, we have D̄ = Dρ âρ normal to the surface. This then gives
I
Q = D̄s · dS̄
∫
cylinder
∫ ∫
= Dρ · dS̄ + 0dS̄ + 0dS̄
∫
sides top bottom
= Dρ 1dS̄
sides
∫ L ∫ 2π
= Dρ ρdϕdz
0 0
= Dρ 2πρL.
65
Therefore,
Q
Ds = Dρ = .
2πρL
ρl
(line charge) Gaussian
surface
L
Figure 37. Calculating the electric field due to a line charge using Gauss’s law.
2.6. Divergence
H
From the previous section we have that Gauss’s law states that S D̄ · dS̄ = Qenc . To illustrate
H
divergence, consider the illustration in Fig. 38. To evaluate S D̄ · dS̄ in Fig. 38 we have
I ∫ ∫ ∫ ∫ ∫ ∫
D̄ · dS̄ = + + + + + .
S f ront back lef t right top bottom
∫
Since the surface area is very small we have f ront
= D̄f ront · ∆S̄f ront = Dx,f ront ∆y∆z. Since P is
⇒
∆x ∂Dx
Dx,f ront = + Dx0 .
2 ∂x
⇒
∫ []
∆x ∂Dx
Dx,f ront dS = + Dx0 ∆y∆z.
f ront 2 ∂x
66
Similarly for the back of the cube in Fig. 38 we have
∫ ( )
∆x ∂Dx
= − Dx0 + ∆y∆z.
back 2 ∂x
⇒
∫ ∫
∂Dx
+ = ∆x∆y∆z,
f ront back ∂x
∫ ∫
∂Dy
+ = ∆x∆y∆z
lef t right ∂y
and
z
Flux density in
the cube
D=Dx0 ax+Dy0 ay+Dz0 az
P(x,y,z) (middle of
∆z the cube)
∆x
∆y
x
Figure 38. Illustration of divergence.
∫ ∫
∂Dz
+ = ∆x∆y∆z.
top bottom ∂z
I ( )
∂Dx ∂Dy ∂Dz
D̄ · dS̄ = + + ∆x∆y∆z
∂x ∂y ∂z
S
( )
∂Dx ∂Dy ∂Dz
= + + ∆V
∂x ∂y ∂z
= Qenc . (2.22)
67
( )
What (2.22) says is that the charged enclosed in a unit volume ∆V is equal to ∂Dx
∂x
+ ∂D
∂y
y
+ ∂D
∂z
z
∆V .
Now let ∆V → 0 and divide both sides of (2.22) by ∆V gives
( )
Q ∂Dx ∂Dy ∂Dz
lim = + + = ρV .
∆V →0 ∆V ∂x ∂y ∂z
∂Dy
This leads to the following notation: ∇ · D̄ = div D̄ = ∂Dx
∂x
+ ∂y
+ ∂Dz
∂z
which gives
∇ · D̄ = ρV . (2.23)
Equation (2.23) is known as one of Maxwell’s equations. Equation (2.23) is very useful because it
is a nice relation between the charge in a region and the fields being radiated from that charge in
the region. As long as (2.23) is satisfied, then the solution to (2.23) is unique in that region.
everywhere.
2.6.2. Solving for volume charge in a region
2 2 2 y(−2) 2
Find ρV for D̄ = 4xy
z
âx + 2xz ây − 2xz2 y âz . Solving gives ∇ · D̄ = 4y
z
+ 0 − 2x z3
= 4y
z
+ 4xz3 y =
4y
z3
(z 2 + x2 ).
I ∫
D̄ · dS̄ = ∇ · D̄dV. (2.24)
S V
2.8.1. Work
Suppose we want to move a charge Q a distance dL in an Ē-field. We know that Ē = F̄ /Q
68
⇒ F̄ = ĒQ where F̄ is the force on Q arising from the electric field. Thus, to move the charge we
have the applied force Fapplied = −QĒ · âL . We know our work is force times distance = energy
expended. ⇒ the differential work dW is dW = −QĒ · dL̄ = −QĒ · âL dL. ⇒
∫ f inal
W = −Q Ē · dL̄. (2.25)
initial
Let Ē = 1
z2
(8xyzâx + 4x2 zây − 4x2 yâz ) V/m. Find the differential amount of work done to
move a 6 nC charge a distance of 2 µm from P(2,-2,3) in the âL = − 67 âx + 37 ây + 27 âz direction.
−32
Solution: The electric field at P is Ē(2, −2, 3) = 3
âx + 16
3 y
â + 32
â .
9 z
⇒ Ē · âL = 64
7
+ 16
7
+ 64
63
=
112
9
= 12.44. ⇒ dW = −QĒ · âL · dL = −6nC · 12.44 · 2µm = −149 fJ.
2.8.3. Line integral example
A line integral is used if the path is continuous. For this example let the electric field be
Ē = yâx + xây + 2âz . Find the work to carry 2 C from B(1,0,1) to A(.8,.6,1) along the arc of the
circle x2 + y 2 = 1 and z = 1.
Solution: We will work in rectangular coordinates for the problem in Fig. 39. Note that dL̄ =
∫ A
W = −Q Ē · dL̄
B
∫ A
= −2 (yâx + xây + 2âz )(dxâx + dyây + dzâz )
B
∫ .8 ∫ .6 ∫ 1
= −2 ydx − 2 xdy − 4 dz.
1 0 1
√ √
Now we need to move along a circular path. This then gives y = 1 − x2 and x = 1 − y 2 .
This then gives
∫ .8 √ ∫ .6 √
W = −2 1 − x dx − 2
2 1 − y 2 dy − 0 = −.96J.
1 0
69
y
A(.8, .6, 1)
B(1, 0, 1)
∫ .8 ∫ .6 ∫ 1
W = −2 ydx − 2
xdy − 4 dz
1 0 1
∫ .8 ∫ .6
= −2 (−3x + 3)dx − 2 (1 − y/3)dy − 0
1 0
= −.96J.
This shows that the work done is independent of the path taken in any electrostatic field.
2.8.4. Potential
Now define the potential difference V as the work done in moving a unit positive charge from
∫ f inal
V =− Ē · dL̄. (2.26)
initial
70
Solution: From before, we know the electric field from a point charge is
Q
Ē = Er âr = âr .
4πε0 r2
∫ A ∫ rA ( )
Q Q 1 1
VAB = − Ē · dL̄ = − dr = − .
B rB 4πε0 r2 4πε0 rA rB
Note that when we use the potential notation, we have a reference point. Usually “ground” or
“infinity”.
2.8.6. Potential field of a line charge example
From before, we know that the potential from a point charge is
∫ A ∫ rA ( )
Q Q 1 1
VAB = − Ē · dL̄ = − dr = − = VA − VB .
B rB 4πε0 r2 4πε0 rA rB
We know work is independent of the path between the two points. This implies VAB is independent
Q
V = .
4πε0 r
∫
ρL (r̄′ )
V (r̄) = dL′ . (2.28)
l 4πε0 |r̄ − r̄′ |
71
and
∫
ρv (r̄′ )
V (r̄) = dV ′ . (2.30)
V 4πε0 |r̄ − r̄ |
′
the x-y plane and centered about the z-axis. Since this is a line charge, we need to evaluate a line
integral over the charge. This then gives the following:
∫ 2π
ρL a′
V (r̄) = dϕ′
0 4πε0 |r̄ − r̄′ |
√
where dL′ = a′ dϕ, r̄ = zâz and r̄′ = a′ âρ . This then gives |r̄ − r̄′ | = a′2 + z 2 . ⇒
∫ 2π
ρ a′ ρ a′
V (r̄) = √L dϕ′ = √L .
0 4πε0 a′2 + z 2 2ε0 a′2 + z 2
This is the potential at infinity. For a potential difference in general we have VAB = VA − VB =
∫A
− B Ē · dL̄. Now say A → B then VAB → 0. ⇒
∫ A
− Ē · dL̄ → 0
B
I
Ē · dL̄ = 0 (2.31)
72
z
(0,0,z)
y
a’
ρl
φ’
dL’ = a’ dφ’
∆L
73
the max rate of change of V . To illustrate this consider the equipotential surfaces in Fig. 42. We
want the electric field at point P. The direction of max change is in the direction of âN pointed
towards the higher potentials. ⇒
dV dV
Ē = − âN = − âN
dL max dN
dV dV dV dV
− âN = − âx,N − ây,N − âz,N = −∇V
dN dx dy dz
Ē = −∇V. (2.32)
+60 +50
+70 +40
+80 +30
Direction of
max change
74
2.8.11. The dipole
The electric dipole has the geometry shown in Fig. 43 where d << r. Now, the potential at
( )
R2 −R1
P w.r.t. ∞ is V = 4πε0 R1 − 4πε0 R2 = 4πε0 R1 R2 . We want to simplify this expression. To do
Q Q Q
this, choose P such that R1 and R2 are much much greater than d. This is shown in Fig. 44. This
results in the following approximations: R1 ≈ R2 ≈ r and R2 − R1 = d cos θ. This then simplifies
( )
Q d cos θ Q d cos θ
V = = .
4πε0 rr 4πε0 r2
P
R1
+Q
θ r
d R2
x
-Q
Figure 43. A dipole.
Ē = −∇V
( )
∂ Q d cos θ
= − âr
∂r 4πε0 r2
( )
1 ∂ Q d cos θ
− âθ
r ∂θ 4πε0 r2
( )
1 ∂ Q d cos θ
− âϕ
r sin θ ∂ϕ 4πε0 r2
Q d cos θ Q d sin θ
= 3
âr + âθ + 0.
2πε0 r 4πε0 r3
75
⇒
Qd
Ē = [2 cos θâr + sin θâθ ].
4πε0 r3
θ R1
P
+Q
r
d
y
R2
{
x R2 - R1 = d cos θ
-Q
Figure 44. Far-field approximation of a dipole.
z
0.4
0.6
0.8
1
0 0
-1
-0.8
-0.6
-0.4
Figure 45. Equipotential surfaces of a dipole.
76
CHAPTER 3. PROPERTIES OF DIELECTRICS AND
CONDUCTORS
3.1. Current
3.1.1. Definition
To introduce current we first consider the illustration in Fig. 46. Current is defined as the
rate of movement of charge passing a given reference point of 1 C/s. ⇒
dQ
I= . (3.1)
dt
The unit for current is Amperes or A. The current density J¯ describes the current distribution over
the reference plane. The unit for current density is A/m2 and ∆I = J¯ · ∆S̄. Thus, in general we
have
∫
I= J¯ · dS̄. (3.2)
S
∆S
Q Q
Q -Q
Reference plane
Figure 46. Current definition.
∆Q ∆x
∆I = = ρV ∆S .
∆t ∆t
77
⇒
∆x
lim ρV ∆S = ρV ∆SVx = ∆I
∆t→0 ∆t
J¯ = ρV V̄ (3.3)
where V̄ is the general velocity vector. Note that (3.3) shows that charge in motion results in a
current called the convection current.
∆Q = ρv∆V z z
∆Q = ρv∆V
∆S
∆S
y y
∆L ∆x
∆L
x x
Figure 47. Convection current.
Consider a region bounded by a closed surface S. The current through the surface S is
I
I= J¯ · dS̄
S
where this outward flow of positive charge must be balanced by a decrease of positive charge. ⇒
I
dQ
I= J¯ · dS̄ = −
S dt
where Q is the charge inside a closed surface and the negative sign represents an outward traveling
78
current. Using the divergence theorem
I ∫
D̄ · dS̄ = ∇ · D̄dV
S V
gives
I ∫ ∫
J¯ · dS̄ = ¯ = − dQ = −
∇ · JdV
∂ρv
dV. (3.4)
S V dt V ∂t
This then gives the current continuity equation in point form as:
∂ρv
∇ · J¯ = − . (3.5)
∂t
where Q = −e. The charges then collide with each other and the conductor and eventually results
in a drift velocity vd . ⇒ vd = −µe Ē where −µe is the mobility of an electron. Substituting this
into (3.3) gives J¯ = −ρv µe Ē. This then gives (Fig. 48)
J¯ = σ Ē (3.6)
J
Area = S
L
Figure 48. Current in metallic conductors.
Now, assume the current is uniform throughout the cylindrical region in Fig. 48. We know
79
∫ ∫a
that I = J¯ · dS̄ = JS and Vab = −
S V b
Ē · d¯l = −ELba = ELab or V = EL. ⇒ J = I
S
= σE = σ VL
⇒V = L
σS
I = RI where
L
R= . (3.7)
σS
Solution: We know that the diameter = 0.05808 in 2.54 cm/1 in = 1.4752 mm. ⇒ S = πr2 =
1.7092µm2 . We also have L = 5280 ft = 1.6093 km. ⇒
1.6093km
Req = = 16.23Ω
5.8 × 107 S/m1.7092µm2
static charges. To do this, consider the problem defined in Fig. 49. We have Dt = 0 otherwise
the flux would change the charge at the surface and then the problem would not be static. For the
normal component, consider the pill box in Fig. 49. For an incremental ∆S we have from Gauss’s
law:
I
Qenc = D̄ · dS̄.
S
⇒
∫ ∫ ∫
+ + = Qenc .
top bottom sides
80
H
For the Ē-field we consider the loop in Fig. 49 and we know: Ē · dL̄ = 0. ⇒
∫ b ∫ c ∫ d ∫ a
+ + + = 0.
a b c d
∫d
In the conductor we have Ē = 0 which gives c
= 0. Now let ∆W → 0 and ∆h → 0. ⇒
Et ∆W − EN 12 ∆h + EN 12 ∆h = 0. ⇒ Et ∆W = 0. ⇒
Et = 0.
Dt = E t = 0 (3.8)
and
DN = ε 0 E N = ρ s . (3.9)
Or in general we have
âN × Ē = 0 (3.10)
and
âN · D̄ = ρs . (3.11)
I
Ē · dL̄ → Etan1 ∆W − Etan2 ∆W = 0.
81
Free-space
EN E
a N
∆S ∆W b
D a
DN
{ ∆h ∆h ∆h
c
Et
Dt d ∆W
ρenc Loop
Pill box
(Gaussian surface)
Conductor
Figure 49. Boundary conditions on a conductor.
Dielectric for
region 1
ε1
EN E
aN
∆S ∆W b
D a
DN
{ ∆h ∆h ∆h
c
Et
Dt d ∆W
ρenc Loop
82
This then gives
Etan1 = Etan2 . (3.12)
I
Qenc = εĒ · dS̄ → DN 1 ∆S − DN 2 ∆S = ∆Q = ∆Sρs .
Or in general we have
âN × (Ē1 − Ē2 ) = 0 (3.14)
and
is the absolute potential of a charge. Therefore, at the boundary of the conducting plane on the
+Q −Q +Q −Q
right side of Fig. 51 we have V = 4πε0 r+
+ 4πε0 r−
. But r+ = r− = r giving V = 4πε0 r
+ 4πε0 r
=0V
along the boundary. Other examples of image theory are shown in Fig. 52.
z z
+Q +Q
d d r+
(Image Theory)
y V=0 y
Conducting plane
with V = 0.
-d r
-
-Q
Figure 51. Illustration of image theory.
83
z z
ρ L ρ L
d d
(Image Theory)
V=0 V=0
y y
-d
−ρ
L
z z
IL IL
IL’ IL’
d d
(Image Theory)
V=0 V=0
y y
-d
IL’
IL
Figure 52. Other examples of image theory.
84
3.5. Capacitance
3.5.1. Definition
In this section we introduce capacitance. To do this, we first consider the problem defined in
Fig. 53. Two conductors are embedded in a uniform dielectric. Conductor M2 carries a total positive
charge +Q and conductor M1 caries an equal negative charge −Q. The electric flux travels from
M2 to M1 . ⇒ there exists a potential Vo between the two conductors. Now define the capacitance
between the conductors as
Q
C= , (3.16)
Vo
or in general H
εi Ē · dS̄
C= ∫+ . (3.17)
− − Ē · dL̄
We can see in equations (3.16) and (3.17) that C is independent of Vo and Q. C is only dependent
+ +
+ +
+ M2 +
+
D + +
+
+ ++
- - -
- - -
- M1 -
- - - - - ε
- -
(dielectric)
which is the B.C. on the lower plane. On the upper plane we have the B.C. DN = −Dz . This then
85
gives
∫ lower ∫ 0
ρs ρs
Vo = − Ē · dL̄ = − dz = d.
upper d εi εi
If the surface area of the plane is finite with area S then Q = ρs S. This then gives
Q εS
C= = . (3.18)
V d
Note that equation (3.18) neglects fringing. The energy stored in a capacitor is
∫
1
WE = εE 2 dV
2 vol
1
= CV 2
2 o
1
= QVo
2
1 Q2
= . (3.19)
2C
εi
E
Bottom plate
z=0 y
+ρs
Figure 54. Capacitance between two parallel plates.
3.5.3. Capacitance between two parallel plates with two dielectrics example
In this example, we want to determine C for the structure shown in Fig. 55. We know that
Q
C= Vo
. Assume a potential of Vo between the plates and that E1 and E2 are uniform. This then
gives Vo = E1 d1 + E2 d2 . At the dielectric interface we know DN 1 = DN 2 (ρs = 0). This then gives
86
ε1 E1 = ε2 E2 . ⇒ Vo = E1 d1 + ε1 E1 d2
ε2
= E1 [d1 + ε1
d ].
ε2 2
This then gives
Vo
E1 = .
d1 + εε12 d2
Vo Vo
ρs1 = ε1 E1 = ε1 = .
d1 + εε12 d2 d1
ε1
+ d2
ε2
Q ρs S S Vo 1 1 C1 C2
C= = = d1
= = 1 1 = .
Vo Vo Vo ε1
+ dε22 d1
ε1 S
+ d2
ε2 S C1
+ C2
C1 + C2
z
Upper plate
Surface area S
d2 ε2 Dielectric between
+ E2 d
Vo the plates
- ε1
E1 d1
Lower plate
Figure 55. Capacitance between two parallel plates with two dielectrics.
87
This then gives
ρv
∇ · ∇V = −
ε
∂ 2V ∂2V ∂ 2V
= + +
∂x2 ∂y 2 ∂z 2
= ∇2 V.
Equation (3.20) is called Poisson’s equation. Thus, if ρv = 0 then we have Laplace’s equation:
∇2 V = 0. (3.21)
In this section we show that the solutions to (3.20) and (3.21) are unique. Assume that V1
and V2 are solutions to Laplace’s equation or are a solution to (3.21). This then gives ∇2 V1 = 0 and
and V2 satisfy the B.C. of a problem. Therefore, there exists only one unique solution to the space.
3.8. Parallel Plate Capacitance Example Using Laplace’s Equation
For this example, we will consider the problem defined in Fig. 56. We have two parallel plates
each with a surface area S separated by a distance d. For this example, assume that the plates
are close enough such that if a voltage is applied between the top and bottom plate the resulting
electric field is in the z-direction only. To compute the equivalent capacitance we need to evaluate
C = Q/Vo where Vo is the voltage across the plates and Q is the resulting charge on each plate from
this voltage. For the dielectric region between the plates, we need to satisfy Laplace’s equation.
This then means we have a solution in the following form: V (z) = Az + B. We will determine the
coefficients A and B by enforcing the B.C. on each plate. At z = 0 the voltage is V = 0 and at z =
d the voltage is V = Vo . Therefore, V (0) = B = 0 and V (d) = Ad = Vo . This then gives A = V o/d.
88
Vo
This then gives the potential between the plates as V (z) = d
z. Now we know the potential in
the dielectric. Taking the gradient of V (z) gives Ē = −∇V = − ∂z
∂ Vo
d
zâz = − Vdo âz , which is the
electric field in the dielectric. ⇒ D̄ = εĒ = −ε Vdo âz = DN . At the plate at z = 0 we have the B.C.
∫
ρs = DN = −ε Vdo âz . This then gives Q = S − εVdo dS = −ε VodS . Therefore, C = VQo = |Q|
Vo
= εS
d
.
z
Upper plate
Surface area S
ε
+
Vo d Dielectric between
-
the plates
Lower plate
Figure 56. Computing the capacitance between two parallel plates of finite size.
In this problem we calculate the electric field between the two plates shown in Fig. 57 in the
( )
1 ∂ ∂V 1 ∂2V ∂ 2V 1 ∂ 2V
ρ + 2 2 + = 0 + + 0 = 0.
ρ ∂ρ ∂ρ ρ ∂ϕ ∂z 2 ρ2 ∂ϕ2
For this second partial differential equation we have the following solution V = Aϕ + B. Applying
the B.C. we have V = 0 at ϕ = 0 ⇒ B = 0. Also, at ϕ = α we have Aα = Vo . ⇒ A = Vo
α
. This then
gives the potential between the plates as V = Vo
α
ϕ. Therefore, Ē = −∇V = − ρ1 ∂V â = − ρ1 Vαo âϕ .
∂ϕ ϕ
89
V = Vo
φ=α V=0
φ=0
Insulating
gap
φ
α
z Infinite plates
(both)
Figure 57. Computing the electric field between two plates of finite size in cylindrical coordinates.
90
CHAPTER 4. MAGNETOSTATIC FIELDS
In this section, we introduce Magnetostatic (or magnetic) fields. To illustrate magnetic fields
we need to consider the differential DC current element in free-space shown in Fig. 58. To relate
the magnetic field to the DC current we start with the Biot-Savart law. The Biot-Savart law is
written in the following manner:
I1 dL̄1 × âR12
dH̄2 = 2
(4.1)
4πR12
where H̄ is the magnetic field intensity with units A/m and R12 = |R̄12 |. Thus, in general and in
integral form we have
I
I1 dL̄1 × âR
H̄ = . (4.2)
4πR2
Point 1
{
dL1
aR 12
I1 R12
P
(Point 2)
Figure 58. Illustration of the Biot-Savart law.
Surface
Js
{
We can also describe the Biot-Savart law in terms of current density J¯s . The total current
91
flowing on the surface in Fig. 59 is IdL̄ = J¯s dS = J¯v dV . ⇒
∫
J¯s × âR
H̄ = dS (4.3)
S 4πR2
and
∫
J¯v × âR
H̄ = dV. (4.4)
vol 4πR2
ρâρ − z ′ âz
âR12 = √ .
ρ2 + z ′ 2
∫ ∞
Idz ′ âz × (ρâρ − z ′ âz )
H̄2 =
−∞ 4π(ρ2 + z ′ 2 )3/2
∫ ∞
I ρâϕ
= dz ′
4π −∞ (ρ + z ′ 2 )3/2
2
∫
Iρâϕ ∞ dz ′
=
4π −∞ (ρ2 + z ′ 2 )3/2
I
= âϕ . (4.5)
2πρ
Equation (4.5) is the magnetic field due to an infinite line of current segment.
4.2. Magnetic Flux Density
Next, we define the magnetic flux density as
B̄ = µ0 H̄ (4.6)
where µ0 = 4π × 10−7 H/m. µ0 is referred to as the permeability of free-space. The units for B̄ are
92
z
Infinite
current I Differential current element
Point 1 { dL
aR
R12
z’ az
y
ρ aρ
P
x (Point 2)
W b/m2 or T (Tesla).
The problem is a current segment with infinite length. The current has a magnitude Ienc . Ampere’s
law states the following:
I
Ienc = H̄ · dL̄. (4.7)
l
Or in words, the integral around a closed path of the magnetic field is equal to the total current
enclosed by the closed path of integration.
path of integration
with radius a.
I enc
93
Fig. 61 is on the z-axis and define our path of integration to be along a circle enclosing the current
with radius a, then we have the following expression for the magnetic field from the current:
I ∫ 2π ∫ 2π
H̄ · dL̄ = Hϕ ρdϕ = Hϕ ρ dϕ = Hϕ 2πρ = I.
0 0
Therefore,
I
Hϕ = .
2πρ
In this example, we calculate the magnetic field in and around the coax in Fig. 62. For this
example, we assume that I is uniformly distributed in the conductors. +I is on the center conductor
and −I is on the outer conductor. Notice H is not a function of ϕ or z. For a ≤ ρ < b we have
2 2 2
Hϕ = I
2πρ
. For 0 < ρ ≤ a we have Ienc = I aρ2 = I πa
πρ ρ
2 and 2πρHϕ = I a2 where 2πρ is the path of
( ) )(
πρ2 − πb2 ρ2 − b2
2πρHϕ = I − I =I −I 2 .
πc2 − πb2 c − b2
Finally, for ρ >C we have Hϕ = 0 because Ienc = 0. A plot of the magnetic field is shown in Fig.
63.
4.4. Curl
In this section we define the curl operation. To do this consider the illustration in Fig. 64.
H
For this, we want to evaluate H̄ · dL̄ = I = Jz ∆x∆y. On a side we have (H̄ · ∆L̄)1−2 = Hy,1−2 ∆y.
We want to write this expression in terms of Hyo . We know that
94
c a
I Infinitely long
I/2πa
I/4πa
0 a 2a 3a = b 4a = c
Figure 63. Plot of the magnetic field along the coax.
95
This then implies
( )
∂Hy 1
Hy,1−2 = Hyo + ∆x . (4.8)
∂x 2
I ( )
∂Hy ∂Hx
H̄ · dL̄ = − ∆x∆y = Jz ∆x∆y.
∂x ∂y
z
H = Ho = Hxo ax + Hyo ay + Hzo az
4 3
∆x
1 ∆y 2
x
Figure 64. Curl of the magnetic field in a region.
H
H̄ · dL̄ ∂Hy ∂Hx
lim = − = Jz .
∆x→0,∆y→0 ∆x∆y ∂x ∂y
H
H̄ · dL̄ ∂Hz ∂Hy
lim = − = Jx
∆y→0,∆z→0 ∆y∆z ∂y ∂z
96
and H
H̄ · dL̄ ∂Hx ∂Hz
lim = − = Jy .
∆z→0,∆x→0 ∆z∆x ∂z ∂x
H
H̄ · dL̄
(CurlH̄)N = lim = ∇ × H̄. (4.10)
SN →0 ∆SN
∇ × H̄ = J.
¯ (4.11)
H
Also, Ē · dL̄ = 0 or
∇ × Ē = 0. (4.12)
Finally,
I ∫ ∫
H̄ · dL̄ = (∇ × H̄) · dS̄. (4.13)
l S
From before we defined the magnetic flux density as B̄ = µH̄. The magnetic flux is defined
as
∫
Φ= B̄ · dS̄. (4.14)
S
I
B̄ · dS̄ = 0.
S
97
Using the divergence theorem we have
∇ · B̄ = 0. (4.15)
∇ × Ē = 0 (4.16)
∇ × H̄ = J¯ (4.17)
∇ · D̄ = ρv (4.18)
∇ · B̄ = 0. (4.19)
Equations (4.16)-(4.19) are referred to as Maxwell’s static equations. The integral form of Maxwell’s
static equations are written in the following manner:
I
Ē · dL̄ = 0 (4.20)
I l
∫
H̄ · dL̄ = I = J¯ · dS̄ (4.21)
I l
∫
S
D̄ · dS̄ = Q = ρV dV (4.22)
IS V
B̄ · dS̄ = 0. (4.23)
S
If the radius is a = 1.0 mm, the conductor is on the z-axis and Iz = 20 A, find:
a) Hϕ at ρ = 0.5 mm
b) Bϕ at ρ = 0.8 mm
c) Total magnetic flux per unit length inside the wire
d) Total flux for ρ < 0.5 mm
98
Ienc
a) From before we had Hϕ = 2πρ
. This then gives
(0.5)2 20
Hϕ = = 1591.5A/m.
(1.0) 2π ∗ 0.5mm
2
Ienc
b) Again from before we had Hϕ = 2πρ
. This then gives
(0.8)2 20
Hϕ = = 2546.4A/m.
(1.0)2 2π ∗ 0.8mm
Computing the magnetic flux density then gives Bϕ = µ0 1Hϕ = 3.2 mT.
c) In general the magnetic flux density in the wire can be written as
ρ2 I µ0 Iρ
Bϕ = 2
µ0 = .
a 2πρ 2πa2
∫ 1m ∫ 1mm
µ0 Iρ µ0 20 1 20µ0 5µ0
Φ= 2
dρdz = 2
(1mm)2 = = = 2 µW b/m.
0 0 2πa 2π(1mm) 2 2π2 π
∫ 1m ∫ .5mm
µ0 Iρ 20µ0 1 2 5µ0 (.5mm)2
Φ= dρdz = (.5mm) = = 0.5 µW b/m.
0 0 2πa2 2π(1mm)2 2 (1mm)2
Next, to compute the TOTAL flux, multiply the previous result by 1 m. This then gives Φ =
0.5 µW b.
e) Computing the total flux for ρ > 1 mm leads to the following result:
∫ ∫ ρ′
ρ′
1.0m
10µ0 10µ0
Φ = lim dρdz = lim ln(ρ) = ∞.
′ ρ →∞ πρ ρ′ →∞ π
0 1.0mm 1mm
In a similar manner to the electric scalar potentials, we define scalar vector and magnetic
99
potentials. To do this, assume there exists a scalar magnetic potential denoted as Vm where
H̄ = −∇Vm . (4.24)
H̄ = −∇Vm (4.25)
∇ · B̄ = µ0 ∇ · H̄
= 0
= µ0 ∇ · (−∇Vm ).
∇2 Vm = 0 (4.26)
where again J¯ = 0.
I
= −∇Vm
2πρ ϕ
1 ∂Vm
= − .
ρ ∂ϕ
100
Therefore, when ϕ = 0 then Vm = 0. However, if ϕ = 2π then Vm = −I. This illustrates that
Vm is not a single valued function. This is because as ϕ rotates around the z-axis a different value
of current is enclosed after each rotation. In summary, if Ienc = 0 then we have a single valued
function. This then gives
∫ a
Vm,ab = − H̄ · dL̄ (4.27)
b
P(ρ,π/4,0)
φ
Iout x
ρ=a
ρ=b
ρ=c
Figure 65. Coax for the scalar potential example.
1
H̄ = ∇ × Ā (4.29)
µ0
and
1
∇ × J¯ = ∇ × ∇ × Ā. (4.30)
µ0
101
Equation (4.31) can be used in general to compute the magnetic vector potential from any current
distribution. This expression results in a purely mathematical result and cannot be measured.
However, once Ā is computed, the magnetic field can then be computed using (4.29). Finally, for
surface and volume currents we have the following:
∫
µ0 J¯s dS
Ā = (4.32)
S 4πR
and
∫
µ0 J¯v dV
Ā = . (4.33)
V 4πR
if the charge Q is moving with velocity v̄ in a region with flux density B̄, then it has been
experimentally determined that the force on the charge Q can be computed as F̄ = Q(v̄ × B̄).
Consider the charge Q = 18 nC traveling with velocity v̄ = 5 × 106 (0.6âx + 0.75ây + 0.3âz )
m/s in the presence of a magnetic flux density B̄ = −3.0âx + 4.0ây + 6.0âz mT. Compute |F̄ |.
Solution:
Using (4.34) with Ē = 0 we first compute v̄ × B̄ in the following manner:
0.6 0.75 0.3
v̄ × B̄ = ∗ 5 × 106 ∗ 10−3
−3.0 4.0 6.0
= 16.5 × 103 âx + 22.5 × 103 ây + 23.25 × 103 âz .
102
This then gives
|F̄ | = Q|v̄ × B̄| = 653µN.
flux density is derived. To compute this, the force equation (4.34) with Ē = 0 is written in
differential form in the following manner: dF̄ = dQ(v̄×B̄). Next, we define J¯ = ρv v̄ and dQ = ρv dV .
This then gives
= (J¯ × B̄)dV.
and
∫
F̄ = J¯v × B̄dV, (4.35)
V
∫
F̄ = J¯s × B̄dS (4.36)
S
and
I
F̄ = IdL̄ × B̄ (4.37)
or
I
F̄ = −I B̄ × dL̄ (4.38)
103
Therefore, for a straight conductor we have
F̄ = I L̄ × B̄. (4.39)
For this example the force from the current on the conducting square loop in Fig. 66 is
computed. Using (4.5), the magnetic field from the current on the wire can be computed as
H̄ = 15/2πxâz A/m. This then gives B̄ = µ0 H̄ = 3/xâz µT . Next, using (4.38) we get
I
F̄ = −I B̄ × dL̄
[∫ ∫ 2 ∫ 1 ∫ 0 ]
3
1 1 1 1
= −2 × 10−3 ∗ 3 × 10−6 âz × dxâx + âz × dyây + âz × dxâx + âz × dyây
1 x 0 3 3 x 2 1
[ 3 2 1 0 ]
−9 y
= −6 × 10 ln(x)ây + (−âx ) + ln(x)ây + y(−âx )
3
1 0 3 2
= −8âx nN. (4.40)
z
Free space
15 A y
(1,0,0) (1,2,0)
(3,0,0)
2 mA
(3,2,0)
x
Figure 66. Force on a square loop example.
104
is written using (4.1) in the following manner:
IdL̄1 × âR12
dH̄2 = 2
. (4.41)
4πR12
Also, the differential force on a differential current element is dF̄ = IdL̄ × B̄. Now let dB̄2 be
the differential flux density at point 2 caused by the current element at point 1. This then gives
IdL̄ = I2 dL̄2 . This then gives
The last term in (4.42) represents the force between two differential current segments. Therefore,
in general if we have two long parallel wires with equal constant currents, the integration of (4.42)
F F
I
I
In this section, the magnetic boundary conditions illustrated in Fig. 68 are derived. First
start by applying Gauss’s law to the surface. This then gives
I
B̄ · dS̄ = 0.
S
105
⇒
∆sBN 1 − ∆sBN 2 = 0.
⇒
BN 1 = BN 2 .
I
H̄ · dL̄ = Ienc .
L
Ht1 ∆L − Ht2 ∆L = Js ∆L = I.
⇒
Ht1 − Ht2 = Js .
Dielectric for
region 1
µ1
HN H
BN 1
∆S Ht 1
a b
B
BN
{ ∆h ∆h
∆L
c
∆h
Ht
Bt d Ht2
Loop
Pill box
(Gaussian surface) Dielectric for
region 2
a
N12 BN 2
µ2
Figure 68. Magnetic boundary conditions.
106
4.17. Inductance
In this section we define the inductance (or self-inductance) as the ratio of the total flux
linkages to the current which they link as
NΦ
L= (4.46)
I
where N is the number of turns in the coil and I is the current flowing in the coil.
4.18. Inductance of a Coax Example
In this example, we calculate the per-unit inductance of the coaxial cable in Fig. 69. First,
using the expression in (4.5) the field from the current on the inner conductor can be written as
I
H̄ = Hϕ = .
2πρ
Next, to compute the total flux, the surface S (also shown in Fig. 69) must be defined. Then, using
∫ d ∫ b
I
Φ = µ0 dρdz
2πρ
0 a
( )
µ0 Id b
= ln .
2π a
(top view)
dielectric Z=d
b (σ,ε,µ)
inner H
conductor
Z=0
(σc) I
Surface S
Figure 69. Coax for computing the per-unit inductance.
107
Therefore,
Φ
L =
I ( )
µ0 d b
= ln .
2π a
Or in per-unit notation:
( )
µ0 b
L= ln H/m. (4.47)
2π a
108
CHAPTER 5. TIME-VARYING FIELDS
This law is called Faraday’s law and it states that a time varying flux passing through a closed loop
induces a voltage on that loop. In differential form this can be written as
Binc(t)
Bind(t)
closed path
Figure 70. Illustration of Faraday’s law.
dΦ
Vemf = − (5.1)
dt
dΦ
Vemf = −N . (5.2)
dt
I
Vemf = Ē · dL̄. (5.3)
l
Note that (5.3) is closed and path specific. Next, the induced emf can be written as
I ∫
d
Vemf = Ē · dL̄ = − B̄ · dS̄. (5.4)
l dt S
109
Now applying Stokes Theorem
I ∫
H̄ · dL̄ = (∇ × H̄) · dS̄ (5.5)
l S
to (5.4) gives
∫ ∫
dB̄
Vemf = (∇ × Ē) · dS̄ = − · dS̄. (5.6)
S S dt
∂ B̄
∇ × Ē = − . (5.7)
∂t
to the port and B̄ = Bo (2ây + 3âz ) sin(ωt) is defined in the region around the antenna. For this
example, compute
a) The total flux Φ for N = 1.
b) Vemf for N = 10, Bo = 0.2 T, a = 10 cm and ω = 103 rad/s.
c) Polarity of Vemf for t = 0.
d) I for R = 1 KΩ.
N-turns z B
I
1
+
R a y
V
-
emf
110
Solution:
a) The total flux can be computed as
∫
Φ = 1 B̄ · dS̄
∫ S
= Bo (2ây + 3âz ) sin ωtâz dS
S
dΦ
Vemf = −N
dt
= −3πωN a2 Bo cos ωt
c) At time t = 0, dΦ
dt
> 0 with amplitude Vemf = −188.5 V. Since Φ is increasing , I is from 2 to 1
V2 − V1
I =
R
188.5
= cos(103 t)
1000
= 0.19 cos(103 t).
111
we add a Ḡ term to both sides. This results in
∇ × H̄ = J¯ + Ḡ.
∇ · J¯ + ∇ · Ḡ = 0.
∂ρv
∇ · Ḡ =
∂t
∂
= (∇ · D̄)
∂t ( )
∂ D̄
= ∇· .
∂t
Therefore,
∂ D̄
Ḡ = .
∂t
∂ D̄
∇ × H̄ = J¯ + . (5.8)
∂t
∂ D̄
Equation (5.8) is called Maxwell’s second equation for time-varying fields and the ∂t
term is denoted
as the displacement current density.
5.4. Maxwell’s Equations for Time-Varying Fields
In this section we summarize Maxwell’s equations for time-varying fields in both the differential
and integral forms. First the differential form:
∂ B̄
∇ × Ē = − (5.9)
∂t
∂ D̄
∇ × H̄ = J¯ + (5.10)
∂t
112
∇ · D̄ = ρv (5.11)
∇ · B̄ = 0. (5.12)
We also have
J¯ = σ Ē. (5.13)
I ∫
∂ B̄
Ē · dL̄ = − · dS̄ (5.14)
L S ∂t
I ∫
∂ D̄
H̄ · dL̄ = I + · dS̄ (5.15)
L S ∂t
I ∫
D̄ · dS̄ = ρv dV (5.16)
S V
I
B̄ · dS̄ = 0. (5.17)
S
A few comments can be made about the differential form of Maxwell’s equations. Equation
(5.9) states that if the H̄-field is changing with respect to time at some point, then the Ē-field has
a curl at that point. Next, equation (5.10) states that a time-varying Ē-field generates an H̄-field.
∂Ex
∇ × Ē = ây
∂z
∂ H̄
= −µ0
∂t
∂Hy
= −µ0 ây .
∂t
113
Notice that the previous result is a magnetic field with a y-component. Using Maxwell’s second
equation then gives
∂Hy
∇ × H̄ = âx
∂z
∂ Ē
= ε0
∂t
∂Ex
= ε0 âx .
∂t
Notice in this case that the previous result is an electric field with an x-component. Comparing the
previous two expressions results in the following
∂Ex ∂Hy
= −µ0 (5.18)
∂z ∂t
and
∂Hy ∂Ex
= −ε0 . (5.19)
∂z ∂t
The expressions in (5.18) and (5.19) have two unknowns (Ex and Hy ). In other words, we have two
equations with two unknowns. Next, if we differentiate (5.18) with respect to z we get
∂ 2 Ex ∂ 2 Hy
= −µ0 . (5.20)
∂z 2 ∂t∂z
∂ 2 Hy ∂ 2 Ex
= −ε0 2 . (5.21)
∂z∂t ∂t
∂ 2 Ex ∂ 2 Ex
= µ ε
0 0 . (5.22)
∂z 2 ∂t2
Equation (5.22) is referred to as the wave equation and is one equation with one unknown Ex . Next,
114
we define the propagation velocity as
1
ν=√ ≈ 3 × 108 = c. (5.23)
µ0 ε0
∂ 2 Hy ∂ 2 Hy
= µ ε
0 0 . (5.24)
∂z 2 ∂t2
Note that the expressions for the electric and magnetic fields in the a region MUST satisfy equations
(5.22) and (5.24), respectively. This will result in a solution to the electromagnetic problem in that
region. Next, to solve the wave equations we assume a solution in the following form: Ex (z, t) =
f1 (t − z/ν) + f2 (t + z/ν). It can be shown that
′
Ex (z, t) = |Exo | cos(ωt − ko z + ϕ1 ) + |Exo | cos(ωt + ko z + ϕ2 ). (5.25)
The first term in (5.25) represents the forward traveling wave and the second term represents the
backward traveling (or reflected) wave. We can also define the wave number in free-space as
ω
ko = rad/m. (5.26)
c
Next, the wavelength is the distance in space that ko z shifts by a value of 2π. This then gives
ko z = ko λ = 2π. ⇒
2π
λ= . (5.27)
ko
λ is called the free-space wavelength. We also have the phasor form of the propagating wave (5.25):
Note that (5.28) is written in the time domain. Next, We can also write Maxwell’s equations in the
115
frequency domain in the following manner:
e = jωε0 E
∇×H e (5.29)
e = −jωµ0 H
∇×E e (5.30)
e = 0
∇·E (5.31)
e = 0.
∇·H (5.32)
Equations (5.29)-(5.32) assume a source free region and a sinusoidal steady state source. This then
leads to the generalization of (5.22) which is the wave equation in the frequency domain:
e = −k 2 E
∇2 E e (5.33)
o
√ ex = −ko E
ex .
where k = ω/c = ω µ0 ε0 . Next, to simplify (5.33) we consider the x-component: ∇2 E
e = −jωµ0 H
Now, using Maxwell’s equation ∇ × E e we can obtain H:
e
[ ]
e 1 −jko z ′ jko z
Hy = − (−jko )Exo e + (jko )Exo e
jωµ0
√ √
ε0 −jko z ′ ε0 jko z
= Exo e − Exo e
µ0 µ0
= Hyo e−jko z + Hyo
′ jko z
e .
116
Next, comparing Hyo to Exo gives
√
µ0
Exo = Hyo = η0 Hyo (5.36)
ε0
and
√
′ µ0 ′ ′
Exo =− Hyo = −η0 Hyo (5.37)
ε0
√
where η0 = µ0 /ε0 = 377Ω ≈ 120πΩ. Equation (5.36) represents the forward traveling wave and
(5.37) represents the reflected wave. Therefore, in the far-field we have
e
E
η0 = . (5.38)
e
H
Note: be careful with vector components in (5.38). Also in instantaneous form we have
√ √
ε0 ′ ε0
Hy (z, t) = Exo cos(ωt − ko z) − Exo cos(ωt + ko z). (5.39)
µ0 µ0
In this section we derive the expressions for a propagating wave in a lossy dielectric. To do
this we define the dielectric medium to have a constant µr and εr . This then gives
e = −k 2 E
∇2 E e (5.40)
√ √
where k (the wave number) in a more general form is k = ω µε = ko εr µr , ε = εr ε0 and µ = µr µ0 .
Next, considering the x-component of the electric field gives
ex
d2 E ex .
= −k 2 E (5.41)
dz 2
The solution to (5.41) allows for a complex solution. This then gives jk = α + jβ where α is the
117
attenuation constant and β is the phase constant. This then gives a solution of the following form
Next, the loss in a dielectric can be represented with complex values of ε where ε = ε′ − jε′′ =
ε0 (ε′r − jε′′r ). Similarly, we have µ = µ′ − jµ′′ = µ0 (µ′r − jµ′′r ). This then results in the following
wave number: √
√ √ jε′′
k = ω µ(ε′ − jε′′ ) = ω µε′ 1− .
ε′
√ [√ ( )2 ]1/2
µε′ ε′′
α = Re[jk] = ω 1+ −1 (5.44)
2 ε′
and √ [√ ( )2 ]1/2
µε′ ε′′
β = Im[jk] = ω 1+ +1 . (5.45)
2 ε′
Note that (5.46) changes with material properties. This then implies λ = 2π/β also changes with
material properties. Next, the wave impedance can ge written as
√ √
µ µ 1
η= = √ ( ) (5.47)
ε − jε′′
′ ε ′
′′
1 − jεε′
which is a complex value. Because of this, the Electric and Magnetic fields are no longer in phase.
118
Also notice that
2π 2π 1 c λ0
λ= = √ ′ = √ ′ = √ ′ = (5.48)
β ω µε f µε f µr εr µr ε′r
where λ0 is the free-space wavelength. As long as µr ε′r > 1 then λ < λ0 . Therefore, waves in
dielectrics have a shorter wavelength. This is useful for 1) printed transmission lines, 2) printed
antennas and 3) printed filters.
119
CHAPTER 6. TOPICS ON ELECTROMAGNETIC
COMPATIBILITY
RL RFE
length L
RNE
Rs
Vs
Ground plane
Next, we can derive an equivalent circuit of the coupled TLs. By looking at one small section
of the coupled TLs we will be able to describe the interaction between the two TLs. The equivalent
circuit for one section of coupled TLs is shown in Fig. 73. We want to simplify the analysis of the
coupling between the two TLs in Fig. 73 to the following two cases:
120
Lg/2 Rg/2 Rg/2 Lg/2
Cgr
Lgr/2 Lgr/2
For this analysis, we assume that we have high impedance loads on the coupled TLs in Fig.
72. If this is the case, then the coupling between the TLs is mostly capacitive. This reduces the
equivalent circuit in Fig. 73 to the equivalent circuits shown in Figs. 74 a) and b).
Receptor conductor
Generator conductor
RL RFE
Cgr = (cgr)L +
+ Cgr
cap cap
Cg + Cr Rs Cg RL Cr R = RNE||RFE VNE = VFE
Vin
Rs cap
RNE Vs
VNE
-
-
Vs -
Ground plane
a) b)
Figure 74. a) Capacitive coupling between the coupled transmission lines and b) the equivalent circuit of
the capacitively coupled transmission lines.
Next, we want to write the voltages induced on RN E and RF E by Vs . Evaluating the equivalent
circuit in Fig. 74 b) gives
( 1 )
R|| jωC
VNcap
E = VFcap
E = Vin 1
r
1
R|| jωC r
+ jωCgr
( jω CgrC+C
gr )
r
= Vin 1 . (6.1)
jω + R(Cgr +Cr )
121
Next, for low frequencies we have the following assumption:
1
ω << . (6.2)
R(Cgr + Cr )
E = VF E ≈ jωCgr Vin R.
VNcap cap
(6.3)
Note again that (6.3) is for high impedance circuits. We can also write
( )
Zin
Vin = Vs (6.4)
Zin + Rs
and
( ) [ ( )]
1 1 1
Zin = ||RL || + R|| . (6.5)
jωCg jωCgr jωCr
But, for low frequencies we can approximate a capacitor as an open. This then gives
( )
RL
Vin ≈ Vg(DC) ≈ Vs . (6.6)
RL + Rs
In this section we want to look at the low frequency inductive coupling between the TLs. For
this case we simply open all capacitors in Fig. 73. This then results in the equivalent circuit shown
Next, KVL around the generator gives −Vs + Rs Ig + jωLg Ig + Ig RL − jωIr Lgr = 0. ⇒
Vs jωLgr
Ig = + Ir . (6.8)
(Rs + RL ) + jωLg (Rs + RL ) + jωLg
122
Receptor conductor
Generator conductor
+ ind
RL VFE
RFE
Lg Lr
-
Lgr
Ig Ir
+
Rs ind RNE
V NE
Vs -
Ground plane
Figure 75. The equivalent circuit of the inductively coupled transmission lines.
Next, if we assume that ωLr << RN E + RF E , ωLg << (Rs + RL ) and a weak coupling condition of
Vs
Ig ≈ = Ig (DC) (6.9)
Rs + RL
and
jωLgr
Ir ≈ Ig (DC). (6.10)
RN E + RF E
Then
−RN E
E ≈ jωLgr
VNind Ig (DC) (6.11)
RN E + RF E
and
RF E
E ≈ +jωLgr
VFind Ig (DC). (6.12)
RN E + RF E
123
6.1.4. Equations for both inductive and capacitive coupling
Noting that both V ind and V cap are proportional to jω and Vs we can write the results as
( )
VN E ind cap
= jω MN E + MN E (6.13)
Vs
and
( )
VF E ind cap
= jω MF E + MF E (6.14)
Vs
where
RN E Lgr
MNind
E = , (6.15)
RN E + RF E Rs + RL
RF E Lgr
MFind
E = (6.16)
RN E + RF E Rs + RL
and
RL Cgr
E = MF E = (RN E ||RF E )
MNcap cap
. (6.17)
Rs + RL
RL
Vg (DC) = Vs (6.18)
RL + Rs
and
Vs
Ig (DC) = . (6.19)
Rs + RL
jωLgrIg(DC)
+ +
+
-
- -
Figure 76. The equivalent inductive/capacitive model of the coupled transmission lines.
124
Note that the equivalent circuit model shown in Fig. 76 breaks down for TL lengths longer
than ≈ 0.2λ. Therefore, if the frequency is high enough then a smaller segment of the TL should
be analyzed or a more accurate distributive model of the coupled TL should be used.
Next, we consider the problem in Fig. 72 for the case when the ground plane has finite
conductivity. We represent the loss with a lumped resistor Ro = ro l where l is the length of the
TL and ro is the per-unit loss of the ground plane. The problem with Ro defined is shown in Fig.
77. Ro is called the common impedance and Vo is called the common impedance voltage where
Vo = ro Iref . ⇒ Iref = Ir + Ig ≈ Ig (DC) at low frequencies. ⇒
( )
Ro
Vo ≈ Vs . (6.20)
Rs + RL
This voltage appears on the receptor circuit also. Therefore, this contributes to the near- and
VNciE
= MNciE (6.21)
Vs
and
VFciE
= MFciE (6.22)
Vs
where
RN E Ro
MNciE = (6.23)
RN E + RF E Rs + RL
and
RF E Ro
MFciE = − . (6.24)
RN E + RF E Rs + RL
With this notation, we have the following final coupling equations that include the common impedance
noise voltages:
( )
VN E cap
= jω MN E + MN E + MNciE
ind
(6.25)
Vs
125
and
( )
VF E cap
= jω MF E + MF E + MFciE .
ind
(6.26)
Vs
Rs Ig
+ + Ir
Vs
VNE RNE VFE RFE RL
Vo
+
-
- -
Iref
Ro
Figure 77. Coupled transmission lines with a lossy reference conductor.
jωLgrIg(DC) RoIg(DC)
+ +
+
+
-
-
- -
Figure 78. The equivalent inductive/capacitive model of the coupled transmission lines with a lossy
reference conductor.
6.2. Shielding
6.2.1. Reducing capacitive coupling
In this section we investigate the use of shielding to reduce the coupling between the TLs in
Fig. 72. From (6.3) we can see that we can reduce the capacitive coupling by :
1) reducing frequency
2) reducing RN E and/or RF E
3) reducing Vg (DC) by reducing Vs
4) reducing Cgr by increasing the spacing between the traces and shielding
For this section we will focus on using a shield around the receptor to reduce the capacitive
126
coupling. This problem is defined in Fig. 79 where Vs is the generator (source) open circuit voltage,
Vs,o is the output voltage of the source, Cg is the capacitance from the generator to the reference
conductor, Cgs is the capacitance from the generator to the shield, Cgr is the capacitance from
the generator to the reference conductor, Cs is the capacitance from the shield to the reference
conductor, Crs is the capacitance from the receptor conductor to the shield, Cr is the capacitance
from the receptor to the reference conductor and Vsh is the voltage between the shield and reference
conductor. If we ground the shield on the receptor conductor (at both ends), then we have the
equivalent circuit shown in Fig. 80.
Receptor conductor
Shield +
Generator conductor
RL RFE VFE
Cgs + Cs
-
Cgr Crs Vsh
+ Cg +
Cr
-
Rs Vs,o RNE
VNE
Vs - -
Ground plane
Figure 79. Capacitively coupled transmission lines with a shielded receptor.
+
+ Cgr
cap cap
Rs Cg RL Cgs Cr Crs R = RNE||RFE VNE = VFE
Vs,o
Vg
-
-
Figure 80. The equivalent circuit of the coupled transmission lines with a shielded receptor.
Using circuit analysis, the following expression can be evaluated from the circuit in Fig. 80:
( ( ) )
R|| jω(Crs1 +Cr )
VNcap
E = Vs,o
( ) (6.27)
R|| jω(Crs1 +Cr ) + jωC1 gr
or
127
( ( Cgr ) )
jω Cgr +Crs +Cr
VNcap
E = Vs,o 1 . (6.28)
jω + R(Cgr +Crs +Cr )
At low frequencies we have Vs,o ≈ Vs RL /(RL + Rs ) = Vg (DC) and ω << 1/(R(Cgr + Cr + Crs )).
This then simplifies (6.28) down to
E = VF E ≈ jωRCgr Vg (DC).
VNcap cap
(6.29)
Equation (6.29) is similar to (6.3), except in (6.29) Cgr is much less than Cgr in (6.3).
When looking at shielding for inductive coupling, we need to first consider the expressions in
(6.11) and (6.12). We can see that we can reduce the inductive coupling by :
1) reducing frequency
2) reducing Ig (DC)
3) reducing Lgr by
-
Figure 81. Inductively coupled transmission lines with a shielded receptor.
To understand the best method of shielding for inductive coupling we need to consider the
problem in Fig. 81. We want to write the voltages across the near- and far-end resistors in
128
terms of the inductive coupling and the source voltage. First, KVL around the shield loop gives:
Ish Rsh + (jωLg,sh )Ig − (jωLsh )Ish − (jωLsh,r )Ir = 0. Solving for Ish gives
+jω(Lg,sh Ig + Lsh,r Ir )
Ish = . (6.30)
Rsh + jωLsh
Next, KVL around the receptor loop gives: Ir RN E −(jωLgr )Ig −(jωLsh,r )Ish +(jωLr )Ir +Ir RF E = 0
[ ]
jω Lgr Ig (Rsh + jωLsh ) − jωLsh,r (Lg,sh Ig + Lsh,r Ir )
Ir = . (6.32)
RN E + RF E + jωLr Rsh + jωLsh
[ ]
jω Ig (Lgr Rsh + jωLsh Lgr − jωLsh,r Lg,sh ) + jω(Lsh,r )2 Ir
Ir = . (6.33)
RN E + RF E + jωLr Rsh + jωLsh
Now assume we have a uniform current distribution along the conductors in Fig. 81. This then
gives Lgr ≈ Lg,sh and Lsh ≈ Lsh,r . Rearranging (6.33) and substituting in the above assumptions
gives:
[ ] [ ]
(ωLsh,r )2 jωLgr Rsh
Ir 1 + = Ig . (6.34)
[RN E + RF E + jωLr ][Rsh + jωLsh ] [RN E + RF E + jωLr ][Rsh + jωLsh ]
Next, if we denote the denominator of (6.35) at DEN = [RN E +RF E +jωLr ][Rsh +jωLsh ]+(ωLsh,r )2
129
and assume that Lr ≈ Lsh ≈ Lsh,r then DEN simplifies to DEN = (RN E +RF E )Rsh +jωLsh [RN E +
RF E + Rsh ]. Next, if we assume that RN E + RF E >> Rsh then (6.35) simplifies to
Using Ig ≈ Ig (DC) = Vs
Rs +RL
. This then gives the following expressions for the near- and far-end
voltages in Fig. 81:
[ ][ ]
RN E Rsh
VNind
E = −RN E Ir ≈ (jωLgr Ig (DC)) (6.37)
RN E + RF E Rsh + jωLsh
and
[ ][ ]
−RF E Rsh
VFind
E = RF E Ir ≈ (jωLgr Ig (DC)) . (6.38)
RN E + RF E Rsh + jωLsh
The
[ ]
Rsh
(6.39)
Rsh + jωLsh
term in (6.38) is referred to as the Shield Factor (SF). This term is the same in both (6.37) and
What we see from (6.41) is that the shield has no effect on the coupling caused by the mutual
inductance at low frequencies. The voltage w.r.t. frequency caused by mutual inductance is
plotted in Fig. 82. We can see that the shield for the inductive coupling starts to reduce V ind at
approximately Rsh /Lsh . Therefore, the shield grounded at one end does not improve the inductive
crosstalk. To do this we need to ground the shield at both ends. However, shields grounded at both
ends are susceptible to ground loop problems.
130
|V IND|
Low frequency
High frequency
20 dB/decade
ω
ω ≅ Rsh/Lsh
Figure 82. Induced voltage caused by inductive coupling.
E = VF E ≈ jωCgr Vin R.
VNcap cap
(6.42)
−RN E
E ≈ jωLgr
VNind Ig (DC) (6.43)
RN E + RF E
and
RF E
E ≈ jωLgr
VFind Ig (DC) (6.44)
RN E + RF E
131
and
( )
VF E ind cap
= jω MF E + MF E (6.46)
Vs
where
RN E Lgr
MNind
E = , (6.47)
RN E + RF E Rs + RL
RF E Lgr
MFind
E = (6.48)
RN E + RF E Rs + RL
and
RL Cgr
E = MF E = (RN E ||RF E )
MNcap cap
. (6.49)
Rs + RL
Next, if we include the common impedance noise voltage (these expressions are valid with or without
a shield) we have:
( )
VN E cap
= jω MN E + MN E + MNciE
ind
(6.50)
Vs
and
( )
VF E cap
= jω MF E + MF E + MFciE
ind
(6.51)
Vs
where
RN E Ro
MNciE = (6.52)
RN E + RF E Rs + RL
and
RF E Ro
MFciE = − . (6.53)
RN E + RF E Rs + RL
With a shield
For capacitive coupling with a shield we have the following equations:
E = VF E ≈ jωRCgr Vg (DC).
VNcap cap
(6.54)
[ ]
−RN E
VNind
E = jωLgr Ig (DC) SF (6.55)
RN E + RF E
132
and
[ ]
RF E
VFind
E = jωLgr Ig (DC) SF (6.56)
RN E + RF E
where
Rsh 1
SF = = Lsh
. (6.57)
Rsh + jωLsh 1 + jω R sh
Then to compute both the inductive and capacitive coupling between TLs with a shield we use the
following expressions:
( )
VN E ind cap
= jω MN E + MN E (6.59)
Vs
and
( )
VF E ind cap
= jω MF E + MF E (6.60)
Vs
where
−RN E Lgr
MNind
E = SF, (6.61)
RN E + RF E Rs + RL
RF E Lgr
MFind
E = SF (6.62)
RN E + RF E Rs + RL
and
RL Cgr
E = MF E = (RN E ||RF E )
MNcap cap
. (6.63)
Rs + RL
133
Coupling
magnetic
fields
IG
LHT
+ - + -
E1 = jωlgr1LHTIg E1 = jωlgr1LHTIg
RNE RFE
E2 = jωlgr2LHTIg E2 = jωlgr2LHTIg
+ - + -
Ground reference
Figure 84. Equivalent circuit of the twisted pair coupling problem.
The ind./cap. model used to investigate the crosstalk in the twisted pair will depend on how
the ground is connected (i.e., if a balanced or unbalanced system is used). First, if we assume
an unbalanced system we get the equivalent circuits in Figs. 85 and 86. The model in Fig. 85
represents the inductive coupling and the model in Fig. 86 represents the capacitive coupling.
134
we have a zero-sum value for the KVL analysis of each twist. This indicates that the inductive
coupling between the single wire and the twisted pair in Fig. 83 is almost zero or negligible. This
then tells us that twisted pair is useful when it is desired to reduce the crosstalk due to inductive
coupling.
+ - + -
+
E1 = jωlgr1LHTIg E2 = jωlgr2LHTIg
IND
RNE V NE RFE
E2 = jωlgr2LHTIg E1 = jωlgr1LHTIg
-
+ - + -
Figure 85. Equivalent circuit for the inductive coupling in the twisted pair coupling problem.
that the twisted pair is tight, then we can assume the the following: cgr1 ≈ cgr2 . This assumption
states that the capacitance between the generator in Fig. 83 and both of the conductors in the
twisted pair are the same. This then reduces the circuit in Fig. 86 to the circuit in Fig. 87. We
VNCAP
E ≈ jωcgr1 LHT NT Vg (DC)RN E ||RF E (6.64)
or
VNCAP
E ≈ jωCgr1
T
Vg (DC)RN E ||RF E (6.65)
T
where LHT is the length of the full-twist, NT is the number of twists and Cgr1 = cgr1 LHT NT .
Therefore, using a twisted pair with unbalanced grounding reduces V IN D but not V CAP . Finally,
in general we have
RN E [ ]
VN E ≈ jω(lgr1 − lgr2 )LHT I(DC) + jωCgr1
T
Vg (DC)RN E ||RF E . (6.66)
RN E + RF E
135
+
CAP
RNE VNE RFE
I2 = jωcgr2LHTVg
I1 = jωcgr1LHTVg I1 = jωcgr1LHTVg
I2 = jωcgr2LHTVg
Figure 86. Equivalent circuit for the capacitive coupling in the twisted pair coupling problem.
CAP
RNE VNE 2I1 RFE
Figure 87. Reduced equivalent circuit for the capacitive coupling in the twisted pair coupling problem.
coupling for this case is shown in Fig. 89. The inductive coupling is the same as for the unbalanced
case. We can see that the near-end voltage is zero for an even number of twists and close to zero
for an odd number of twists. Therefore, for the balanced case the capacitive coupling is essentially
zero. In summary, if possible it is best to use the twisted pair in a balanced case because both the
inductive and capacitive coupling is minimized.
RNE/2 RFE/2
RNE/2 RFE/2
Figure 88. The twisted pair coupling problem with a balanced load.
136
{
0 for an even number of twists
+ -
+
RNE/2 RFE/2
CAP
VNE
RNE/2 RFE/2
-
N(I1+II2) N(I1+II2)
Figure 89. The equivalent circuit for the capacitive coupling in the twisted pair coupling problem with a
balanced load.
137