今天教大家快速入门matlab的
PDE
工具箱。在学习数学物理方程的时候,小伙伴们肯定被各式各样的偏微分方程
弄得十分头疼。matlab提供了一个非常厉害的工具箱,能够在一个二维有界
区域求解三类
偏微分方程的数值解,不仅为偏微分方程的学习提供了直观物理感受,也在数学建模中发挥巨大作用。
1 问题选择
在一个二维的有界区域Ω上,matlab的PDE工具箱主要使用``有限元`方法解决下面四类问题:
在求解一个偏微分方程的时候,我们首先要做的是确定方程所属的类型。
2 边界条件选择
•在菜单栏中选择矩形框;•然后在二维平面画一个矩形,可以看到出现了四条带小箭头的线段;•鼠标左键双击其中其中一条小线段,就会跳出选择狄利克雷边界条件和诺伊曼边界条件参数的图形框;
•边界条件形式如下,只需根据需要更改参数。
3 菜单按钮和简单使用
3.1 matlab pde工具箱打开方式
有两种:1.在命令行输入:pdetool
2.在APP工具栏数学、统计和优化中选择Pde Modeler
3.2 Pde Modeler工具栏使用
求解整个pde的过程可以概括为从左到右依次使用如下工具栏。
•1) 利用绘图(draw)模式,定义需要求解的问题的空间范围(domain)Ω 。•2) 利用 boundary 模式,指定边界条件。•3) 利用 PDE 模式,指定 PDE 系数,即输入 c,a,f 和 d 等 PDE 模式中的系数。•4) 在 mesh 模式下,产生 mesh 点,以便将原问题离散化。•5) 在 solve 模式下,求解。•6) 最后,在 Plot 模式下,显示答案。
3.3 简单的小例子入手
假定R0 = 0.25,求解的时间t =0.025
•命令行输入pdetool打开工具箱•使用矩形工具(第一个)拖动绘制矩形。双击画出的矩形填写Left=0,Bottom=0,Width=1,Height =1
•点击∂Ω,边界以红色显示(红蓝绿依次表示狄利克雷、诺依曼和混合边界条件)。使用shift和鼠标点击,选中四条边界,设置为诺依曼边界条件,根据问题,填写q = 0 , g = 0,点选OK。
•点击PDE按钮,选择抛物方程,设置参数如下:
•点击三角形按钮,以及后面的加密按钮,对区域进行三角形剖分。
•选择Solve下拉菜单中的参数,根据问题设置参数如下:
•点击等号进行求解
•点击带图案的那个按钮,选中Color选项,以及 Height(3-D plot)选项和Animation选项,其他想勾的也可以勾,比如show mesh。不同组合可以呈现出不同的图,根据需要选择。
更多偏微分方程的自学资料请后台回复
PDE
作者:li | 编辑:图图