毕业论文文献综述信息与计算科学椭圆型偏微分方程的求解及其应用一、 前言部分微积分产生以后，人们就开始把力学中的一些问题，归结为偏微分方程进行研究。

早在18世纪初，人们已经将弦线振动的问题归结为弦振动方程，并开始探讨了它的解法。

随后，人们又陆续了解了流体的运动、弹性体的平衡和振动、热传导、电磁相互作用、原子核和电子的相互作用、化学反应过程等等自然现象的基本规律，把它们写成偏微分方程的形式，并且求出了典型问题的解答，从而能通过实践，验证这些基本规律的正确性，显示了数学物理方程对于认识自然界基本规律的重要性。

有了基本规律，人们还要利用这些基本规律来研究复杂的自然现象和解决复杂的工程技术问题，这就需要求出数学物理方程中的许多特定问题的解答。

随着电子计算机的出现及计算技术的发展，即使是相当复杂的问题，也有可能计算出解得足够精确的数值来，这对于预测自然现象的变化（如天气预报）和进行各种工程设计（如机械强度的计算）都有着很重要的作用[1]。

许多复杂的自然现象，其运动规律、过程和状态都是通过微分方程这种数学形式来描述的。

当我们研究只有一个自变量的运动过程时出现的微分方程称为常微分方程。

当一个微分方程除了含有几个自变量和未知数外，还含有未知数的偏导数时，称为偏微分方程[2]-[6]。

在偏微分方程中，偏导数自然是不可缺少的。

例如： ()(),,u ua x y f x y x y∂∂+=∂∂ （1.1.1） 拉普拉斯方程22232220u u uu x y z∂∂∂∆=++=∂∂∂（1.1.2） 热传导方程()222,,u u a f x t u t x ∂∂=+∂∂（1.1.3） 波动方程()2222,,u a u f t x y t∂=∆+∂（1.1.4）等都是偏微分方程。

其中，u 为未知数，a 为常数，(),a x y 、f 为已知函数。

偏微分方程的一般形式为()112,,,,,,,,0n n x x F x x x u u u ⋅⋅⋅⋅⋅⋅⋅⋅⋅= （1.1.5） 其中：F 为已知函数；12,,,n x x x ⋅⋅⋅为自变量；u 是关于这些自变量的未知数。

应注意F 中必须含有未知函数u 的偏导数。

偏导数方程（1.1.5）中所含有偏导数的最高阶数为该偏微分方程的阶。

如（1.1.1）是一阶偏微分方程，方程（1.1.2）~（1.1.4）是二阶偏微分方程。

如果一个偏微分方程对于未知函数及其所有偏导数都是线性的，则称之为线性偏微分方程，否则称为非线性方程。

如（1.1.1）、（1.1.2）、（1.1.4）都是线性方程。

我们将主要研究二阶线性偏微分方程，因为它们在物理、力学和其它自然科学以及工程技术中经常出现，常称为数学物理方程。

n 个自变量的二阶线性偏微分方程的一般形式为,,11i j i nni jx x i x i j i au u b u cu f ==++=∑∑ （1.1.6）不失一般性，可以假设ij ji a a =，且ij a ，i b ，c 及f 是空间nR 中某区域Ω内的函数，如果方程（1.1.6）中的自由项0f ≡，则称方程为齐次方程，否则称为非齐次方程。

设方程（1.1.5）的阶数为m ，函数()12,,,n u u x x x =⋅⋅⋅在区域nR Ω⊂中具有m 阶连续偏导数，且代入方程（1.1.5）后成为恒等式，则称u 为区域Ω内方程(1.1.5)的一个解。

容易验证函数()2u x y =+，()sin v x y =-都是方程222220u uu x y∂∂∆=+=∂∂ （1.1.7）的一个解。

称方程（1.1.7）为二维拉普拉斯（Laplace ）方程或二维调和方程。

由复变函数理论知，任何一个解析函数()f z 的实部和虚部都是方程（1.1.7）的解。

考察两个自变量的二阶线性偏微分方程 111222122xx xy yy x y a u a u a u b u b u cu f+++++=，（1.1.8）其中ij a ，i b ，c ，f 都是x ，y 的连续可微实值函数，并且11a ，12a ，22a 不同时为零。

在你一点00(,)x y ∈Ω的一个领域内考察自变量变换(,)x y ξξ=，(,)x y ηη=. （1.1.9）假设它的Jacobi 行列式(,)(,)0(,)x yx x x y D J D x y ξξξηηη==≠，由隐函数存在定理知该变换是可逆的，即存在逆变换(,)x x ξη=，(,)y y ξη=。

直接计算，有x x x u u u ξηξη=+，y y y u u u ξηξη=+，…将其代入方程(1.1.8)，得 *******111222122a u a u a u b u b u c f ξξξηηηξη+++++=，（1.2.0）其中*c c =，*f f =，*ij a ，*i b 可以分别用ij a ，i b 以及ξ和η的各阶偏导数表示。

特别地*22111112222x y x y a a a a ξξξξ=++，*22221112222x yx y a a a a ηηηη=++.（1.2.1）希望选取一个变换（1.1.9），使方程（1.2.0）有比方程（1.1.8）更简单的形式。

注意到（1.2.1）式中的*11a 与*22a 有相同的形式，如果我们能够解出方程2211122220x x y y a a a ϕϕϕϕ++= （1.2.2）的两个线性无关的解1(,)x y ϕ，2(,)x y ϕ，那么取1(,)x y ξϕ=，2(,)x y ηϕ=，就能保证**11220a a =≡。

