用PDE Toolbox求解三类典型方程求解抛物型方程.docVIP

实验二 用PDE Toolbox求解三类典型方程 ----求解抛物型方程 一、实验目的 本次实验在掌握PDE Toolbox的功能和PDE图形用户界面的基础上，求解抛物型方程。 二、实验内容 1、解抛物型方程的一个例子 热传导方程：金属板的导热问题 考虑一个带有矩形孔的金属板上的热传导问题。板的左边保持在100℃，板的右边热量从板向环境空气定常流动，其他边及内孔边界保持绝缘。初始时板的温度为0℃ 的外边界定点坐标为（-0.5，-0.8），（0.5，-0.8），（0.5，0.8），（-0.5，0.8），内边界顶点坐标为（-00.5，-0.4），（00.5，-0.4），（00.5，0.4），（-00.5，0.4）。 使用GUI求解这一问题。在PDE Toolbox窗口的工具栏中选择Generic Scalar. 第一步：区域设置 单击工具，在窗口拖拉出一个矩形，双击矩形区域，在对话框中输入Left为-0.5，Bottom为-0.8，Width为1，Height为1.6，单击OK按钮，显示矩形区域R1。用同样方法做内孔R2，只要设置Left为-00.5，Bottom为-0.4，Width为0.1，Height为0.8即可。然后在栏中键入R1-R2。 第二步：设置边界条件 单击工具，图形边界变红，逐段双击边界，打开Boundary Condition对话框，输入边界条件。对于同一类型的边界，可按shift键，将多个边界同时选择，统一设置边界条件。本题在左边界上选择Dirichlet条件，输入h 为1，r为100；右边界上，选择Neumann条件，输入g为-1，q为0；其他边界上，选择Neumann条件，输入g为0，q为0 。 第三步：设置方程类型 单击PDE菜单中PDE Specification…选项，打开PDE Specification对话框，选择方程类型。本题单击Parabolic(抛物型)，输入d为1,c为1，a为0，f为0，然后单击OK按钮。 第四步：网格剖分 单击工具，或者单击Mesh菜单中Initialize Mesh选项，可进行初始网格剖分，如果需要网格加密，再单击，或者单击Mesh菜单中Refine Mesh选项。 第五步：初值和误差设置 单击Solve菜单中Parameters…选项，打开Solve Parameters…对话框，输入Time为0：5，u(t0)为0 ，Relative tolerance为0.01，Absolute tolerance为0.001，然后单击OK按钮。 第六步：解的图形 单击Plot菜单中Parameters…选项，出现Plot Selection对话框，从中可以选择Color, Contour, Arrows, Deformed mesh, Height(3-D plot)，还可以设置等值线的数目等。本例中选择Color, Contour和Arrows三项，然后单击Plot按钮，窗口显示出时的解的彩色图形。 第七步：数值解的输出 单击Solve菜单中Export Solution选项，在打开的Export对话框中输入u，单击OK按钮确定。再在MATLAB命令窗口键入u，按回车键确定，即显示按节点编号排列的数值解 u u = 0 79.6291 95.3774 98.0380 98.5706 98.6395 0 79.6292 95.3774 98.0380 98.5705 98.6395 100.0000 100.0000 100.0000 100.0000 100.0000 100.0000 0 83.3577 96.3858 98.5869 99.0275 99.0845 0 83.3653 96.3874 98.5875 99.0279 99.0849 0 89.8188 97.8111 99.1613 99.4315 99.4664 0 89.8118 97.8096 99.1607 99.4310 99.4660 100.0000 100.0000 100.0000 100.0000 100.0000 100.0000 0 79.3324 95.3139 98.0139 98.5544 98.6243 0 78.5592 95.1462 97.9486 98.5096 98.5