Complex Exponential Representation of Oscillations

Solve enough oscillation problems with sines and cosines and you start to feel the friction: every derivative shuffles cos⁡\cos into −sin⁡-\sin and back, every superposition demands a sum-to-product identity you half-remember, every phase shift spawns a page of trigonometric bookkeeping. Now watch a physicist one year further along handle the same problem. She writes x=Re⁡ A~eiωtx = \operatorname{Re}\,\tilde{A}e^{i\omega t}, differentiates by multiplying by iωi\omega, divides two complex numbers, and reads off the amplitude and the phase lag from a single result — no identities, no shuffling, three lines. The machine that makes this possible is Euler's formula, eiθ=cos⁡θ+isin⁡θe^{i\theta} = \cos\theta + i\sin\theta, which Richard Feynman called "the most remarkable formula in mathematics." The trick is not a notational fad. The exponential eiωte^{i\omega t} is an eigenfunction of the derivative operator — differentiation merely rescales it — so any linear differential equation with constant coefficients collapses into a polynomial equation the moment you feed it exponentials. Trigonometric functions are what you get when you insist on splitting that one clean object into two real shadows. Mastering this representation buys you the entire frequency-domain worldview: phasors, impedance, transfer functions, Fourier analysis, and eventually the e−iEt/ℏe^{-iEt/\hbar} phases of quantum mechanics all speak this language. Damped and driven oscillators — the next articles in this track — become almost embarrassingly easy once you own it. This page builds the machine, proves it is legitimate, and shows the one place it bites.

Key takeaways
  • Linear ODEs with constant coefficients are solved most cleanly by replacing real sinusoids with complex exponentials eiωte^{i\omega t} and taking the real part at the end.
  • The derivative of eiωte^{i\omega t} is just iω eiωti\omega\,e^{i\omega t}: differentiation becomes multiplication, so differential equations become algebra.
  • A real oscillation Acos⁡(ωt+ϕ)A\cos(\omega t + \phi) is encoded by a single complex number — the phasor A~=Aeiϕ\tilde{A} = Ae^{i\phi} — which bundles amplitude and phase together.
  • Because the equations are linear with real coefficients, the real and imaginary parts never mix; you may carry the full complex solution and project at the end.
  • Superposing sinusoids of one frequency becomes vector addition of arrows in the complex plane — the law of cosines replaces every trig identity.
  • The single trap: real parts do not survive multiplication. For energies, powers, and any nonlinear quantity, project first, multiply second.
  • Key idea: differentiation of eiωte^{i\omega t} is multiplication by iωi\omega, so solve the complexified linear ODE algebraically and take the real part — amplitude and phase emerge automatically as the modulus and argument of one complex amplitude.

Where this lives today

Open any wireless chip shipped in 2026 — the 5G modem in a phone, a Wi-Fi 7 radio, a Bluetooth earbud — and the heart of it is a pair of signal paths labeled I and Q: in-phase and quadrature, the real and imaginary parts of a complex baseband signal. Every symbol the radio transmits is a phasor A~=Aeiϕ\tilde{A} = Ae^{i\phi}, a point in the complex plane chosen from a constellation diagram; every received symbol is demodulated by multiplying against e−iωte^{-i\omega t} and low-pass filtering — a hardware implementation of "factor out the rotation and freeze the phasor." Equation (3) is executed billions of times per second in your pocket. The same complex arithmetic runs the electric grid, where power engineers dispatch gigawatts by managing phasor angles between cities (with jj, following Steinmetz, since ii was taken by current), and monitor grid stability through GPS-synchronized phasor measurement units streaming A~\tilde{A} thirty times a second. The deepest tenants, though, are quantum. The wavefunction of a stationary state evolves as e−iEt/ℏe^{-iEt/\hbar} — not as a bookkeeping convenience but as physical law; interference of those complex phases is quantum mechanics. Lock-in amplifiers pull femtovolt signals out of noise a million times larger by multiplying against a reference eiωte^{i\omega t} and time-averaging — equation (10)'s constant term, weaponized. When a gravitational-wave observatory reports a chirp, the strain data was heterodyned, filtered, and matched — every step phasor arithmetic. Euler's spinning arrow turned out to be one of the load-bearing structures of the technological world.

The physical picture

Picture a clock hand of length AA spinning counterclockwise in a plane, sweeping out angle at a steady ω\omega radians per second. Shine a light from above and watch its shadow on the horizontal axis: the shadow oscillates back and forth, compressing near the ends, sweeping fastest through the middle — exactly simple harmonic motion. That plane is the complex plane, the spinning hand is A~eiωt\tilde{A}e^{i\omega t}, and the shadow is Re⁡(A~eiωt)=Acos⁡(ωt+ϕ)\operatorname{Re}(\tilde{A}e^{i\omega t}) = A\cos(\omega t + \phi). Every sinusoid you will ever meet is the shadow of uniform circular motion — which is why simple harmonic motion and circular motion share one ω\omega, and why the geometry of circles keeps solving the algebra of oscillations.

Figure
The spinning arrow and its shadow — a real oscillation is the projection of uniform rotation in the complex plane
Re Im A e^(iωt) ωt+φ ω the shadow t x(t) x(t) = A cos(ωt+φ) one spinning arrow carries both amplitude (its length) and phase (its angle); the cosine is its shadow, unrolled in time

