Question

matlab中dblquad函数求一给定参数二元整天分布的积分，数值怎么大于1

d=1;
a=0;
b=0;

q=dblquad(@(x,y)(2.*pi).^-0.5.*d.^(-1).*exp(-0.5.*d.^(-2).*((x-a).^2+(y-b).^2)), -1000, 1000, -1000,1000)

q =

2.5066

高斯分布积分不管怎么积也不能大于1啊，求大神看看哪里错了

Answer 1

你有点像当然了。

 

多变量正态分布的概率密度函数定义为

    

其中Σ为半正定的协方差矩阵（对σ²的推广）。

 

需要注意的是，指数项前面系数的分母里，2π并不是开平方，而是N/2次方。

 

对于二元的情况，N=2：

>> q=dblquad(@(x,y)(2.*pi).^-1.*d.^(-1).*exp(-0.5.*d.^(-2).*((x-a).^2+(y-b).^2)), -1000, 1000, -1000,1000)

q =

    1.0000