Electronics Guide

Laplace Transform

The Laplace transform converts a function of time into a function of a complex frequency variable, s. In electronics it turns the linear differential equations of circuits with capacitors and inductors into algebraic equations, which are solved in the s-domain and converted back into functions of time, usually with a table and a partial fraction expansion. The same framework yields the transfer function, the pole-zero description of a circuit, and a direct test for stability.

The transform extends two tools that engineers usually meet first. Phasor analysis handles sinusoidal steady state but ignores the transients that follow switching. The Fourier transform has no built-in place for initial conditions and does not exist for signals that grow. The unilateral Laplace transform accepts initial capacitor voltages and inductor currents as part of the problem, and it converges for many signals whose Fourier transforms do not exist.

The transform carries the name of Pierre-Simon Laplace, who used a similar transform in his work on probability theory; integrals of this kind appear to have first attracted his attention in 1782. Its route into engineering ran through the operational calculus of Oliver Heaviside, developed between 1880 and 1887, which replaced the differential operator with an algebraic variable and recovered solutions from tables. The method was highly successful but unproved until Thomas Bromwich put it on a rigorous basis by treating the operators as contour integrals. Gustav Doetsch published a text on the Laplace transform and its engineering applications in 1937, which the MacTutor History of Mathematics archive describes as the first such text, and Murray F. Gardner and John L. Barnes published Transients in Linear Systems Studied by the Laplace Transformation in 1942. During and soon after World War II, the Laplace transform replaced Heaviside's operational calculus in engineering practice.

Definition and Complex Frequency

The unilateral Laplace transform of a function f(t) is

F(s) = ∫0− f(t) e−st dt

where s = σ + jω is a complex variable. The operation is written L{f(t)} = F(s), and its inverse L−1{F(s)} = f(t); textbooks set the operator as a script capital L. Time functions are written in lowercase and their transforms in uppercase, so a voltage v(t) has the transform V(s). Because the unilateral transform looks only at t ≥ 0, time functions are treated as zero before the origin, as if multiplied by the unit step u(t), which is 0 for t < 0 and 1 for t > 0.

The variable s is called complex frequency. In the kernel e−st = e−σte−jωt, the imaginary part ω is an ordinary angular frequency in radians per second, and the real part σ is an exponential rate in reciprocal seconds, so the transform describes a signal in terms of exponentially weighted sinusoids rather than the constant-amplitude sinusoids of Fourier analysis. Because st must be dimensionless, s has units of reciprocal seconds, and F(s) carries the units of f(t) multiplied by seconds: the transform of a voltage is measured in volt-seconds.

The Lower Limit of Integration

The integral starts at t = 0, an instant just before zero, and the choice matters whenever something happens exactly at the origin, such as a switch closing or an impulse arriving. With the 0 limit, an impulse at t = 0 lies inside the range of integration, so the transform of the unit impulse δ(t) is 1, and the initial conditions in the transformed equations are the values just before the event, such as the capacitor voltage before a switch closes. Those are the values an engineer actually knows.

A transform defined from 0+ excludes an impulse at the origin and calls for initial values just after the event, which are often the very unknowns being sought. Many textbooks leave the lower limit unstated or mix the two conventions. Kent Lundberg, Haynes Miller, and David Trumper of the Massachusetts Institute of Technology examined the resulting confusion in "Initial Conditions, Generalized Functions, and the Laplace Transform: Troubles at the Origin" (IEEE Control Systems Magazine, February 2007) and advocated the 0 form with its matching derivative rule, which this article uses throughout.

Existence and the Region of Convergence

The defining integral does not converge for every function or every value of s. The unit step shows the mechanism: its transform, ∫0− e−st dt = 1/s, exists only when Re(s) > 0, because only then does e−st die away as t grows. The growing exponential eat gives 1/(s − a) for Re(s) > a, so the transform exists even for a signal that grows without bound, provided s lies far enough to the right.

A sufficient condition for existence is that f(t) be piecewise continuous on every finite interval and of exponential order, meaning that constants M and σ0 exist such that |f(t)| ≤ M exp(σ0t) for all sufficiently large t; the weighting e−σt then overpowers the growth of f(t) whenever σ > σ0. Steps, ramps, sinusoids, polynomials, and exponentials all qualify, and impulses are handled separately as generalized functions. A signal that grows faster than any exponential, such as exp(t2), has no Laplace transform, and neither does 1/t, which cannot be integrated near the origin.

The Region of Convergence

The set of values of s for which the integral converges is the region of convergence. For the unilateral transform it is a half-plane to the right of a vertical line, Re(s) > σc, where σc is the abscissa of convergence; for a signal of finite duration it is the entire plane. For a rational F(s), the boundary passes through the rightmost pole.

Outside the region the integral diverges, but the expression for F(s) extends beyond it by analytic continuation, and its poles, which lie outside the region or on its boundary, carry most of the information the method uses. Because the unilateral transform assumes that f(t) is zero for t < 0, F(s) alone identifies f(t): by Lerch's theorem, two piecewise continuous functions with the same transform can differ only at isolated points of discontinuity. Many unilateral tables therefore omit the region of convergence, which the bilateral transform, described near the end of this article, cannot do.

