matched_filter_snr — ONED `matched` op

• Data kinds: signal × signal × pairs → signal

• Call: import fullseye as fs; fs.ledger.matched_filter_snr(data, template, psd, rate: 'float' = 4096.0, *, f_low_hz: 'float | None' = 20.0, f_high_hz: 'float | None' = None, template_ref_index: 'int | None' = None, taper_s: 'float' = 0.0) -> 'np.ndarray' (to call the implementation directly, import gwsearch; gwsearch.matched_filter_snr(data, template, psd, rate: 'float' = 4096.0, *, f_low_hz: 'float | None' = 20.0, f_high_hz: 'float | None' = None, template_ref_index: 'int | None' = None, taper_s: 'float' = 0.0) -> 'np.ndarray'; from the registry, ops1d.get("matched_filter_snr"))

Usage

Matched-filter SNR time series |ρ(t)|, maximised over phase, as a signal.

> The detailed description below is the original text — the summary and the headings are translated.

Parameters

data : array_like

観測時系列(1-D)。

template : array_like

テンプレート(data 以下の長さ。末尾を 0 で埋めて同じ長さにする)。

psd : array_like

雑音の片側 PSD の pairs(psd_welch)。

rate : float

標本化周波数 [Hz]。

f_low_hz, f_high_hz : float or None

内積の帯域。

template_ref_index : int or None

テンプレート内の「事象の時刻」の標本番号(None なら argmax |h|)。返り値はこの分だけ

回してあり、|ρ| の山の標本番号 = data 内の事象の時刻になる。

taper_s : float

data の両端を余弦で 0 へ落とす長さ [s] (Tukey 窓の端)。実データの端の不連続が

FFT の循環で全体に漏れるのを抑える(GWOSC の例題は 32 s に対し 2 s)。

Returns

numpy.ndarray

data と同じ長さの非負配列。ρ(t) = ⟨d|h_t⟩/σ、⟨a|b⟩ = 4Re∫ÃB̃*/S df、

σ = optimal_snr。位相は解析信号で最大化する(|z|)。テンプレートが data の外へ

はみ出す時刻(先頭の ref 標本と末尾の len(template) − ref 標本)は循環の折り返し

なので 0 にしてある(山の探索に混じらない)。

Notes

白色ガウス雑音だけなら |ρ| の典型値は ~1.25(レイリー)、閾値 8 が重力波探索の慣例。

References (sample data, literature)

• Sample-data catalog (download URLs / licences) — 2-D uses skimage.data (BSD/public domain) plus synthetic images; 3-D lists download URLs for real data sources (Stanford, PDS, …).

• Operator provenance and references — the sources of the research/methods this op family came from.

• The canonical algorithm (author, year) and its uses are named in the family usage guide above.

Runnable examples (verified samples that actually call this op)

• poc_gw_detector_pseudo_simulation — py -3.11 examples/poc_gw_detector_pseudo_simulation.py

Ops the type connects to (they accept `signal` as input)

create_funct_1d_array · create_funct_1d_pairs · smooth_funct_1d_gauss · smooth_funct_1d_mean · derivate_funct_1d · integrate_funct_1d · zero_crossings_funct_1d · local_min_max_funct_1d

Same category (`matched`)

psd_welch · whiten_by_psd · chirp_newtonian · chirp_mass_from_track · optimal_snr · coincidence_delay

*Provenance: gwsearch.py — ONED operator registry. This per-op note is generated by tools/opdocs.py md (do not hand-edit).*

matched_filter_snr — ONED matched op