偏微分方程解法.ppt

《偏微分方程解法.ppt》由会员分享，可在线阅读，更多相关《偏微分方程解法.ppt（26页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、关于偏微分方程的解法第一张，PPT共二十六页，创作于2022年6月 1、设置PDE的定解问题即设置二维定解区域、边界条件以及方程的形式；和系数 2、用有限元法（FEM）求解PDE即网格的生成、方程的离散以及求出数值解；3、解的可视化主要讲述如何用主要讲述如何用MATLABMATLAB实现对偏微分方程的仿真求解实现对偏微分方程的仿真求解MATLABMATLAB的偏微分方程工具箱（的偏微分方程工具箱（PDE ToolboxPDE Toolbox）的出现，）的出现，为偏微分方程的求解以及定性研究提供了捷径主要步为偏微分方程的求解以及定性研究提供了捷径主要步骤为：骤为：第二张，PPT共二十六页，创作于

2、2022年6月PDEToolbox注意事项只能解决二维模型，一维的扩成二维，三维的只能解决二维模型，一维的扩成二维，三维的缩成二维，时间维不计算在内缩成二维，时间维不计算在内公式类型公式类型，只能解决部分偏微分方程，由公式类，只能解决部分偏微分方程，由公式类型决定型决定边界条件两种，边界条件两种，DirichletDirichlet和和NeumannNeumann初始条件初始条件第三张，PPT共二十六页，创作于2022年6月先确定方程大类第四张，PPT共二十六页，创作于2022年6月Draw Mode画图模式，先将处理的区域画出来，画图模式，先将处理的区域画出来，二维二维，方，方形，圆形，支持

3、多边形，可以手动更改坐标，形，圆形，支持多边形，可以手动更改坐标，旋转旋转rotaterotate例如，对于细杆导热，虽然是一维问题，可以例如，对于细杆导热，虽然是一维问题，可以将宽度将宽度y y虚拟出来，对应于虚拟出来，对应于y y的边界条件和初始的边界条件和初始条件按照题意制定条件按照题意制定第五张，PPT共二十六页，创作于2022年6月Boundary Mode第六张，PPT共二十六页，创作于2022年6月PDE ModePDE Specification，确定偏微分方程类型共有四种：第七张，PPT共二十六页，创作于2022年6月椭圆形Elliptic第八张，PPT共二十六页，创作于20

4、22年6月抛物型Parabolic第九张，PPT共二十六页，创作于2022年6月双曲型Hyperbolic第十张，PPT共二十六页，创作于2022年6月Mesh Mode网格划分，细化网格划分，细化第十一张，PPT共二十六页，创作于2022年6月Solve，Plot如果有初始条件（与如果有初始条件（与t t有关），则在有关），则在SolveSolve的的ParametersParameters里有其设定，如果没有初始条件里有其设定，如果没有初始条件（与（与t t无关），则不必设定无关），则不必设定PlotPlot只是确定画图只是确定画图的参数，包括是否动画，是否的参数，包括是否动画，是否3D3

5、D，是否画出，是否画出等温线，是否有箭头。等温线，是否有箭头。第十二张，PPT共二十六页，创作于2022年6月Save As保存成保存成M-fileM-file，自动生成，自动生成第十三张，PPT共二十六页，创作于2022年6月 例：解热传导方程 边界条件是齐次类型，定解区域自定。【解解】第一步：启动MATLAB，键入命令pdetool并回车，就进入GUI在Options菜单下选择Gid命令，打开栅格，栅格使用户容易确定所绘图形的大小第二步：选定定解区域本题为自定区域：自拟定解区域如图22 1所示：E1-E2+R1-E3具体用快捷工具分别画椭圆E1、圆E2、矩形R1、圆E3然后在Set for

6、mula栏中进行编辑并用算术运算符将图形对象名称连接起来(或删去默认的表达式，直接键入E1-E2+R1-E3)第十四张，PPT共二十六页，创作于2022年6月第十五张，PPT共二十六页，创作于2022年6月第三步：选取边界第三步：选取边界 首先选择Boundary菜单中Boundary Mode命令，进入边界模式然后单击Boundary菜单中Remove All Subdomain Borders选项。从而去掉子域边界，如图22 2单击Boundary菜单中Specify Boundary Conditions选项，打开Boundary Conditions对话框，输入边界件本例取默认条件，即

7、将全部边界设为齐次Dirichlet条件，边界显示为红色如果想将几何与边界信息存储，可选Boundary菜单中的Export Decomposed GeometrvBoundary Conds命令，将它们分别存储在g、b变量中，并通过MATLAB形成M文件第十六张，PPT共二十六页，创作于2022年6月第十七张，PPT共二十六页，创作于2022年6月第四步：设置方程类型第四步：设置方程类型选择PDE菜单中PDE Mode命令，进入PDE模式，再单击PDE菜单中PDE Secification选项，打开PDE Secification对话框，设置方程类型本例取抛物型方程故参数c,a,f,d,分别

8、是l，0，10，1第五步：第五步：选择Mesh菜单中Initialize Mesh命令，进行网格剖分，选择Mesh菜单中Refine Mesh命令，使网格密集化，如图22.3第十八张，PPT共二十六页，创作于2022年6月第十九张，PPT共二十六页，创作于2022年6月第六步：第六步：解偏微分方程并显示图形解解偏微分方程并显示图形解选择Solve菜单中Solve PDE命令，解偏微分方程并显示图形解，如图 2.4 所示第二十张，PPT共二十六页，创作于2022年6月第二十一张，PPT共二十六页，创作于2022年6月 第七步：第七步：单击Plot菜单中Parameter选项，打开Plot Sel

9、ection对话框，选中Color,Height(3D plot)和Show mesh三项.再单击Polt按钮，显示三维图形解，如图22.5所示.第二十二张，PPT共二十六页，创作于2022年6月第八步第八步：若要画等值线图和矢量场图，单击plot菜单中parameter 选项，在plot selection对话框中选中contour 和arrow两选项。然后单击plot按钮，可显示解的等值线图和矢量场图，如图2.6所示。图2.6 解的等值线图和矢量场图 第二十三张，PPT共二十六页，创作于2022年6月求解椭圆型方程的例子第二十四张，PPT共二十六页，创作于2022年6月求解抛物型方程的例子考虑一个带有矩形孔的金属板上的热传导问题。板的左边保持在100 C，板的右边热量从板向环境空气定常流动，其他边及内孔边界保持绝缘。初始 是板的温度为0 C，于是概括为如下定解问题;区域的边界顶点坐标为(-0.5,-0.8),(0.5,-0.8),(-0.5,0.8),(0.5,0.8)。内边界顶点坐标(-0.05,-0.4),(-0.05,0.4),(0.05,-0.4),(0.05,0.4)。第二十五张，PPT共二十六页，创作于2022年6月感感谢谢大大家家观观看看第二十六张，PPT共二十六页，创作于2022年6月2022/10/16