线性算子方程求解方法的初步探究.docVIP

线性算子方程求解方法的初步探究 工业设计 张静 1014201056 随着科学技术的发展，越来越多的问题出现在工程计算领域中。许多学者用各种各样的方程来描述这些问题，形成了人们熟知的带有各式各样复杂定解条件的微分方程、积分方程和积分-微分方程。为了便于研究这些方程，人们在不同的空间中将这些方程写成统一的形式：算子方程。 1. 线性算子方程简介 设X,Y为Banach空间，A:→为有界线性算子，考虑线性算子方程 (1-1) 的数值求解大致分为两类： 第一，有限维逼近，包括投影法(伽略金(Galerkin) 型方法，最小二乘，配置法等)，差分法和有限元法。有限维逼近的核心是用有限维空间的近似方程来代替原方程，从而用近似方程的解来逼近原方程的解。对于方程(1-1)，我们通常按某种方法选取,的有限维子空间序列和，以及到的映射和到的映射，以 （1-2） 作为(1-1) 的近似方程。因为,,,构造方式多种多样，所以我们就得到了各种不同的数值方法。对于大多数近似求解方法，和是线性算子。 如果和是线性投影算子，就得到了一般的投影法。 当，，时，就是通常的伽略金法。特别的，在Hilbert空间中取为→的正交投影算子，就得到了正交伽略金法。 当A为微分算子或积分算子，所用的投影空间为样条函数空间时，就是有限元法。有限元法是求解微分方程的一种有效的数值方法，它是在古典的发的基础上采用分片多项式插值的技巧而形成的一种新的适应性很强的数值方法。 配置法。设，为某区域上的函数空间，，为的基底，：→为插值投影算子 则近似方程（1-2）为 它等价于 令，则上式变为： （1-3） 于是，为求解方程（1-2），只须由方程（1-3）解出，这就是配置法。简单的说，配置法就是取为插值投影算子，从而将求解算子方程（1-2）转化为求解线性方程组（1-3）这样一种求解方法。 若，为Hilbert空间，和分别为和的完全坐标序列，，为正交投影算子，，，则近似方程（1-2）等价于 （1-4） 若为线性算子，，则为最小二乘法，因为此时（1-4）与求使得极小的问题等价。 查分法采用各种各样的方法把方程定解问题离散化。它包含两个主要问题：一是差分格式的选取，不同差分格式的收敛性、收敛速度、稳定性和精度都不相同。如何构造高精度、稳定、收敛速度快的格式是差分法的一个重要课题；二是区域的离散化问题，不同的区域离散化方法对方程的求解，解的精度，稳定性和收敛性都有不同的影响。 第二，逐次逼近法，各种类型的迭代法属于其中。迭代法的常见形式是 迭代格式的收敛性是指：，这里为方程（1-1）的真解。用迭代 法求方程解的关键问题是：如何正确选取迭代格式才能保证迭代的收敛性和 收敛速度。 2. 线性算子方程求解方法 求解线性算于方程 常要求解唯一，或存在且连续等。可是，这样的条件在实际问题求解过程中经常无法得到验证。在此探讨一种求解方法，解决上述问题。 2.1 精确解 由空间的可分性及格兰姆—施米特正交化方法，给出空间新的标准正交基，且。 引理1.1 可分Hilbert空间H上的连续线性算子方程（1-1）等价于 （1-1） 证明：令是上完全正交基，其中，于是有 （1.1）成立，充分性得证。有（1.1），有 。 由的完全性，可知方程（1-1）成立。等价性得证。 令，，表示空间的正交补空间，,,，是从H到上的投影算子。则方程（1-1）也可以表示为一般的形式： （1-2） 表示，对于任意。显然，方程组（1-1）等价于（1-2）式。 引理1.2 若，则。 证明：因为，则，由于 因此。引理1.2成立。 定理1.1 线性算子方程（1）解存在的充要条件是。且当时，令，方程（1）的解可以表示为 （1-3） 证明：若方程（1）存在解，由引理1.1及（1-2）式，有 （1-4） （1-5） 假若，由（1.5）式，有 与（1-4）式矛盾，因此，充分性得证。 若，则。由引理1.2，，这样令 （1-6） 于是有 对于任意的，有 再由引理1.1，可知是方程（1-1）的解。从而定理1.1得证。 2.2 近似解 记,表示的正交补，，，是到上的投影。 推论