Skip to main content

β: This English translation is in beta — the Traditional-Chinese original is the authoritative version.

Stochastic Noise Basics

Prerequisites: lti_vs_ltv · Notation | Next: white_noise_to_phase_noise

This page answers a prerequisite question: before we can talk about the ISF, we have to pin down "noise" as a stochastic process. Phase jitter in an oscillator originates from the random currents in(t)i_n(t) of transistors and resistors. Before feeding it into the phase integral of [P1] Eq.(11), we need to know: how to quantify the "strength" of this random current (PSD), how to recover time-domain power from the PSD (Parseval), whether it is correlated across time instants (autocorrelation), whether a single measurement can represent the whole ensemble (ergodicity), and — most critical for oscillators — whether the noise intensity itself varies periodically with the oscillation phase (cyclostationary).

That last point is the core of effective_isf later on: the transistor leaks noise only during the short window when it conducts, and where that window falls in the waveform's phase determines which segment of the ISF weights it.

Physical intuition (big picture first): think of noise as "an invisible hand randomly kicking the oscillator at every instant." We need two independent pieces of information: (1) how hard the hand kicks on average (intensity, measured by the PSD); (2) when the hand kicks hardest (time modulation, described by cyclostationarity / the NMF). The ISF is the third item: "at which phase of the waveform a kick changes the phase the most." Multiply the three together and integrate — that is phase noise.

1. White noise

White noise (a stochastic process whose power spectrum is flat at all frequencies). The name comes from white light — all frequency components with equal strength. Its single-sided PSD (power spectral density) is constant:

Si(f)=in2Δf=const.S_i(f)=\frac{\overline{i_n^2}}{\Delta f}=\text{const}.
  • Units: current-noise PSD is in A2/Hz\text{A}^2/\text{Hz} (mean-square current per Hz of bandwidth). Typical thermal-noise expressions: in2/Δf=4kT/R\overline{i_n^2}/\Delta f=4kT/R (resistor) or 4kTγgm4kT\gamma g_m (MOS channel).
  • Dimension check: in2\overline{i_n^2} is A2\text{A}^2; dividing by the bandwidth Δf\Delta f (Hz) gives A2/Hz\text{A}^2/\text{Hz} ✓.
  • Physical origin: thermal noise (random thermal motion of carriers) and shot noise (discreteness of carriers crossing a barrier) can both be treated as white over the frequency range we care about.
  • Canonical value: this site's white-noise source uses Si=1024 A2/HzS_i=10^{-24}\ \text{A}^2/\text{Hz} (see canonical example B). This corresponds to a current-noise density of Si=1012 A/Hz=1 pA/Hz\sqrt{S_i}=10^{-12}\ \text{A}/\sqrt{\text{Hz}}=1\ \text{pA}/\sqrt{\text{Hz}}.

Why oscillators care about white noise: white noise has equal-strength components near every nω0n\omega_0; each ISF harmonic cnc_n "downconverts" the white noise of the corresponding band to near the carrier, forming the 1/f² (20-20 dB/dec) phase-noise skirt. This is exactly the content of [P1] Eq.(21) — see white_noise_to_phase_noise.

2. Flicker (1/f) noise

Flicker noise (also called 1/f noise: low-frequency noise whose power rises inversely with frequency). In MOS devices it mainly comes from carriers being randomly captured/released by interface traps. Its PSD is not flat:

Si,1/f(f)=in2ω1/fΔω(Δω<ω1/f),S_{i,1/f}(f)=\overline{i_n^2}\cdot\frac{\omega_{1/f}}{\Delta\omega}\qquad(\Delta\omega < \omega_{1/f}),