Transforms of Common Signals

A few direct integrations produce almost every transform needed in circuit analysis. The sifting property of the unit impulse gives ∫0− δ(t)e−st dt = e0 = 1. The sinusoids follow from Euler's formula: because cos ωt = (ejωt + e−jωt)/2, its transform is

½[1/(s − jω) + 1/(s + jω)] = s/(s2 + ω2)

and the same steps give ω/(s2 + ω2) for sin ωt. Repeated integration by parts gives n!/sn+1 for tn. In the table, a and ω are real constants, and T is a positive delay.

Common unilateral Laplace transform pairs
Signal f(t), t ≥ 0 Transform F(s) Region of convergence
Unit impulse δ(t)1All s
Unit step u(t)1/sRe(s) > 0
Ramp t1/s2Re(s) > 0
tn, n = 1, 2, 3, …n!/sn+1Re(s) > 0
e−at1/(s + a)Re(s) > −a
te−at1/(s + a)2Re(s) > −a
tne−atn!/(s + a)n+1Re(s) > −a
1 − e−at, a > 0a/[s(s + a)]Re(s) > 0
sin ωtω/(s2 + ω2)Re(s) > 0
cos ωts/(s2 + ω2)Re(s) > 0
e−at sin ωtω/[(s + a)2 + ω2]Re(s) > −a
e−at cos ωt(s + a)/[(s + a)2 + ω2]Re(s) > −a
t sin ωt2ωs/(s2 + ω2)2Re(s) > 0
t cos ωt(s2 − ω2)/(s2 + ω2)2Re(s) > 0
Delayed impulse δ(t − T)e−sTAll s
Delayed step u(t − T)e−sT/sRe(s) > 0
Rectangular pulse u(t) − u(t − T)(1 − e−sT)/sAll s

The pattern in the denominators runs through the rest of the subject: a factor (s + a) goes with e−at, a factor s2 + ω2 with a sinusoid of frequency ω, and a squared factor with a multiplying t. The roots of the denominators are the poles of F(s), and their locations decide what the signal does.

Properties of the Transform

Each property follows from the defining integral, and together they build nearly every transform met in circuit work from a few basic pairs. Below, f(t) and g(t) have transforms F(s) and G(s) and are zero for t < 0.

Linearity and Time Scaling

The transform of a weighted sum is the weighted sum of the transforms, L{af(t) + bg(t)} = aF(s) + bG(s), which lets a waveform be split into table entries and carries superposition into s-domain circuit analysis. Compressing or stretching the time axis rescales s in the opposite sense:

L{f(ct)} = (1/c) F(s/c), for c > 0

Filter designers apply the relationship in reverse: replacing s with s/ωc in a prototype designed for a cutoff of 1 rad/s multiplies every pole and zero by ωc, moves the cutoff to ωc, and compresses the impulse response in time by the same factor.

Time Delay and Frequency Shift

Delaying a signal by T multiplies its transform by e−sT:

L{f(t − T) u(t − T)} = e−sT F(s), for T > 0

The unit-step factor matters, because the property applies only to a signal shifted as a whole, start included: sin ω(t − T) u(t), a sinusoid with a phase offset that starts at t = 0, has a different transform from the delayed sinusoid sin ω(t − T) u(t − T). Multiplying by an exponential has the dual effect, a shift in s:

L{e−at f(t)} = F(s + a)

For a > 0, every pole and zero moves left by a, which turns the sine and cosine entries of the table into their damped forms and is why damping appears in the s-plane as a horizontal displacement.

Differentiation and Initial Conditions

The derivative property makes the transform a tool for differential equations. Integrating the defining integral by parts gives

L{df/dt} = sF(s) − f(0)

and applying the rule twice gives

L{d2f/dt2} = s2F(s) − s f(0) − f′(0)

In general, the nth derivative transforms to snF(s) minus a polynomial in s built from the initial values of f and its first n − 1 derivatives. Differentiation becomes multiplication by s, and initial conditions enter as additive terms, with no separate step to fit constants.

The 0 limit handles jumps correctly. The unit step jumps from 0 to 1 at the origin, so its derivative is the unit impulse, and the rule gives s(1/s) − u(0) = 1, the transform of δ(t). A rule written with 0+ would subtract u(0+) = 1 and return 0, losing the impulse.

Integration and Multiplication by t

Integration divides by s:

L{∫0−t f(τ) dτ} = F(s)/s

Applied to the capacitor law v(t) = v(0) + (1/C)∫0−t i(τ) dτ, the rule gives V(s) = v(0)/s + I(s)/(sC), the series capacitor model used later. Multiplying a signal by t differentiates its transform, L{t f(t)} = −dF(s)/ds; the table entries te−at, t sin ωt, and t cos ωt follow from this rule, and so does the link between repeated poles and powers of t in a response.

Convolution

The convolution of two signals that start at t = 0 is

(f * g)(t) = ∫0−t f(τ) g(t − τ) dτ

