| 查看: 2206 | 回复: 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/179
香港中文大学(深圳)靳羽华教授交叉实验室招募2026年材料学博士生(光致变色)
+2/134
北京-89175-事业单位-诚征女友
+1/90
中国电气装备集团研究院招收AI+材料等相关岗位实习生
+1/85
【全额奖学金+英国博士学位】国科大杭高院Nathan课题组招收2026级联培博士生1名!
+1/49
武汉高校接受研究生调剂,材料、化学、纺织方向
+2/44
天津城建大学招收物理、化学或者生物调剂生一名,前往天津大学联合培养
+1/33
春眠不觉晓
+1/21
同济大学环境学院 肖倩特聘研究员课题组 招聘硕士/博士(长期有效)
+1/19
双一流大学湘潭大学“化工过程模拟与强化”国家地方联合工程研究中心招收各类博士生
+1/18
中国地质大学(武汉)清洁能源与资源循环团队诚招科研助理/2026级硕士、博士生/博士后
+1/15
---大龄的未婚男找结婚女对象
+1/14
双一流湘潭大学化工学院招收化工、能源动力、工程热物理相关研究生调剂
+1/14
中南大学喻桂朋教授团队新增2026年博士研究生
+1/14
宁波东方理工大学余鹏课题组2026年上海交大联培博士生、博士后、科研助理招聘
+1/8
多伦多大学 & Synetic Labs 招博士后/研究员
+1/6
南医大基础医学院招收调剂生(2生信1生物))
+1/5
宁波东方理工大学余鹏课题组2026年上海交大联培博士生、博士后、科研助理招聘
+1/3
环境工程专业 2026年硕士研究生调剂名额(有学硕1名,专硕2名)
+1/3
“机械、材料”招生
+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













回复此楼