MATLAB求解PDE问题分析.docVIP

MATLAB求解PDE问题（1）——概述、例子(转) (2011-07-20 16:48:45) oolbox 提供利用有限元方法求解偏微分方程的GUI以及相应的命令行函数。利用该工具箱可以求解椭圆型方程、抛物型方程、双曲型方程、特征值方程以及非线性方程。PDE Toolbox的功能非常强大，网上有许多利用PDE Toolbox解决各种物理问题的论文，还有专门介绍工具箱的参考书。 网上的例子虽然很多，但是大部分是介绍PDE工具箱自带的一些例子，这些例子中解的区域，边界条件是PDE工具箱已经编写好的，直接调用就可以。对于该如何自己设定求解区域及边界条件，却很少有人涉及。网上搜索发现只有刘平在博客中详细介绍过求解区域的设定。下面以一个椭圆型方程的例子来详细说明求解的各个步骤，希望对大家能有所帮助。 设要求如下形式的椭圆方程的解： 按照PDE的要求，将方程化为标准形式 求解后的图像如下，第一幅图是解的图像，第二幅是计算误差。从第二幅图可以看到，计算的最大误差是10-3方量级。 通过这个例子我们可以基本掌握PDE求解偏微分方程的步骤和方法，后面我将详细介绍如何设置区域及边界条件。掌握了区域和边界条件的设定，就可以轻松求解遇到的偏微分方程了。图后是附带的matlab命令以及注释，并提供m文件附件下载，下载后解压即可。希望能对大家有所帮助。 下面是编写的求解上述方程的matlab