---
title: Instantaneous Spectrum Analysis
url: https://www.emergentmind.com/topics/instantaneous-spectrum-is
type: topic
---

# Instantaneous Spectrum Analysis

Instantaneous Spectrum (IS) denotes a class of time-local frequency representations for nonstationary signals, and, in some measurement systems, a single-shot physical mapping from frequency to space. In AM–FM signal models, IS assigns amplitude or energy to instantaneous-frequency trajectories; in synchrosqueezing, it is a concentrated and invertible reassignment of a time–frequency transform; in sparse and nonperiodic formulations, it is a sharp spectrum assembled from locally estimated modes; and in analog microwave analyzers, it is realized as an optical image of resonant response in a magnetic-field gradient [1504.07554][1404.7550][2508.05380][1509.01395].

## 1. Core definitions and formal variants

For a nonstationary signal written as
$$
f(t)=\sum_{k=1}^K A_k(t)e^{2\pi i\phi_k(t)},
$$
with slowly varying amplitudes $A_k(t)>0$ and strictly increasing phases $\phi_k(t)$, the instantaneous frequency (IF) of the $k$th mode is $\phi_k'(t)$. In this setting, the instantaneous spectrum at time $t$ is a distribution $S_f(t,\omega)$ that concentrates the energy of $f$ near the curves $\omega=\phi_k'(t)$; operationally, it answers the question “at time $t$, which frequencies $\omega$ carry energy in $f$?” [1404.7550].

A different but closely related formulation arises in the Hilbert-spectrum framework. Starting from a multicomponent latent signal
$$
z(t)=\sum_k A_k(t)e^{j\theta_k(t)}
=\sum_k A_k(t)e^{j\int^t\omega_k(\tau)\,d\tau},
$$
the instantaneous spectrum is defined by depositing the amplitude $A_k(t)$ at the instantaneous frequency $\omega_k(t)$,
$$
H(\omega,t)=\sum_{k=0}^{K-1}A_k(t)\,\delta\!\bigl(\omega-\omega_k(t)\bigr).
$$
This representation is described as exact and free of time–frequency smoothing once a particular AM–FM decomposition has been fixed [1504.07554].

A third formulation is atom-based. In the quadratic-chirplet framework, one defines a local transform
$$
V_x(\tau,\alpha,\beta)=\langle x(\cdot),\psi_{\gamma,\alpha,\beta}(\cdot;\tau)\rangle,
$$
and then sets
$$
IS_x^{(\gamma,\beta)}(\tau,f)=|V_x(\tau,\alpha=f,\beta)|^2.
$$
Here the IS is the squared magnitude of a local analysis coefficient indexed by time $\tau$, frequency $f$, chirp-rate parameter $\beta$, and window-spread parameter $\gamma$ [2508.05380].

| Framework | IS expression | Distinguishing feature |
|---|---|---|
| Hilbert spectrum | $H(\omega,t)=\sum_k A_k(t)\delta(\omega-\omega_k(t))$ | Exact AM–FM deposition |
| Synchrosqueezing | $T_f(\omega,b)$ | Reassigned, sparse, invertible |
| Sparse Frequency Analysis | $\mathrm{IS}(f,n)=\sum_k A_k[n]\delta(f-f_k[n])$ | TV-sparse amplitudes on a discrete grid |
| Quadratic-chirplet framework | $IS_x^{(\gamma,\beta)}(\tau,f)=|V_x(\tau,f,\beta)|^2$ | Unifies common analyses |

Taken together, these definitions suggest that IS is not a single universal object but a family of representations whose exact form depends on the signal model, transform, or measurement principle.

## 2. Exact AM–FM and Hilbert-spectrum formulation

In the Hilbert-spectrum theory, one begins with a real signal $x(t)\in\mathbb{R}$ and seeks a complex extension
$$
x_a(t)=z(t)=\sum_{k=0}^{K-1}\psi_k(t),\qquad x(t)=\Re\{z(t)\}.
$$
Each component is written as
$$
\psi_k(t)=A_k(t)\exp\!\bigl[j\,\theta_k(t)\bigr]=s_k(t)+j\,\sigma_k(t),
$$
where $A_k(t)$ is the instantaneous amplitude, $\theta_k(t)$ is the instantaneous phase, the quadrature $\sigma_k(t)$ is not fixed a priori, and the instantaneous frequency is assumed nonnegative:
$$
\omega_k(t)=\theta_k'(t)\ge 0.
$$
The instantaneous parameters are then defined exactly by
$$
A_k(t)=|\psi_k(t)|,\qquad
\theta_k(t)=\arg\{\psi_k(t)\},\qquad
\omega_k(t)=\frac{d}{dt}\theta_k(t).
$$
Under this construction, the Hilbert or instantaneous spectrum is
$$
H(\omega,t)=\sum_{k=0}^{K-1}A_k(t)\delta(\omega-\omega_k(t)).
$$
The central claim of this framework is that, with appropriate assumptions on the multicomponent AM–FM model, one obtains exact instantaneous amplitudes and instantaneous frequencies rather than a smeared time–frequency display [1504.07554].