Freeze the arrow at t=0t = 0 and you hold a single complex number A~=Aeiϕ\tilde{A} = Ae^{i\phi}: its length is the amplitude, its angle is the starting phase. This frozen snapshot is called a phasor, and it is the profound bookkeeping upgrade of this whole subject — the two real numbers you used to carry separately (AA and ϕ\phi, or equivalently the coefficients of cos⁡ωt\cos\omega t and sin⁡ωt\sin\omega t) fuse into one algebraic object. Multiplying by eiωte^{i\omega t} sets it spinning; multiplying by any other complex number stretches it and turns it. All the algebra of oscillations becomes the geometry of rotating arrows. Why is this so much easier than trigonometry? Because cos⁡\cos and sin⁡\sin are two functions that the derivative shuffles into each other, while eiωte^{i\omega t} is a single function that the derivative merely rescales. In the language of linear algebra: the exponential is an eigenfunction of d/dt\mathrm{d}/\mathrm{d}t with eigenvalue iωi\omega. Working in the eigenbasis of the operator that runs the whole problem diagonalizes everything — the same reason you diagonalize matrices, choose normal modes, or Fourier-transform a signal. And phase, that perennial source of sign errors, becomes something you can see: the angle of an arrow. When a driven oscillator's response lags its drive, the response phasor literally hangs at an angle behind the drive phasor, and sweeping the drive frequency rotates it smoothly from 0°0° through 90°90° at resonance to 180°180° — a rotation that is invisible in real-trig algebra and obvious in the complex plane.

Euler's formula and the eigenfunction property

Everything rests on one identity. It deserves a slow look before we put it to work, because it is the bridge between the exponential function — the natural language of growth and decay — and the circular functions — the natural language of oscillation. The claim is that a purely imaginary exponent converts growth into rotation:

eiθ=cos⁡θ+isin⁡θ(1)e^{i\theta} = \cos\theta + i\sin\theta \qquad (1)

where: ii — the imaginary unit, i2=−1i^2 = -1; θ\theta — a real angle in radians. Read it geometrically: eiθe^{i\theta} is the point on the unit circle at angle θ\theta, with the cosine as its horizontal coordinate and the sine as its vertical coordinate. Setting θ=π\theta = \pi gives the famous special case eiπ=−1e^{i\pi} = -1: half a turn lands you at minus one.

Derivation AEuler's formula from the power series

We are deriving Euler's formula (1), the identity on which the entire complex-exponential method stands. The strategy is the cleanest one available: expand eiθe^{i\theta} in its Maclaurin series — which converges for every complex argument, so no analytic caution is needed — and watch the powers of ii sort the series into two real sub-series. The thing to track is the four-step cycle i,−1,−i,+1i, -1, -i, +1 of the powers of ii: it is precisely this period-4 rhythm that carves the exponential into the period-4 derivatives of cosine and sine. No differential equations, no integrals — just term-by-term bookkeeping.

Step 1. Expand the exponential in its Maclaurin series, valid for all complex arguments:

eiθ=∑n=0∞(iθ)nn!=1+iθ+(iθ)22!+(iθ)33!+(iθ)44!+⋯e^{i\theta} = \sum_{n=0}^{\infty}\frac{(i\theta)^n}{n!} = 1 + i\theta + \frac{(i\theta)^2}{2!} + \frac{(i\theta)^3}{3!} + \frac{(i\theta)^4}{4!} + \cdots

Why this is safe: the exponential series has infinite radius of convergence, so it defines eze^{z} for any complex zz — this series is the definition being used.

Step 2. Evaluate the powers of ii and sort by parity. Even powers give i2k=(−1)ki^{2k} = (-1)^k, odd powers give i2k+1=i(−1)ki^{2k+1} = i(-1)^k:

eiθ=∑k=0∞(−1)kθ2k(2k)!+i∑k=0∞(−1)kθ2k+1(2k+1)!e^{i\theta} = \sum_{k=0}^{\infty}\frac{(-1)^k\theta^{2k}}{(2k)!} + i\sum_{k=0}^{\infty}\frac{(-1)^k\theta^{2k+1}}{(2k+1)!}

Why the split is legal: absolutely convergent series may be rearranged freely, so grouping even and odd terms changes nothing.

Step 3. Recognize the two real series: the first is exactly the Maclaurin series of cos⁡θ\cos\theta, the second of sin⁡θ\sin\theta. Therefore eiθ=cos⁡θ+isin⁡θe^{i\theta} = \cos\theta + i\sin\theta. Verified: at θ=0\theta = 0 both sides give 1; at θ=π\theta = \pi the right side gives −1+0i-1 + 0i, reproducing eiπ=−1e^{i\pi} = -1; differentiating both sides with respect to θ\theta and multiplying by −i-i returns the original identity, confirming self-consistency.

With Euler's formula in hand, the property that powers the whole method is one differentiation away. Ask what the derivative operator does to a complex exponential of time — the answer is the reason this page exists:

ddt eiωt=iω eiωt(2)\frac{\mathrm{d}}{\mathrm{d}t}\,e^{i\omega t} = i\omega\,e^{i\omega t} \qquad (2)

where: ω\omega — angular frequency (rad/s); tt — time. The function comes back unchanged in shape, merely multiplied by iωi\omega: it is an eigenfunction of d/dt\mathrm{d}/\mathrm{d}t. Iterating, dneiωt/dtn=(iω)neiωt\mathrm{d}^n e^{i\omega t}/\mathrm{d}t^n = (i\omega)^n e^{i\omega t}, so an nn-th order linear ODE acting on eiωte^{i\omega t} becomes a degree-nn polynomial in iωi\omega — calculus has been demoted to algebra. Note the geometric meaning of the eigenvalue: multiplying by ii rotates a complex number by 90°90°, so the velocity phasor leads the position phasor by a quarter turn and is ω\omega times longer, exactly the familiar SHM relation vmax⁡=ωAv_{\max} = \omega A with its quarter-period lead. The dictionary between the complex bookkeeping and the physical motion is the projection. A real oscillation is the real part of a spinning phasor:

x(t)=Acos⁡(ωt+ϕ)=Re⁡(A~ eiωt),A~=Aeiϕ(3)x(t) = A\cos(\omega t + \phi) = \operatorname{Re}\big(\tilde{A}\,e^{i\omega t}\big), \qquad \tilde{A} = Ae^{i\phi} \qquad (3)

