| 查看: 105 | 回复: 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 ] |
回复此楼» 猜你喜欢
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复