| 查看: 255 | 回复: 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 ] |
回复此楼» 猜你喜欢
三无产品还有机会吗
已经有6人回复
-
垃圾破二本职称评审标准
已经有7人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有7人回复
-
博士申请都是内定的吗？
已经有14人回复
-
谈谈两天一夜的“延安行”
已经有13人回复
-
氨基封端PDMS和HDI反应快速固化
已经有11人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有11人回复
-
论文投稿求助
已经有4人回复
-
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛
已经有3人回复
2楼2009-02-01 12:27:39
3楼2009-05-30 09:24:08