PDA

查看完整版本 : [求助]Matlab解一维偏微分方程(对流方程)


qingxin12
2008-12-18, 22:02
因为在写毕业论文,里面涉及到一个对流方程:du/dt+a*(du/dx)=0,要用数值的方法求他的近似解。现在问题是用matlab程序编写出来,然后作出图像,分析等等。

再次之前我从未接触过matlab,搜了很多的帖子还是不知其解,所以恳请这里的高人能帮助一下,非常特别的感谢。

问题: 一维对流方程 du/dt+a*(du/dx)=0
初值: u(0,x)=1, if x|<0.1
else u(0,x)=0
边值: u(t,0)=1, u(t,1)=0
x 的取值范围在0到1之间,其中时间步长取0.005,空间步长取0.01
求当t=0.5时的u随x的分布图。

请高手们一定要帮一下忙,我已经毫无头绪忙了2天了,因为在国外,所看的书也极其有限,再次感谢大家了!!