and its transform is the product F(s)G(s). The output of a linear time-invariant system at rest is the convolution of its input with its impulse response, so in the s-domain the output is a product, and cascaded stages that do not load one another multiply their transfer functions. The converse does not hold: the transform of a product f(t)g(t) is not F(s)G(s), a common error in the analysis of modulated or gated signals.

Summary of Properties

A signal that repeats with period T for t ≥ 0 has the transform F1(s)/(1 − e−sT), where F1(s) is the transform of its first period, because the delayed copies of that period sum as a geometric series. The table collects the properties.

Properties of the unilateral Laplace transform
Property Time function Transform
Linearityaf(t) + bg(t)aF(s) + bG(s)
Time scaling, c > 0f(ct)(1/c) F(s/c)
Time delay, T > 0f(t − T) u(t − T)e−sT F(s)
Frequency shifte−at f(t)F(s + a)
First derivativedf/dtsF(s) − f(0)
Second derivatived2f/dt2s2F(s) − s f(0) − f′(0)
Integration0−t f(τ) dτF(s)/s
Multiplication by tt f(t)−dF(s)/ds
Convolution(f * g)(t)F(s) G(s)
Periodic signal, period Tf(t + T) = f(t) for t ≥ 0F1(s)/(1 − e−sT)

Initial- and Final-Value Theorems

Two limit theorems extract the start and the end of f(t) directly from F(s), without inverting the transform. Both are widely used to check results, and both give wrong answers without warning when their conditions are ignored.

Statements

The initial value theorem gives the value just after the origin, and the final value theorem gives the value at which the signal settles:

f(0+) = lims→∞ sF(s)

limt→∞ f(t) = lims→0 sF(s)

Both follow from the derivative property, ∫0− f′(t)e−st dt = sF(s) − f(0). As s approaches zero, the weighting disappears and the integral tends to the total change f(∞) − f(0), provided that a final value exists. As s grows without bound along the positive real axis, the weighting suppresses everything except the instant at the origin, where only the jump f(0+) − f(0) survives. The initial value theorem therefore returns the value at 0+ even though the transform starts at 0.

Conditions of Validity

The initial value theorem requires that f(t) contain no impulse, or derivative of an impulse, at the origin. For a rational F(s), that means F(s) must be strictly proper, with a numerator of lower degree than the denominator; otherwise sF(s) grows without bound. Dividing first removes the problem: (s + 3)/(s + 1) = 1 + 2/(s + 1) is the transform of δ(t) + 2e−t, and the theorem applied to the remainder gives f(0+) = 2.

The final value theorem requires that f(t) actually settle. For a rational F(s), every pole of sF(s) must lie in the open left half-plane; equivalently, F(s) may have at most one simple pole at the origin, and all its other poles must have negative real parts. For f(t) = sin ωt, sF(s) = sω/(s2 + ω2) tends to 0 as s approaches 0, yet the sine oscillates forever, because its poles at ±jω sit on the imaginary axis. For f(t) = e2t, sF(s) = s/(s − 2) also tends to 0, yet the signal grows without bound, because its pole at s = 2 lies in the right half-plane.

Uses

Within their conditions, the theorems are quick and reliable. The final value of the unit-step response of a stable transfer function H(s) is lims→0 s × H(s)/s = H(0), the DC gain, and the same argument shows that an integrator inside a stable feedback loop drives the steady-state error for a step command to zero. Applied to df/dt, the initial value theorem gives the initial slope, f′(0+) = lims→∞ s[sF(s) − f(0)], valid when f(t) does not jump at the origin; the worked RLC example uses it to check the initial rate of rise of current.

Inverse Transform by Partial Fractions

Returning from F(s) to f(t) turns an algebraic answer into a waveform. The general inversion formula is the Bromwich integral

f(t) = [1/(2πj)] ∫c−j∞c+j∞ F(s) est ds

taken along a vertical line Re(s) = c to the right of every singularity of F(s), but circuit work rarely evaluates it. Lumped linear circuits driven by the usual test inputs produce rational transforms, ratios of polynomials N(s)/D(s) with real coefficients, and partial fraction expansion splits them into table entries. The roots of D(s) are the poles, and their pattern sets the form of the expansion. The procedures below assume that N(s) has lower degree than D(s); the last subsection covers the other case.

Distinct Poles

When D(s) has n distinct roots p1, p2, …, pn, the expansion is

F(s) = K1/(s − p1) + K2/(s − p2) + … + Kn/(s − pn)

Each coefficient, the residue at its pole, comes from the Heaviside cover-up method: multiply F(s) by (s − pi), cancel the common factor, and evaluate what remains at s = pi. Each term inverts to Ki exp(pit). For F(s) = 6(s + 2)/[s(s + 1)(s + 3)]:

  1. At s = 0, cover the factor s: K1 = 6(2)/[(1)(3)] = 4.
  2. At s = −1, cover (s + 1): K2 = 6(1)/[(−1)(2)] = −3.
  3. At s = −3, cover (s + 3): K3 = 6(−1)/[(−3)(−2)] = −1.
  4. Assemble and invert: F(s) = 4/s − 3/(s + 1) − 1/(s + 3), so f(t) = 4 − 3e−t − e−3t for t ≥ 0.