which is exactly the device flicker model of [P1] Eq.(22), p.185 (site spec Section 3, formula 13).

  • Units: still A2/Hz\text{A}^2/\text{Hz}. ω1/f\omega_{1/f} is the device's 1/f corner (angular frequency, rad/s) — the frequency at which the 1/f component's strength equals the white-noise floor.
  • Dimension check: in2\overline{i_n^2} (A2/Hz\text{A}^2/\text{Hz}) × ω1/f/Δω\times\ \omega_{1/f}/\Delta\omega (dimensionless ratio) =A2/Hz=\text{A}^2/\text{Hz} ✓.
  • Key pitfall: the device 1/f corner ω1/f\omega_{1/f} is not equal to the phase-noise 1/f³ corner. The latter is Δω1/f3=ω1/fc02/(2Γrms2)\Delta\omega_{1/f^3}=\omega_{1/f}\cdot c_0^2/(2\Gamma_{rms}^2) ([P1] Eq.(24)), upconverted only through the ISF's DC coefficient c0c_0. A symmetric waveform (small c0c_0) can push the 1/f³ corner far below the device corner — see symmetry (claim C5).
Noise typePSD shapePhysical originBecomes, in phase noise
whiteflat (const A2/Hz\text{A}^2/\text{Hz})thermal / shot1/f² (20-20 dB/dec)
flicker (1/f)1/f\propto 1/ftrap capture/release1/f³ (30-30 dB/dec), via c0c_0 only

3. PSD and Parseval: frequency-domain intensity ↔ time-domain power

The PSD (power spectral density) tells you "how much power per unit bandwidth." Integrating it over frequency gives the total time-domain power (variance). This is Parseval's theorem (time-domain energy = frequency-domain energy) in its stochastic-process form:

in2=0Si(f)df(單邊 PSD).\overline{i_n^2}=\int_0^{\infty}S_i(f)\,df\qquad(\text{單邊 PSD}).
  • Units: left side A2\text{A}^2; right side (A2/Hz)×Hz=A2(\text{A}^2/\text{Hz})\times\text{Hz}=\text{A}^2 ✓.
  • Intuition: the PSD is a "density of power," and integrating means "multiply density by bandwidth and add it up." For white noise this integral diverges — so any real system has a finite bandwidth, and white noise is only an in-band approximation.
  • Why it matters: phase noise plays the same trick — integrate the phase PSD Sϕ(f)S_\phi(f) over offset frequency to get the phase variance σϕ2\sigma_\phi^2, then convert to rms jitter. Full derivation in psd_phase_noise_jitter.

Convention: this site uses the single-sided PSD (positive frequencies only, 0f<0\le f<\infty). If you switch to the double-sided (<f<-\infty<f<\infty) PSD, the value at each frequency is divided by 2. This factor-of-2 convention is the same class of bookkeeping issue as the SSB phase-noise relation L12Sϕ\mathcal{L}\approx\frac12 S_\phi, discussed in detail in white_noise_to_phase_noise.

4. Autocorrelation and the real definition of "white"

Autocorrelation (how correlated one stochastic process is between two time instants):

Ri(τ)=in(t)in(t+τ).R_i(\tau)=\overline{i_n(t)\,i_n(t+\tau)}.
  • Units: A2\text{A}^2 (the expectation of a product of two currents).
  • Wiener–Khinchin theorem: autocorrelation and PSD are a Fourier pair. White noise has a flat PSD, corresponding to a delta-function autocorrelation: Ri(τ)=in2Δfδ(τ)R_i(\tau)=\dfrac{\overline{i_n^2}}{\Delta f}\,\delta(\tau).
  • Intuition: white noise's "value now" is completely uncorrelated with its "value the next instant" (zero memory). This is exactly why, in the phase integral of [P1] Eq.(11), noise contributions at different instants can be superposed independently — it turns the phase-variance calculation into "square and add term by term" with no cross-correlation terms to handle.
  • Flicker counterexample: 1/f noise has a long autocorrelation tail (long memory, strong correlation), which is also why it is hard to integrate and hard to simulate (see the aliasing caveats on the DSP page).

5. Ergodicity

