| 查看: 313 | 回复: 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 ] |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有10人回复
-
什么时候能放榜呀？
已经有3人回复
-
一个有机合成实验室都需要哪些设备？
已经有8人回复
-
UV压敏胶开发
已经有5人回复
-
面上提前没消息，有中的吗
已经有16人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有25人回复
-
2026年国自然面上资助率
已经有22人回复
-
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有12人回复
-
微信指数没变化，科研之友没阅读
已经有19人回复
2楼2009-02-01 12:27:39
3楼2009-05-30 09:24:08