As checks, the residues sum to zero, matching f(0+) = lims→∞ sF(s) = 0, and the final value, 4, matches lims→0 sF(s) = 12/3.

Repeated Poles

A pole p of multiplicity m contributes m terms:

A1/(s − p) + A2/(s − p)2 + … + Am/(s − p)m

The cover-up method gives the highest-order coefficient: Am is (s − p)mF(s) evaluated at s = p. The others need derivatives: Am−k is 1/k! times the kth derivative of (s − p)mF(s), evaluated at s = p. Each term Ak/(s − p)k inverts to Aktk−1ept/(k − 1)!, which is why repeated poles put factors of t, t2, and so on in front of exponentials. For F(s) = (s + 3)/[(s + 1)2(s + 2)] = A1/(s + 1) + A2/(s + 1)2 + B/(s + 2):

  1. Cover (s + 2) and set s = −2: B = (1)/(−1)2 = 1.
  2. Cover (s + 1)2 and set s = −1: A2 = (2)/(1) = 2.
  3. Differentiate the covered expression, (s + 3)/(s + 2), to get −1/(s + 2)2, and set s = −1: A1 = −1.
  4. Invert: f(t) = −e−t + 2te−t + e−2t for t ≥ 0.

Evaluating both sides at a convenient point can replace the derivative step: at s = 0, F(0) = 3/2 must equal A1 + A2 + B/2 = A1 + 2.5, which again gives A1 = −1.

Complex-Conjugate Poles

Real coefficients make complex poles arrive in conjugate pairs, whose terms combine into damped sinusoids. One method completes the square and matches the damped sine and cosine entries of the table. For F(s) = (s + 7)/(s2 + 4s + 13), the denominator is (s + 2)2 + 32, with poles at s = −2 ± j3, and rewriting the numerator as (s + 2) + 5 gives

F(s) = (s + 2)/[(s + 2)2 + 32] + (5/3) × 3/[(s + 2)2 + 32]

so f(t) = e−2t[cos 3t + (5/3) sin 3t], or 1.944e−2t cos(3t − 59.04°).

The other method treats each complex pole like a distinct real pole. A pair at −α ± jβ whose upper pole has the residue K inverts to 2|K|e−αt cos(βt + θ), where θ is the angle of K; the lower pole has the conjugate residue. Here the cover-up method at p = −2 + j3 gives K = (p + 7)/(p + 2 + j3) = (5 + j3)/(j6) = 0.5 − j0.833, so 2|K| = 1.944 and θ = −59.04°, as before. The real part of the pole sets the rate of decay, and the imaginary part sets the frequency of oscillation.

Improper Functions and Delays

When the numerator degree equals or exceeds the denominator degree, divide first; the polynomial part of the quotient inverts to impulses and their derivatives, as in the example under the initial value theorem. In a circuit, an impulse usually marks an idealization, such as an ideal voltage step applied directly across a capacitor. When F(s) contains factors such as e−sT, expand the rational function that multiplies each exponential, invert it, and apply the time-delay property.

MATLAB and SciPy both provide a residue function for numerical partial fraction expansion, and SymPy inverts common transforms symbolically. Numerical expansion becomes ill-conditioned when poles are repeated or closely spaced, so hand results remain a useful check.

Solving Linear Differential Equations

The same machinery solves any linear differential equation with constant coefficients, whether it describes a circuit, a servo motor, or a thermal model:

  1. Transform both sides term by term, using the derivative property so that every initial condition appears explicitly.
  2. Solve the algebraic equation for the transform of the unknown.
  3. Expand that transform in partial fractions.
  4. Invert each term with the table, and state the result for t ≥ 0.
  5. Check the answer against the initial conditions, the final value, and the original equation.

For example, take

d2y/dt2 + 3 dy/dt + 2y = 4u(t), with y(0) = 1 and y′(0) = 0

Transforming term by term gives [s2Y(s) − s] + 3[sY(s) − 1] + 2Y(s) = 4/s, so

Y(s) = 4/[s(s + 1)(s + 2)] + (s + 3)/[(s + 1)(s + 2)]

The cover-up method expands the first term to 2/s − 4/(s + 1) + 2/(s + 2) and the second to 2/(s + 1) − 1/(s + 2). Collecting terms and inverting gives

y(t) = 2 − 2e−t + e−2t, for t ≥ 0

At t = 0 the solution equals 1 and its derivative, 2e−t − 2e−2t, equals 0, and the final value, 2, is the constant solution of 2y = 4.

Zero-Input and Zero-State Responses

The first part of Y(s), 4/[s(s + 1)(s + 2)], is the zero-state response, produced by the input with the system at rest. The second, (s + 3)/[(s + 1)(s + 2)], is the zero-input response, produced by the initial conditions alone. Linearity makes the complete response their sum. Both parts share the characteristic polynomial s2 + 3s + 2, whose roots, −1 and −2, are the natural frequencies of the system and produce the modes e−t and e−2t whatever the input. The input contributes its own pole, here at s = 0, which generates the forced response, the constant 2.

