Question

MATLAB 里二重积分integral2里带矩阵怎么算啊？？

fun = @(x,y) x+y;
q = integral2(@(x,y) fun(x,y),0,1,0,1)；
没有问题，但是加入一个矩阵就一直在报错：
m=[1,2,3];
fun = @(x,y,m) x+y+m;
q = integral2(@(x,y) fun(x,y,m),0,1,0,1,'ArrayValued',true)

这么回事啊！求大神解答！！！！

Answer 1

可以直接调用MATLAB向量化积分函数quadv或integral，这两个函数是MATLAB提供的仅有的两个支持向量化积分的函数。具体做法是：将被积函数f(x,t) 定义为可接受向量化参数t的关于x的函数句柄，然后直接调用quadv或integral，其中，integral函数要求将'ArrayValued' 参数设置为 true。
integral(@(x) f(x,t), 0, 1, 'ArrayValued', true)

出自MATLAB论坛帖子《MATLAB高效编程之向量化积分 》--winner245
我预测integral2同理。

Answer 2

x^2+y^2+z^2+ax+by+cz+d=0当a^2+b^2+c^2-4d>0时表示球。
x^2/a^2+y^2/b^2+z^2/c^2=1当a>=b>=c>0时表示椭球面。
x^2/a^2+y^2/b^2-z^2/c^2=1当a>=b>0,c>0时表示单叶双曲面。
z^2/c^2-x^2/a^2-y^2/b^2=1当a>=b>0,c>0时表示双叶双曲面。
x^2/a^2+y^2/b^2-z^2/c^2=0当a>=b>0,c>0时表示锥面。
x^2/a^2+y^2/b^2=z当a>=b>0时表示椭圆抛物面。（你的方程就是这一类型）
x^2/a^2-y^2/b^2=z当a>=b>0时表示双曲抛物面。

Answer 3

integral2 函数没有 'ArrayValued' 选项，应该是不能向量化积分的。不知道你解决了吗，最近也遇到这种计算，用循环的话，感觉计算效率有点不能接受。