偏微分方程讲义.pptVIP

双曲型方程的差分解法 可以使用下面的差分进行近似， 上面的格式都是显格式。当然，也可以采用不同的格式组合得到其他的格式。这里不再讲解。 2.3偏微分方程的Matlab求解工具 PDEPE 命令的使用方法及例子 PDEPE（ initial-boundary value problems for parabolic-elliptic PDEs in 1-D ）是求解抛物线和椭圆型一维偏微分方程的命令。为求解这些类型的各种各样的问题，在求解方程的形式，初始条件以及边界条件上，PDEPE使用下面的统一结构。方程形式 函数f用来描述“通量”，函数s用来描述“源”。 初始条件为 边界条件为（可以通过取不同的函数，对应全部的三类边界条件。） Matlab的命令格式为pdepe。 M是参数，可以取0，1，2. pdefun是上面的函数f，icfun描述初始条件，bcfun描述边界条件，xmesh和tspan是x和t方向上的划分，options是选择项。下面的两个例子是matlab帮助文件中自带的例子，这里做一个详细的讲解。 例子5，求偏微分方程 其中x在[0,1]之间，且满足以下初边值条件 步骤1： 将欲求解的偏微分方程改写成如下的标准形式。 步骤 2: 编写偏微分方程的系数向量函数。 function [c,f,s]=ex20_1pdefun(x,t,u,dudx) c=pi^2;f=dudx;s=0; 步骤 3: 编写起始值条件。 function u0=ex20_1ic(x) u0=sin(pi*x); 步骤 4: 编写边界条件。在编写之前，先将边界条件改写成标准形式，找出相对应的函数pa()和pb()，然后写出MATLAB 的边界条件函数，例如，原边界条件可写成 因而，边界条件函数可编写成 function [pa,qa,pb,qb]=ex20_1bc(xa,ua,xb,ub,t) pa=ua;qa=0; pb=pi*exp(-t);qb=1; 步骤 5: 取点。例如 x=linspace(0,1,20); %x 取 20 点 t=linspace(0,2,5); %时间取5 点输出 步骤 6: 利用 pdepe求解。 m=0; %依步骤1 之结果 sol=pdepe(m,@ex20_1pdefun,@ex20_1ic,@ex20_1bc,x,t); %这里 sol实际上是二维矩阵 步骤 7 显示结果。 surf(x,t,sol) title(pde 数值解),xlabel(位置),ylabel(时间 ),zlabel(u) 例6, 试解以下联立的偏微分方程系统 步骤1：改写偏微分方程为标准形式 其中m=0和 初始条件:在x=0处 在x=1处 步骤 2：编写偏微分方程的系数向量函数。 function [c,f,s]=ex20_2pdefun(x,t,u,dudx) c=[1 1];f=[0.024 0.170].*dudx; y=u(1)-u(2);F=exp(5.73*y)-exp(-11.47*y); s=[-F F]; 步骤 3：编写初始条件函数 function u0=ex20_2ic(x) u0=[1 0]; 步骤 4：编写边界条件函数 function [pa,qa,pb,qb]=ex20_2bc(xa,ua,xb,ub,t) pa=[0 ua(2)];qa=[1 0]; pb=[ub(1)-1 0];qb=[0 1]; 步骤 5：取点 m=0; x=[0 0.005 0.01 0.05 0.1 0.2 0.5 0.7 0.9 0.95 0.99 0.995 1]; t=[0 0.005 0.01 0.05 0.1 0.5 1 1.5 2]; %************************************* %利用pdepe求解 %************************************* -382- sol=pdepe(m,@ex20_2pdefun,@ex20_2ic,@ex20_2bc,x,t); u1=sol(:,:,1); %第一个状态之数值解输出 u2=sol(:,:,2); %第二个状态之数值解输出 %************************************* %绘图输出 %************************************* figure(1), surf(x,t,u1) title(u1 之数值解), xlabel(x), ylabel(t) % f