时间戳

首页 > 解决方案 > 为什么我的高斯 pdf 没有整合到一个？

python - 为什么我的高斯 pdf 没有整合到一个？

问题描述

普通pdf：

import numpy as np
import scipy

def gaussian(x, mu = 0, sigma = 1):
  return 1/(np.sqrt(2*np.pi*sigma)) * np.exp(-(x-mu)**2 / (0.5*sigma)**2)

整合整个支持：

scipy.integrate.quad(gaussian, -np.inf, np.inf)

返回

(0.3535533905932738, 1.4635936470160148e-11)

我知道我在 pdf 的某个地方搞砸了，但我已经开始研究了一个小时，但我看不到它

标签： pythonnumpyscipy

解决方案

你的高斯函数不正确。它应该是：

def gaussian(x, mu = 0, sigma = 1):
  return (1/(sigma * np.sqrt(2 * np.pi))) * np.exp((-(x-mu)**2) / (2 * sigma ** 2))

推荐阅读