24小时热门版块排行榜

>论坛更新日志 (2974)
>虫友互识 (344)
>休闲灌水 (170)
>导师招生 (124)
>文献求助 (112)
>论文道贺祈福 (53)
>博后之家 (38)
>论文投稿 (32)
>硕博家园 (26)
>招聘信息布告栏 (24)
>考博 (24)
>考研 (21)
>基金申请 (18)
>找工作 (14)
>公派出国 (14)
>教师之家 (12)

返回列表

当前主题已经存档。

autumnlyy

铜虫 (小有名气)

应助: 1 (幼儿园)
金币: 118.6
帖子: 63
在线: 10.1小时
虫号: 687587
注册: 2009-01-04
性别: GG
专业: 光学

[交流] 【求助】请教matlab 重积分的问题

一个简单的问题：
假如在一个二维平面上，x1，x2，……xn
                              y1，y2，……yn
上面每一个点的数值我都知道
我想把这些点的值组成的一个函数z=f（xm，yn）（并没有解析表达式）做一个积分
相当于，对z=f（x，y）从x1->xn, y1->yn作一个重积分

我利用matlab的辛普森公式做以下的操作，是否能够得出我想要的结果？

for jj=1:n
      IIxtmp(1,jj)=(dh/3)*(z(jj,1)+4*sum(z(jj,2:2:n-1))+2*sum(z(jj,3:2:n-2))+z(jj,n));  %对x方向先积分
      IIytmp(1,jj)=(dh/3)*(z(1,jj)+4*sum(z(2:2:n-1,jj))+2*sum(z(3:2:n-2,jj))+z(n,jj));  %对y方向先积分
end

IIx=(dh/3)*(IIxtmp(1)+4*sum(IIxtmp(2:2:n-1))+2*sum(IIxtmp(3:2:n-2))+IIxtmp(n)); %对y方向后积分
      IIy=(dh/3)*(IIytmp(1)+4*sum(IIytmp(2:2:n-1))+2*sum(IIytmp(3:2:n-2))+IIytmp(n)); %对x方向后积分

II=0.5*(IIx+IIy);  %积分的结果
我采用上面的方法，能够的到我想要的结果么？

其中h为步长

[ Last edited by ddx-k on 2009-2-12 at 19:50 ]

回复此楼