Compared with the Classical Method

The classical method adds a guessed particular solution to a homogeneous solution and then fits constants to the initial conditions. The Laplace method folds that last step into the first, handles steps, pulses, and delays without special cases, and turns coupled mesh or node equations into simultaneous algebraic equations in s. It does require constant coefficients: t d2y/dt2 + y = 0 transforms into another differential equation, because multiplication by t becomes differentiation with respect to s.

Circuit Elements in the s-Domain

Engineers usually skip the differential equation and transform the circuit itself. Kirchhoff's laws are linear, so they hold for transformed currents and voltages, and each element is replaced by an s-domain model that includes its stored energy.

Element Models

A resistor gives V(s) = RI(s). The inductor law v = L di/dt and the capacitor law i = C dv/dt become

V(s) = sL I(s) − L i(0)

I(s) = sC V(s) − C v(0)

Each law can be drawn as a series form, suited to mesh analysis, or a parallel form, suited to nodal analysis. The inductor's series form is an impedance sL in series with a voltage source L i(0), polarized to drive current in the direction of the initial current; dividing by sL gives the parallel form, sL in parallel with a current source i(0)/s directed the same way. The capacitor's parallel form is an impedance 1/(sC) in parallel with a current source C v(0) that drives current out of the positive terminal, as the charged capacitor itself would; dividing by sC gives the series form, 1/(sC) in series with a voltage source v(0)/s of the same polarity as the initial voltage. A constant in the s-domain is an impulse in time, so the sources C v(0) and L i(0) deliver the stored charge and flux linkage to the circuit at t = 0.

s-domain element models with the passive sign convention
Element Time domain Series form Parallel form
Resistorv = RiV(s) = R I(s)I(s) = V(s)/R
Inductorv = L di/dtV(s) = sL I(s) − L i(0)I(s) = V(s)/(sL) + i(0)/s
Capacitori = C dv/dtV(s) = I(s)/(sC) + v(0)/sI(s) = sC V(s) − C v(0)

Impedance and the Familiar Methods

With the initial conditions set to zero, V(s)/I(s) is the impedance Z(s): R, sL, or 1/(sC). The transformed circuit then behaves like a resistive network, and series and parallel combination, voltage and current division, source transformation, superposition, nodal and mesh analysis, and Thevenin and Norton equivalents all apply unchanged. A resistor R in series with a capacitor C, for example, gives the capacitor voltage VC(s) = Vin(s) × [1/(sC)]/[R + 1/(sC)] = Vin(s)/(1 + sRC).

Controlled sources, ideal operational amplifiers, and ideal transformers keep their algebraic relationships, which contain no derivatives. Coupled coils add mutual terms: the voltage across one coil gains sM I2(s) from the current in the other, with an initial-condition term M i2(0) of the same form as the self-inductance term.

Switched Circuits

Switching problems follow a fixed pattern. First, analyze the circuit before switching, usually in DC steady state with capacitors open and inductors shorted, to find every capacitor voltage and inductor current at t = 0. Next, redraw the circuit as it stands after switching, with s-domain models and their initial-condition sources. Finally, solve for the wanted transform and invert it. Superposition then separates the zero-input response to stored energy from the zero-state response to external inputs.

Transfer Functions, Poles, and Zeros

For a linear time-invariant system with input x(t) and output y(t), the transfer function is the ratio of their transforms with every initial condition set to zero:

H(s) = Y(s)/X(s)

The transfer function belongs to the system rather than to any input. Because the transform of a unit impulse is 1, H(s) is also the transform of the impulse response h(t), and the response to any input follows from Y(s) = H(s)X(s).

Rational Form

A lumped circuit of resistors, capacitors, inductors, and controlled sources has a rational transfer function, written as a ratio of polynomials or in factored form:

H(s) = (bmsm + … + b1s + b0)/(ansn + … + a1s + a0) = K(s − z1)…(s − zm)/[(s − p1)…(s − pn)]

The roots zi of the numerator are the zeros, the roots pi of the denominator are the poles, and K is a gain constant. The denominator degree n is the order of the system, which for a circuit cannot exceed the number of independent energy-storage elements. Real coefficients force complex poles and zeros into conjugate pairs. Transfer functions of real circuits are proper, with m ≤ n, because an improper one would demand gain that grows without limit as frequency rises. Pole-zero plots mark poles with crosses and zeros with circles.

Examples from Simple Circuits

The RC divider above is a low-pass filter, H(s) = 1/(1 + sRC), with one pole at s = −1/(RC) and no finite zeros. With R = 1 kΩ and C = 100 nF, the pole sits at −10,000 s−1, the time constant is 100 µs, and the −3 dB frequency is 10,000 rad/s, about 1.59 kHz. Taking the output across the resistor instead gives the high-pass filter H(s) = sRC/(1 + sRC), with the same pole and a zero at the origin that blocks DC.

