偏微分方程边值问题的数值解法.docVIP

PAGE PAGE 14 MACROBUTTON MTEditEquationSection2 SEQ MTEqn \r \h \* MERGEFORMAT SEQ MTSec \r 1 \h \* MERGEFORMAT SEQ MTChap \r 1 \h \* MERGEFORMAT 求解偏微分方程的边值问题 本实验学习使用MATLAB的图形用户命令pdetool来求解偏微分方程的边值问题。这个工具是用有限元方法来求解的，而且采用三角元。我们用内个例题来说明它的用法。 一、MATLAB支持的偏微分方程类型 考虑平面有界区域D上的二阶椭圆型PDE边值问题： MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 1. SEQ MTEqn \c \* Arabic \* MERGEFORMAT 1) 其中 未知函数为。它的边界条件分为三类： （1）Direchlet条件： MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 1. SEQ MTEqn \c \* Arabic \* MERGEFORMAT 2) （2）Neumann条件： MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 1. SEQ MTEqn \c \* Arabic \* MERGEFORMAT 3) （3）混合边界条件：在边界上部分为Direchlet条件，另外部分为Neumann条件。 其中是定义在边界的已知函数，另外也可以是一个2*2的函数矩阵，是沿边界的外法线的单位向量。 在使用pdetool时要向它提供这些已知参数。 二、例题 例题1 用pdetool求解 MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 1. SEQ MTEqn \c \* Arabic \* MERGEFORMAT 4) 解 ：首先在MATLAB 的工作命令行中键入pdetool ，按回牟键确定，于是出现PDE Toolbox 窗口，选Genenic Scalar模式． ( l ）画区域圆 单击椭圆工具按钮，大致在（0，0）位置单击鼠标右键，拖拉鼠标到适当位置松开。为了保证所绘制的圆是标准的单位园，在所绘园上双击，打开 Object Dialog 对话框，精确地输入圆心坐标X-center 为0 、Y-center 为0 及半径Radius 为l ，然后单击OK 按钮，这样单位画已画好． ( 2 ）设置边界条件单击工具边界模式按钮 ，图形边界变红，逐段双击边界，打开Boundary condition 对话框．输入边界条件．对于同一类型的边界，可以按Shift键，将多个边界同时选择，统一设边界条件．本题选择Dirichlet 条件，输入h 为1 , r 为0。，然后单击OK 按钮．也可以单击Boundary菜单中Spocify Boundary Condition …选项，打开Boundary Condition 对话框输入边界条件． ( 3 ）设置方程单击偏微分方程按钮 ，打开PDE Specification 对话框，选择方程类型· 本题选Ellintic（椭圆型）,输入c为1 , a 为O , f 为1 ，然后单击OK 按钮． ( 4 ）网格剖分 单击网格工具，或者单击Mesh 菜单中Initialize Mesh项，可进行初始网格剖分．这时在PDE Toolbox 窗口下方的状态栏内显示出初始网格的节点数和三角形单元数．本题节点数为144 个，三角形单元数为254 个(图?? )。如果要细化网格，单击细化工具，或者单击Mesh 菜单中Refine Mesh 选项，节点数成为541 个，三角形单元数为1016 个。 ( 5 ）解方程 单击解方程工具，或者单击S olve菜单中Solve PDE 选项，可求得方程数值解并用彩色图形显示。单击作图工具，或者单击Plot 菜单中Parameter…选项，出现Plot selection 对话框．从中选择于Height ( 3-D plo