li_huiyuan 发表于 2009-3-17 19:15

三重积分vpa问题

做论文需要三重积分，用matlab中的int函数做的，积分结果为：f1 =
int(int(80669135367/537500000000*pi^(1/2)*14190^(1/2)*exp(-73189/10000*y+221341123/20000000-exp(-73189/10000*y+606809999/100000000)-5412/625*z-exp(-5412/625*z+7810869/1562500))-80669135367/537500000000*erf(47/625*14190^(1/2)-1/25*14190^(1/2)*y-1/25*14190^(1/2)*z)*pi^(1/2)*14190^(1/2)*exp(-73189/10000*y+221341123/20000000-exp(-73189/10000*y+606809999/100000000)-5412/625*z-exp(-5412/625*z+7810869/1562500)),y = -inf .. inf),z = -inf .. inf)
接下来用vpa(f1,10),可总显示错误，请问各位大侠，这是怎么回事啊？本人初学，请大家帮帮忙。。。

[ 本帖最后由 sigma665 于 2009-3-18 08:51 编辑 ]

friendchj 发表于 2009-3-17 19:33

回复 楼主 li_huiyuan 的帖子

估计这个积分没有解析式，所以matlab肯本没有给出解析的结果，自然不能用vpa了

li_huiyuan 发表于 2009-3-17 19:38

错误提示显示：??? Error using ==> sym/maple
Error, (in evalf/int) function does not evaluate to numeric

Error in ==> C:\MATLAB6p5\toolbox\symbolic\@sym\vpa.m
On line 32==> r = maple('evalf',s,d);
我把y的积分限改为-3，+3；z的积分限改为-3，+3；
然后用vpa（f1，4）却可以得到结果：.3402e-1

sigma665 发表于 2009-3-18 08:51

回复 板凳 li_huiyuan 的帖子

把无穷大换个大数试试

li_huiyuan 发表于 2009-3-18 14:56

换过了，也不好使。哪位大侠把f1的表达式拷到matlab程序里试一试，然后再回答我啊.

sigma665 发表于 2009-3-18 16:42

我matlab出问题，算不了符号

[ 本帖最后由 sigma665 于 2009-3-18 16:46 编辑 ]
页: [1]
查看完整版本: 三重积分vpa问题