| 查看: 107 | 回复: 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 ] |
回复此楼» 猜你喜欢
材料调剂
已经有10人回复
-
一志愿哈尔滨工业大学085600英一数二337分求调剂
已经有5人回复
-
288求调剂 一志愿哈工大 材料与化工
已经有11人回复
-
求调剂，一志愿南京航空航天大学 ，080500材料科学与工程学硕
已经有9人回复
-
本科211，293分请求调剂
已经有4人回复
-
求调剂
已经有3人回复
-
085400电子信息319求调剂（接受跨专业调剂）
已经有5人回复
-
268求调剂
已经有9人回复
-
一志愿双非085502，267分，过四级求调剂
已经有3人回复
-
一志愿085404，总分291，四级已过，求调剂
已经有4人回复
