Функция отображения логарифма PDF в Matlab

Question

Мне нужно рассчитать плотность 50 в стандартном нормальном распределении, на самом деле результат вернется к 0 pdf('norm',50,0,1), затем я пытаюсь отобразить логарифм этого, он не работает с log(pdf('norm',50,0,1)), В R мы можем сделать как dnorm(50,log = TRUE)Любая подобная функция в MATLAB? Я не нашел это. Спасибо

Answer 1

Если вы вычислите pdf, оцененный как 50, вы получите 0 из-за конечной точности (в частности, результат меньше realmin). Вы должны вычислить логарифм напрямую. Из определения, логарифм стандартного нормального PDF, оцененный в 50

>> format long %// to show more decimals
>> -50^2/2 - .5*log(2*pi)
ans =
   -1.250918938533205e+003