| 查看: 315 | 回复: 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 ] |
回复此楼» 猜你喜欢
面上函评意见出来了，像什么等级？
已经有15人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有11人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有7人回复
找到一些相关的精华帖子，希望有用哦~
计算扩散力或能量的方程~~求大神进来看看，小弟新虫并无金币麻烦了
已经有6人回复
- 麻烦有机大神帮我看看这个怎么读？
已经有3人回复
- matlab程序问题，总是出错，麻烦大家看看问题在那里？
已经有6人回复
- matlab作业哪位大神麻烦给做一下,十万火急,周五之间要交
已经有3人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
5
