偏微分方程的解法课件.pptVIP

偏微分方程的解法 第一页，编辑于星期四：十三点 三十一分。 1、设置PDE的定解问题．即设置二维定解区域、边界条件以及方程的形式；和系数 2、用有限元法（FEM）求解PDE．即网格的生成、方程的离散以及求出数值解； 3、解的可视化． 主要讲述如何用MATLAB实现对偏微分方程的仿真求解．MATLAB的偏微分方程工具箱（PDE Toolbox）的出现，为偏微分方程的求解以及定性研究提供了捷径．主要步骤为： 第二页，编辑于星期四：十三点 三十一分。 PDEToolbox注意事项 只能解决二维模型，一维的扩成二维，三维的缩成二维，时间维不计算在内 公式类型，只能解决部分偏微分方程，由公式类型决定 边界条件两种，Dirichlet和Neumann 初始条件 第三页，编辑于星期四：十三点 三十一分。 先确定方程大类 第四页，编辑于星期四：十三点 三十一分。 Draw Mode 画图模式，先将处理的区域画出来，二维，方形，圆形，支持多边形，可以手动更改坐标，旋转rotate 例如，对于细杆导热，虽然是一维问题，可以将宽度y虚拟出来，对应于y的边界条件和初始条件按照题意制定 第五页，编辑于星期四：十三点 三十一分。 Boundary Mode ? 第六页，编辑于星期四：十三点 三十一分。 PDE Mode PDE Specification，确定偏微分方程类型共有四种： 第七页，编辑于星期四：十三点 三十一分。 椭圆形Elliptic ? 第八页，编辑于星期四：十三点 三十一分。 抛物型Parabolic ? 第九页，编辑于星期四：十三点 三十一分。 双曲型Hyperbolic ? 第十页，编辑于星期四：十三点 三十一分。 Mesh Mode 网格划分，细化 第十一页，编辑于星期四：十三点 三十一分。 Solve，Plot 如果有初始条件（与t有关），则在Solve的Parameters里有其设定，如果没有初始条件（与t无关），则不必设定Plot只是确定画图的参数，包括是否动画，是否3D，是否画出等温线，是否有箭头。。。 第十二页，编辑于星期四：十三点 三十一分。 Save As 保存成M-file，自动生成 第十三页，编辑于星期四：十三点 三十一分。 例： 解热传导方程 边界条件是齐次类型，定解区域自定。 【解】 第一步：启动MATLAB，键入命令pdetool并回车，就进入GUI．在Options菜单下选择Gid命令，打开栅格，栅格使用户容易确定所绘图形的大小． 第二步：选定定解区域本题为自定区域：自拟定解区域如图22 1所示：E1-E2+R1-E3．具体用快捷工具分别画椭圆E1、圆E2、矩形R1、圆E3．然后在Set formula栏中进行编辑并用算术运算符将图形对象名称连接起来(或删去默认的表达式，直接键入E1-E2+R1-E3) 第十四页，编辑于星期四：十三点 三十一分。 第十五页，编辑于星期四：十三点 三十一分。 第三步：选取边界 首先选择Boundary菜单中Boundary Mode命令，进入边界模式．然后单击Boundary菜单中Remove All Subdomain Borders选项。从而去掉子域边界，如图22 2．单击Boundary菜单中Specify Boundary Conditions选项，打开Boundary Conditions对话框，输入边界件．本例取默认条件，即将全部边界设为齐次Dirichlet条件，边界显示为红色． 如果想将几何与边界信息存储，可选Boundary菜单中的Export Decomposed Geometrv．Boundary Conds命令，将它们分别存储在g、b变量中，并通过MATLAB形成M文件． 第十六页，编辑于星期四：十三点 三十一分。 第十七页，编辑于星期四：十三点 三十一分。 第四步：设置方程类型 选择PDE菜单中PDE Mode命令，进入PDE模式， 再单击PDE菜单中PDE Secification选项，打开 PDE Secification对话框，设置方程类型． 本例取抛物型方程 故参数c,a,f,d,分别是l，0，10，1． 第五步：选择Mesh菜单中Initialize Mesh命令， 进行网格剖分， 选择Mesh菜单中Refine Mesh命令，使网格密集化， 如图22.3． 第十八页，编辑于星期四：十三点 三十一分。