偏微分方程的MATLAB解法1.docVIP

偏微分方程的MATLAB解法 　　第37卷第23期·学术第37卷第3No.2期Vol.37No.3 　　张向利：户县拖拉机及驾驶员年度检审现状及对策HUNANAGRICULTURALMACHINERY 　　湖南农机 　　2010年35月 　　May.2010ar.2010 　　偏微分方程的ＭＡＴＬＡＢ解法 　　李 　　明 　　（长江大学信息与数学学院，湖北荆州434023） 　　摘要：一般来说，在实际中偏微分方程的通解是不容易求出的，用定解条件确定函数更是比较困难。本文在运 　　列出两种方法：pdepe函数和PDE工具箱，并应用实例展现出两种方法的实现过程。结用MATLAB解偏微分方程时， 　　MATLAB对解偏微分方程带来极大的方便，并且在此基础上可以解决更多更复杂的问题。果表明：关键词：偏微分方程；MATLAB；PEDToolbox中图分类号：TE121文献标识码：A文章编号：1007-8320（2010）03-0089-02 　　MATLABSolutionpartialdifferentialequations 　　LIMing 　　（SchoolofMathematics，YangtzeUniversity，Jingzhou，Hubei434023，China） 　　Abstract：Generallyspeaking，inpractice，thegeneralsolutionofpartialdifferentialequationsisnoteasytofind，andusetheboundaryconditionstodeterminethefunctionismoredifficult.ThisarticlusingMATLABforsolvingpartialdifferentialequations，twomethodsarelisted：pdepefunctionandthePDEToolbox.Onthebasisofdescribingthetwomethodsindetail，thisarticleusesexamplestoshowtherealizationintwoways.Theresultsshowedthat:MATLABbringsgreatconvenienceonthesolutionsofpartialdifferentialequations，onthisbasis，itcansolvemorecomplexissues. 　　Keywords：partialdifferentialequations；MATLAB；PEDToolbox 　　1前言 　　偏微分方程定解问题有着广泛的应用背景。人们用偏微分方程来描述、解释或者预见各种自然现象，并用于科学和工对于广大应用工作者来说，从偏程技术的各个领域［1］。然而， 　　微分方程模型出发，使用有限元法或有限差分法求解都要耗费很大的工作量，才能得到数值解。现在，MATLABPDE准确的求解过程。Toolbox已实现对于空间二维问题高速、 　　输出参数：sol：是一个三维数组，sol（：，：，i）表示ui的解，换句话说uk对应x（i）和（tj）时的解为sol（i，j，k），通过sol，我（）直接计算某个点的函数值。们可以使用pdeval 　　实例讲解 　　例：试求解下面的偏微分 　　?? 　　???????????????? 　　坠u=0.024坠2u-F（u1-u2）坠t坠x2坠u2=0.170坠2u2-F（u1-u2）（1） 　　2Matlab解偏微分方程 　　MATLAB提供了两种方法［2］解决PDE问题，一是pdepe函数［3］，它可以求解一般的PDEs，具有较大的通用性，但只支持命令行形式调用。二是PDE工具箱，可以求解特殊PDE 　　［4］ 　　F（x）=e5.73x-e-11.46x，且满足初始条件u（0）＝1，u（其中，1x，2x，坠u21，0）＝0及边界条件坠u10，t）=0，u（t）=0，u（t）=1t）20，11， 　　=0。 　　解：（1）对照给出的偏微分方程，根据标注形式，则原方程可以改写为 　　1u 　　·*坠2=坠22u1 　　12 　　2 　　2222222222 　　问题，但有较大的局限性，比如只能求解二阶PDE问题，并且不能解决偏微分方程组。它提供了GUI界面，可以从繁杂的编程中解脱出来，同时还可以通过FileSaveAs直接生成M代码。2.1函数解法 　　pdepe函数介绍 　　它的调用格式为sol=pdepe（m，@pdefun，@pdeic，@pde－bc，x，t） 　　输入参数：@pdefun：是PDE的问题描述函数 　　@pdebc：是PDE的边界条件描述函数@pdeic：是PDE的初值条