一般化正規分布を数値的に統合しようとしています。その下の面積は 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)
一般化正規分布を数値的に統合しようとしています。その下の面積は 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)
機能は実際には
trapz(X, Y)
だから試してみてください
trapz(x, densityGG)
そして、物事を整理する必要があります。
Trapzのヘルプ ファイルを読んでください。
書く必要があります: trapz(x,densityGG)