where: A=∣A~∣≥0A = |\tilde{A}| \ge 0 — the real amplitude; ϕ=arg⁡A~\phi = \arg\tilde{A} — the phase offset; A~\tilde{A} — the complex amplitude, or phasor. One complex number now stores what used to take two reals. In particular the initial conditions live entirely in A~\tilde{A}: the real part of A~\tilde{A} is x(0)x(0), and its imaginary part is −x˙(0)/ω-\dot{x}(0)/\omega. There is a second, more symmetric way to write the projection, and it repays attention because it is the seed of Fourier analysis. Using Re⁡(z)=12(z+z∗)\operatorname{Re}(z) = \tfrac12(z + z^{*}) for any complex number,

Re⁡(A~eiωt)=12 A~ eiωt+12 A~∗e−iωt(4)\operatorname{Re}\big(\tilde{A}e^{i\omega t}\big) = \tfrac12\,\tilde{A}\,e^{i\omega t} + \tfrac12\,\tilde{A}^{*}e^{-i\omega t} \qquad (4)

where: A~∗\tilde{A}^{*} — the complex conjugate of A~\tilde{A}. A real oscillation is the sum of two counter-rotating phasors: one spinning at +ω+\omega, its conjugate twin at −ω-\omega. The "negative frequency" that puzzles every student meeting a two-sided Fourier spectrum is nothing exotic — it is the clockwise-spinning half that any purely real signal must carry so that the imaginary parts cancel at every instant.

Solving equations by algebra

Now put the machine to work on the equation that owns the rest of this track. A mass mm on a spring of stiffness kk with a drag force proportional to velocity (bx˙b\dot{x}) obeys mz¨+bz˙+kz=0m\ddot{z} + b\dot{z} + kz = 0. Rather than guessing combinations of sines, cosines, and decaying envelopes, feed it the one function that differentiation cannot deform — a trial exponential z=estz = e^{st} with ss complex and unknown:

z(t)=c1es1t+c2es2t,ms2+bs+k=0(5)z(t) = c_1 e^{s_1 t} + c_2 e^{s_2 t}, \qquad ms^2 + bs + k = 0 \qquad (5)

where: s1,2s_{1,2} — the two roots of the characteristic equation ms2+bs+k=0ms^2 + bs + k = 0; c1,c2c_1, c_2 — complex constants fixed by the initial position and velocity. Substituting este^{st} turns each time-derivative into a factor of ss, leaving (ms2+bs+k)est=0(ms^2 + bs + k)e^{st} = 0; since an exponential never vanishes, the polynomial must. The quadratic formula then delivers the entire physics of damping in one line: s=(−b±b2−4mk)/2ms = \big(-b \pm \sqrt{b^2 - 4mk}\big)/2m. When b2<4mkb^2 < 4mk the roots are complex, s=−β±iω02−β2s = -\beta \pm i\sqrt{\omega_0^2 - \beta^2} with β=b/2m\beta = b/2m and ω0=k/m\omega_0 = \sqrt{k/m} — a spiral in the complex plane: rotation (oscillation) at the reduced frequency, times shrinkage (decay) at rate β\beta. When b2>4mkb^2 > 4mk both roots are real and negative — pure decay, no rotation, no oscillation. The full story of damping regimes is the story of where these two roots sit.

Figure
Where the roots live — the characteristic roots in the complex plane classify the motion at a glance
Re s Im s growth — unstable decay — stable β = 0 : s = ±iω₀ pure oscillation underdamped: s = −β ± iω₁ roots ride the circle of radius ω₀ critical: s = −ω₀ twice overdamped: two real roots split along the axis distance from the imaginary axis = decay rate; height above the real axis = oscillation frequency

The same trial-exponential move solves the driven problem, and this is where the method pays its full fare. The physical equation has a cosine on the right-hand side; we complexify it, solve by division, and project back.

Go deeper · Derivation B · 2 minThe driven damped oscillator solved by complexification
Derivation BThe driven damped oscillator solved by complexification

The goal is the steady-state response of x¨+2βx˙+ω02x=(F0/m)cos⁡ωt\ddot{x} + 2\beta\dot{x} + \omega_0^2 x = (F_0/m)\cos\omega t — the workhorse equation of resonance physics — and we want it with its amplitude and phase lag in closed form. The strategy: replace the cosine drive by the complex drive (F0/m)eiωt(F_0/m)e^{i\omega t} whose real part it is, solve the complexified equation with a one-term ansatz, and take the real part at the end; linearity with real coefficients guarantees the projection is legitimate. What used to be a two-page exercise in matching coefficients of cos⁡ωt\cos\omega t and sin⁡ωt\sin\omega t becomes a single complex division. The step to watch is the last one: the phase lag must be extracted from the argument of a complex number, and its quadrant matters — the lag runs from 00 to 180°180°, not −90°-90° to 90°90°.

Step 1. Complexify the equation. Introduce z(t)z(t) obeying

z¨+2βz˙+ω02z=F0m eiωt\ddot{z} + 2\beta\dot{z} + \omega_0^2 z = \frac{F_0}{m}\,e^{i\omega t}

Why this is allowed: cos⁡ωt=Re⁡ eiωt\cos\omega t = \operatorname{Re}\,e^{i\omega t}, and the differential operator has real coefficients, so taking the real part of this entire equation reproduces the physical one — meaning x=Re⁡ zx = \operatorname{Re}\,z solves the original problem. The imaginary part quietly solves the sine-driven twin and never contaminates our answer.

Step 2. Insert the steady-state ansatz z(t)=A~ eiωtz(t) = \tilde{A}\,e^{i\omega t}. Why one term suffices: after transients die (they decay as e−βte^{-\beta t}, the homogeneous solutions of (5)), a linear system driven at frequency ω\omega can only respond at frequency ω\omega — all that is left to determine is one complex constant A~\tilde{A}.

