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

时间:2019-02-03 21:20:39

标签: bilby

是否可以自己注入没有彩色高斯噪声的信号?

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

1 个答案:

答案 0 :(得分:1)

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

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起正确