| 查看: 314 | 回复: 0 | |||
[交流]
各位大神麻烦看看，究竟我的程序是对的还是，书上写的是是对的
|
|
本人菜鸟，正在学习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 ] |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
找到一些相关的精华帖子，希望有用哦~
计算扩散力或能量的方程~~求大神进来看看，小弟新虫并无金币麻烦了
已经有6人回复
- 麻烦有机大神帮我看看这个怎么读？
已经有3人回复
- matlab程序问题，总是出错，麻烦大家看看问题在那里？
已经有6人回复
- matlab作业哪位大神麻烦给做一下,十万火急,周五之间要交
已经有3人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
30