A series RLC circuit driven by a voltage source gives second-order forms. The output across the resistor is a band-pass function, H(s) = (R/L)s/[s2 + (R/L)s + 1/(LC)], and the output across the capacitor is the standard low-pass form

H(s) = ω02/(s2 + 2ζω0s + ω02)

with ω0 = 1/√(LC) and damping ratio ζ = (R/2)√(C/L). For ζ < 1 the poles are s = −ζω0 ± jω0√(1 − ζ2). Each lies at distance ω0 from the origin, and ζ is the cosine of the angle between the negative real axis and the line from the origin to the pole.

What Poles and Zeros Do

Each simple pole p contributes a term proportional to ept to the impulse response, and the same natural modes appear in the response to any input. Zeros create no modes of their own; they change the residues and so set how strongly each mode is excited. A zero at s = z blocks transmission of an exponential input ezt, which is why zeros on the imaginary axis produce the nulls of notch filters.

A right-half-plane zero, found for example in the control-to-output response of a boost converter and in a Miller-compensated amplifier stage, raises gain like an ordinary zero but adds phase lag like a pole, and it can make a step response start in the wrong direction. Such nonminimum-phase zeros limit how fast a feedback loop around the stage can be made.

When a zero coincides with a pole, both vanish from H(s), but the mode still exists inside the system. Canceling a right-half-plane pole therefore hides an instability rather than removing it, and component tolerances make exact cancellation impossible in hardware anyway.

The s-Plane and Stability

Plotting poles on the s-plane, with σ on the horizontal axis and jω on the vertical axis, turns a transfer function into a map of behavior.

Pole locations and the modes they produce
Pole location Term in h(t) Behavior
Negative real axis, s = −ae−atDecays with time constant 1/a
Origin, s = 0ConstantNeither decays nor grows
Positive real axis, s = aeatGrows exponentially
Left half-plane pair, s = −α ± jβe−αt cos(βt + θ)Oscillates at β rad/s inside a decaying envelope
Imaginary-axis pair, s = ±jβcos(βt + θ)Oscillates with constant amplitude
Right half-plane pair, s = α ± jβeαt cos(βt + θ)Oscillates with growing amplitude
Pole of multiplicity k at ptk−1eptAs for a simple pole at p, with a polynomial factor; grows if p lies on the imaginary axis

In the table, a, α, and β are positive. Poles farther to the left decay faster, and pairs farther from the real axis oscillate faster. Poles on a circle about the origin share an undamped natural frequency, and poles on a ray from the origin share a damping ratio. A pole or pair much closer to the imaginary axis than the rest decays slowest and dominates the response once the faster modes have died away.

Bounded-Input, Bounded-Output Stability

A system is bounded-input, bounded-output (BIBO) stable when every bounded input produces a bounded output. For a linear time-invariant system this holds exactly when the impulse response is absolutely integrable, and for a proper rational transfer function it holds exactly when every pole has a negative real part. A single pole to the right of the imaginary axis makes the system unstable, however small its residue.

Simple poles on the axis, with none to the right, give marginal stability: the impulse response stays bounded but never decays. Such a system is still not BIBO stable, because an input at the pole frequency produces an unbounded output. An ideal integrator, H(s) = 1/s, turns a step into a ramp, and a lossless resonator, H(s) = ω/(s2 + ω2), driven by sin ωt produces the transform ω2/(s2 + ω2)2, whose inverse contains a growing t cos ωt term. Oscillators exploit this boundary: their poles start slightly in the right half-plane so that oscillation builds up, and an amplitude-limiting nonlinearity then holds them effectively on the axis.

Testing Without Finding Roots

Stability depends only on the roots of the characteristic polynomial, the denominator of H(s). A necessary condition for every root to have a negative real part is that all coefficients be nonzero and share the same sign, and for first- and second-order polynomials that condition is also sufficient. A third-order polynomial a3s3 + a2s2 + a1s + a0 with positive coefficients has all its roots in the left half-plane exactly when a2a1 > a3a0. Thus s3 + 2s2 + 3s + 5 is stable, since 6 > 5, while s3 + 2s2 + 3s + 10 has a pair of right-half-plane roots, since 6 < 10. The Routh-Hurwitz criterion extends the test to any order through a tabular array built from the coefficients.

Feedback Moves Poles

An amplifier with forward gain A(s) and feedback factor β(s) has closed-loop gain A(s)/[1 + A(s)β(s)], whose poles are the roots of 1 + A(s)β(s) = 0. Raising the loop gain moves those roots and, with enough phase lag, pushes a pair across the imaginary axis into oscillation. Root locus plots trace that movement in the s-plane, while Bode and Nyquist methods predict it from the loop gain evaluated at s = jω. Reducing a larger feedback structure to one such transfer function is the subject of Control System Modeling and Block Diagrams.

Frequency Response, Phasors, and the Fourier Transform

The two frequency-domain tools that most engineers learn first are both recovered by restricting s to the imaginary axis, s = jω, and each carries a condition that the Laplace framework makes explicit.

Complex Exponentials and Frequency Response

If a complex exponential est has been applied to a linear time-invariant system since the distant past, the output is the convolution of this input with the impulse response:

