matlab2011a怎样计算重积分啊？

牛行天下

怎样用matlab实现内积分限是外积分限函数的重积分？其中a，b是常数，谢谢

bainhome

2012b版本之后，应能使用新函数:

1. 
   
2. f=@(x,y) sin(x+y).*cos(x+2*y);
   
3. g=@(y) sin(y);
   
4. q = integral2(f,0,1,0,g)
   
5. q =
   
6.     0.0530

复制代码

老版本参看http://forum.vibunion.com/forum. ... =1&extra=#pid429038和http://forum.vibunion.com/forum. ... p;page=1&extra=。
不清楚2011a能不能用integral2这个函数，或可参见quad2d函数，这个一般区域积分的函数自2009a之后貌似就有了。

牛行天下

bainhome 发表于 2013-8-3 22:20

                               
登录/注册后可看大图

2012b版本之后，应能使用新函数:
老版本参看这里和这里。
不清楚2011a能不能用integral2这个函数，或可参 ...

嗯，谢谢了哈，刚才我试了试，integral2函数好像在matlab2011a不能用

米斯兰达

假设楼主的f(x,y)是x+y，而φ(y)=y^2吧 那么这个重积分可用下面的方法进行计算

s=quadl(@(x)arrayfun(@(xx)quadl(@(y)xx+y,0,xx.^2),x),a,b)

注意如果f(x,y)和φ(y)里面有*和/，一定要是.*和./

牛行天下

米斯兰达 发表于 2013-8-4 17:58

                               
登录/注册后可看大图

假设楼主的f(x,y)是x+y，而φ(y)=y^2吧 那么这个重积分可用下面的方法进行计算

s=quadl(@(x)arrayfun(@( ...

嗯，谢谢了哈，现在能想到的函数基本上都用了，不大行啊，可能要自己编程了

bainhome

就目前所给出这点有限的信息，没觉得哪里用quad2d解决不了，调用方式一样：

1. q = quad2d(f,0,1,0,g)
   
2. q =
   
3. 0.3500

复制代码

正由于quadl、dblquad等函数系列的算法有局限，例如默认调用方式只能用于矩形区域、振荡积分的易奇异等等，09版本之后对这个积分系列连续做了几次调整，现在已经相对成熟，自编函数如果不是对MATLAB底层函数有所了解，很难比之更好——因为我自己很早前写过一些。
个人看法，仅供参考

米斯兰达

牛行天下 发表于 2013-8-4 20:46

                               
登录/注册后可看大图

嗯，谢谢了哈，现在能想到的函数基本上都用了，不大行啊，可能要自己编程了

是在不行就用蒙特卡洛法吧 点数取多一点

牛行天下

bainhome 发表于 2013-8-4 22:04

                               
登录/注册后可看大图

就目前所给出这点有限的信息，没觉得哪里用quad2d解决不了，调用方式一样：
正由于quadl、dblquad等函数系 ...

嗯，非常感谢您的指点，也可能是我对这些函数不够了解，导致积分时老报错，

牛行天下

米斯兰达 发表于 2013-8-5 07:35

                               
登录/注册后可看大图

是在不行就用蒙特卡洛法吧 点数取多一点

嗯，谢谢您的提示，可能的话就尝试一下