Step 3. Differentiate by multiplication. Each dot becomes a factor iωi\omega:

(−ω2+2iβω+ω02) A~ eiωt=F0m eiωt\big(-\omega^2 + 2i\beta\omega + \omega_0^2\big)\,\tilde{A}\,e^{i\omega t} = \frac{F_0}{m}\,e^{i\omega t}

Why the exponential cancels: it is common to every term and never zero — this cancellation is the moment the differential equation officially becomes algebra.

Step 4. Solve for the phasor by one complex division:

A~=F0/mω02−ω2+2iβω\tilde{A} = \frac{F_0/m}{\omega_0^2 - \omega^2 + 2i\beta\omega}

Why this object is called a transfer function: it is the ratio of response to drive, and it depends only on the system (ω0,β\omega_0, \beta) and the probing frequency ω\omega — the entire input-output behavior of the oscillator in one complex-valued function of frequency.

Step 5. Extract modulus and argument. Multiplying numerator and denominator by the conjugate of the denominator,

∣A~∣=F0/m(ω02−ω2)2+4β2ω2,tan⁡δ=2βωω02−ω2|\tilde{A}| = \frac{F_0/m}{\sqrt{(\omega_0^2 - \omega^2)^2 + 4\beta^2\omega^2}}, \qquad \tan\delta = \frac{2\beta\omega}{\omega_0^2 - \omega^2}

Why the quadrant convention: δ\delta is defined as the lag of the response behind the drive, A~=∣A~∣e−iδ\tilde{A} = |\tilde{A}|e^{-i\delta}. Below resonance the denominator's real part is positive and δ\delta sits in the first quadrant; above resonance the real part flips sign and δ\delta climbs into the second quadrant, approaching 180°180° — take the arctangent in the correct quadrant or the high-frequency physics comes out wrong.

Step 6. Project back to the physical motion:

xss(t)=Re⁡(A~eiωt)=F0/m(ω02−ω2)2+4β2ω2 cos⁡(ωt−δ)x_{\mathrm{ss}}(t) = \operatorname{Re}\big(\tilde{A}e^{i\omega t}\big) = \frac{F_0/m}{\sqrt{(\omega_0^2 - \omega^2)^2 + 4\beta^2\omega^2}}\,\cos(\omega t - \delta)

Why the lag is physical: the oscillator needs time to respond to the push; at resonance (ω=ω0\omega = \omega_0) the denominator is purely imaginary, δ=90°\delta = 90° exactly, and the drive force is aligned with the velocity — optimal power delivery, the fingerprint of resonance. *Verified: matches Taylor ch. 5 and Goldstein ch. 6; as β→0\beta \to 0 with ω→ω0\omega \to \omega_0 the amplitude diverges, recovering the undamped resonance catastrophe; at ω→0\omega \to 0 it reduces to the static stretch F0/mω02=F0/kF_0/m\omega_0^2 = F_0/k.*

For steady-state work the two payoff formulas deserve their own display. The response phasor of the driven oscillator, and its polar decomposition, are

A~=F0/mω02−ω2+2iβω(6)\tilde{A} = \frac{F_0/m}{\omega_0^2 - \omega^2 + 2i\beta\omega} \qquad (6)

where: F0F_0 — drive amplitude; ω0=k/m\omega_0 = \sqrt{k/m} — natural frequency; β=b/2m\beta = b/2m — damping rate; ω\omega — drive frequency. One complex ratio holds the complete steady-state response at every drive frequency, and

∣A~∣=F0/m(ω02−ω2)2+4β2ω2,tan⁡δ=2βωω02−ω2(7)|\tilde{A}| = \frac{F_0/m}{\sqrt{(\omega_0^2 - \omega^2)^2 + 4\beta^2\omega^2}}, \qquad \tan\delta = \frac{2\beta\omega}{\omega_0^2 - \omega^2} \qquad (7)

where: δ\delta — the phase lag of response behind drive. Amplitude is a modulus, phase is an argument: the resonance curve and the phase curve that fill the next articles of this track are just ∣A~(ω)∣|\tilde{A}(\omega)| and arg⁡A~(ω)\arg\tilde{A}(\omega) read off equation (6).

Which equation answers which question

you wantusenotes
turn trig into algebraeiθ=cos⁡θ+isin⁡θe^{i\theta} = \cos\theta + i\sin\theta, eq. (1)the master identity
differentiate a sinusoidd/dt→iω\mathrm{d}/\mathrm{d}t \to i\omega, eq. (2)eigenfunction property
encode amplitude and phaseA~=Aeiϕ\tilde{A} = Ae^{i\phi}, eq. (3)one phasor = two reals
explain negative frequencyconjugate pair, eq. (4)real signal = two counter-rotating phasors
free damped motioncharacteristic roots, eq. (5)root location = damping regime
driven steady statetransfer function, eqs. (6)-(7)modulus = amplitude, argument = phase lag
add same-frequency sinusoidsphasor sum, eq. (8)law of cosines, no identities
energies and powersproduct rule, eq. (9)project FIRST, then multiply

Adding oscillations: phasor arithmetic

Superposition is where the representation stops being a convenience and becomes a superpower. Add two oscillations of the same frequency with different amplitudes and phases — the kind of sum that appears in interference, in coupled oscillators, in AC circuit nodes — and trigonometry demands the sum-to-product identities. Phasors demand only that you add two arrows:

A1cos⁡(ωt+ϕ1)+A2cos⁡(ωt+ϕ2)=Re⁡(A~ eiωt),A~=A1eiϕ1+A2eiϕ2(8)A_1\cos(\omega t + \phi_1) + A_2\cos(\omega t + \phi_2) = \operatorname{Re}\big(\tilde{A}\,e^{i\omega t}\big), \qquad \tilde{A} = A_1e^{i\phi_1} + A_2e^{i\phi_2} \qquad (8)

