В любом случае вам нужно сохранить непрерывную случайную величину параметризованным способом, если вы не работаете с приближением. В этом случае вы обычно работаете с объектами распределения; а для известных распределений существуют формулы для дифференциальной энтропии в терминах параметров.
Scipy соответственно предоставляет entropy
метод для rv_continuous
, который вычисляет дифференциальную энтропию, где это возможно:
In [5]: import scipy.stats as st
In [6]: rv = st.beta(0.5, 0.5)
In [7]: rv.entropy()
Out[7]: array(-0.24156448)