A major conceptual feature of this theory is its treatment of uncertainty. Classical time–frequency methods such as STFT and wavelets are described as obeying a Heisenberg trade-off, whereas the AM–FM model shifts uncertainty into the model assumptions and the quadrature choice. In particular, the paper states that if one does not assume harmonic correspondence, there exist infinitely many choices of quadrature $\sigma(t)\neq\mathcal{H}\{x(t)\}$ for which the complex extension remains analytic in the complex-time sense; correspondingly, no universal linear operator can recover the “true” quadrature for all possible latent signals. Exactness is therefore conditional on adopting an appropriate component model rather than on applying a universal transform [1504.07554].

The practical computation described in this literature is the HSA–IMF algorithm. It proceeds by decomposing $x(t)$ into Intrinsic Mode Functions via Empirical Mode Decomposition, enforcing carried assumptions to obtain a unique amplitude estimate for each IMF, removing the amplitude to obtain a real frequency-modulated signal, normalizing to unit amplitude and recovering a quadrature, differentiating the recovered phase to obtain $\hat\omega_k(t)$, and finally assembling
$$
H(\omega,t)=\sum_k \hat A_k(t)\,\delta(\omega-\hat\omega_k(t)).
$$
The paper characterizes this as returning the exact instantaneous spectrum under the AM–FM component assumptions [1504.07554].

## 3. Transform-based, sparse, and nonperiodic computational constructions

Synchrosqueezing provides a rigorous transform-based route to IS. Starting from the continuous wavelet transform
$$
W_f(a,b)=a^{-1/2}\int_{\mathbb{R}} f(u)\,\overline{\psi((u-b)/a)}\,du,
$$
one defines a local frequency estimate wherever $W_f(a,b)\neq 0$:
$$
\omega(a,b)=\frac{1}{2\pi}\,\Im\!\left\{\frac{\partial_b W_f(a,b)}{W_f(a,b)}\right\}
= -\,\frac{i\,\partial_b W_f(a,b)}{2\pi\,W_f(a,b)}.
$$
The synchrosqueezing transform then reallocates each wavelet coefficient in the frequency direction toward $\omega=\omega(a,b)$, yielding
$$
T_f(\omega,b)=\lim_{\delta\to 0}\int_{a\in[M^{-1},M],\,|W_f(a,b)|>\epsilon}
W_f(a,b)\,\delta^{-1}h\!\left(\frac{\omega-\omega(a,b)}{\delta}\right)a^{-3/2}\,da.
$$
This IS is sparse, invertible, and stable: ridge-tracking in $|T_f(\cdot,b)|$ permits recovery of the constituent modes, and the paper reports robustness under small deterministic perturbations and Gaussian white noise [1404.7550].

Sparse Frequency Analysis (SFA) approaches IS from a convex variational standpoint. A discrete-time real signal is modeled as
$$
x[n]=\sum_{k=0}^{K-1}\bigl(a_k[n]\,c(n,k)+b_k[n]\,s(n,k)\bigr),
$$
with narrow-band sinusoids whose amplitudes vary in time. The instantaneous parameters are
$$
A_k[n]=\sqrt{a_k[n]^2+b_k[n]^2},\qquad
\phi_k[n]=\mathrm{atan2}(b_k[n],a_k[n]),\qquad
f_k[n]=\frac{1}{2\pi}\Delta_n\phi_k[n].
$$
The amplitudes are estimated by minimizing total variation in $a_k$ and $b_k$ together with a spectral sparsity term, either under perfect reconstruction or under a quadratic data-fidelity penalty in noise. The optimization is solved by ADMM with variable splitting and a nested MM step for the coupled $\ell_2$ sparsity term. Once the amplitudes are estimated, the instantaneous spectrum is defined on the discrete grid by
$$
\mathrm{IS}(f,n)=\sum_{k=0}^{K-1}A_k[n]\delta\bigl(f-f_k[n]\bigr)\approx \sum_{k=0}^{K-1}A_k[n]\delta(f-f_k),
$$
which preserves abrupt amplitude and phase changes more sharply than classical linear filtering [1302.6523].