这样，（1.2.0）式就较（1.1.8）式大为化简。

现在考察这种选取的可能性[7]。

我们知道关于ϕ的一阶偏微分方程（1.2.2）的求解问题可以化为求下述常微分方程在(),x y 平面上的积分曲线问题：()()2211122220a dy a dxdy a dx -+=.（1.2.3）设1(,)x y c ϕ=是方程（1.2.3）的一族积分曲线，则1(,)z x y ϕ=就是方程（1.2.2）的一个解。

称方程（1.2.3）的积分曲线为方程（1.1.8）的特征线，方程（1.2.3）有时亦称为特征方程。

偏微分方程可根据它的数学特征分为三大类型，即抛物型、双曲型、椭圆型。

这三类偏微分方程描述了不同本质的物理现象，其应用是极其广泛的[8]。

我们可以看到，两个自变量的二阶线性方程通过自变量的可逆变换能够化成哪种标准形，要看二次型()22111222,2Q l m a l a lm a m =++的代数性质如何让来定，或者说，由于,l m 平面上的二次曲线(),1Q l m =的性质而定。

由于这个曲线可以是一个椭圆、一个双曲线或者一个抛物线，故我们相应地定义方程在一点的类型如下：若方程（1.1.8）中二阶偏导数项的系数111222,,a a a 在区域Ω中某点()00,x y 满足 12211220a a a ∆=->， 则称方程在()00,x y 为双曲型的；若在点()00,x y 满足 12211220a a a ∆=-=， 则称方程在点()00,x y 为抛物型；若在点()00,x y 满足 12211220a a a ∆=-<， 则称方程在点()00,x y 为椭圆型的[1]。

给定一个常微分方程，有通解和特解的概念。

对于偏微分方程也一样。

换句话说，为了完全确定一个物理状态，只有相应的偏微分方程是不够的，必须给出它的初始状态和边界状态，即给出外加的特定条件，这种特定条件称为定解条件。

描述初始时刻物理状态的定解条件称为初值条件或初始条件，描述边界上物理状态的条件称为边界条件或边值条件。

一个方程配上定解条件就构成定解问题[7]。

那么我们如何求解偏微分方程的定解问题呢？数学物理方程中有许多是线性方程，与其对应的已经给出很多求准确解的方法，如特征线法、分离变量法、格林函数法、积分变换法及复变函数法等[1]。

求解有界区域上的线性偏微分方程定解问题的基本方法——分离变量法。

分离变量法的理论基础是Fourier级数展开（一个函数按照某个具体的完备正交函数系展开）。

对于包含无限区域或半无限区域的偏微分方程的定解问题，经常采用积分变换法。

这种方法是通过函数变换，减少泛定方程中自变量的个数，把偏微分方程问题转化为常微分方程问题，使计算大为化简。

差分法是求解偏微分方程常用的数值解法。

它的基本原理是：首先，将问题离散化，用差商代替微商，将微分方程和定解条件都用代数方程来代替；然后，解这些代数方程构成的方程组，得到定界问题的近似解[7]-[10]。

二、主题部分众所周知，17世纪微积分创立后，常微分方程理论立刻就发展起来。

即应用常微分方程于几何与力学问题的全新的计算。

结果是在天体力学中不仅能得到并解释早先已经知晓的那些事实，而且得到了新的发现（例如，海王星的发现就是在对微分方程的分析的基础上作出的）。

开始研究偏微分方程要晚得多。

对在物理学中碰到的偏微分方程的研究在18世纪中叶导致了分析学的一个新分支——数学物理方程的建立。

J.达浪贝尔（1717-1783）、L.欧拉（1707-1783）、D.伯努利（1700-1782）、J.拉格朗日（1736-1813）、P.拉普拉斯（1749-1827）、S.泊松（1781-1840）、J.傅里叶（1768-1830）等人的工作为这一科学分支奠定了基础。

他们在考察具体地数学物理问题中，所提出的思想与方法，竟适用于众多类型的微分方程，成为19世纪末偏微分方程一般理论发展的基础[5]。

早期建立的数学物理方程有根据牛顿引力理论而推导出的描述引力势的拉普拉斯方程和泊松方程。

对于建立的数学物理方程，需要作出各种附有具体条件而构成典型问题的解，然后根据实际测量结果来检验和修正相应的物理理论。

通过求解数理方程，使人们对自然现象获得更深刻的认识，并能预见新的现象。

随着现代科学和技术的进步，将会不断涌现新的数学物理方程，而其产生和应用的范围已经并且更多地超出了传统的物理学、力学、天文学等领域。

例如，在化学、生命科学、经济学等自然科学和社会科学各个领域，以及在资源勘探与开发、大型建筑与水利工程、金属冶炼工程、通信工程、新能源开发、大气物理、气象预报、航天工程、医疗诊断与材料无损探伤、遗传工程等广泛的工程技术各个领域都涉及到数学物理方程的理论及其重要应[4][5]。

数值天气预报、大型水坝应力分析等许多例子，说明数值求解偏微分方程在各门学科和工程中的应用，解偏微分方程已经成为科学与工程计算的核心内容，包括一些大型的计算和很多已经成为常规的计算。

为什么它在当代能发挥这样大的作用呢？第一是计算机本身有了很大的发展；第二是数值求解方程的计算方法也有了很大的发展，这两者对人们计算能力的发展都是十分重要的[9]。