偏微分方程Matlab数值解法（补充4）.docVIP

PAGE PAGE 5 偏微分方程Matlab数值解法（补充4） Matlab可以求解一般的偏微分方程，也可以利用偏微分方程工具箱中给出的函数求解一些偏微分方程。 1 偏微分方程组求解 Matlab语言提供了pdepe()函数，可以直接求解偏微分方程 （4.1） 这样，偏微分方程可以编写为以下函数的描述，其入口为 其中：pdefun为函数名。由给定输入变量可计算出这三个函数。 边界条件可以用下面的函数描述 （4.2） 这样的边值函数可以写为Matlab函数 初始条件数学描述为，编写一个简单的函数即可 还可以选择和的向量，再加上描述这些函数，就可以用pdepe（）函数求解次偏微分方程，需要用如下格式求解 试求下列偏微分方程 其中：，且满足初始条件，及边界条件 解：对照给出的偏微分方程和（4.1），可将原方程改写为 可知，且 编写下面的Matlab函数 function [c,f,s]=c7mpde(x,t,u,du) c=[1;1];y=u(1)-u(2);F=exp(5.73*y)-exp(-11.46*y);s=F*[-1;1]; f=[0.024*du(1);0.17*du(2)]; 套用（4.2）中的边界条件，可以写出如下的边值方程 左边界，右边界 编写下面的Matlab函数 function [pa,qa,pb,qb]=c7mpbc(xa,ua,xb,ub,t) pa=[0;ua(2)];qa=[1;0];pb=[ub(1)-1;0];qb=[0,1]; 另外，描述初值的函数 function u0=c7mpic(x) u0=[1;0]; 有了这三个函数，选定和向量，则可以由下面的程序直接求此微分方程，得出解和。 x=0:0.05:1;t=0:0.05:2;m=0; sol=pdepe(m,@c7mpde,@c7mpic,@c7mpbc,x,t); surf(x,t,sol(:,:,1)) figure;surf(x,t,sol(:,:,2)) 2 二阶偏微分方程的属性描述 2.1 椭圆型偏微分方程 椭圆型偏微分方程的一般形式为 其中：若，为的梯度，则其定义为 散度的定义为 这样，可以更明确地表示为 若为常数，则进一步化简为 其中，又称为Laplace算子。这样椭圆型偏微分方程可以简单地写为 2.2 抛物型偏微分方程 抛物型偏微分方程的一般形式为 根据上面叙述，若为常数，则该方程可以更简单地写为 2.3双曲型偏微分方程 双曲型偏微分方程的一般形式为 若为常数，则可以将该方程简化为 三类方程的直接的区别在于对的导数的阶次。 若对没有求导，可以理解为其值为常数，故称为椭圆型的。 若取对时间的一阶导数，则与对的二阶导数直接构成了抛物线关系，故称为抛物型偏微分方程。 若取对时间的二阶导数，称其为双曲型偏微分方程。 2.4 特征值型偏微分方程 特征值型偏微分方程为 对常数该方程可以化简为 该方程是椭圆型偏微分方程的特例。 3 偏微分方程求解界面应用举例 3.1 偏微分方程求解程序概述（已讲） 3.2 偏微分方程求解区域绘制（已讲） 3.3 偏微分方程边界条件描述 求解边界在偏微分方程界面下用。在偏微分方程工具箱支持包括Dirichlet条件和Neumann条件。 Dirichlet条件。 一般描述为 其中：表示求解区域的边界。假设在边界上满足该方程，则只需要给出和函数即可，这两个参数可以为常数，也可以为的函数，甚至可以是的函数，为了方便起见，一般可以令。 Neumann条件。 其扩展形式为 其中：为向量法向量的偏导数。 选择Boundary→Specify Boundary Condition菜单，设置即可。 3.4 偏微分方程求解举例 【10.2】试求双曲型偏微分方程 解：（双曲型方程的标准形式为 因此，。 区域绘制：（R1+E1+E2)-E3 网格绘制（网格加密） Boundary→Specify Boundary Condition菜单，设置。 PDE→PDE Specification→Parabolic选项，输入 也可修改边界条件为（边界上所有的值为5）。 Plot→Plot Solution Paramters 3.5 时变动画显示 偏微分方程默认的时间向量，以所给的是最终时刻时的解。从双曲型偏微分方程看，方程的解应该是时间t的函数，所以应该用动画来显示出来。 （以例2为例） 可由Solve→Parameters菜单引出的对话框设置时间向量，在time栏内填入，这样求解即为这段时间的解。 选择其中的动画选项：Animation 单击Options设置动画的播放速度，默认为6帧/秒，这样就可获得动画。 可以用Plot→Export Movie菜单将动画输出到Matlab工作空间，例