我想使用SciPy的去卷积函数来找到两个高斯分布的未知分布.在SciPy中没有与此功能相关的文档,所以我只是在寻找一个关于如何在我的情况下使用此功能的示例.例如,给定两个正态分布N(100,1),N(300,2),我想了解我怎样才能找到去卷积n的分布(200,1).
>>> sample1 = np.round(scipy.around(scipy.stats.norm(100, 1).rvs(size=1000)))
>>> sample2 = np.round(scipy.stats.norm(300, 2).rvs(size=2000))
>>> signal.deconvolve(sample1, sample2)
Run Code Online (Sandbox Code Playgroud)
上面的代码给出了负值,这似乎是错误的.如何从此反卷积中恢复分布N(200,1)?特别是,我认为我的问题是我不明白如何得到除数.
我真正想要看的是如何使用SciPy的反卷积从这些样本中恢复~N(200,1)的示例.