matched_filter_snr — ONED `matched` op

• 数据种类:signal × signal × pairs → signal

• 调用: 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'(要直接调用实现,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';从台账取用则 ops1d.get("matched_filter_snr"))

用法

匹配滤波的信噪比时间序列 |ρ(t)|（对相位取最大），以 signal 返回。

> 以下的详细说明为原文 —— 摘要与标题已翻译。

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 が重力波探索の慣例。

参考(示例数据、文献)

• 示例数据目录(下载 URL / 许可证) —— 2-D 用 skimage.data(BSD/公有领域)加合成图,3-D 给出真实数据源(Stanford/PDS 等)的下载 URL。

• 算子来历与参考文献 —— 该算子族所依据的研究/方法出处。

• 算法的正典(作者、年份)与用途见上面的族使用指南。

可运行的示例(实际调用该算子并已验证的样例)

• poc_gw_detector_pseudo_simulation — py -3.11 examples/poc_gw_detector_pseudo_simulation.py

类型可衔接的下一个算子(可接受 `signal` 作为输入)

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

同类别(`matched`)

psd_welch · whiten_by_psd · chirp_newtonian · chirp_mass_from_track · optimal_snr · coincidence_delay

*Provenance: gwsearch.py — ONED 算子登记表。本条目由 tools/opdocs.py md 自动生成(请勿手工编辑)。*

matched_filter_snr — ONED matched op

用法