Skip to content

Hilbert Transform

The Hilbert transform constructs a quadrature companion to a signal. Its principal uses are analytic-signal formation, envelope and phase analysis, and single-sideband (SSB) generation. The transform is a singular convolution in time and a sign-dependent phase rotation in frequency. This section uses the Fourier convention already established in Section 4, X(ω)=∫−∞∞x(t)e−jωtdtX(\omega)=\int_{-\infty}^{\infty}x(t)e^{-\mathrm{j}\omega t}\mathrm{d}t, without repeating the broader Fourier theory.

The Hilbert transform of x(t)x(t) is

x^(t)=H{x(t)}=1πPV⁡∫−∞∞x(τ)t−τ dτ.\hat{x}(t)=\mathcal H\{x(t)\} =\frac{1}{\pi}\operatorname{PV} \int_{-\infty}^{\infty}\frac{x(\tau)}{t-\tau}\,\mathrm{d}\tau.

Here PV⁡\operatorname{PV} denotes the Cauchy principal value; the ordinary integral is generally undefined because its kernel is singular at τ=t\tau=t.

For a sufficiently well-behaved signal, the principal value means the symmetric limiting process

PV⁡∫−∞∞x(τ)t−τ dτ=lim⁡R→∞lim⁡ϵ↓0[∫t−Rt−ϵx(τ)t−τ dτ+∫t+ϵt+Rx(τ)t−τ dτ].\operatorname{PV}\int_{-\infty}^{\infty} \frac{x(\tau)}{t-\tau}\,\mathrm{d}\tau =\lim_{R\to\infty}\lim_{\epsilon\downarrow0} \left[ \int_{t-R}^{t-\epsilon}\frac{x(\tau)}{t-\tau}\,\mathrm{d}\tau +\int_{t+\epsilon}^{t+R}\frac{x(\tau)}{t-\tau}\,\mathrm{d}\tau \right].

The equal exclusion on the two sides allows the singular contributions to cancel. Replacing the principal value by an ordinary improper integral is a fundamental error.

The same definition can be written as the convolution

The kernel is understood as a distribution. Pointwise existence follows under standard smoothness and decay conditions. More generally, H\mathcal H is bounded on Lp(R)L^p(\mathbb R) for 1<p<∞1<p<\infty; in particular, every finite-energy signal in L2L^2 has a Hilbert transform in the mean-square sense even when the pointwise integral requires care.

If x(t)↔X(ω)x(t)\leftrightarrow X(\omega) and x^(t)↔X^(ω)\hat{x}(t)\leftrightarrow \hat{X}(\omega), then

The adopted DC convention is sgn⁡(0)=0\operatorname{sgn}(0)=0, so

