| 查看: 2264 | 回复: 19 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】feval的使用
|
|||
|
现在在看《数值分析导论》韩渭敏译，379页中的边值问题， 在对这个偏微分方程的求解过程中用到一条语句 feval(f,(1:n-1)*h,(1:n-1)*h) 注释中说明f方程的右边的函数值，我不知道如何写这个f的m函数，有人可以指点一下啊  |
回复此楼» 猜你喜欢
材料工程085601，270求调剂
已经有42人回复
-
279求调剂
已经有11人回复
-
271求调剂
已经有24人回复
-
297工科调剂?
已经有6人回复
-
材料工程281还有调剂机会吗
已经有35人回复
-
0856专硕求调剂 希望是a区院校
已经有27人回复
-
考研调剂
已经有10人回复
-
求助调剂，跨调
已经有14人回复
-
085600材料与化工，求调剂
已经有14人回复
-
求调剂
已经有12人回复
» 抢金币啦！回帖就可以得到:
-
安徽农业大学李亚辉教授招收 化学、农药学、植保类2026级博士研究生
+3/192
-
双一流南京医科大学招计算机、生物信息等方向研究生调剂
+1/138
-
物理学0702 调剂
+1/89
-
武汉工程大学联合培养招收硕士（食品、环境和化学相关专业）
+1/88
-
烟台大学环境与材料工程学院先进储能材料与器件课题组招收材料、化学、环境方向硕士
+1/86
-
齐鲁工业大学轻工学部曹珊副教授接收硕士调剂1人
+1/38
-
岭南师范学院化学化工学院招收生物与医药制药工程专业硕士调剂
+1/37
-
山东航空学院安全工程专业085702第三轮接收研究生调剂，欢迎咨询
+1/37
-
大湾区大学李红庚课题组招聘研究助理
+1/32
-
【香港城大】范战西教授课题组现招1名博士生（全奖, 2026年秋季 或 2027年春季 入学）
+1/18
-
省重点实验室招收调剂生（一志愿07开头专业），海洋生物专业名额充足
+1/15
-
贵州铜仁学院研究生招生
+1/11
-
新出名额：欢迎药学、生物、人工智能、计算机等有班长团支书经历加入-上海
+1/10
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
+1/9
-
浙江农林大学 林业与生物技术学院 生物学(生物化学与分子生物学) 第二轮调剂生招生
+1/8
-
重庆大学附属涪陵医院2026年度全职博士后招聘简章（年薪30-65万）
+1/8
-
活性环氧稀释剂
+1/4
-
26申博，ai+材料方向，求课题组推荐
+1/3
-
二维材料方向博士后招聘
+1/2
-
汕头大学化学系陈广慧教授课题组招收计算化学方向考研调剂生
+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
