Skip to content

Harmonic Balance

Approximate a nonlinear system's periodic steady state by a truncated Fourier series and solve for coefficients whose retained residual harmonics balance to zero.

Version
v3 · 2026-09-06 · History
Domain-specific #
1978
Origin domain
applied mathematics
Subdomain
nonlinear periodic systems
Aliases
Harmonic balance method, Frequency-domain harmonic balance

Core Idea

Harmonic balance turns the search for a periodic steady state of a nonlinear differential or differential-algebraic system into a finite nonlinear algebra problem. It assumes that the state is represented by a truncated Fourier series, substitutes that representation into the governing equations, and requires the Fourier coefficients of the retained residual harmonics to vanish. Nakhla and Vlach established an influential circuit formulation in which frequency-domain linear behavior is balanced against nonlinear device currents.

The answer is approximate because nonlinearities generate harmonics outside the retained band. Increasing the harmonic order can improve accuracy, but only when nonlinear evaluation, sampling, aliasing control, conditioning, and nonlinear-solver convergence are also adequate.

Scope of Application

Harmonic balance is used for forced and self-sustained oscillations, radio-frequency and microwave circuits, power electronics, rotating machinery, structural vibration, and other systems whose long-run response is periodic or multiperiodic. Review treatments distinguish the representation and device-evaluation choices that make it practical for nonlinear circuits.

Clarity

State whether the system is autonomous or forced; the assumed fundamental frequencies; retained positive, negative, and zero modes; residual projection; nonlinear evaluation grid; oversampling or dealiasing rule; phase condition; nonlinear solver; stopping rule; and independent error check. A reported residual without those choices is not reproducible.

Manages Complexity

The method avoids simulating many decaying transients when only periodic steady state matters. In circuit applications, linear subnetworks can remain in the frequency domain while nonlinear devices are evaluated on a time grid through FFT transforms. Krylov methods can avoid explicitly factoring the full dense Jacobian for large systems.

Abstract Reasoning

  1. Choose a period or a vector of fundamental frequencies.
  2. Truncate the Fourier representation to a declared mode set.
  3. Express derivatives algebraically by multiplying each mode by its frequency.
  4. Evaluate nonlinear terms without contaminating retained modes through uncontrolled aliasing.
  5. Project the governing residual onto each retained basis function.
  6. Add a phase condition when time-translation symmetry makes an autonomous orbit nonunique.
  7. Solve the resulting nonlinear coefficient equations.
  8. validate against added harmonics, a time-domain residual, or an independent computation.

Knowledge Transfer

The portable pattern is choose a structured finite basis, enforce balance only in that basis, and test the information discarded by truncation. It transfers to spectral discretization and periodic-orbit computation. The proposed immediate parent is Approximation.

Relationships to Other Abstractions

Local relationship map for Harmonic BalanceParents appear above the current abstraction, mutual partners to the right, and children below. Node labels state whether each abstraction is prime or domain-specific; colors identify relation types.Harmonic BalanceDOMAINPrime abstraction: Approximation — is a kind ofApproximationPRIME

Current abstraction Harmonic Balance Domain-specific

Parents (1) — more general patterns this builds on

  • Harmonic Balance is a kind of Approximation Prime

    Approximation is the proposed immediate parent.

Hierarchy path (1) — routes to 1 parentless root

Neighborhood in Abstraction Space

Harmonic Balance sits in a sparse region of the domain-specific corpus (96th percentile for distinctiveness): few abstractions share its structure, so a faithful description tends to retrieve it precisely.

Family — Iterative Numerical Methods & Stability (7 abstractions)

Nearest neighbors

Computed from structural-signature embeddings · 2026-09-08