| 查看: 134 | 回复: 0 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】网格点积分。
|
|||
|
在二维矩形区域网格节点处计算正常的数值积分，（详见附件），结果部分节点计算出现警告： Warning: Infinite or Not-a-Number function value encountered. 导致少量节点无计算结果，绘图从而出现“缺块”。做了几天，但不知如何解决问题，求助各位。 ＝＝＝＝＝＝＝＝＝＝＝＝＝ clear all clc P=1.2e-3;gamma=0.177e6;L=2e-3;w=8;Ls=22;a=30; coef=P/(pi*gamma*L*w);coef1=pi/(2*Ls); global coef1; Xmin=-a; Xmax=a; Ymin=0; Ymax=Ls; dx=1; dy=1; [x,y]=meshgrid(Xmin:dx:Xmax,Ymin:dy:Ymax); % 在区域 上生成XY平面网格， for yy=1:2*a+1 for xx=1:Ls+1 T(xx,yy)=coef*quad(@func,-w/2,w/2,[],[],x(xx,yy),y(xx,yy)); end end TT=T*1e6; surf(x,y,TT) view(-35,6) ＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝ function f=func(t,x,y) global coef1; f=atanh(sin(coef1*y)./cosh(coef1*(t+x))); [ Last edited by sunxiao on 2009-3-8 at 13:20 ] |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有7人回复
-
面上意见出来了
已经有23人回复
-
能否申诉？
已经有8人回复
-
梦想
已经有5人回复
30