| 查看: 106 | 回复: 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 ] |
回复此楼» 猜你喜欢
情人节自我反思：在爱情中有过遗憾吗？
已经有11人回复
-
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有12人回复
-
过年走亲戚时感受到了所开私家车的鄙视链
已经有9人回复
-
今年春晚有几个节目很不错，点赞！
已经有10人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有5人回复