A different line of work replaces the periodic boundary condition (PBC) underlying STFT implementations with the Linear eXtrapolation Condition (LXC). In that formulation, a real signal is locally modeled as a sum of complex AM–FM modes
$$
S(t)\approx \sum_{m=1}^M c_m(t_k)\,e^{[2\pi i f_m(t_k)+\lambda_m(t_k)](t-t_k)}
$$
near each analysis time $t_k$, where the poles
$$
h_m=H_m'(t_k)=2\pi i f_m(t_k)+\lambda_m(t_k)
$$
are extracted via a small Toeplitz system under the no-wrap-around LXC constraint. Amplitudes are then obtained by least squares, and each mode contributes a Lorentzian-shaped local spectrum
$$
F_m(f,t_k)=\left|\frac{c_m\,\lambda_m}{\lambda_m+2\pi i(|f_m|-f)}\right|.
$$
Summing over modes yields a framewise instantaneous spectrum suitable for pulse-series signals [2602.03680].

The quadratic-chirplet framework of instantaneous time-frequency atoms places several of these constructions into a single two-parameter continuum. With
$$
\psi_{\gamma,\alpha,\beta}(s;\tau)
=(\gamma/\pi)^{1/4}\exp\!\left[-\frac{\gamma}{2}(s-\tau)^2+j2\pi\alpha(s-\tau)+j\pi\beta(s-\tau)^2\right],
$$
the IS is $|V_x(\tau,\alpha=f,\beta)|^2$, and time-domain analysis, frequency-domain analysis, fractional Fourier analysis, synchrosqueezed STFT, and synchrosqueezed STFrFT arise as limiting or specialized cases in the $(\beta,\gamma)$ plane. This framework explicitly treats those analyses as decompositions into AM–FM components using specialized or limiting forms of a quadratic chirplet [2508.05380].

## 4. Multivariate instantaneous spectrum and physical moments

In multivariate analysis, IS is generalized from scalar amplitude–frequency traces to matrix-valued instantaneous moments. For a two- or three-component real-valued time series $x(t)\in\mathbb{R}^2$ or $\mathbb{R}^3$, the analytic signal is
$$
x_+(t)=x(t)+i\,H[x](t).
$$
The instantaneous amplitude and instantaneous frequency are defined by
$$
A(t)\equiv \kappa(t)=\sqrt{\tfrac12\,x_+(t)\cdot x_+(t)},
\qquad
\omega(t)=\Im\{x_+(t)\cdot x_+'(t)\}/\|x_+(t)\|^2.
$$
One then introduces the one-sided spectral matrix and its global moments, together with instantaneous moment matrices
$$
S_0(t)=x_+(t)x_+(t)^{\mathrm{H}},\qquad
S_1(t)=\mathrm{Herm}\{(-i)x_+'(t)x_+(t)^{\mathrm{H}}\},\qquad
S_2(t)=x_+'(t)x_+'(t)^{\mathrm{H}}.
$$
Their traces recover instantaneous power, instantaneous frequency, and bandwidth [1201.5916].

The distinctive result in this line of work is that these instantaneous spectral moments coincide exactly with the physical moments of a canonical time-varying ellipse traced by the signal. The ellipse is parameterized by semi-axes, orbital phase, and, in three dimensions, a rotation. Scalar, vector, and tensor moments such as variance, normal vector, moment of inertia, circulation, and angular momentum are defined for the corresponding ring of particles. The paper proves exact identities linking the matrix-valued spectral moments to those physical moments, and, taking the trace of $S_1(t)/2$, obtains the formula
$$
\Gamma(t)=A^2(t)\,\omega(t).
$$
Thus circulation equals instantaneous power times instantaneous frequency [1201.5916].

This correspondence gives IS a geometric and mechanical interpretation. The amplitude becomes the root-mean-square radius of the ellipse, while the frequency decomposes into orbital and precessional contributions. In this setting, the instantaneous spectrum is not merely a sharper spectrogram; it is a local moment decomposition of an evolving geometric object. The paper further notes that multivariate wavelet-ridge analysis can be used in noisy data to estimate $x_+(t)$ and hence the instantaneous moment structure [1201.5916].

## 5. Physical realizations in RF and microwave spectrum analyzers

The term instantaneous spectrum also appears

Source: https://www.emergentmind.com/topics/instantaneous-spectrum-is