where: A~\tilde{A} — the resultant phasor, obtained by ordinary vector addition in the complex plane. The resultant amplitude follows from the law of cosines, A=A12+A22+2A1A2cos⁡(ϕ1−ϕ2)A = \sqrt{A_1^2 + A_2^2 + 2A_1A_2\cos(\phi_1 - \phi_2)}: constructive when the arrows align, destructive when they oppose, Pythagoras when they are perpendicular.

Figure
Adding oscillations is adding arrows — two phasors and their resultant, drawn with the numbers of worked example 1
A₁ = 3 cm, φ₁ = 0° A₂ = 4 cm, φ₂ = 90° A = 5 cm φ = 53° tip-to-tail in the complex plane resultant length: law of cosines at 90° the cross term vanishes, leaving Pythagoras: 3-4-5 aligned arrows: A = 7 (constructive) opposed arrows: A = 1 (destructive) no sum-to-product identities were harmed: adding sinusoids of one frequency is vector addition of their phasors
Go deeper · Derivation C · 1 minTwo cosines into one by adding phasors
Derivation CTwo cosines into one by adding phasors

We are deriving the superposition rule (8): two sinusoids of the same frequency combine into a single sinusoid, with the resultant amplitude and phase given by vector addition of phasors. The traditional route expands both cosines, collects the coefficients of cos⁡ωt\cos\omega t and sin⁡ωt\sin\omega t, and re-assembles them with a half-page of identities; the phasor route is four lines, and — more valuably — it explains why the answer is the law of cosines: because we literally are adding two arrows. Watch how the shared factor eiωte^{i\omega t} carries all the time dependence, leaving the geometry frozen at t=0t = 0.

Step 1. Write each oscillation as the real part of a spinning phasor sharing the same eiωte^{i\omega t}:

A1cos⁡(ωt+ϕ1)+A2cos⁡(ωt+ϕ2)=Re⁡(A1eiϕ1eiωt)+Re⁡(A2eiϕ2eiωt)A_1\cos(\omega t + \phi_1) + A_2\cos(\omega t + \phi_2) = \operatorname{Re}\big(A_1e^{i\phi_1}e^{i\omega t}\big) + \operatorname{Re}\big(A_2e^{i\phi_2}e^{i\omega t}\big)

Why the same frequency matters: only then can the two terms share one exponential factor — phasor addition is a same-frequency privilege.

Step 2. Use linearity of the real part to combine under one Re⁡\operatorname{Re}, and factor out the common rotation:

=Re⁡[(A1eiϕ1+A2eiϕ2)eiωt]=Re⁡(A~ eiωt)= \operatorname{Re}\Big[\big(A_1e^{i\phi_1} + A_2e^{i\phi_2}\big)e^{i\omega t}\Big] = \operatorname{Re}\big(\tilde{A}\,e^{i\omega t}\big)

Why this is the whole trick: the bracket is a single complex constant A~\tilde{A} — the time dependence never touched the addition.

Step 3. Compute the resultant's magnitude from ∣A~∣2=A~A~∗|\tilde{A}|^2 = \tilde{A}\tilde{A}^{*}:

A=∣A~∣=A12+A22+2A1A2cos⁡(ϕ1−ϕ2)A = |\tilde{A}| = \sqrt{A_1^2 + A_2^2 + 2A_1A_2\cos(\phi_1 - \phi_2)}

Why the law of cosines appears: the cross terms A1A2(ei(ϕ1−ϕ2)+e−i(ϕ1−ϕ2))=2A1A2cos⁡(ϕ1−ϕ2)A_1A_2(e^{i(\phi_1-\phi_2)} + e^{-i(\phi_1-\phi_2)}) = 2A_1A_2\cos(\phi_1-\phi_2) — the phasor diagram is a triangle, and this is its third side.

Step 4. Compute the resultant's phase from the real and imaginary parts of A~\tilde{A}:

ϕ=arg⁡A~=arctan⁡A1sin⁡ϕ1+A2sin⁡ϕ2A1cos⁡ϕ1+A2cos⁡ϕ2\phi = \arg\tilde{A} = \arctan\frac{A_1\sin\phi_1 + A_2\sin\phi_2}{A_1\cos\phi_1 + A_2\cos\phi_2}

Verified: setting ϕ1=ϕ2\phi_1 = \phi_2 gives A=A1+A2A = A_1 + A_2 (constructive); ϕ1−ϕ2=180°\phi_1 - \phi_2 = 180° gives A=∣A1−A2∣A = |A_1 - A_2| (destructive); the general case reproduces the classical sum-to-product result, cross-checked in Morin ch. 4.

What if the two frequencies are close but not equal? Then the two phasors rotate at slightly different rates, and their relative angle drifts slowly: the arrows alternately align and oppose, so the resultant amplitude swells and shrinks. Factoring the sum of two equal-amplitude cosines makes this audible:

cos⁡ω1t+cos⁡ω2t=2cos⁡ ⁣(ω1−ω22 t)cos⁡ ⁣(ω1+ω22 t)(9)\cos\omega_1 t + \cos\omega_2 t = 2\cos\!\Big(\frac{\omega_1 - \omega_2}{2}\,t\Big)\cos\!\Big(\frac{\omega_1 + \omega_2}{2}\,t\Big) \qquad (9)

where: the second factor oscillates at the average frequency, while the first is a slow envelope at half the difference frequency. Your ear hears one loudness maximum per envelope half-cycle, so the perceived beat frequency is the full difference ∣f1−f2∣|f_1 - f_2| — two guitar strings at 440440 and 442 Hz442\ \mathrm{Hz} throb twice per second, and a tuner's job is to slow that throb to zero. One last identity carries the method's only serious danger, and it earns its own display precisely because everything up to here has encouraged you to forget the projection and compute with zz directly. That habit is safe for anything linear — and fatal for anything quadratic:

