嵌套二重积分的插值解法和RBF网络解法

pasuka

23# wanglu

lz计算shunfly示例所用积分式

解析结果

219531

数值积分

219532
messenger 发表于 2009-10-12 10:20

贴一个利用ode45计算积分的方法

1. 
   
2. %%
   
3. clear all;clc
   
4. num = 5;
   
5. result = zeros(num,2);
   
6. %%
   
7. tic,
   
8. % 定义匿名函数，增加 0*y 项
   
9. fc = @(k) @(x,y) besselj(1,k.*x).*exp(-x)+0*y;
   
10. % 利用龙格库塔法计算积分
    
11. for k = 1:num
    
12.     [~,yy2] = ode45(fc(k),[0,100],0);
    
13.     result(k,1) = yy2(end);
    
14. end
    
15. toc;
    
16. %%
    
17. tic,
    
18. % 定义匿名函数，此处略有不同，没有点乘
    
19. fd = @(k) @(x,y) besselj(1,k*x)*exp(-x)+0*y;
    
20. % 使用数组函数
    
21. temp = arrayfun(@(k) ode45(fd(k),[0,100],0),1:num);
    
22. % 导出结果
    
23. for k = 1:num
    
24.     result(k,2) = temp(k).y(end);
    
25. end
    
26. toc;
    

复制代码

计算结果

1. 
   
2. result =
   
3.    0.292912105440771   0.292912105440771
   
4.    0.276513031388514   0.276513031388514
   
5.    0.228445753183351   0.228445753183351
   
6.    0.188887531708083   0.188887531708083
   
7.    0.160933838415485   0.160933838415485
   

复制代码