Bartlett's Method
Understanding Bartlett's Method
The periodogram (|FFT|²/N) is an unbiased but inconsistent estimator of the power spectral density. "Inconsistent" means that its variance does not decrease as more data is collected. Longer records produce more spectral bins but each bin still fluctuates by approximately ±100% of its true value. Bartlett's insight was that averaging independent periodograms reduces this variance.
The trade-off is fundamental: shorter segments mean fewer frequency bins (coarser resolution) but averaging K of them reduces variance by K. With N = 10,000 samples and K = 10, each segment has M = 1,000 samples. Frequency resolution degrades by 10x, but the spectral estimate is 10x smoother (10 dB lower variance). This is the classical bias-variance trade-off in spectral estimation.
Bartlett's Method Formulas
1. Divide x[0..N−1] into K segments of M = N/K
2. Pk(f) = (1/M)|∑ xk[n] e−j2πfn/M|²
3. PBartlett(f) = (1/K) ∑ Pk(f)
Variance Reduction:
Var[PBartlett] = Var[Pperiodogram] / K
Resolution vs. Variance:
Δf = K/N Hz (K× worse than periodogram)
Variance: 1/K of periodogram
Resolution × Variance = constant (for fixed N)
Spectral Estimation Methods
| Method | Overlap | Window | Variance | Leakage |
|---|---|---|---|---|
| Periodogram | None | Rectangular | High | High |
| Bartlett | 0% | Rectangular | σ²/K | High |
| Welch | 50% | Hann/Hamming | ~σ²/K | Low |
| Blackman-Tukey | N/A | Correlation | Low | Tunable |
Frequently Asked Questions
How does Bartlett's method work?
Divide N samples into K segments of M=N/K. Periodogram each segment. Average K periodograms. Variance reduced by K. Resolution degraded by K. Fundamental variance-resolution trade-off for fixed N.
Bartlett vs. Welch?
Welch: overlapping + windowed segments. 50% overlap: ~82% of Bartlett's variance reduction. Windows reduce spectral leakage. Welch universally preferred in practice. Bartlett is the theoretical foundation.
RF applications?
Spectrum analyzer trace averaging (video average). Noise PSD measurement. Phase noise estimation. Channel power. Interference characterization. Digital receivers: hardware-accelerated FFT averaging for real-time spectral monitoring.