MATLAB三重积分问题

时间:2019-07-25 16:00:07

标签: matlab convolution gaussian integral

我要评估的方程式:https://i.stack.imgur.com/W0HxL.jpg

我试图评估这个积分(和等式),其中“ a”和“ x”只是输入,而x′是(x,y,z)。还:x将具有[a,b,c]的形式。在下面的“ int_func”中,我采用了转置和原始点的点积,以产生下面可以观察到的^ 2幂,例如在(X(:,1)-x)。^ 2

到目前为止,我的代码是:

int_func = @(x,y,z) func([x,y,z]).*exp((-3/(2*a^2)).*((X(:,1)-x).^2 + (X(:,2)-y).^2 + (X(:,3)-z).^2));

val = integral3(int_func,-inf,inf,-inf,inf,-inf,inf);
disp(val * (3/(2*pi*a^2))^(3/2));

其中:

X = [0, 0, 0];

a = 1;

使用此代码,我应该得到的答案是 0.0370,但是我得到的是0.5556。

有人可以建议另一种评估此积分的方式还是找出我列出的内容有什么问题?

0 个答案:

没有答案