从写完上一篇常微分方程的求解到现在已经很长时间了,这周也一直忙于报到的各种事宜,无暇坐下来写些东西,趁着这个周末,终于完成了这个姊妹篇。
对于偏微分方程的求解,Matlab提供了两种工具。第一种是pdepe()函数,它的特点是通用性好,不受求解阶次的限制,不足之处是只支持命令行的格式;第二种是PDE工具箱,它的特点是提供了一个GUI界面,简洁易懂可视,可以从枯燥的编程中解脱出来,不足之处是使用有限制,只能求解二阶的PDE,且不支持偏微分方程组的求解。
(1)、首先,我们来介绍pdepe()函数的使用。
pdepe()函数的调用格式为:
sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t)
其中,@pdefun是PDE方程的函数描述,它必须写成下面这种固定的格式:
这样,偏微分方程可以编写下面的函数描述,其入口为
[c,f,s]=pdefun(x,t,u,ux)
其中,pdefun为函数名;m,x,t就是对应于标准格式中的相关参数。
@pdebc是PDE的边界条件描述函数,必须先化成下面的标准形式: