Question

matlab问题，求大神指教

matlab问题，求大神指教这个integral3哪里用错了

Answer 1

问题出在自定义函数上，点运算符不完整，而导致错误使用  *（内部矩阵维度必须一致）。 

完整的自定义函数为

fun3=@(x,y,z)z.^2.*log(x.^2+y.^2+z.^2+1)./(x.^2+y.^2+z.^2+1)

另外还存在的问题是

1、自定义函数fun3，后面附带约束条件是不正确的，可以去掉。

2、变量z的积分区间，应该为约束函数，即

zmin = 0；

zmax = @(x,y) sqrt(1 - x.^2 - y.^2);

3、处理后的代码

fun3=@(x,y,z)z.^2.*log(x.^2+y.^2+z.^2+1)./(x.^2+y.^2+z.^2+1) %.*(z>=0 & z<=sqrt(1-x.^2-y.^2));

xmin = -1;xmax = 1;

ymin = -1;ymax = 1;

zmin = 0;zmax = @(x,y) sqrt(1 - x.^2 - y.^2);

I=integral3(fun3,xmin,xmax,ymin,ymax,zmin,zmax)

4、最后，你还得检查一下，ymin，ymax ；zmin ，zmax 的约束函数是否正确

追问

我的书上除了integral3用的是triplequad

我只是换了triplequad按道理应该对的啊

这又错在哪呢

追答

这又错在哪呢，log前面*少了一个点。

用triplequad是可以的。应这样来写。

追问

但是按我书上的做法加点就不对了

追答

书上的不一定是对，可能有错