python - 如何将高斯函数分割为等体积部分

Question

我正在尝试使用 Python 将高斯形曲线分割为 K 个等体积段，以实现信号过滤目的。

我正在寻找伪代码、总体思路或执行它的库。 任何帮助将不胜感激。

谢谢!

例如下图中:K=6。卷 s1 = s2 = ... = s6:

Answer 1

您需要确定分布的百分位数。您可以使用此 scipy.stats.norm 类及其 .ppf() 方法。

import numpy as np
import scipy.stats as sps
import matplotlib.pyplot as plt

mu = 25
sigma = 4
splits = 8

# define the normal distribution and PDF
dist = sps.norm(loc=mu, scale=sigma)
x = np.linspace(dist.ppf(.001), dist.ppf(.999))
y = dist.pdf(x)

# calculate PPFs
step = 1 / splits
quantiles = np.arange(step, 1.0 - step / 2, step)
ppfs = dist.ppf(quantiles)  # boundaries

# plot results
fig, ax = plt.subplots(figsize=(10, 4))
ax.plot(x, y, color='k')
for i, ppf in enumerate(ppfs):
    ax.axvline(ppf, color=f'C{i}', label=f'{quantiles[i]:.3f}: {ppf:.1f}')
ax.legend()
plt.show()

这是基于this answer