Re⁡(z1) Re⁡(z2)=12Re⁡(z1z2)+12Re⁡(z1z2∗)(10)\operatorname{Re}(z_1)\,\operatorname{Re}(z_2) = \tfrac12\operatorname{Re}\big(z_1 z_2\big) + \tfrac12\operatorname{Re}\big(z_1 z_2^{*}\big) \qquad (10)

where: z1,z2z_1, z_2 — any two complex signals. The product of real parts is not the real part of the product: writing each real part as a conjugate average 12(z+z∗)\tfrac12(z + z^{*}) and multiplying out produces cross terms z1z2∗z_1z_2^{*} that Re⁡(z1z2)\operatorname{Re}(z_1z_2) alone misses. For two phasors at the same frequency, z1z2z_1z_2 spins at 2ω2\omega and time-averages to zero, while z1z2∗z_1z_2^{*} is constant — which yields the single most used corollary in physics: the time-averaged product of two same-frequency oscillations is ⟨x1x2⟩=12Re⁡(A~1A~2∗)\langle x_1x_2\rangle = \tfrac12\operatorname{Re}(\tilde{A}_1\tilde{A}_2^{*}), and in particular ⟨x2⟩=12∣A~∣2\langle x^2\rangle = \tfrac12|\tilde{A}|^2.

Limiting cases

CaseResultMeaning
undamped free oscillationz=A~eiω0tz = \tilde{A}e^{i\omega_0 t}, s=±iω0s = \pm i\omega_0roots on the imaginary axis: the phasor spins on a circle of fixed radius, forever
underdampeds=−β±iω02−β2s = -\beta \pm i\sqrt{\omega_0^2 - \beta^2}complex roots: rotation plus shrinkage — a spiral; modulus decays as e−βte^{-\beta t}
overdampeds=−β±β2−ω02s = -\beta \pm \sqrt{\beta^2 - \omega_0^2}, both realno imaginary part, no rotation: oscillation degenerates into two pure decays
driven, ω→0\omega \to 0A~→F0/mω02=F0/k\tilde{A} \to F_0/m\omega_0^2 = F_0/k, δ→0\delta \to 0quasi-static limit: the mass tracks the force in phase, stretched by the spring law
driven, ω=ω0\omega = \omega_0A~=F0/(2iβω0m)\tilde{A} = F_0/(2i\beta\omega_0 m), δ=90°\delta = 90°purely imaginary denominator: response lags by a quarter period, drive feeds velocity
driven, ω→∞\omega \to \infty∣A~∣→F0/mω2|\tilde{A}| \to F_0/m\omega^2, δ→180°\delta \to 180°mass can't keep up: response opposite the drive, amplitude dying as 1/ω21/\omega^2
equal amplitudes, nearby frequencieseq. (9)beats: envelope at 12∣ω1−ω2∣\tfrac12\vert\omega_1 - \omega_2\vert, heard as ∣f1−f2∣\vert f_1 - f_2\vert throbs per second

Worked examples

1 · Two speakers, one microphone (confidence builder)

Two signals of the same frequency arrive at a microphone: one with amplitude A1=3 cmA_1 = 3\ \mathrm{cm} (in suitable units) and phase ϕ1=0°\phi_1 = 0°, the other with A2=4 cmA_2 = 4\ \mathrm{cm} and phase ϕ2=90°\phi_2 = 90°. Find the amplitude and phase of the combined oscillation.

Step 1. Encode each signal as a phasor: A~1=3ei0=3\tilde{A}_1 = 3e^{i0} = 3, and A~2=4eiπ/2=4i\tilde{A}_2 = 4e^{i\pi/2} = 4i — one arrow along the real axis, one along the imaginary axis.

Step 2. Add the arrows: A~=3+4i\tilde{A} = 3 + 4i. No identities, no angle-chasing — complex addition is the entire calculation.

Step 3. Read off the amplitude as the modulus: A=∣A~∣=32+42=25=5 cmA = |\tilde{A}| = \sqrt{3^2 + 4^2} = \sqrt{25} = \boxed{5\ \mathrm{cm}}. The law-of-cosines cross term 2A1A2cos⁡(ϕ1−ϕ2)2A_1A_2\cos(\phi_1 - \phi_2) vanished because the arrows are perpendicular — this is Pythagoras in phasor clothing, the 3-4-5 triangle of the figure above.

Step 4. Read off the phase as the argument: ϕ=arctan⁡(4/3)=53°\phi = \arctan(4/3) = \boxed{53°} (first quadrant, since both components are positive). The combined signal is x(t)=5cos⁡(ωt+53°) cmx(t) = 5\cos(\omega t + 53°)\ \mathrm{cm}.

2 · A driven cart, three lines instead of two pages (exam level)

A cart of mass m=0.5 kgm = 0.5\ \mathrm{kg} on a spring k=200 N/mk = 200\ \mathrm{N/m} moves through oil giving a drag coefficient b=4 kg/sb = 4\ \mathrm{kg/s}. It is driven by F(t)=10cos⁡(19 t) NF(t) = 10\cos(19\,t)\ \mathrm{N} (SI units). Find the steady-state amplitude and phase lag, and compare the amplitude with the static stretch produced by a constant 10 N force.

Step 1. Assemble the system constants: ω0=k/m=400=20 rad/s\omega_0 = \sqrt{k/m} = \sqrt{400} = 20\ \mathrm{rad/s}; β=b/2m=4/(2×0.5)=4 s−1\beta = b/2m = 4/(2 \times 0.5) = 4\ \mathrm{s^{-1}}; drive ω=19 rad/s\omega = 19\ \mathrm{rad/s}, just below resonance; F0/m=20 N/kgF_0/m = 20\ \mathrm{N/kg}.

