| 查看: 131 | 回复: 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人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有3人回复
-
国自然结果
已经有5人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有12人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有8人回复
-
据悉今年马上要出结果了
已经有8人回复
-
关于代码变化问题，想知道的进来
已经有23人回复
-
FileCode能看出啥？
已经有6人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有20人回复
-
2026国自然放榜时间
已经有9人回复
20