y(t) = ∫0− h(τ) es(t − τ) dτ = est0− h(τ) e−sτ dτ = H(s) est

provided s lies in the region of convergence of H(s). The exponential passes through unchanged in form, scaled by the complex number H(s), which is the precise sense in which s is a frequency. Setting s = jω gives the frequency response H(jω): a stable system driven by cos ωt settles to |H(jω)| cos(ωt + φ), where φ is the angle of H(jω).

The transform also shows how that steady state emerges from a sinusoid switched on at t = 0. The input cos ωt u(t) has the transform s/(s2 + ω2), with poles at ±jω. In the expansion of Y(s) = H(s)s/(s2 + ω2), the residues at those poles are H(jω)/2 and its conjugate, which together invert to |H(jω)| cos(ωt + φ). The remaining terms come from the poles of H(s); they decay if the system is stable, and otherwise no steady state exists. For H(s) = 1/(s + 1) driven by cos 2t, the complete response is y(t) = 0.447 cos(2t − 63.4°) − 0.2e−t.

Phasors as a Special Case

Phasor impedances are s-domain impedances evaluated at s = jω: jωL is sL, and 1/(jωC) is 1/(sC). Phasor analysis keeps only the terms from the input poles at ±jω and discards those from the circuit's own poles, which is legitimate only when the circuit is stable and those terms have died away. The Laplace method keeps both, which is the difference between a steady-state answer and a complete response.

Relationship to the Fourier Transform

The Fourier transform of a signal is

F(jω) = ∫−∞ f(t) e−jωt dt

For a signal that is zero before t = 0, the Laplace transform at s = σ + jω is the Fourier transform of f(t)e−σt. The factor e−σt tames signals the Fourier integral cannot handle, which is why the Laplace transform exists for growing exponentials and for the responses of unstable systems.

When the region of convergence includes the imaginary axis, as it does for any absolutely integrable signal such as the impulse response of a stable system, substituting s = jω into F(s) gives the Fourier transform. Otherwise the substitution fails. The unit step has the Laplace transform 1/s, with region Re(s) > 0, but its Fourier transform is πδ(ω) + 1/(jω), and etu(t) has no Fourier transform at all. Periodic signals, whose Fourier transforms reduce to impulses at the harmonic frequencies, are analyzed term by term in Fourier Series.

Reading Frequency Response from Poles and Zeros

In factored form, |H(jω)| equals |K| times the product of the distances from the zeros to the point jω, divided by the product of the distances from the poles, and for positive K the phase is the sum of the angles of the vectors from the zeros to jω minus the sum of the angles of the vectors from the poles. As ω sweeps past a pole close to the imaginary axis, the distance to that pole shrinks and the magnitude peaks, which is resonance. On a Bode plot, a real pole at s = −a contributes an asymptote falling at 20 dB per decade above ω = a, and at ω = a its exact contribution is −3.01 dB and −45°.

Worked Example: Series RLC Step Response

A DC source of VS = 10 V is switched at t = 0 into a series loop of R = 120 Ω, L = 10 mH, and C = 1 µF. The capacitor starts discharged and no current flows, so vC(0) = 0 and i(0) = 0. The task is to find the loop current i(t) and the capacitor voltage vC(t), as a demonstration of the method rather than of the physics of damping.

Solution

  1. Transform the circuit. The switched source becomes 10/s, the resistor remains 120 Ω, the inductor becomes the impedance sL = 0.01s, and the capacitor becomes 1/(sC) = 106/s. Both initial-condition sources are zero.
  2. Apply Kirchhoff's voltage law around the loop: 10/s = I(s)[120 + 0.01s + 106/s]. Multiplying through by s and dividing by 0.01 gives I(s) = 1000/(s2 + 12,000s + 108).
  3. Find the capacitor voltage: VC(s) = I(s)/(sC) = 109/[s(s2 + 12,000s + 108)].
  4. Locate the poles. The quadratic has roots s = −6000 ± j8000, and VC(s) has one more pole at s = 0, from the step input.
  5. Expand in partial fractions. Write VC(s) = k1/s + (k2s + k3)/(s2 + 12,000s + 108). The cover-up method gives k1 = 109/108 = 10. Matching the coefficients of s2 and s in 109 = 10(s2 + 12,000s + 108) + (k2s + k3)s gives k2 = −10 and k3 = −120,000.
  6. Complete the square to match the table. Because s2 + 12,000s + 108 = (s + 6000)2 + 80002 and 10s + 120,000 = 10(s + 6000) + 60,000, the transform becomes VC(s) = 10/s − 10(s + 6000)/[(s + 6000)2 + 80002] − 7.5 × 8000/[(s + 6000)2 + 80002].
  7. Invert term by term. For t ≥ 0, with t in seconds, vC(t) = 10 − 10e−6000t cos 8000t − 7.5e−6000t sin 8000t volts. The current follows in the same way from I(s) = 0.125 × 8000/[(s + 6000)2 + 80002], so i(t) = 125e−6000t sin 8000t milliamperes.

