Разница в Scipy.stats между нормой и многомерной нормой ⇐ Python

[Anonymous]

Функция scipy.stats «multivariate_normal» не дает желаемого распределения. Я использую многомерную нормаль, чтобы обеспечить общность в рамках физического проекта, но распределение z-компонента вектора не такое, как ожидалось, даже если матрица стандартного отклонения имеет только z-компонент, как показано ниже: /p>

Код: Выделить всё

from scipy import stats
import numpy as np

mu = np.array([0,0,2000])
sigma = np.diag([0,0,50])
vals = stats.multivariate_normal(mu, sigma, allow_singular=True).rvs(10000)

не дает желаемого распределения в значениях[2], как ожидалось, нормального распределения со стандартным отклонением 50 и средним значением 2000. Вместо этого среднее значение правильное, но стандартное отклонение намного меньше, чем должно быть.
Что я здесь делаю не так?

Подробнее здесь: https://stackoverflow.com/questions/792 ... riate-norm