| 查看: 2205 | 回复: 19 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】feval的使用
|
|||
|
现在在看《数值分析导论》韩渭敏译，379页中的边值问题， 在对这个偏微分方程的求解过程中用到一条语句 feval(f,(1:n-1)*h,(1:n-1)*h) 注释中说明f方程的右边的函数值，我不知道如何写这个f的m函数，有人可以指点一下啊  |
回复此楼» 猜你喜欢
284求调剂
已经有10人回复
-
一志愿山东大学药学学硕求调剂
已经有4人回复
-
07化学280分求调剂
已经有4人回复
-
298-一志愿中国农业大学-求调剂
已经有12人回复
-
求材料，环境专业调剂
已经有3人回复
-
335求调剂
已经有5人回复
-
求调剂
已经有7人回复
-
一志愿吉大化学322求调剂
已经有4人回复
-
环境学硕288求调剂
已经有8人回复
-
341求调剂(一志愿湖南大学070300)
已经有6人回复
» 抢金币啦！回帖就可以得到:
-
中国科学院杭州医学研究-常皓研究员课题组博士后招聘
+1/276
-
福建师范大学叶应祥课题组招收2026级博士研究生
+2/132
-
中山大学环境材料团队招聘
+1/83
-
广东工业大学-木质纤维素高值化利用团队招博士研究生
+1/81
-
浙江大学侯阳教授团队诚聘博士后---高性能固态电解质与固态电池器件设计
+1/78
-
中山大学生医工学院课题组招硕士生（高分子/纳米材料/分析传感/荧光免疫检测背景）
+2/52
-
招收荧光探针/有机合成/生物学 博士后--薪金: 30万-50万
+1/31
-
湖北工业大学 生命科学与健康学院 课题组招收2026级食品/生物方向硕士
+1/28
-
巴黎西岱大学（Université Paris Cité）Fayna MAMMERI 导师寻找CSC公派博士生
+1/22
-
丽水学院 材料与能源工程学院 陈江飞老师课题组硕士招生----材料与化工专业
+1/16
-
双一流湘潭大学化工学院招收化工、能源动力、工程热物理相关研究生调剂
+1/14
-
宁波东方理工大学余鹏课题组2026年上海交大联培博士生、博士后、科研助理招聘
+1/6
-
211/双一流---石河子大学---有机化学方向招调剂生
+1/5
-
211/双一流石河子大学化学化工学院--电催化，CO2转化
+1/3
-
085602调剂
+1/3
-
中科院赣江院黄小磊实验室招聘特别研究助理—材料模拟与计算领域
+1/3
-
环境工程专业 2026年硕士研究生调剂名额（有学硕1名，专硕2名）
+1/3
-
丽水学院 资源与环境专业 硕士招生/调剂
+1/3
-
26申博自荐
+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
