research34 min studyLesson 840 in Physics

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

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.

a~(f)=∫−∞∞a(t)e−2πift dt,(a∣b)=4Re⁡∫0∞a~(f)b~∗(f)Sn(f) df,ρopt2=(h∣h)\widetilde a(f)=\int_{-\infty}^{\infty}a(t)e^{-2\pi ift}\,dt,\qquad (a|b)=4\operatorname{Re}\int_0^\infty\frac{\widetilde a(f)\widetilde b^*(f)}{S_n(f)}\,df,\qquad \rho_{\mathrm{opt}}^2=(h|h)
Sn is a one-sided strain PSD in Hz⁻¹, while Fourier amplitudes have seconds as units. ρopt is the expected optimal matched signal-to-noise ratio for the specified signal and ideal noise model. Calibrated response, integration band and analysis-window conventions remain part of a practical calculation.↔ Scroll wide equations horizontally.

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.

ln⁡L(ϑ)=−12(d−h(ϑ)∣d−h(ϑ))+const,O=(h1∣h2)(h1∣h1)(h2∣h2)\ln\mathcal L(\vartheta)=-\frac12(d-h(\vartheta)|d-h(\vartheta))+\mathrm{const},\qquad \mathcal O=\frac{(h_1|h_2)}{\sqrt{(h_1|h_1)(h_2|h_2)}}
The likelihood assumes the stated Gaussian noise covariance and known PSD. ϑ denotes source and nuisance parameters. The real normalized overlap lies between −1 and 1 before any time/phase maximization; a fitting factor additionally optimizes the intended template family. Neither expression removes prior dependence or detector-calibration uncertainty.↔ Scroll wide equations horizontally.

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.

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.

  1. 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.

  2. 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.

  3. 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 result

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.

Explore the model

Interactive lab

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.

Noise weighting rescales data directionsThe 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. Current readings: Area scale · det A: 1.4; Singular values · σmax, σmin: 1.597, 0.876; Condition number · κ₂(A): 1.823. All parameters are adjustable with the labeled controls below.Matrix A with first row 1.4, 0.6 and second row zero, one; display interpolates between identity and A while final determinant remains 1.4.−4−2024−2−1012x · equal geometric scaleyA = [ 1.4 0.6 ; 0 1 ]Display blend: 0% · final values below
Interpolated geometryOriginal unit squareSelected final transform
Area scale · det A
1.4
Singular values · σmax, σmin
1.597, 0.876
Condition number · κ₂(A)
1.823
t = 0 s

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.

Check your understanding

Which conclusion follows from a high best noise-weighted template overlap?

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.

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.