The cosine and sine terms combine into a single sinusoid, so the capacitor voltage can also be written vC(t) = 10 − 12.5e−6000t sin(8000t + 53.13°) V.

Checking the Result

  • Initial values: sVC(s) = 109/(s2 + 12,000s + 108) and sI(s) = 1000s/(s2 + 12,000s + 108) both tend to 0 as s grows, so vC(0+) = 0 and i(0+) = 0, consistent with the continuity of capacitor voltage and inductor current.
  • Initial slope: s × sI(s) = 1000s2/(s2 + 12,000s + 108) tends to 1000 A/s, which equals VS/L, since the full 10 V appears across the inductor at t = 0+.
  • Final value: the poles of sVC(s) lie in the open left half-plane, so the final value theorem gives vC(∞) = 109/108 = 10 V, and the current settles to zero.
  • Consistency: C dvC/dt reproduces i(t) exactly.

Reading the Poles

The real part of the pole pair sets the envelope e−6000t, with a time constant of 1/6000 s, about 167 µs, and its imaginary part sets the oscillation frequency, 8000 rad/s, about 1.27 kHz. Its distance from the origin, 10,000 rad/s, equals ω0 = 1/√(LC), and the ratio 6000/10,000 gives ζ = 0.6, the value that (R/2)√(C/L) predicts. How that damping ratio sets overshoot, ringing, and settling time is worked out in the time domain in Transient Response of RC, RL, and RLC Circuits.

A Precharged Capacitor

Suppose instead that the capacitor holds 4 V before the switch closes, with the polarity the source would give it. Its series model adds a source 4/s that opposes the 10/s source around the loop, so (10 − 4)/s = I(s)[120 + 0.01s + 106/s], I(s) = 600/[(s + 6000)2 + 80002], and i(t) = 75e−6000t sin 8000t mA. The capacitor voltage is the voltage across both parts of its model, VC(s) = I(s)/(sC) + 4/s, which inverts to

vC(t) = 10 − 6e−6000t cos 8000t − 4.5e−6000t sin 8000t V, for t ≥ 0

The result starts at 4 V and still settles at 10 V. The poles did not move: initial conditions change the residues, and with them the amplitudes of the modes, but never the natural frequencies.

Limits of the Method

The circuit methods built on the transform rest on two assumptions: the system is linear, and its parameters do not change with time.

Nonlinear and Time-Varying Circuits

Diodes, transistors driven over large swings, and saturating magnetic cores have no s-domain models, because the transform of a nonlinear expression such as v2 cannot be written algebraically in terms of V(s). Small-signal analysis linearizes the circuit about its operating point, and piecewise-linear analysis solves each linear interval in turn, carrying the initial conditions across every switching event as a new t = 0. Switched-mode power converters are routinely analyzed both ways, interval by interval for their waveforms and through averaged small-signal models for control loop design. SPICE-class simulators handle the general case by integrating the circuit equations numerically in time.

A circuit whose element values or connections change with time is time-varying even when it is linear. A mixer, which multiplies the signal by a local-oscillator waveform, is an example, and no single transfer function describes it.

Distributed and Irrational Systems

Transmission lines and other distributed elements produce transfer functions that are not rational. A lossless line terminated in its characteristic impedance delivers its input to the far end after a transit time T, so its transfer function is e−sT, and skin-effect loss introduces terms in √s. Partial fractions do not apply directly. Engineers either substitute rational approximations, such as the first-order Padé approximation e−sT ≈ (1 − sT/2)/(1 + sT/2), or invert the transform numerically with algorithms such as Talbot's method or the Gaver-Stehfest algorithm.

The Bilateral Transform

The bilateral, or two-sided, Laplace transform integrates over all time:

FB(s) = ∫−∞ f(t) e−st dt

It accepts signals that are nonzero before t = 0, and its region of convergence is in general a vertical strip, σ1 < Re(s) < σ2, so the expression alone no longer identifies the signal. With the region Re(s) > −a, 1/(s + a) is the transform of e−atu(t); with the region Re(s) < −a, the same expression is the transform of −e−atu(−t), a signal that exists only before the origin. A bilateral transform must therefore always be quoted with its region of convergence.

The bilateral form has no place for initial conditions, so circuit analysis, which starts from a known state at a switching instant, uses the unilateral transform; the two agree for any signal that is zero for t < 0. Sampled-data systems have their own counterpart, the z-transform, linked to the Laplace variable through z = esT for a sampling period T.

Summary

The unilateral Laplace transform, defined from t = 0, turns linear constant-coefficient differential equations into algebra and brings initial conditions in automatically. A short table and a few properties cover nearly every circuit problem, partial fractions invert the rational results, and the limit theorems check answers when their conditions hold. Impedances sL and 1/(sC), with initial-condition sources, let the methods of resistive circuit analysis solve transient problems. The poles of the transfer function fix the natural modes and must lie in the open left half-plane for stability. On the imaginary axis the transform reduces to frequency response, phasors, and, when the region of convergence allows, the Fourier transform. Its reach ends where linearity, time invariance, or lumped-element behavior ends.

Related Topics