Ergodicity (the property that "time average = ensemble average").

  • Plain language: you have only one oscillator and measure one long waveform (time average); the theory speaks of the statistics of infinitely many identical oscillators at the same instant (ensemble average). Ergodicity guarantees the two are equal — otherwise "measuring one device" could never validate "the ensemble statistics predicted by theory."
  • Applicability: the process must be stationary (statistics do not drift with time) before ergodicity is even on the table. Thermal noise qualifies.
  • Failure warning: the oscillator's excess phase ϕ(t)\phi(t) itself is not stationary — it is a random walk (variance grows linearly with time; see accumulated jitter σΔt=κΔt\sigma_{\Delta t}=\kappa\sqrt{\Delta t}, [P2] Eq.(8)). So we do statistics on "the time derivative of phase (frequency)" or on "phase differences," not on the absolute value of ϕ\phi. This subtlety is retold in dsp_view_of_phase_noise from the "1/(jω) integrator turns white noise into 1/f²" viewpoint.

6. Cyclostationary noise — this page's focus, straight into the effective ISF

Cyclostationary noise (a stochastic process whose statistics vary "periodically" with time). [P1] Sec. III-D, p.186 says it plainly:

"the channel noise of a MOS device in an oscillator is cyclostationary because the noise power is modulated by the gate-source overdrive which varies with time periodically." ([P1], p.186, cyclostationary paragraph)

In other words: the transistor does not leak the same amount of noise at every instant. It is noisy only during the short window when it conducts (large current, large overdrive); the rest of the time it is nearly silent. And "conducting" is locked to the oscillation waveform — repeating once per period TT. That is "cyclostationarity."

[P1] decomposes a white cyclostationary current as ([P1] Eq.(25), p.186):

in(t)=in0(t)α(ω0t),i_n(t)=i_{n0}(t)\cdot\alpha(\omega_0 t),

where in0(t)i_{n0}(t) is a white stationary process (fixed intensity) and α(ω0t)\alpha(\omega_0 t) is a deterministic, periodic, normalized (max = 1) function describing "when the noise amplitude gets turned up or down." This site calls it the NMF — noise-modulating function, symbol α(ω0t)\alpha(\omega_0 t), 0α10\le\alpha\le1, dimensionless (see the notation page).

  • Units: α\alpha is dimensionless (it is a normalized envelope); in0i_{n0} and ini_n are both in A.
  • Dimension check: A=A×(無因次)\text{A}=\text{A}\times(\text{無因次}) ✓.
  • Meaning of α2\alpha^2: because the PSD deals in power (i2i^2), the noise's instantaneous mean-square power is proportional to α2(ω0t)\alpha^2(\omega_0 t). The maximum instantaneous PSD =in2/Δf=\overline{i_n^2}/\Delta f; at all other times it is suppressed by α21\alpha^2\le1.

Why this leads straight to the ISF: the effective ISF

Substitute the cyclostationary decomposition into the phase integral [P1] Eq.(11): the noise contribution at instant τ\tau is weighted by Γ(ω0τ)\Gamma(\omega_0\tau), while its own intensity is modulated by α(ω0τ)\alpha(\omega_0\tau). The product of the two periodic functions can be packaged into a single effective ISF ([P1] Eq.(27), p.186):

 Γeff(ω0τ)=Γ(ω0τ)α(ω0τ) \boxed{\ \Gamma_{eff}(\omega_0\tau)=\Gamma(\omega_0\tau)\cdot\alpha(\omega_0\tau)\ }

In [P1]'s own words:

"the cyclostationary noise can be treated as a stationary noise applied to a system with an effective ISF." ([P1], p.186)

Meaning: just replace Γ\Gamma with Γeff=Γα\Gamma_{eff}=\Gamma\cdot\alpha, and every downstream formula (cnc_n, Γrms\Gamma_{rms}, Eq.(21), Eq.(23)) carries over unchanged. Cyclostationarity does not need a new theory; it only changes which curve you take the Fourier coefficients of.

  • Units: Γeff\Gamma_{eff} is still dimensionless (dimensionless Γ\Gamma × dimensionless α\alpha) ✓.
  • Teaching point: when computing cnc_n and Γrms\Gamma_{rms}, use Γeff\Gamma_{eff} (i.e., as [P1] puts it, "the effective ISF should be used in all subsequent calculations"), not the bare Γ\Gamma.