Step 2. Build the complex denominator of equation (6): real part ω02−ω2=400−361=39 s−2\omega_0^2 - \omega^2 = 400 - 361 = 39\ \mathrm{s^{-2}}; imaginary part 2βω=2×4×19=152 s−22\beta\omega = 2 \times 4 \times 19 = 152\ \mathrm{s^{-2}}. Note the imaginary part dominates — we are deep in the resonance region.

Step 3. Amplitude = modulus: ∣A~∣=20392+1522=201521+23104=20156.9=0.127 m≈12.7 cm|\tilde{A}| = \dfrac{20}{\sqrt{39^2 + 152^2}} = \dfrac{20}{\sqrt{1521 + 23104}} = \dfrac{20}{156.9} = \boxed{0.127\ \mathrm{m} \approx 12.7\ \mathrm{cm}}.

Step 4. Phase lag = argument: δ=arctan⁡(152/39)=75.6°\delta = \arctan(152/39) = \boxed{75.6°} — the response trails the drive by just over a fifth of a cycle, on its way to the 90°90° lag it would show exactly at ω0\omega_0.

Step 5. Compare with statics: a constant 10 N force stretches the spring by F0/k=10/200=0.05 m=5 cmF_0/k = 10/200 = 0.05\ \mathrm{m} = 5\ \mathrm{cm}. Driving near resonance amplifies this to 12.7 cm12.7\ \mathrm{cm} — a factor of 2.52.5, which is no accident: the quality factor of this oscillator is Q=ω0/2β=20/8=2.5Q = \omega_0/2\beta = 20/8 = 2.5, and QQ is precisely the resonant amplification of the static response. The complex division knew all of resonance physics before we asked.

Origins · 1 minLeonhard Euler and how this was discovered
LE
Leonhard Euler
1707–1783 · Basel, St. Petersburg and Berlin academies

Euler is the most prolific mathematician who ever lived — his collected works run past seventy volumes — and an astonishing share of the symbols on this page are his: ee, ii, f(x)f(x), π\pi in its modern role, even sin⁡\sin and cos⁡\cos treated as functions rather than chord lengths. The identity eiθ=cos⁡θ+isin⁡θe^{i\theta} = \cos\theta + i\sin\theta appears in his 1748 masterpiece Introductio in analysin infinitorum, derived exactly as in this article, from the series. Euler computed on regardless of circumstance: he lost sight in one eye in his thirties and in both by 1771, then increased his output, dictating calculations — including a complete lunar theory — to assistants from memory. He worked at the St. Petersburg and Berlin academies, returned to Russia in 1766, and was calculating the orbit of the newly discovered Uranus on the afternoon he died. The formula joining growth, rotation, and oscillation into one exponential is his monument; physics has been cashing it ever since.

History
1714
— Roger Cotes publishes iθ=ln⁡(cos⁡θ+isin⁡θ)i\theta = \ln(\cos\theta + i\sin\theta), Euler's formula in logarithmic disguise, decades early and largely unnoticed.
1748
— Euler's Introductio in analysin infinitorum derives eiθ=cos⁡θ+isin⁡θe^{i\theta} = \cos\theta + i\sin\theta from the power series and displays the special case eiπ=−1e^{i\pi} = -1.
1806
— Jean-Robert Argand publishes the geometric picture of complex numbers as points and arrows in a plane (Caspar Wessel had it in 1799, in a Danish memoir nobody read), turning Euler's identity into rotating-vector geometry.
1822
— Joseph Fourier's Théorie analytique de la chaleur shows arbitrary signals decompose into sinusoids — complex exponentials become the universal basis of linear physics.
1893
— Charles Proteus Steinmetz presents the phasor method to the International Electrical Congress in Chicago; AC circuit analysis collapses from differential equations into complex arithmetic, and electrical engineering adopts jj overnight.

Going further

References. Morin, Introduction to Classical Mechanics, ch. 4 (Oscillations) · Goldstein, Poole & Safko, Classical Mechanics, 3rd ed., ch. 6 · Landau & Lifshitz, Mechanics, 3rd ed., §21–§26 · Taylor, Classical Mechanics, ch. 5 · Feynman, Lectures on Physics, vol. I, chs. 22–24 (algebra, resonance, transients).

Provenance. Euler's formula and series expansion verifiable from any analysis text. Driven steady-state amplitude/phase match Taylor Classical Mechanics ch. 5 and Goldstein ch. 6; damped characteristic roots standard. Phasor superposition reduces to the law of cosines. A verifier should recompute the complex division for A~\tilde A, confirm ⟨cos⁡2⟩=1/2\langle\cos^2\rangle=1/2, and check the phase limits δ(0)=0, δ(ω0)=90∘, δ(∞)=180∘\delta(0)=0,\ \delta(\omega_0)=90^\circ,\ \delta(\infty)=180^\circ. Worked-example numbers on this page (A=5 cm, ϕ=53°A = 5\ \mathrm{cm},\ \phi = 53° for the phasor sum; ∣A~∣=12.7 cm, δ=75.6°, Q=2.5|\tilde{A}| = 12.7\ \mathrm{cm},\ \delta = 75.6°,\ Q = 2.5 for the driven cart) are illustrative arithmetic recomputed from the boxed formulas.

Builds on: Simple Harmonic Motion · Amplitude, Phase, and Angular Frequency · Small Oscillations About Equilibrium. Leads to: Damped Harmonic Oscillator · Driven Oscillations · Resonance · Transient and Steady-State Response.

Is the imaginary part of the solution physically real? What is it doing there?

It is a perfectly real second solution, riding along unused. If z=x+iyz = x + iy solves the complexified equation with drive Re⁡+iIm⁡\operatorname{Re} + i\operatorname{Im}, then xx solves the cosine-driven problem and yy solves the sine-driven one — two independent physical problems processed in one calculation. Nothing about yy contaminates xx, because the equation is linear with real coefficients. Geometrically, the imaginary part is the quarter-period-shifted twin of your oscillation: the shadow of the same spinning arrow on the vertical axis instead of the horizontal one. You solve two problems and keep the one you asked for.

