matlab怎么求解偏微分方程
跪求MATLAB解抛物型偏微分方程的程序
如何用Matlab解偏微分方程组该方程组由两个抛物型偏微分方程组成
如何用Matlab解偏微分方程组该方程组由两个抛物型偏微分方程组成 这个没有自带的函数，需要把插分格式写出来以后自己编程。
急求！！！ 大学数学，用matlab解决问题，题目是一维抛物型偏微分方程差分解法 显式前向欧拉法源程序：function[u，x，t]=EF_Euler(A，xf，T，it0，bx0，bxf，M，N)解方程 A u_xx=u_t，0，0初值：u(x，0)=it0(x)边界条件：u(0，t)=bx0(t)，u(xf，t)=bxf(t)M：x 轴的等分段数N：t 轴的等分段数dx=xf/M；x=[0：M]*dx；dt=T/N；t=[0：N]'*dt；for i=1：M+1u(i，1)=it0(x(i))；endfor j=1：N+1u([1 M+1]，j)=[bx0(t(j))；bxf(t(j))]；endr=A*dt/dx/dx，r1=1-2*r；if(r>；0.5)disp('r>；0.5，unstability')；endfor j=1：Nfor i=2：Mu(i，j+1)=r*(u(i+1，j)+u(i-1，j))+r1*u(i，j)；(9.2.3)endendu=u'；在MATLAB中编写脚本文件：A=0.5；方程系数it0=inline('sin(pi*x)'，'x')；初始条件bx0=inline('0')；bxf=inline('0')；边界条件xf=2；M=80；T=0.1；N=100；[u1，x，t]=EF_Euler(A，xf，T，it0，bx0，bxf，M，N)；figure(1)，clf，mesh(u1)xlabel('x')ylabel('t')zlabel('U')title('r>；0.5')M=50；[u1，x，t]=EF_Euler(A，xf，T，it0，bx0，bxf，M，N)；figure(2)，clf，mesh(u1)xlabel('x')ylabel('t')zlabel('U')title('r)隐式后向欧拉法源程序：function[u，x，t]=IB_Euler(A，xf，T，it0，bx0，bxf，M，N)解方程 A1 u_xx=u_t，0，0初值：u(x，0)=it0(x)边界条件：u(0，t)=bx0(t)，u(xf，t)=bxf(t)M：x 轴的。
求解抛物线型偏微分方程matlab程序 MATLAB提供两种解决PDE问题：pdepe()函数求解般PDEs据用较通用性支持命令行形式调用二PDE工具箱求解特e69da5e6ba903231313335323631343130323136353331333365666232殊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()直接计算某点函数值