LC vs ring: the cyclostationary impact differs enormously ([P1]'s beautiful conclusion)

[P1], p.186 gives two contrasting examples of great pedagogical value:

OscillatorLarge current (large α\alpha) occurs atΓ\Gamma thereCyclostationary impact
Colpitts LCtank-voltage minimum (trough)small (ISF is small at peaks/troughs)Γeff\Gamma_{eff} differs a lot from Γ\Gamma; cannot be ignored
ring (single-ended)transitions (zero crossings, max slope)large (ISF is large at transitions)ΓeffΓ\Gamma_{eff}\approx\Gamma; small impact

[P1]'s reading: the Colpitts transistor "turns on only where the ISF is small" — nature effectively injects its noise at insensitive phases, so cyclostationarity helps the LC. The ring is unlucky: the maximum current occurs exactly at the most sensitive transitions, so noise and sensitivity peak simultaneously. [P1] calls it out directly:

"This unfortunate coincidence is one of the reasons why ring oscillators in general have inferior phase noise performance compared to a Colpitts LC oscillator." ([P1], p.186)

This is one of the two big reasons "LC is cleaner than ring" (the other: a ring burns all of its stored energy every period). Quantitative discussion in [P2] and lc_vs_ring.

Toy-model reminder: later pages use toy shapes for α(ω0t)\alpha(\omega_0 t) (e.g., setting the NMF to a conduction window) to demonstrate how Γeff=Γα\Gamma_{eff}=\Gamma\cdot\alpha changes c0c_0 and Γrms\Gamma_{rms}; these are pedagogical toy models, not transistor-level. The real α\alpha must be extracted from the device noise model and operating point ([P1] states "α\alpha can be derived easily from device noise characteristics and operating point").

Connecting this page to the ISF in one sentence

The full phase-noise recipe is the product of three periodic/stochastic quantities, then integrated:

The two "×" signs in the middle are the origin of the effective ISF Γeff=Γα\Gamma_{eff}=\Gamma\cdot\alpha; the final integrator is what turns white noise into 1/f² (see the DSP viewpoint linked in the next section).

Applicability and failure conditions

ConditionWhen it holdsWhen it fails
stationary (thermal noise)ergodicity holds; measuring one device represents the ensembledrift/aging → re-measure
white approximationflat PSD, delta autocorrelation, contributions superpose independentlyinside the 1/f region, switch to the flicker model
cyclostationarity packaged as Γeff=Γα\Gamma_{eff}=\Gamma\cdot\alphasmall perturbation; α\alpha periodic and normalizedstrong nonlinearity coupling α\alpha to the signal calls for more care
phase described as a random walkdo statistics on phase differences/frequencytime-averaging the absolute value of ϕ\phi directly diverges

Key takeaways

  • White noise: flat PSD (A2/Hz\text{A}^2/\text{Hz}), delta autocorrelation, memoryless; becomes 1/f² under the ISF.
  • Flicker (1/f) noise: S1/fS\propto1/f ([P1] Eq.(22)); upconverts to 1/f³ only via c0c_0; device corner \neq 1/f³ corner.
  • PSD ↔ Parseval: integrating the PSD gives time-domain power (variance); units self-consistent A2/Hz×Hz=A2\text{A}^2/\text{Hz}\times\text{Hz}=\text{A}^2.
  • Autocorrelation: white = delta (memoryless); 1/f = long tail (strongly correlated).
  • Ergodicity: time average = ensemble average; but ϕ\phi is a random walk, non-stationary — do statistics on phase differences/frequency.
  • Cyclostationary (this page's protagonist): noise intensity varies periodically with the oscillation phase (the transistor leaks noise only while conducting); described by the NMF α(ω0t)\alpha(\omega_0 t) and packaged as Γeff=Γα\Gamma_{eff}=\Gamma\cdot\alpha ([P1] Eq.(25)–(27)); LC benefits, ring suffers.
  • Sources: [P1] Sec. III-D (cyclostationary) Eqs.(25)–(27), p.186; device flicker Eq.(22), p.185; claims C9.

Further reading