偏微分方程实验.docVIP

预备实验 偏微分方程工具箱的功能和使用 一、实验目的 PDE Toolbox提供了研究和求解空间二维偏微分方程问题的一个强大而又灵活实用的环境。本次实验了解PDE Toolbox的功能和PDE图形用户界面。 PDE Toolbox的功能包括：设置PDE 定解问题，即设置二维定解区域、边界条件以及方程的形式和系数；用有限元法求解PDE，即网格的生成、方程的离散以及求出数值解；解得可视化。 二、实验内容 1、解偏微分方程的一个例子 解Poisson方程：，边界条件为齐次dirichlet类型。 第一步：启动MATLAB，键入pdetool，按回车键确定便可启动GUI(图形用户界面)，然后在options菜单下选择grid命令，打开栅格。栅格的使用，能使用户容易确定所绘图形的大小。 第二步：分布完成平面几何造型：R1-C1-E1+R2+C2。用菜单或快捷工具分别画矩形R1，矩形R2，椭圆E1，圆C1，圆C2。画圆时，首先选中椭圆工具，按鼠标右键并拖动即可，或者在按Ctrl的同时，拖动鼠标也可绘制图。然后在set formula栏，进行编辑并用算术运算符将图形对象名称连接起来，或删除默认的表达式直接键入R1-C1-E1+R2+C2。 选择Boundary菜单中Boundary Mode 命令，进入边界模式。单击Boundary菜单中Remove All Subdomain Borders选项，去掉子域边界，如果想将几何信息和边界信息进行存储，应选择Boundary菜单中的 Export Decomposed Geometry, Boundary Cond’s…命令，将它们分别储存于g,b变量中，通过MATLAB形成 M文件。 第三步：选取边界，单击Boundary菜单中Specify Boundary Conditions…选项，打开Boundary Conditions对话框，输入边界条件，本例取缺省条件，即将全部边界设为齐次dirichlet条件，边界颜色显示为红色。 第四步：选择PDE菜单中PDE Mode命令，进入PDE模式。单击PDE菜单中PDE Specification…选项，打开PDE Specification对话框，设置方程类型。本例取缺省设置，类型为椭圆型，参数c,a,f分别为1，0，10 第五步：选择Mesh菜单中Initialize Mesh命令，进行网格剖分。 第六步：选择Mesh菜单中Refine Mesh命令，对网格加密。 第七步：选择Solve菜单中Solve PDE命令，解偏微分方程并显示图形解。 第八步：单击Plot菜单中Parameters选项，打开Plot Selection对话框，选中Color，Height(3-D plot)和Show mesh三项，然后单击Plot按钮，显示三维图形解。 第九步：如果要画等值线图和矢量场图，单击Plot菜单中Parameters选项，打开Plot Selection对话框，选中Contour和Arrows两项，然后单击Plot按钮，可显示解的等值线图和矢量场图。 2、PDE Toolbox菜单 File菜单 New ----------------新建一个几何结构实体模型，默认文件名为“Untitled”。 Open… ---------------从硬盘装载M文件。 Save --------------将在GUI内完成的成果储存到一个M文件中。 Save As… ---------将在GUI内完成的成果储存到另外一个M文件中。 Print… ----------将PDE工具箱完成的图形送到打印机进行拷。 Exit ------------退出PDE工具图形用户界面。 Edit菜单 Undo -------------在绘制多边形时退回到上一步操作 Cut ----------------将已选实体剪切到剪贴板上 Copy --------------将已选实体拷贝到剪贴板上 Paste --------------将剪贴板上的实体粘贴到当前几何结构实体模型中 Clear ----------------- 删除已选的实体 Select All ----------- 选择当前几何结构实体造型中的所有实体及其边界和子域 Options菜单 Grid ----------------------绘图时打开或关闭栅格 Grid Spacing… -----------调整栅格的大小 Snap - ----------------------打开或关闭捕捉栅格功能 Axes Limits…