Matched filtering and waveform systematics
Define a noise-weighted waveform comparison with a one-sided PSD and distinguish signal recovery, parameter bias and model validation.
On this page
The question for today
What you’ll understand
- Compute a discretized optimal signal-to-noise ratio with consistent Fourier and PSD units.
- Explain the Gaussian-noise likelihood behind matched filtering.
- Quantify a simple overlap loss while recognizing waveform and parameter degeneracies.
A weak gravitational-wave signal is compared with data through a model of both the waveform and detector noise. Matched filtering weights frequencies according to their noise power rather than treating every visible oscillation equally. Its optimality statement has explicit assumptions about noise and signal shape. Recovering a large filter value, estimating source parameters, and demonstrating the accuracy of a waveform model are related tasks with different standards. A transparent finite example can illustrate the normalization without pretending to reproduce a real search or inference pipeline.
The noise spectrum defines the inner product
Write calibrated detector strain as signal plus noise. For an ideal stationary zero-mean Gaussian noise model, its one-sided power spectral density defines the frequency weighting. Adopt a Fourier transform with a negative two-pi-i-f-t exponent. Real time series then have redundant negative-frequency information, so an integral over positive frequencies uses the conventional factor four in the real noise-weighted inner product. A two-sided spectrum would require a corresponding change of normalization.
Strain is dimensionless, its Fourier transform has seconds as units, and a one-sided strain-noise PSD has inverse-hertz units. Dividing two Fourier amplitudes by the PSD and integrating over hertz therefore gives a dimensionless inner product. This unit check helps expose mistakes involving amplitude spectral density, which is the square root of the PSD and cannot simply replace it in the denominator.
Likelihood and overlap answer different questions
Under the stated Gaussian model the log likelihood is minus one half the squared noise-weighted residual, up to a parameter-independent constant. Maximizing over an amplitude for a fixed template shape yields matched filtering. Searches also consider arrival time, phase and template parameters, and assess a statistic's background under the actual noise process. Nonstationary glitches or uncertain spectra can violate the simple assumptions, so an optimal formula is not itself a false-alarm probability.
Normalize waveforms by their noise-weighted norms to define an overlap, often maximized over relative time and phase. A best overlap below one reduces the signal-to-noise ratio recoverable with that template family. It measures similarity under the chosen PSD and band, not a universal pointwise waveform error. Optimizing more physical parameters can improve recovery while shifting their inferred values, so a high match does not guarantee negligible parameter bias.
Systematics compete with statistical information
Waveform errors can arise from truncated PN phasing, uncertain merger calibration, omitted precession or eccentricity, inaccurate higher modes and an inappropriate ringdown attachment. Detector calibration and PSD uncertainty affect the same comparison. A model should be validated against independent calculations over its intended parameter region. Agreement with its own fitted training cases does not establish accuracy everywhere, especially when extrapolating to mass ratios or spins not covered by those checks.
Chirp mass, mass ratio, spins, redshift, inclination, distance and antenna projection leave partially degenerate imprints. A stronger signal can narrow statistical errors while making a fixed modeling error more important. Real inference therefore compares calibrated detector data with propagated waveforms, priors and a noise model, and studies systematic robustness. The elementary filter geometry teaches how information is weighted; it cannot supply trustworthy source posteriors, discovery claims or a complete error budget from a toy signal alone.
Work it through
A two-bin optimal norm
In a deliberately finite positive-frequency approximation, two bins each have width 1 Hz. The waveform values are htilde1=2×10⁻²³ s and htilde2=1×10⁻²³ s, and both PSD values are 4×10⁻⁴⁶ Hz⁻¹. Compute ρopt² and ρopt. Then double both waveform amplitudes at fixed PSD.
Replace the integral by the declared two-bin sum. The squared Fourier magnitudes are 4×10⁻⁴⁶ s² and 1×10⁻⁴⁶ s²; adding them before division is valid because the two PSD values and bin widths are equal.
The norm is 4×(1 Hz)×[5×10⁻⁴⁶ s²/(4×10⁻⁴⁶ Hz⁻¹)]=5. Thus the optimal signal-to-noise ratio is √5≈2.236. Fourier seconds and the spectral units cancel to leave a dimensionless norm.
Doubling the waveform multiplies each squared magnitude by four, so the norm becomes 20 and ρopt becomes 2√5≈4.472. The change is linear in strain amplitude, not quadratic in signal-to-noise ratio itself.
The finite example gives ρopt²=5 and ρopt≈2.236; doubling strain doubles ρopt. These values illustrate a known-signal norm, not a real detection threshold or an analysis of measured detector noise.
The living notebook
Explore the model
Interactive lab · matrix
Noise weighting rescales data directions
View two real coordinates only as an analogy for two frequency components. Reset and set both controls to Diagonal scale a=1 and Shear b=0; predict selected determinant 1 and singular values 1,1. For a fresh comparison, Reset and set both controls to a=0.5 and b=0; predict determinant 0.5 and singular values 1,0.5. Geometry starts at identity after Reset and interpolates during playback; readouts describe the selected final matrix.
- Area scale · det A
- 1.4
- Singular values · σmax, σmin
- 1.597, 0.876
- Condition number · κ₂(A)
- 1.823
The selected matrix is A = [[a,b],[0,1]]. Playback interpolates geometry between the identity and A; readouts always describe A. Its determinant is a and its singular values come from AᵀA. Positive a makes every supported matrix invertible. This two-dimensional family is not a general eigensolver or a learned network. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forA relative whitening weight of one half corresponds algebraically to four times the noise power in that coordinate. This positive triangular matrix lab does not estimate a PSD, generate Gaussian noise, whiten actual Fourier data or compute a likelihood.
Pause & reason
Check your understanding
Your turn
Try it yourself
A source has optimal SNR 20 and a template family's best normalized overlap is 0.98. Find the recovered SNR. In a simple Euclidean uniform-source model with fixed threshold and all other factors fixed, estimate the retained detection volume if distance reach scales with recovered SNR.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
The SNR ratio is the fitting factor. Euclidean volume scales as the cube of the distance ratio; this is a population approximation with stated assumptions.
Compare with the solution
The recovered SNR is 0.98×20=19.6. The reach ratio is 0.98 and the retained volume fraction is 0.98³=0.941192, a loss of about 5.88 percent. This simplified scaling omits cosmological rates, orientation selection and search-background changes. It says nothing by itself about the bias in masses, spins or distance: parameter robustness needs additional tests beyond signal-recovery overlap.
Keep the idea
Take it with you
- A one-sided PSD and Fourier convention fix the dimensionless noise-weighted waveform norm.
- Matched-filter optimality and a Gaussian likelihood rely on explicit noise and signal assumptions.
- High overlap, accurate parameter inference and a validated waveform error budget are distinct achievements.
Read further
Original sources and deeper explanations to follow at your own pace.
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.