求救：百思不解的一个符号积分结果

szldh2005

有一个符号积分，在不同的机上竟然会有不同的计算结果。百思不解。请高手解惑！是matlab的设置问题还是matlab的本身的问题？还是同时安装了matlab2008a,2010a的干扰问题？谢谢！

1. clc;clear;
   
2. syms x y u b
   
3. u=10^-54
   
4. 
   
5. 
   
6. %%%%%%% B 符号积分->数值计算 %%%%%%
   
7. 
   
8. % 一、原无限积分
   
9. disp('%%%%%% B 符号积分->数值计算 %%%%%%%%')
   
10. disp('一、原无限积分')
    
11. b=-log(10^-15);
    
12. f1=@(x) exp(-x)./x;
    
13. w1k=@(u) vpa(int(f1(x),x,u,inf),20); % 符号积分-数值计算(原无限积分)
    
14. w1d=@(u) vpa(int(f1(x),x,u,b),20); % 符号积分-数值计算(原无限截断积分)
    
15. w1k=w1k(u)
    
16. w1d=w1d(u)
    
17. 
    
18. 
    
19. % 二、等价无限积分
    
20. disp('二、等价无限积分')
    
21. b=log(-log(10^-15));
    
22. f2=@(y) exp(-exp(y));
    
23. w2k=@(u) vpa(int(f2(y),y,log(u),inf),20); % 符号积分-数值计算(等价无限积分)
    
24. w2d=@(u) vpa(int(f2(y),y,log(u),b),20); % 符号积分-数值计算(等价截断积分)
    
25. w2k=w2k(u)
    
26. w2d=w2d(u)
    
27. 
    
28. 
    
29. % 三、等价有限积分1
    
30. disp('三、等价有限积分1')
    
31. f3=@(y) 1./log(y);
    
32. w3k=@(u) -vpa(int(f3(y),y,0,exp(-u)),20); % 符号积分-数值计算(等价有限积分1)
    
33. w3k=w3k(u)
    
34. 
    
35. 
    
36. % 四、等价有限积分2
    
37. disp('四、等价有限积分2')
    
38. f4=@(y) exp(-u)./(u-log(y));
    
39. w4k=@(u) vpa(int(f4(y),y,0,1),20); % 符号积分-数值计算(等价有限积分2)
    
40. w4k=w4k(u)
    
41. 
    
42. 
    
43. % 五、等价有限积分3
    
44. disp('五、等价有限积分3')
    
45. f5=@(y) exp(-u./y)./y;
    
46. w5k=@(u) vpa(int(f5(y),y,0,1),20); % 符号积分-数值计算(等价有限积分3)
    
47. w5k=w5k(u)
    

复制代码

nwcwww

08a和10a的结果不同很正常。08a的时候symbolic工具箱内核用的是maple, 08b开始变成了mupad。
部分maple能解决的问题mupad搞不定，而且mupad总体来说稍慢。但是mupad的界面做的不错，而且由于掌握在mathworks手中，从产品的兼容和统一性来说应该更好。而且mupad在最近的版本中也在不断改进。

szldh2005

在winxp3，2010a版本中，不同的机器得到不同的结果（而且结果差得太大！），这才是困扰我很久的问题。

szldh2005

已搞清楚。主要原因是matlab2010a在winxp2中的符号计算有bug。由于我笔记本的winxp3是由在winxp2中安装应用软件升级到winxp3的，其应用软件性能可能继承了winxp2的性能。在独立安装winxp3后再安装matlab2010a，则符号计算正常。