| 查看: 302 | 回复: 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 ] |
回复此楼» 猜你喜欢
听说今天filecode变了
已经有30人回复
-
filecode
已经有14人回复
-
固定端突然变了，今天
已经有7人回复
-
filecode
已经有11人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有6人回复
-
基金中了
已经有4人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
-
求各位大神看下
已经有33人回复
-
咨询面上基金
已经有7人回复
-
影响面上的因素
已经有11人回复
找到一些相关的精华帖子，希望有用哦~
计算扩散力或能量的方程~~求大神进来看看，小弟新虫并无金币麻烦了
已经有6人回复
- 麻烦有机大神帮我看看这个怎么读？
已经有3人回复
- matlab程序问题，总是出错，麻烦大家看看问题在那里？
已经有6人回复
- matlab作业哪位大神麻烦给做一下,十万火急,周五之间要交
已经有3人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
50