| 查看: 2226 | 回复: 19 | |||
[交流]
【求助】feval的使用
|
|
现在在看《数值分析导论》韩渭敏译,379页中的边值问题, 在对这个偏微分方程的求解过程中用到一条语句 feval(f,(1:n-1)*h,(1:n-1)*h) 注释中说明f方程的右边的函数值,我不知道如何写这个f的m函数,有人可以指点一下啊 |
» 猜你喜欢
0703化学求调剂,各位老师看看我!!!
已经有3人回复
339求调剂
已经有5人回复
287求调剂
已经有8人回复
一志愿北京化工大学 070300 学硕 336分 求调剂
已经有8人回复
085602 化工专硕 338分 求调剂
已经有5人回复
0856,材料与化工321分求调剂
已经有4人回复
一志愿北京化工大学材料与化工 264分各科过A区国家线
已经有10人回复
266求调剂
已经有8人回复
265求调剂
已经有3人回复
286求调剂
已经有6人回复
» 抢金币啦!回帖就可以得到:
湖北大学陈苗苗教授食品安全研究团队诚招博士生
+1/184
资源与环境专硕
+1/73
青岛科技大学生物工程学院招调剂硕士研究生
+1/71
坐标广州,征女友
+2/58
2026 年山东建筑大学金属增材制造方向硕士研究生调剂公告
+1/40
浙江省人民医院博士后岗位招聘
+1/39
环境化学、地学、环境毒理方向,擅长污染物迁移转化降解、分析检测、计算模拟者优先
+1/35
哈尔滨工业大学(深圳)-何自开教授团队诚招化学工程与技术专业博士研究生
+1/30
福建师范大学环境与资源学院环境电化学课题组招收2026硕士(含调剂)
+1/29
北京理工大学-化学与化工学院-招收2026级博士生 [申请-考核制]
+1/27
赣南师范大学智能制造与未来能源学院2026年硕士调剂信息
+1/16
大连理工大学气体监测与环境催化团队诚聘副教授、助理教授及博士后3-4名
+1/6
同济大学环境学院 肖倩特聘研究员课题组 招聘硕士/博士(长期有效)
+1/6
湖北大学材料学招26级硕士研究生若干名
+1/5
双一流天津工业大学电信学院李鸿强教授招收2026年申请审核制博士
+1/5
【博士招生】广东工业大学国家优青课题组招收2026年环境方向博士生
+1/5
金属材料先进连接技术
+1/3
085602调剂
+1/2
上海交通大学化学化工学院张智涛课题组诚聘博士后
+1/1
上海交通大学环境智能感知实验室招聘博士后
+1/1
2楼2011-01-07 20:04:38
3楼2011-01-07 21:47:10
4楼2011-01-08 10:12:37
5楼2011-01-08 10:28:14
|
>> Poisson('f','g',16,1e-8,10000) Warning: Could not find an exact (case-sensitive) match for 'Poisson'. C:\Documents and Settings\lenovo\My Documents\MATLAB\poisson.m is a case-insensitive match and will be used instead. You can improve the performance of your code by using exact name matches and we therefore recommend that you update your usage accordingly. Alternatively, you can disable this warning using warning('off','MATLAB:dispatcher:InexactCaseMatch'). This warning will become an error in future releases. ??? Error using ==> mtimes Inner matrix dimensions must agree. Error in ==> f at 2 f=-2*pi^2*sin(pi*x)*sin(pi*y); Error in ==> poisson at 13 Fr(2:n,2:n)=h^2*feval(f,(1:n-1)*h,(1:n-1)*h); 错误的提示 |
6楼2011-01-08 10:51:59
7楼2011-01-08 10:57:32
8楼2011-01-08 11:04:41
|
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
10楼2011-01-08 11:06:14
11楼2011-01-08 11:06:47
12楼2011-01-08 11:09:33
13楼2011-01-08 11:34:09
14楼2011-01-08 11:48:22
15楼2011-01-08 11:53:22
16楼2011-01-08 12:03:07
17楼2011-01-08 13:37:10
18楼2011-01-08 14:13:18
19楼2011-01-08 14:32:43
20楼2011-01-08 15:21:26














回复此楼