如何注入零噪声信号紧凑二进制合并信号



是否可以在没有彩色高斯噪声的情况下自行注入信号?

Arunava Mukherjee通过电子邮件提出的问题

是的。有两种简单的方法可以做到这一点。

1( 使用现有的辅助函数

生成干涉仪对象时,bilby提供了几个用 bilby.gw.detector.get_interferometer_with... 表示的辅助例程。在这种情况下,您需要使用此函数(我已经截断了文档(

bilby.gw.detector.get_interferometer_with_fake_noise_and_injection(
    name, injection_parameters, injection_polarizations=None, 
    waveform_generator=None, sampling_frequency=4096, duration=4, 
    start_time=None, outdir='outdir', label=None, plot=True, save=True, 
    zero_noise=False)
Docstring:
Helper function to obtain an Interferometer instance with appropriate
power spectral density and data, given an center_time.
Note: by default this generates an Interferometer with a power spectral
density based on advanced LIGO.
Parameters
----------
    name: str
        Detector name, e.g., 'H1'.
...
    zero_noise: bool
        If true, set noise to zero.

因此,您只需将标志传入,它将创建一个仅包含注入信号的干涉仪(然后您需要在传入可能性的干涉仪列表中为您想要的每个干涉仪创建一个。

2( 使用低水平集应变数据方法

或者,您可能希望使用低级方法本身。作为一般经验法则,您始终可以查看通用帮助程序函数的源代码,以确定应如何执行此操作。在这里,我们创建了一个H1干涉仪,以零噪声设置应变数据并注入信号:


interferometer = get_empty_interferometer("H1")
interferometer.power_spectral_density = PowerSpectralDensity.from_aligo()
interferometer.set_strain_data_from_zero_noise(
    sampling_frequency=sampling_frequency, duration=duration,
    start_time=start_time)
injection_polarizations = interferometer.inject_signal(
    parameters=injection_parameters,
    waveform_generator=waveform_generator)

截至 v.0.3.5 时正确的信息

最新更新