1

我正在尝试对Generalized Normal Distribution进行数值积分。我希望它下面的面积是1。我得到的是-1。为什么会这样？

mu=0;
alpha=0.5;
beta=2;

x=-2:0.01:2;
densityGG = beta/(2*alpha*gamma(1/beta)) * exp(-power(abs(x-mu)/alpha,beta));

plot(x,densityGG);
trapz(densityGG,x)
4

2 回答 2

4

该功能实际上是

trapz(X, Y)

所以试试

trapz(x, densityGG)

它应该解决问题。

于 2011-08-04T15:47:49.723 回答
2

阅读Trapz的帮助文件......它就在那里

你需要写：trapz(x,densityGG)

于 2011-08-04T15:49:54.497 回答