24小时热门版块排行榜

返回列表

chuckxian

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 33
帖子: 24
在线: 13.3小时
虫号: 2421002
注册: 2013-04-13
专业: 偏微分方程

[交流] 各位大神麻烦看看，究竟我的程序是对的还是，书上写的是是对的

本人菜鸟，正在学习matlab编程，初试编了个程序，始终与与书中所写的答案，在初值条件编写上存在差异，我始终觉得是书上写错了，但我的图像好像又不对，麻烦大家看看

其中u1=sint ,u2=cost, l=1,a^2=0.001,φ(x)=x,x的和t范围都是[0,1],都分为100份
Δx=Δt=0.01；题目和最后求解的差分格式都在图片中

我的编程如下
u=zeros(100,100);
for i=1:100
u(1,i)=sin(i/100); u(100,i)=cos(i/100); u(i,1)=i/100;我写的初值条件
书上写的初值条件为u(i,1)=sin(i/100); u(1100i)=cos(i/100); u(1,i)=i/100;
下面都程序都一样·····
end
a=0.001;
delax=0.01; delat=0.01;
r=a*a*delat/delax/delax;
for j=1:100
u(2:99,j+1)=(1-2*r)*u(2:99,j)+r*(u(1:98,j)+u(3:100,j));
pause(0.01);
end
surf(u)

IMG_20rrrr130420_164207.png

IMG_20rrrr130420_164207.png

IMG_201ee30420_170824.png

[ Last edited by chuckxian on 2013-4-20 at 17:11 ]

回复此楼

» 猜你喜欢

1楼 2013-04-20 17:05:33

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

相关版块跳转我要订阅楼主 chuckxian 的主题更新

返回列表