偏微分方程求解课件.ppt

求解实例 第三十一页，编辑于星期四：十三点 三十分。 【解】由给定的PDE，可以得出d=1,c=1,a=2,f=10 第三十二页，编辑于星期四：十三点 三十分。 step1:点击工具栏的【PDE】按钮，如下输入PDE的参数，注意选择Hyperbolic 第三十三页，编辑于星期四：十三点 三十分。 step2:绘制求解域对坐标轴的操作可以在【Options】主菜单中操作，包括设置网格、坐标系范围等 (1)【Options】-Axis Limits设置如下 第三十四页，编辑于星期四：十三点 三十分。 第三十五页，编辑于星期四：十三点 三十分。 偏微分方程求解课件 第一页，编辑于星期四：十三点 三十分。 基础知识预习 微分方程的ＭＡＴＬＡＢ求解包含 １：常微分方程的ＭＡＴＬＡＢ求解（上　　　　　　　　　 节课已经讲过）这里不再赘述。 ２：偏微分方程的ＭＡＴＬＡＢ求解（本　　　 次教学内容） 第二页，编辑于星期四：十三点 三十分。 偏微分方程概念 偏微分方程（Partial Differential Equation，简称PDE）指含有未知函数及其偏导数的方程。描述自变量、未知函数及其偏导数之间的关系。 偏微分方程分为①线性偏微分方程式与②非线性偏微分方程式，常常有几个解而且涉及额外的边界条件。 常微分方程：在微分方程中，若自变量的个数只有一个的微分方程。 偏微分方程：自变量的个数有两个或两个以上的微分方程。 第三页，编辑于星期四：十三点 三十分。 求解偏微分方程的方法 求解偏微分方程的数值方法： 1. 有限元法（Finite Element Method, FEM）--- hp-FEM 2. 有限体积法（Finite Volume Method, FVM） 3. 有限差分法（Finite Difference Method, FDM）。 其它：广义有限元法（Generalized Finite Element Method, FFEM）、扩展有限元法（eXtended Finite Element Method, XFEM）、无网格有限元法（Meshfree Finite Element Method）、离散迦辽金有限元法（Discontinuous Galerkin Finite Element Method, DGFEM）等。 第四页，编辑于星期四：十三点 三十分。 MATLAB解偏微分方程 MATLAB提供了两种方法解决PDE 问题：①pdepe()函数，它可以求解一般的PDEs，具有较大的通用性，但只支持命令行形式调用。 ②PDE 工具箱，可以求解特殊PDE 问题，PDEtool 有较大的局限性，比如只能求解二阶PDE 问题，并且不能解决偏微分方程组，但是它提供了GUI界面，从繁杂的编程中解脱出来了，同时还可以通过File-Save As直接生成M代码 使用pdeval()直接计算某个点的函数值？？？ 第五页，编辑于星期四：十三点 三十分。 一般偏微分方程组(PDEs)的MATLAB求解 直接求解一般偏微分方程(组)，它的调用格式为sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t) 问题描述函数 初值条件 边界条件 输出参数 自变量 参数 第六页，编辑于星期四：十三点 三十分。 【输入参数】（1） @pdefun：是PDE 的问题描述函数，它必须换成下面的标准形式 PDE 就可以编写下面的入口函数 [c,f,s]=pdefun(x,t,u,du) m,x,t就是对应于(式1)中相关参数和自变量，du是u的一阶导数，由给定的输入变量即可表示出出c,f,s这三个函数 第七页，编辑于星期四：十三点 三十分。 【输入参数】（2） @pdeic：是PDE 的初值条件，必须化为下面的形式 我们使用下面的简单的函数来描述为u0=pdeic(x) 第八页，编辑于星期四：十三点 三十分。 【输入参数】（3） @pdebc：是PDE的边界条件描述函数，必须先化为下面的形式 于是边值条件可以编写下面函数描述为[pa,qa,pb,qb]=pdebc(x,t,u,du)其中a 表示下边界，b 表示下边界 第九页，编辑于星期四：十三点 三十分。 【输入参数】（4） m：就是对应于(式1)中相关参数 x,t：就是对应于(式1)中自变量 第十页，编辑于星期四：十三点 三十分。 【输出参数】 sol：是一个三维数组，sol(:,:,i)表示ui的解，换句话说uk对应x(i)和t(j)时的解为sol(i,j,k) 第十一页，编辑于星期四：十三点 三十分。 实例讲解（题目） 例： 第十二页，编辑于星期四：十三点 三