| 查看: 2527 | 回复: 19 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】feval的使用
|
|||
|
现在在看《数值分析导论》韩渭敏译，379页中的边值问题， 在对这个偏微分方程的求解过程中用到一条语句 feval(f,(1:n-1)*h,(1:n-1)*h) 注释中说明f方程的右边的函数值，我不知道如何写这个f的m函数，有人可以指点一下啊  |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有8人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
» 抢金币啦！回帖就可以得到:
-
告别流水线爽文！108万字完结仙侠《女娲石之灵石奇缘（仙界篇）》值得细品
+1/9018
-
散金祈福
+5/3970
-
基金祈福
+1/965
-
散金币
+1/933
-
散金祈福
+3/783
-
何不食肉糜
+1/483
-
华西生物制药研究院招聘科研助理
+2/400
-
五大联赛看好哪支球队
+5/185
-
Polymers 期刊特刊 组稿
+1/76
-
那年春暖花开，今朝如梦初醒 —— 北漂十五载，相亲十年，归乡孑然一身
+1/69
-
哈尔滨工业大学（校本部）电子与信息工程学院2027级博士招生
+1/63
-
上海大学王亮课题组诚招博士后（碳材料、光催化、电催化方向）
+2/44
-
瑞典-博后-固态电池-粘结剂-锂离子电池方向
+1/35
-
中山大学化学工程与技术学院欧阳钢锋/吴鑫课题组博士生（2027年入学）/博士后招聘启事
+1/33
-
南方医科大学深圳医院陈烨/余涛教授课题组联合培养博士后招聘
+1/6
-
南京大学 智能驱动与感知材料实验室 诚招申请考核博士生/科研助理/博士后
+1/4
-
清华大学环境课题组招聘环境工作专业客座生研究生2-3名（应用导向）
+1/2
-
2026年国自然
+1/1
-
清华大学环境课题组招聘环境工作专业客座生研究生2-3名（应用导向）
+1/1
-
找女友
+1/1
|
function U=poisson(f,g,n,tol,max_it) if nargin<5 max_it=10000 end if nargin<4 tol=1e-5 end n1=n+1; h=1/n; toln=(h^2)*tol; Fr=zeros(n,n); Fr(2:n,2:n)=h^2*feval(f,(1:n-1)*h,(1:n-1)*h); U=zeros(n1,n1); U(1,1:n1)=feval(g,0,(0:n)*h); U(n1,1:n1)=feval(g,1,(0:n)*h); U(1:n1,1)=feval(g,(0:n)*h,0); U(1:n1,n1)=feval(g,(0:n)*h,1); rel_err=1; itnum=0; while((rel_err>toln)&(itnum<=max_it)) err=0; umax=0; for j=2:n for i=2:n temp=(U(i+1,j)+U(i-1,j)+U(i,j+1)+U(i,j-1))/4-Fr(i,j); dif=abs(temp-u(i,j)); if(err<=diff) err=dif; end U(i,j)=temp; temp=abs(temp); if(umax<=temp) umax=temp; end end end itnum=itnum+1; rel_err=err/umax; end %X=(0:h,:n*h)'; %Y=X; %surf(X,Y,U') %xlabel('x-axis') %ylabel('y-axis') %zlabel('the numerical solution') %title('Plot of the mumerical solution') X=(0:h:n*h)';Y=X; subplot(1,2,1) surf(X,Y,U') xlabel('x-axis') ylabel('y-axis') zlabel('The numerical solution') s1=springtf('h=6.4f',h) title(s1) hold on Err=sin(pi*X)*sin(pi*Y')-U; subplot(1,2,2) surf(X,Y,Err') |
9楼2011-01-08 11:05:13
2楼2011-01-07 20:04:38
3楼2011-01-07 21:47:10
4楼2011-01-08 10:12:37
