Создайте гистограмму распределения Коши с нижним и верхним пределом

Question

Я хотел бы смоделировать следующее распределение Лоренца с гистограммой

^? / _2π
L = ————————
(E - E ₀ ) + 0,25 ? ²

Я нашел scipy.stats.cauchy и хотел бы обрезать распределение по нижнему и верхнему пределам примерно так:

L = cauchy.rvs(size=300, loc = 5, scale =2.5, limits = [0,15] )

Возможно ли это?

Answer 1

Нельзя добавить limits к методу rvs. Насколько я знаю, только truncnorm может сделать это. Что вы можете сделать, это либо обрезать значения с помощью scipy.clip (или numpy.clip), либо отфильтровать значения за пределами ваших пределов с помощью маски.

Первый метод создаст много 0 и 15:

import scipy as sp
L = sp.clip(cauchy.rvs(size=300, loc = 5, scale =2.5), 0, 15)

Второй будет случайным образом распределен в вашем интервале:

import scipy as sp
L = cauchy.rvs(size=10000, loc = 5, scale =2.5), 0, 15)   #create a larger set to filter it out
L = L[sp.logical_and(L<15,L>0)][:300]