声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2196|回复: 6

[综合讨论] 请问这forwdif()等函数是什么意思啊?

[复制链接]
发表于 2007-11-6 14:39 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
matlab程序里的forwdif()、crnich()、trisys()、dirich()等函数,在命令里找不到,不知道含义,请高手指教,谢谢!
回复
分享到:

使用道具 举报

发表于 2007-11-6 14:56 | 显示全部楼层

如果不是 matlab 自带的函数,请把程序下载完整。如果找不到源程序,就没有办法,你这样发问基本上没有人能够帮得上忙
 楼主| 发表于 2007-11-6 15:00 | 显示全部楼层
谢谢先。
看下面程序
function U=forwdif(f,c1,c2,a,b,c,n,m)

%Input - f=u(x,0) as a string 'f'
%      - c1=u(0,t) and c2=u(a,t)
%      - a and b right endpoints of [0,a] and [0,b]
%      - c the constant in the heat equation
%      - n and m number of grid points over [0,a] and [0,b]
%Output - U solution matrix; analogous to Table 10.4

%Initialize parameters and U

h=a/(n-1);
k=b/(m-1);
r=c^2*k/h^2;
s=1-2*r;
U=zeros(n,m);

%Boundary conditions

U(1,1:m)=c1;
U(n,1:m)=c2;

%Generate first row
U(2:n-1,1)=feval(f,h:h:(n-2)*h)';

%Generate remaining rows of U

for j=2:m
   for i=2:n-1
      U(i,j)=s*U(i,j-1)+r*(U(i-1,j-1)+U(i+1,j-1));
   end
end

U=U';
发表于 2007-11-6 15:04 | 显示全部楼层
原帖由 zhouwh 于 2007-11-6 15:00 发表
谢谢先。
看下面程序
function U=forwdif(f,c1,c2,a,b,c,n,m)

%Input - f=u(x,0) as a string 'f'
%      - c1=u(0,t) and c2=u(a,t)
%      - a and b right endpoints of [0,a] and [0,b]
%      -  ...

这是什么程序?啥都不说,换个角度思考,我给你一个莫名其妙的程序,你看得懂不?不要以为我们都是神仙,如果是对你有用的话,你就只能硬着头皮慢慢看,论坛不是养懒人的地方,“把程序一放就什么都不管了”并不适合咱们论坛。建议你还是自己摸索吧,没有人有那么多空余时间帮别人看程序,尤其是这种
发表于 2007-11-6 22:07 | 显示全部楼层

回复 #3 zhouwh 的帖子

自己一句一句读吧,好像没有什么特别难懂的语句
发表于 2007-11-6 22:44 | 显示全部楼层
数值方法  (MATLAB版)  (第三版) 好象有第四版了
第十章的
可以参考一下理论分析过程
以前读这本书的时候看到过该程序

评分

1

查看全部评分

 楼主| 发表于 2007-11-8 22:26 | 显示全部楼层
谢谢大家了!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-28 01:28 , Processed in 0.132766 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表