求解抛物线型偏微分方程matlab程序

MATLAB提供两种解决PDE问题:
pdepe()函数求解般PDEs据用较通用性支持命令行形式调用
二PDE工具箱求解特殊PDE问题PDEtool较局限性比能求解二阶PDE问题并且能解决偏微程组提供GUI界面繁杂编程解脱同通File->Save As直接M代码

MATLAB语言提供pdepe()函数直接求解般偏微程(组)调用格式

sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t)

【输入参数】
@pdefun:PDE问题描述函数必须换面标准形式

PDE编写面入口函数
[c,f,s]=pdefun(x,t,u,du)
m,x,t应于(式1)相关参数duu阶导数由给定输入变量即表示c,f,s三函数

@pdebc:PDE边界条件描述函数必须先化面形式

于边值条件编写面函数描述
[pa,qa,pb,qb]=pdebc(x,t,u,du)
其a表示边界b表示边界

@pdeic:PDE初值条件必须化面形式

股我使用面简单函数描述
u0=pdeic(x)

m,x,t:应于(式1)相关参数

【输参数】
sol:三维数组sol(:,:,i)表示ui解换句说uk应x(i)t(j)解sol(i,j,k)

通sol我使用pdeval()直接计算某点函数值