Artificial Transmission Line Dec 16
Artificial Transmission Line Dec 16
2
Dr. Milica Marković Applied Electromagnetics Laboratory page 3
Figure 1.2: Voltage as a function of time at the generator side z=0 (top) and the load side z=l
(bottom) of the transmission line in Figure 1.4, if the switch closes at t=0 the voltage arrives at
t=l/c=T at the load. These graphs can be obtained by observing the voltage on an oscilloscope at
the load and at the generator side.
Figure 1.3: Voltage along the transmission line in Figure 1.4, for four different time intervals t=0,
switch closes, t=T/4, t=T/2 and t=T. It is assumed that the length of the transmission line is equal
to l=T/c.
Figure 1.4: Electronic Circuit with an emphasis on cables that connect the generator and the load.
l
vBB ‘ (t) = vAA‘ (t − ) (1.2)
c
l
vBB ‘ (t) = Acos(ω(t − )) (1.3)
c
l
vBB ‘ (t) = Acos(ωt − ω ) (1.4)
c
ω
vBB ‘ (t) = Acos(ωt − l) (1.5)
c
Since we know that angular frequency is ω = 2πf
2πf
vBB ‘ (t) = Acos(ωt − l) (1.6)
c
c
The quantity f
is the wavelength λ. We will talk more about wavelength in Section 1.1.5
2π
vBB ‘ (t) = A cos(ωt − l) (1.7)
λ
The quantity 2π
λ
is the propagation constant β
Finally the expression for the voltage at BB end is
l
vBB ‘ (t) = A cos(ωt − 2π ) (1.10)
λ
The signal will experience a phase shift of 2π λl . If this phase shift is small, there will not be much
difference between the phase of the signal at the generator and at the load. This means that we
don’t have to use transmission line theory to account for the effects of the line. If the phase shift is
significant, then we do have to use the transmission line theory. Let’s look at some numbers in the
following example.
1. If λl < 0.01 then the angle 2π λl is of the order of 0.0314 rad or about 20 . In this case, the phase
is obviosly something that we don’t have to worry about. When the length of the transmission
λ
line is much smaller than λ, l << 100 the wave propagation on the line can be ingnored.
2. If λl > 0.01, say λl = 0.1, then the phase is 200 , which is a significant phase shift. In this case it
may be necessary to account for reflected signals, power loss and dispersion on the transmission
line.
Problem 4. Plot the previous three plots on the same graph. Which signal is lagging the sinus
Problem 5. Plot two periods of two sinusoidal signals y1(t) and y2(t) with different phases
Problem 6. Now assume that the sinusoidal signal is generated by a source positioned at x=0
Problem 7. If the voltage on the generator side is y=sin(t) then the voltage on the load sid
Solutions
Problem 1
Using Matlab plot the sinusoidal signal
y=sint as a function of time. Assume that the frequency of the signal is f=1GHz. Then repea
Code:
f=1e6;
t=0:1/100/f:4/f;
y=sin(2*pi*f*t);
plot(t,y,’Color’,’red’);
hold on
Graph of y(t)=sin(t)
In this part of the lab, we were asked to display a simple sine wave with a frequency of 1GHz
Problem 2
Using Matlab plot the sinusoidal signal
Assume that the frequency of the signal is f=1GHz. Then repeat the exercise by plotting the s
Code:
fc=1e6;
t=0:1/100/fc:4/fc;
y=sin((2*pi*fc*t)+((2*pi)/9));
plot(t,y,’Color’,’black’);
Graph of y(t)=sin(t+40)
With the phase change of +40 degrees, we see our graph shift to the left on the x-axis. This
Problem 3
Using Matlab plot the sinusoidal signal y=sin(t-40)
Assume that the frequency of the signal is f=1GHz. Then repeat the exercise by plotting the s
Code:
fc=1e6;
t=0:1/100/fc:4/fc;
y=sin((2*pi*fc*t)-((2*pi)/9));
plot(t,y);
Problem 4
Plot the previous three plots on the same graph. Which signal is lagging the sinusoidal signa
Code:
f=1e6;
t=0:1/100/f:4/f;
y=sin(2*pi*f*t);
plot(t,y,’Color’,’red’);
hold on
fc=1e6;
t=0:1/100/fc:4/fc;
y=sin((2*pi*fc*t)+((2*pi)/9));
plot(t,y,’Color’,’black’);
fc=1e6;
t=0:1/100/fc:4/fc;
y=sin((2*pi*fc*t)-((2*pi)/9));
plot(t,y);
This Code allows us to plot all 3 graphs simultaneously on the same page. (different colors a
Plot of :
y(t) =sin(t) // RED
y(t) =sin(t-40) // BLUE
y(t) =sin(t+40) // BLACK
This graph intricately shows all three sine waves together and their lagging/leading elements
Problem 5
Plot two periods of two sinusoidal signals y1(t) and y2(t) with different phases and at the
Code:
fc=1e6;
t=0:1/100/fc:4/fc;
y=sin(2*pi*fc*t);
plot(t,y);
hold on
fc=1e6;
t=0:1/100/fc:4/fc;
y=sin((2*pi*fc*t)+.025);
plot(t,y);
Plots both sine functions simultaneously with the minimal amount of phase change that display
Problem 6
Now assume that the sinusoidal signal is generated by a source positioned at x=0 and the sou
Code:
c=3e8
t=y./c
y=1:.001:100
plot(t,y);
Problem 7
If the voltage on the generator side is y=sin(t) then the voltage on the load side will lag
Code:
c=3e8;
fc=1e11; %change this number for different frequencies
t=0:1/100/fc:1/fc;
y=sin(2*pi*fc*(t-(1/c)));
plot(t,y);
hold on
F=1000Hz
Dispersion is an effect where different frequencies travel with different speeds on the transmission
line.
Example Find what is the length of the cable at which we need to take into account transmission
line effects if the frequency of operation is 10 GHz.
• Write KCL, KVL for the piece in the time domain (we get differential equations)
• Solve the linear system of equations to get the expression for the voltage and current on the
transmission line as a function of z.
KVL
0 0 ∂i(z, t)
−v(z, t) + R ∆z i(z, t) + L ∆z + v(z + ∆z, t) = 0
∂t
KCL
KVL
∂i(z, t)
−(v(z + ∆z, t) − v(z, t)) = R∆zi(z, t) + L∆z (1.11)
∂t
v(z + ∆z, t) − v(z, t) ∂i(z, t)
− = Ri(z, t) + L (1.12)
∆z ∂t
∂v(z, t) ∂i(z, t)
− = Ri(z, t) + L (1.13)
∂z ∂t
KCL
∂v(z + ∆z, t)
−(i(z + ∆z, t) − i(z, t)) = G∆zv(z + ∆z, t) + C∆z (1.14)
∂t
i(z + ∆z, t) − i(z, t) ∂v(z + ∆z, t)
− = Gv(z + ∆z, t) + C (1.15)
∆z ∂t
∂i(z, t) ∂v(z + ∆z, t)
− = Gv(z + ∆z, t) + C (1.16)
∂z ∂t
We just derived Telegrapher’s equations in time-domain:
∂v(z, t) ∂i(z, t)
− = Ri(z, t) + L
∂z ∂t
∂i(z, t) ∂v(z + ∆z, t)
− = Gv(z + ∆z, t) + C
∂z ∂t
These are two differential equations with two unknowns. It is not impossible to solve them,
however we would prefer to have linear differential equations. So what do we do now?
Express time-domain variables as phasors!
∂V (z)
− = (R + jωL)I(z) (1.17)
∂z
∂I(z)
− = (G + jωC)V (z) (1.18)
∂z
Two equations, two unknowns. To solve these equations, we first integrate both equations over
z,
∂ 2 V (z) ∂I(z)
− 2
= (R + jωL)
∂z ∂z
∂ 2 I(z) ∂V (z)
− = (G + jωC)
∂z 2 ∂z
1 ∂I(z) ∂ 2 V (z)
− = (1.19)
(R + jωL) ∂z ∂z 2
1 ∂V (z) ∂ 2 I(z)
− = (1.20)
(G + jωC) ∂z ∂z 2
now substitute Eq.1.19 into Eq.1.17 and Eq.1.20 into Eq.1.18 and we get
∂ 2 V (z)
− = (G + jωC)(R + jωL)V (z)
∂z 2
∂ 2 I(z)
− = (G + jωC)(R + jωL)I(z)
∂z 2
Or if we rearrange
∂ 2 V (z)
− − (G + jωC)(R + jωL)V (z) = 0 (1.21)
∂z 2
∂ 2 I(z)
− − (G + jωC)(R + jωL)I(z) = 0 (1.22)
∂z 2
The above Eq.1.21 and Eq.1.22 are the equations of the current and voltage wave on a transmission
line. γ = (G + jωC)(R + jωL) is the complex propagation constant. This constant has a real and
an imaginary part.
γ = α + jβ
p
α = Re{ (G + jωC)(R + jωL)}
p
β = Im{ (G + jωC)(R + jωL)}
What is the general solution of the differential equation of the type 1.21 or 1.22?
In this equation V0+ and V0− are the phasors1 of forward and backward going voltage waves, and
I0+ and I0− are the phasors of forward and backward going current waves
1
complex numbers having an amplitude and the phase
The time domain expression for the current and voltage on the transmission line we get
clear all
clc
f = 10^9;
w = 2*pi*f
c=3*10^8;
beta=2*pi*f/c;
lambda=c/f;
t1=0.1*10^(-9)
t2=0.2*10^(-9)
x=0:lambda/20:2*lambda;
Figure 1.8: Forward (top) and reflected (bottom) waves on a transmission line.
y1=sin(w*t1 - beta.*x);
y2=sin(w*t2 - beta.*x);
y3=sin(w*t1 + beta.*x);
y4=sin(w*t2 + beta.*x);
subplot (2,1,1),
plot(x,y1,’r’),...
hold on
plot(x,y2,’--b’),...
hold off
subplot (2,1,2),
plot(x,y3,’r’)
hold on
plot(x,y4,’--b’)
hold off
1.1.8 Relating forward and backward current and voltage waves on the
transmission line -Transmissio Line Impedance Z0
In this equation V0+ and V0− are the phasors of forward and backward going voltage waves, and
I0+ and I0− are the phasors of forward and backward going current waves
We will relate the phasors of forward and backward going voltage and current waves.
When substitute the voltage wave equation into Telegrapher’s Eq. 1.17. The equation is repeated
here Eq.1.30.
∂V (z)
− = (R + jωL)I(z) (1.30)
∂z
γV0+ e−γz − γV0− eγz = (R + jωL)I(z) (1.31)
γ
I(z) = (V0+ e−γz + V0− eγz )
R + jωL
γV0+ −γz γV0− γz
I(z) = e − e (1.32)
R + jωL R + jωL
Now we compare Eq.1.32 with the Eq.1.29.
γV0+
I0+ =
R + jωL
γV0−
I0− = −
R + jωL
We can define the characteristic impedance of a transmission line as the ratio of the voltage to
current amplitude of the forward going wave.
V0+
Z0 =
I0+
R + jωL
Z0 =
γ
s
R + jωL
Z0 =
G + jωC
V (0)
ZL =
I0
Substitute the boundary condition in Eq.1.46
V0+ + V0−
ZL = Z0 (1.35)
V0+ − V0−
We can now solve the above equation for V0−
ZL +
(V0 − V0− ) = V0+ + V0−
Z0
ZL ZL
( − 1)V0+ = ( + 1)V0−
Z0 Z0
ZL
V0− Z0
−1
=
V0+ ZL
Z0
+1
V0− ZL − Z0
+ = (1.36)
V0 ZL + Z0
V−
The quantity V0+ is called voltage reflection coefficient Γ. It relates the reflected to incident
0
voltage phasor. Voltage reflection coefficient isi n general a complex number, it has a magnitude and
a phase.
Examples
p
γ= (R + jωL)(G + jωC)
p
γ = jωL jωC
√
γ = jω LC = jβ
s
R + jωL
Z0 =
G + jωC
s
jωL
Z0 =
jωC
r
L
Z0 =
C
• Phase velocity
ω
v=
β
ω
v= √
ω LC
1
v=√
LC
2
metal resistance is low
3
dielectric conductance is low
• Group velocity
• Wavelength
2π
λ=
β
2π
λ= √
ω LC
2π
λ= √
0 µ0 r
c
λ= √
f r
λ0
λ= √
r
p
γ=
(R + jωL)(G + jωC) (1.37)
s
√
R G RG
γ = jω L C 1 − j + − 2 (1.38)
ωL ωC ω LC
s
√
R G
γ ≈ jω L C 1 − j + (1.39)
ωL ωC
√ q
R G
Taylor’s series for function 1 + x = 1 − j ωL
+ ωC
in Equation 1.39 is shown in Equa-
tions 1.40-1.41.
√
x x 2 x3
1+x=1+ − + − ... f or |x| < 1 (1.40)
s 2 8 16
√ √
R G j R G
γ ≈ jω L C 1 − j + = jω L C 1 − + (1.41)
ωL ωC 2 ωL ωC
4
metal resistance is lower than the inductive impedance
5
dielectric conductance is lower than the capacitive impedance
√
ω LC R G
α= + (1.42)
2 ωL ωC
√
β = ω LC (1.43)
We see that the phase constant β is the same as in the lossless case, and the attenuation
constant α is frequency independent. This means that all frequencies of a modulated signal
are attenuated the same amount, and there is no dispersion on the line. When the phase
constant is a linear function of frequency, β = const ω, then the phase velocity is a constant
1
vp = ωβ = const , and the group velocity is also a constant, and equal to the phase velocity.
In this case, all frequencies of the modulated signal are propagated with the same speed, and
there is no distortion of the signal. This is the case only when the losses in the transmission
line are small.
We usually represent phase and group velocity on ω − β diagrams, shown in Figure 1.1.11(a).
At a frequency ω1 , the ratio of ωβ gives the phase velocity (graphically, this is the slope of the
red line on the graph), whereas the slope of the ω − β curve (blue line on the graph) gives the
group velocity at this freqency. These diagrams are usefull, as they show how phase constant
β varies with frequency, and it also shows how phase and group velocities vary with frequency.
We can see that in this case, both group and phase velocity for this line are positive quantities,
which is a representation of what is called a forward wave. In a forward wave, both signal and
the energy propagate in the forward direction. Backward waves are waves where the signal
propagates forward, however energy propagates backwards. In backward waves, phase and
group velocities have opposite signs.
When the phase constant is a linear function of frequency, then, phase velocity and group
velocity are equal, and do not depend on frequency. Both velocities are equal to the slope of
the line in Figure 1.1.11(b).
• Transmission line impedance Transmission line impedance is the impedance that the forward
voltage and current see as they propagate down the line from the generator to the load.
s
R + jωL
Z0 =
G + jωC
s
jωL
Z0 =
jωC
r
L
Z0 =
C
• Group velocity Group velocity is the velocity of the modulated signal impressed on the carrier
frequency.
1
vg = ∂β
(1.44)
∂ω
1
vg = √ (1.45)
LC
• Wavelength Wavelenght is how far the signal travels on the line while the carrier signal goes
through one whole period in time, see Section??.
2π
λ=
β
2π
λ= √
ω LC
2π
λ= √
0 µ0 r
c
λ= √
f r
λ0
λ= √
r
V0− ZL − Z0
Γ= + = (1.48)
V0 ZL + Z0
In Eq.1.51, Γ is the voltage reflection coefficient, V0+ is the phasor of the forward going wave, z
is the axis in the direction of wave propagation, β is the phase constant7 , Z0 is the impedance of the
transmission line8 . V (z) is a complex number, phasor. We will find the magnitude and phase of the
voltage on the transmission line. √
The magnitude of a complex number can be found as |z| = zz ∗ 9 .
p
|V (z)| = V (z)V (z)∗
q
|V (z)| = V0+ (e−jβz − |Γ|ejβz+Θr )V0+ (ejβz − |Γ|e−(jβz+Θr ) )
q
+
|V (z)| = V0 (e−jβz − |Γ|ejβz+|Θr )(ejβz − |Γ|e−(jβz+Θr ) )
q
+
|V (z)| = V0 1 + |Γ|e−(2jβz+|Θr ) + |Γ|ej2βz+Θr + Γ2 )
q
+
|V (z)| = V0 1 + |Γ|2 + |Γ|(e−(2jβz+Θr ) + e(j2βz+Θr ) )
p
|V (z)| = V0+ 1 + |Γ|2 + 2|Γ|cos(2βz + Θr ) (1.54)
The magnitude of the total voltage on the transmission line is given by Eq.1.54. It seems like a
complicated function.
• Let’s start from a simple case when the voltage reflection coefficient on the tranmission line is
Γ = 0 and draw the magnitude of the total voltage.
HERE PICTURE OF THE FLAT LINE
• Let’s look at another case, Γ = 0.5 and Θr = 0. The equation for the magnitude
r
5
|V (z)| = V0+ + cos2βz (1.55)
4
The function 1.55 is at it’s maximum when cos(2βz) = 1 or z = k2 λ, and the function value is
V (z) = 1.5V0+ . It is at it’s minimum when cos(2βz) = −1 or z = 2k+1
4
λ and the function value
is V (z) = 0.5V0+
It is important to mention here that the function that we see looks like a cosine with an average
value of V0+ , but it is not. The minimums of the function are sharper then the maximums, so
when the reflection coefficient is at it’s maximum of Γ = 1 the function looks like this:
PICTURE OF STANDING WAVE PATTERN WITH SHORT
7
imaginary part of the complex propagation constant
8
defined as the ratio of forward going voltage and current
9
Prove this in Carthesian and Polar coordinate system
• General Case.
In general the voltage maximums will occur when cos(2βz) = 1
p
|V (z)max | = V0+ 1 + kGamma|2 + 2|Γ|
p
|V (z)max | = V0+ (1 + kGamma|)2
|V (z)max | = V0+ (1 + |Γ|) (1.56)
p
|V (z)min | = V0+ 1 + kGamma|2 − 2|Γ|
p
|V (z)min | = V0+ (1 − kGamma|)2
|V (z)min | = V0+ (1 − |Γ|) (1.57)
The ratio of voltage minimum on the line over the voltage maximum is called the Voltage
Standing Wave Ratio (VSWR) or just Standing Wave Ratio (SWR).
V (z)max
SW R =
V (z)min
1 + |Γ|
SW R = (1.58)
1 − |Γ|
cos(2βz) = 1
2βz + Θr = 2nπ
2nπ − Θr
z=
2β
2nπ − Θr
z= (1.59)
4π
Connect to the experimental setup as shown in Figure 2.2(a). The artificial transmission line box
is shown in Figure 2.2(b).
1
BW ≈ N (2.1)
10 τ
This means √
that a 50 section artificial transmissio line with a delay of
Where τ = Ctotal Ltotal is the time delay of the transmission line.
2.2.1 Question 1
How many wavelengths long is the transmission line at f = 8M Hz frequency if L = 3.55H/mile,
C = .1µF/mile, and the length of the whole transmission line is miles?
Each capacitor in the
27
Dr. Milica Marković Applied Electromagnetics Laboratory page 28
2.2.2 Question 2
What is the characteristic impedance of the transmission line?
To calculate the characteristic impedance of the line, use the equations for the lossy transmission
line impedance given in Equation ??. Comment on the value of the imaginary part of the transmission
line impedance. Is this a low-loss line? Check out the section about low-loss lines, and see if you can
2.3.3 Question 4
Are there any losses in the transmission line? Was your line well matched?
2.3.4 Question 5
Calculate the attenuation coefficients α from your measurements, and from the theory. How do they
compare? How will you find the current from the measured voltage along the line?
1. Terminate the line with an impedance equal to 2Z0 , where Z0 in the characteristic impedance
of the lossless line.
4. Measure the voltage along the line in the same manner you did before, and make a new graph
with Matlab.
6. Now select the frequency so that the line is a quarter of a wavelength long and plot the line
voltage on the same graph as in the previous part. As the last part, terminate the lossless line
in the short and open circuit and plot the line voltage for both cases.
7. What would happen to the output voltage if λ/4 transmission line was a perfect one? How do
the currents and voltages compare the case of half and quarter long transmission lines? What
standard circuit element is the quarter wavelength long line remind you off? Explain.
30
Dr. Milica Marković Applied Electromagnetics Laboratory page 31
be calculated using CRF assuming ideally lossless airline in assuming that the operating frequencies
known..
Draw the SWR circle on the Smit chart.. This is the circle whose radius is equal to the appropriate
SWR ratio on the chart. The point of the minimum for short-circuit on the chart is at the intersection
of SWR circle in the negative GR axis.
6. To measure the antenna input impedance versus frequency,, first obtain proxy load positions
beta over the desired frequency range by putting their reference short at the desired antenna load
plane (which is at the end of the cable used to connect the slotted line and the antenna).. Take data
on the proxy load plane locations versus the frequency.. Then with a monopole antenna replacing the
short, take SWR in the minimum location data in each of the same Dial set frequencies.. obtained
the input impedance for each of those frequencies and label them on the Smith charts. Plot also real
and imaginary parts of the antennas impedance versus frequency on using Matlab.
to measure. All connect tutors in this system are GR (general radio) connectors. The slotted line
probe has the diode detector built-in one of the two lands. Connect the send to the HP 415 ESW
art meter. The diode detects the amp to the of the 1 kHz low-frequency signal that modulates the
art signal picked up by the probe. The SWR meter is just a narrow band high gain amplifier for the
1 kHz that the 500 MHz is modulated with. It it’s high gain assures good measurement sensitivity.
The other GR connectors on the probe is used for optimizing the match of the probe to the line.
Connect the adjustable stop tuna to this end. This unit is just a section of quarks with a short bed
the other and so that this impedance can change by physically moving this short.
Number one in this first part of the experiment you will use a short circuit termination. Set the
frequency at 500 MHz.find and also the standing wave pattern.
Question one: why are the nulls of shorted circuit standing wave pattern sharp?
the sharp nulls
will be buried in noise. The best way to find their position is to move towards the zero from both
sides, locate points a and B from figure 42, and then take them zero to be between them. measure
and record the position of several successive zeros along an arbitrary position scale.
question two: what is the wavelength in the line equal to?
Number two: what is the wavelength in the line equal to?
Increase the generator power and repeat the measurement.
Question four: why is this measurement not as good as the previous one? What gets worse: the
minimum or the maximum of the voltage?
Number three in this part of the experiment, you will measure two loads at both 500 and 250
MHz. The first load is the unknown load from part two, and the second one is along to sister of the
same value, mounted in the metal jig from experiment to. First determine the position of the proxy
plane by connecting a short termination. Then connect the load, and measure the distance between
the first the next of the proxy plane in the proxy plane. You can either move towards or away from
the load, but make sure you are consistent on the Smith chart. Repeat this measurement for the
two loads at F1 equals 250 MHz and after equals 500 MHz. Feel a table as the one shown in table
41. The distance fell in the fifth column is the electrical vistas between columns 3 and four. Use this
this this and this may chart to find the normalized load impedance Z from the measures SWR. Then
find the load impedance Z by multiplying by the characteristic impedance of the Coax. include your
Smith chart.
Four at 500 MHz measured the impedance of a car X 50 ohm termination, repeating the procedure
from part three. Include the Smith chart.
Question 5Y this the harder measurement than for badly matched loads?
In the past part of this experiment you will measure the impedance of the monopole antenna
above a metal plane at 500 MHz. The antennae 15 cm long, which is about a quarter wavelength of
this frequency. This is called the quarter wave monopole. They impedance of the monopole with or
without the ground plane is very different.
Question six: this antenna can be viewed as a current flowing through a Short St., Finn Y which
is located about an incident perfectly conducting plane. What basic fear massage among medics
would you use to find electric and magnetic field of such a system? Illustrate your answer with a
simple drawing you do not have time to find the fields.
Question seven: the impedance of an antenna is defined as the voltage across its terminal by
current.
Conclusion
Solve all of the above examples by hand and compare the results with the ones you found
during the lab. In addition, for each of the problems, write what you have learned. Do not write an
essay, just in a few sentences write what you have done in each problem, and specifically what you
have learned from the work done.
Due Date
Submit the printout of the lab work and the conclusion by NEXT FRIDAY. LAB IS DUE BY
noon (stamped) next Friday in the main office. No late labs. Do not bring previous lab writeup for
the next Lab session.