Why the exponential function specifically? What makes it so special for differential equations?

It is the eigenfunction of the derivative: este^{st} is the only function whose derivative is a fixed multiple of itself for all time. Constant-coefficient linear ODEs are built entirely from derivatives, so exponentials pass through them intact, picking up only polynomial factors in ss — the differential equation becomes the characteristic polynomial. Sines and cosines fail this test individually (each differentiates into the other), but the combination cos⁡θ+isin⁡θ\cos\theta + i\sin\theta passes it — and Euler's formula says that combination is an exponential. In linear-algebra language: you are diagonalizing d/dt\mathrm{d}/\mathrm{d}t, exactly as you would diagonalize a matrix to decouple a system.

Can I use the complex trick on a pendulum at large angles?

No — and the reason is the licensing clause. The method requires linearity: Re⁡\operatorname{Re} commutes with derivatives and sums, but not with nonlinear operations, so a term like sin⁡x\sin x in x¨+ω02sin⁡x=0\ddot{x} + \omega_0^2\sin x = 0 destroys the guarantee that the real part of the complex solution solves the real equation (sin⁡(Re⁡z)≠Re⁡(sin⁡z)\sin(\operatorname{Re} z) \neq \operatorname{Re}(\sin z)). At small angles you linearize sin⁡x≈x\sin x \approx x and the machine applies; beyond that you need the genuinely nonlinear toolkit — elliptic integrals, perturbation series — of the anharmonic pendulum article. The complex exponential is the reward for linearity, not a universal solvent.

Why do electrical engineers write j instead of i, and is their method different?

Identical method, colliding notation: ii was already the standard symbol for electric current, so Steinmetz's community adopted jj for the imaginary unit. An engineer's impedance calculation — resistor RR, inductor jωLj\omega L, capacitor 1/(jωC)1/(j\omega C), combined by series-parallel algebra — is precisely this article's transfer-function move: each impedance is the phasor ratio of voltage to current, i.e. the factor a component applies to ejωte^{j\omega t}. Ohm's law with complex numbers is equation (6) wearing a hard hat. When your circuits course and your mechanics course feel eerily parallel, that is because they are the same linear mathematics.

What does a negative frequency even mean in equation 4?

A direction of spin, nothing more mysterious. The phasor e+iωte^{+i\omega t} rotates counterclockwise; its conjugate e−iωte^{-i\omega t} rotates clockwise — "negative frequency" labels the second sense of rotation. A purely real signal must contain both in conjugate-balanced amounts so their imaginary parts cancel at every instant; that is exactly the content of Re⁡(A~eiωt)=12A~eiωt+12A~∗e−iωt\operatorname{Re}(\tilde{A}e^{i\omega t}) = \tfrac12\tilde{A}e^{i\omega t} + \tfrac12\tilde{A}^{*}e^{-i\omega t}. When you later meet a two-sided Fourier spectrum with peaks at ±ω\pm\omega, you are seeing the two counter-rotating halves of each real sinusoid displayed separately. Only complex-valued signals (like the I/Q stream inside a radio) can carry an unbalanced, genuinely one-sided spectrum.

My textbook writes x equals A e to the iωt and never takes a real part. Is that sloppy?

It is a convention, stated once and then suppressed: "the physical quantity is the real part of everything that follows." For any linear manipulation — solving, adding, differentiating — carrying the complex object and projecting at the very end gives identical results, so authors stop writing Re⁡\operatorname{Re} to reduce clutter. The convention is safe until the moment you form a product: energy, power, intensity, anything quadratic must go through the real signal or the 12Re⁡(A~1A~2∗)\tfrac12\operatorname{Re}(\tilde{A}_1\tilde{A}_2^{*}) average rule of equation (10). Good texts flag exactly that exception; good students remember it on the exam. When in doubt, write the Re⁡\operatorname{Re} back in and check what commutes.

Is a phasor just a vector? Why use complex numbers instead of 2D vectors?

For addition, yes — phasors add tip-to-tail exactly like plane vectors, which is why superposition becomes geometry. The upgrade is multiplication: complex numbers multiply by multiplying lengths and adding angles, an operation plane vectors simply do not have. That one extra structure is what lets a single division solve the driven oscillator — the transfer function multiplies the drive phasor to produce the response phasor, rotating it by the phase lag and rescaling it by the amplification, in one stroke. A 2D-vector formalism would need a separate rotation matrix bolted on; the complex plane has rotation built into its arithmetic. That is the entire engineering case for ii.

How do I get real initial conditions into the complex solution?

Through the complex constants, which hold exactly two real numbers each — that is why they are complex. For the free oscillator x=Re⁡(A~eiω0t)x = \operatorname{Re}(\tilde{A}e^{i\omega_0 t}): expand to get x(0)=Re⁡A~x(0) = \operatorname{Re}\tilde{A} and x˙(0)=Re⁡(iω0A~)=−ω0Im⁡A~\dot{x}(0) = \operatorname{Re}(i\omega_0\tilde{A}) = -\omega_0\operatorname{Im}\tilde{A}, so A~=x(0)−i x˙(0)/ω0\tilde{A} = x(0) - i\,\dot{x}(0)/\omega_0. Modulus and argument of that one number then hand you the amplitude A=x(0)2+(x˙(0)/ω0)2A = \sqrt{x(0)^2 + (\dot{x}(0)/\omega_0)^2} and phase in a single step — the same formulas you would grind out by matching cos⁡\cos and sin⁡\sin coefficients, delivered automatically. For the damped case the same matching fixes c1,c2c_1, c_2 in equation (5).

Before this

Where this leads