| 查看: 264 | 回复: 2 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】请教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 ] |
回复此楼» 猜你喜欢
参与限项
已经有4人回复
-
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
2楼2009-02-01 12:27:39
3楼2009-05-30 09:24:08