Saturday, 29 November 2014

Telegrapher's equations

The telegrapher's equations (or just telegraph equations) are a pair of linear differential equations which describe the voltage and current on an electrical transmission line with distance and time. They were developed by Oliver Heaviside who created the transmission line model, and are based on Maxwell's Equations.
Schematic representation of the elementary component of a transmission line.
The transmission line model represents the transmission line as an infinite series of two-port elementary components, each representing an infinitesimally short segment of the transmission line:
  • The distributed resistance R of the conductors is represented by a series resistor (expressed in ohms per unit length).
  • The distributed inductance L (due to the magnetic field around the wires, self-inductance, etc.) is represented by a series inductor (henries per unit length).
  • The capacitance C between the two conductors is represented by a shunt capacitor C (farads per unit length).
  • The conductance G of the dielectric material separating the two conductors is represented by a shunt resistor between the signal wire and the return wire (siemensper unit length).
The model consists of an infinite series of the elements shown in the figure, and that the values of the components are specified per unit length so the picture of the component can be misleading. RLC, and G may also be functions of frequency. An alternative notation is to use R'L'C' and G' to emphasize that the values are derivatives with respect to length. These quantities can also be known as the primary line constants to distinguish from the secondary line constants derived from them, these being the propagation constant,attenuation constant and phase constant.
The line voltage V(x) and the current I(x) can be expressed in the frequency domain as
\frac{\partial V(x)}{\partial x} = -(R + j \omega L)I(x)
\frac{\partial I(x)}{\partial x} = -(G + j \omega C)V(x).
When the elements R and G are negligibly small the transmission line is considered as a lossless structure. In this hypothetical case, the model depends only on the L and C elements which greatly simplifies the analysis. For a lossless transmission line, the second order steady-state Telegrapher's equations are:
\frac{\partial^2V(x)}{\partial x^2}+ \omega^2 LC\cdot V(x)=0
\frac{\partial^2I(x)}{\partial x^2} + \omega^2 LC\cdot I(x)=0.
These are wave equations which have plane waves with equal propagation speed in the forward and reverse directions as solutions. The physical significance of this is that electromagnetic waves propagate down transmission lines and in general, there is a reflected component that interferes with the original signal. These equations are fundamental to transmission line theory.
If R and G are not neglected, the Telegrapher's equations become:
\frac{\partial^2V(x)}{\partial x^2} = \gamma^2 V(x)
\frac{\partial^2I(x)}{\partial x^2} = \gamma^2 I(x)
where γ is the propagation constant
\gamma = \sqrt{(R + j \omega L)(G + j \omega C)}
and the characteristic impedance can be expressed as
Z_0 = \sqrt{\frac{R + j \omega L}{G + j \omega C}}.
The solutions for V(x) and I(x) are:
V(x) = V^+ e^{-\gamma x} + V^- e^{\gamma x} \,
I(x) = \frac{1}{Z_0}(V^+ e^{-\gamma x} - V^- e^{\gamma x}). \,
The constants V^\pm and I^\pm must be determined from boundary conditions. For a voltage pulse V_{\mathrm{in}}(t) \,, starting at x=0 and moving in the positive x-direction, then the transmitted pulse V_{\mathrm{out}}(x,t) \, at position xcan be obtained by computing the Fourier Transform, \tilde{V}(\omega), of V_{\mathrm{in}}(t) \,, attenuating each frequency component by e^{\mathrm{-Re}(\gamma) x} \,, advancing its phase by \mathrm{-Im}(\gamma)x \,, and taking the inverse Fourier Transform. The real and imaginary parts of \gamma can be computed as
\mathrm{Re}(\gamma) = (a^2 + b^2)^{1/4} \cos(\mathrm{atan2}(b,a)/2) \,
\mathrm{Im}(\gamma) = (a^2 + b^2)^{1/4} \sin(\mathrm{atan2}(b,a)/2) \,
where atan2 is the two-parameter arctangent, and
a \equiv \omega^2 LC \left[ \left( \frac{R}{\omega L} \right) \left( \frac{G}{\omega C} \right) - 1 \right]
b \equiv \omega^2 LC \left( \frac{R}{\omega L} + \frac{G}{\omega C} \right).
For small losses and high frequencies, to first order in R / \omega L and G / \omega C one obtains
\mathrm{Re}(\gamma) \approx \frac{\sqrt{LC}}{2} \left( \frac{R}{L} + \frac{G}{C} \right) \,
\mathrm{Im}(\gamma) \approx \omega \sqrt{LC}. \,
Noting that an advance in phase by - \omega \delta is equivalent to a time delay by \deltaV_{out}(t) can be simply computed as
V_{\mathrm{out}}(x,t) \approx V_{\mathrm{in}}(t - \sqrt{LC}x) e^{- \frac{\sqrt{LC}}{2} \left( \frac{R}{L} + \frac{G}{C} \right) x }. \,