我需要在标准正态分布中计算50的密度，实际上，结果将由pdf('norm',50,0,1)返回0，然后我尝试显示它的对数，它不适用于{{ 1}}。在R中，我们可以像log(pdf('norm',50,0,1))那样在MATLAB中执行任何类似的功能吗？我没找到它。感谢
答案 0 :(得分：2)
如果计算以50计算的pdf，则由于精度有限而得到0（特别是，结果小于realmin）。您应该直接计算对数。从definition开始，在50处评估的标准正常pdf的对数是
>> format long %// to show more decimals
>> -50^2/2 - .5*log(2*pi)
ans =
-1.250918938533205e+003