Matlab中求二重积分的问题

hufengbin

大家好，我是初学matlab，主要是来解一个二重积分，积分如下：
请问大家该怎么求解呢？

matlab1234

1759553924
先离散，后求和了

qibbxxt

1. >> R=1,r=0.8,d=1.2,h=0.5,f=@(s1,s2)(cos(s1).*cos(s2)+sin(s1).*sin(s2))./sqrt((R*cos(s1)-r*cos(s2)).^2+(R*sin(s1)-r*sin(s2)-d).^2+h.^2)
   
2. 
   
3. R =
   
4. 
   
5.      1
   
6. 
   
7. 
   
8. r =
   
9. 
   
10.     0.8000
    
11. 
    
12. 
    
13. d =
    
14. 
    
15.     1.2000
    
16. 
    
17. 
    
18. h =
    
19. 
    
20.     0.5000
    
21. 
    
22. 
    
23. f =
    
24. 
    
25.     @(s1,s2)(cos(s1).*cos(s2)+sin(s1).*sin(s2))./sqrt((R*cos(s1)-r*cos(s2)).^2+(R*sin(s1)-r*sin(s2)-d).^2+h.^2)
    
26. 
    
27. >> quad2d(f,0,2*pi,0,2*pi)
    
28. 
    
29. ans =
    
30. 
    
31.     2.1391

复制代码

nwcwww

除了用quad2d函数，在MUPAD里面也可以算。
红色的是输入的代码，蓝色的是MATLAB每步的输出结果。