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