В исследовании «Новая модель распада НАСА для Evolve 4.0» предоставляется реализация для массового распределения фрагментов, образовавшихся в результате столкновений и взрывов.
Распределение определяется следующим образом:
mu определяется как среднее значение,
sigma определяется как стандартное отклонение,
alpha - коэффициент масштабирования,
lambda_ {c} - это log10 (характеристика длины c)
При реализации этого распределения, скажем, в python, будет ли уместно складывать вместе масштабированные функции нормального распределения вероятностей? Или я пробую каждое нормальное распределение для значения, масштабирую их и складываю вместе?
Мой вопрос касается реализации этого в python. Первый способ реализовать это было бы использовать numpy для генерации случайного значения в каждом нормальном распределении с помощью np.random.normal (), поэтому я бы определите выражение как alpha * np.random.normal (...) + (1 - alpha) * np.random.normal (). В качестве альтернативы я мог бы использовать SciPy для вычисления функций плотности вероятности нормальных распределений и добавить их: alpha * norm.pdf (...) + (1 - alpha) * norm.pdf (...). Я не уверен, какой из этих методов приведет к правильной реализации.