HH(ω)={e−jπ/2=−j,ω>0,0,ω=0,e+jπ/2=+j,ω<0.H_{\mathrm H}(\omega)= \begin{cases} e^{-\mathrm{j}\pi/2}=-\mathrm{j}, & \omega>0,\\ 0, & \omega=0,\\ e^{+\mathrm{j}\pi/2}=+\mathrm{j}, & \omega<0. \end{cases}

At DC the response has zero magnitude and no phase, so a pure constant is mapped to zero. For an ordinary L2L^2 spectrum, changing the response at the single point ω=0\omega=0 does not alter the time signal. The convention does matter for a constant or for the discrete DC line of a periodic signal.

For a nonzero-frequency complex exponential,

Consequently, for ω0>0\omega_0>0 and arbitrary phase ϕ\phi,

These signs should be derived rather than guessed: split a sinusoid into its positive- and negative-frequency exponentials and apply the two multipliers separately.

Applying the transform twice gives [−jsgn⁡(ω)]2=−1[-\mathrm{j}\operatorname{sgn}(\omega)]^2=-1 away from DC. If x(t)=xDC+xAC(t)x(t)=x_{\mathrm{DC}}+x_{\mathrm{AC}}(t), then

For a periodic signal, xDCx_{\mathrm{DC}} is its mean. For an L2L^2 signal there is no persistent constant component, so H2{x}=−x\mathcal H^2\{x\}=-x in the mean-square sense and H−1=−H\mathcal H^{-1}=-\mathcal H.

Energy, Orthogonality, and Operator Properties

Section titled “Energy, Orthogonality, and Operator Properties”

For any finite-energy signal, the unit magnitude of the multiplier almost everywhere gives

Ex^=12π∫−∞∞∣X^(ω)∣2 dω=12π∫−∞∞∣X(ω)∣2 dω=Ex.E_{\hat{x}} =\frac{1}{2\pi}\int_{-\infty}^{\infty} |\hat{X}(\omega)|^2\,\mathrm{d}\omega =\frac{1}{2\pi}\int_{-\infty}^{\infty} |X(\omega)|^2\,\mathrm{d}\omega =E_x.

If x(t)x(t) is also real, ∣X(ω)∣2|X(\omega)|^2 is even, and

∫−∞∞x(t)x^(t) dt=−j2π∫−∞∞sgn⁡(ω)∣X(ω)∣2 dω=0.\int_{-\infty}^{\infty}x(t)\hat{x}(t)\,\mathrm{d}t =-\frac{\mathrm{j}}{2\pi}\int_{-\infty}^{\infty} \operatorname{sgn}(\omega)|X(\omega)|^2\,\mathrm{d}\omega=0.

Thus a real finite-energy signal and its Hilbert transform are orthogonal in the global inner-product sense.

The following identities hold whenever the displayed operations and transforms exist; aa in the scaling row is real and nonzero.

PropertyResult
LinearityH{a1x1+a2x2}=a1x^1+a2x^2\mathcal H\{a_1x_1+a_2x_2\}=a_1\hat{x}_1+a_2\hat{x}_2
Time shiftH{x(t−t0)}=x^(t−t0)\mathcal H\{x(t-t_0)\}=\hat{x}(t-t_0)
Time scalingH{x(at)}=sgn⁡(a)x^(at)\mathcal H\{x(at)\}=\operatorname{sgn}(a)\hat{x}(at)
DifferentiationH{dxdt}=ddtH{x}\mathcal H\{\frac{\mathrm{d}x}{\mathrm{d}t}\}=\frac{\mathrm{d}}{\mathrm{d}t}\mathcal H\{x\}
ConvolutionH{x∗y}=x^∗y=x∗y^\mathcal H\{x*y\}=\hat{x}*y=x*\hat{y}
Parity for real signalsEven xx gives odd x^\hat{x}; odd xx gives even x^\hat{x}
Inverse on the AC subspaceH−1=−H\mathcal H^{-1}=-\mathcal H

For a real signal x(t)x(t), its analytic signal is

xa(t)=x(t)+jx^(t).x_a(t)=x(t)+\mathrm{j}\hat{x}(t).

Its real part is the original signal and its imaginary part is the Hilbert transform.

Because X^(ω)=−jsgn⁡(ω)X(ω)\hat{X}(\omega)=-\mathrm{j}\operatorname{sgn}(\omega)X(\omega),

Thus negative frequencies are suppressed, positive frequencies are doubled, and DC is retained once. These factors ensure

x(t)=ℜ{xa(t)}=xa(t)+xa∗(t)2.x(t)=\Re\{x_a(t)\} =\frac{x_a(t)+x_a^*(t)}{2}.

The shorthand Xa(ω)=2u(ω)X(ω)X_a(\omega)=2u(\omega)X(\omega) is correct only with the symmetric convention u(0)=1/2u(0)=1/2; otherwise it doubles or deletes DC incorrectly.

For a real finite-energy signal, the positive and negative frequency halves contain equal energy. Doubling one half gives

Exa=2Ex.E_{x_a}=2E_x.

For a real periodic signal with mean xˉ\bar{x}, the corresponding power result is Pxa=2Px−∣xˉ∣2P_{x_a}=2P_x-|\bar{x}|^2, because its DC line is not doubled.

Envelope, Instantaneous Phase, and Frequency

Section titled “Envelope, Instantaneous Phase, and Frequency”

Where xa(t)≠0x_a(t)\neq0, write the analytic signal in polar form:

The phase must be unwrapped before differentiation. These definitions also require A(t)>0A(t)>0 and sufficient differentiability. More importantly, they have a clear physical interpretation mainly for a monocomponent or suitably narrowband AM–FM signal. For a multicomponent waveform, A(t)A(t) can show beating and ωi(t)\omega_i(t) can become erratic or even negative; computability alone does not guarantee a useful physical interpretation.

Let m(t)m(t) be a real message bandlimited to ∣ω∣≤Ωm|\omega|\leq\Omega_m, let ma(t)=m(t)+jm^(t)m_a(t)=m(t)+\mathrm{j}\hat{m}(t), and choose ωc>Ωm\omega_c>\Omega_m so the translated spectra do not overlap. The direct and Hilbert-transform paths modulate quadrature carriers. Their sum or difference produces

The two product-modulator outputs are I(t)=m(t)cos⁡ωctI(t)=m(t)\cos\omega_ct and Q(t)=m^(t)sin⁡ωctQ(t)=\hat{m}(t)\sin\omega_ct. With the convention used here, I−QI-Q selects USB and I+QI+Q selects LSB. the figure also retains the referenced filter-method alternative for comparison.

SSB generation methods. (a) A balanced modulator generates DSB-SC and a sharp BPF selects one sideband. (b) The message is split into direct and Hilbert-transform paths, quadrature product modulators form the two components, and subtraction selects USB while addition selects LSB.

SSB generation methods. (a) A balanced modulator generates DSB-SC and a sharp BPF selects one sideband. (b) The message is split into direct and Hilbert-transform paths, quadrature product modulators form the two components, and subtraction selects USB while addition selects LSB.

Ideal sideband cancellation requires equal path gain and exact quadrature. If the quadrature path has relative gain gg and phase error ϵ\epsilon, a single-tone model gives desired and residual sideband phasors proportional to 1+gejϵ1+ge^{\mathrm{j}\epsilon} and 1−gejϵ1-ge^{\mathrm{j}\epsilon}. Their power ratio is

Perfect gain and phase balance gives infinite ideal rejection. For g=1g=1 and small ∣ϵ∣|\epsilon| in radians, IRR≈4/ϵ2\mathrm{IRR}\approx4/\epsilon^2; even a small phase error therefore leaves a measurable unwanted sideband. Gain ripple and unequal delay have the same practical consequence.

Ideal Phase Shifter and Practical Realization

Section titled “Ideal Phase Shifter and Practical Realization”

The Hilbert transformer is often called an ideal 90∘90^\circ phase shifter, but this phrase needs two qualifications:

  • positive and negative frequencies receive opposite phase shifts, as required for a real quadrature output;

  • the phase shift is frequency independent, but it is not a fixed time delay, whose phase is −ωtd-\omega t_d.

Its impulse response hH(t)=1/(πt)h_{\mathrm H}(t)=1/(\pi t) is nonzero for t<0t<0 and extends infinitely in both directions. It is singular at the origin and not absolutely integrable. Therefore an exact Hilbert transformer is neither causal nor BIBO-stable in the ordinary LTI impulse-response sense, and it cannot be realized physically over all frequencies. This does not contradict the equation: boundedness and energy preservation on L2L^2 are weaker conditions than BIBO stability.

Analog all-pass networks and digital FIR/IIR Hilbert transformers approximate the required phase over a finite band. A causal digital implementation adds delay, so the direct path must receive the same delay to remain aligned with the quadrature path. Finite transition bands, gain ripple, phase error, and delay mismatch limit envelope accuracy and SSB suppression. The discontinuity at DC is another reason no finite causal network can implement the ideal response from zero to infinite frequency.