| 查看: 2244 | 回复: 19 | |||
[交流]
【求助】feval的使用
|
|
现在在看《数值分析导论》韩渭敏译,379页中的边值问题, 在对这个偏微分方程的求解过程中用到一条语句 feval(f,(1:n-1)*h,(1:n-1)*h) 注释中说明f方程的右边的函数值,我不知道如何写这个f的m函数,有人可以指点一下啊 |
» 猜你喜欢
考研调剂
已经有3人回复
317分 一志愿南理工材料工程 本科湖工大 求调剂
已经有9人回复
297求调剂
已经有10人回复
一志愿武汉理工,总分321,英一数二,求老师收留。
已经有10人回复
343求调剂
已经有5人回复
085601材料工程找调剂
已经有7人回复
求调剂,一志愿 南京航空航天大学 ,080500材料科学与工程学硕,总分289分
已经有6人回复
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研经历
已经有5人回复
340求调剂
已经有5人回复
一志愿北化求调剂
已经有4人回复
» 抢金币啦!回帖就可以得到:
【招聘】温州医科大学/附属医院/瓯江实验室李校堃院士团队诚聘博士后和科研人员
+1/192
汕头大学 海洋科学-生物学-生物医药大量接受调剂
+1/191
三峡大学生态学招硕士调剂
+2/160
中国水产科学研究院黄海水产研究所张旭志课题组招聘 柔性传感器研究 科研助理
+2/138
五邑大学池振国教授团队2026年大量接收化学化工材料背景考研生
+5/95
青岛农业大学海洋学院招调剂硕士研究生
+1/73
电力全国重点实验室双一流A类长江学者团队招2026年全日制博士1名(4月3日截止)
+2/68
坐标广州,征女友
+2/48
长江师范学院 材料与化工专业 招收调剂学生
+1/40
球差相关知识
+1/37
沈阳农业大学生物与医药招收调剂生!!!!!!!!!!
+1/35
澳大利亚昆士兰大学化工学院招收全奖博士生2名
+1/26
【博士招生】广东工业大学国家优青课题组招收2026年环境方向博士生
+1/7
中国地震局地质研究所+地球物理专业+接收调剂
+1/7
教育部长江学者和创新团队发展计划”入选团队招收 材料,化学与化工博士研究生
+1/5
中国地质大学(武汉)清洁能源与资源循环团队诚招科研助理/2026级硕士、博士生/博士后
+1/5
三峡大学材化学院磷硅精细化工课题组招收化学类、材料类、材料与化工类硕士
+1/4
招收从事大气复合污染成因领域研究生
+1/3
求调剂
+1/2
哈工大韩飞课题组招收3D微纳加工、光学工程、纳米化学等方向博士生及博后(长期有效)
+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














回复此楼