二阶常微分方程边值问题的数值解-外文翻译.docVIP

毕业设计（论文）外文翻译 题 目：二阶常微分方程边值问题的数值解 系 ：数理与信息系 学生姓名：施逸斌 学 号：081315123 专　　业：信息与计算科学 班 级：C08信息 指导老师： 黄红英 起止日期：2012.3.17-2011.3.27 2012年 3月 27 日 11.5 Rayleigh-Ritz法 打靶法是将边值问题转化为一对初值来求边值问题的近似求解法。有限差分法将连续微分运算离散化为有限差分运算。Rayleigh-Ritz法是从第三种方法去解决这个问题的变分技术。边值问题首先被重新转述为从满足边界条件的所有充分可微函数中选择函数来使得某个积分达到最小值的问题。这样，可行性函数的集合范围被缩小，从而得到最小化问题的近似解且得到边值问题的近似解。 为了描述Rayleigh-Ritz法，我们考虑来自梁受力分析的线性两点边值问题的近似解。这个边值问题可描述为微分方程： （11.21） 边界条件： （11.22） 这个微分方程描述了长度为1的可变交叉截面(记为)的横梁的挠度。挠度取决于施加的压力和。 在下面的讨论中，我们假设且。进一步我们假设存在常数，使得 且， 这些假设保证了式(11.22)和式(11.23)所给出的边值问题有唯一解(见参考文献[BSW])。 如同许多描述物理现象的边值问题一样，梁方程的解满足变分特性。粱方程的变分原理对于开发Rayleigh-Ritz法是很基本的，同时它把粱方程的解表征为对中的所有函数使得某积分达到最小值的函数，其中函数集是指上那些满足性质的函数的集合。下面的定理描述了该特征。 定理11.4 设，，且 ，， 函数是微分方程 ， （11.23） 的唯一解，当且仅当是上唯一使得积分 （11.24） 达到最小的解。 该定理的具体证明可见参考文献[Shul，pp．88-89]。该证明步骤分三步进行。首先证明式（11.24）的任意解对所有都满足方程 （11.25） 第二步证明是式（11.25）的解当且仅当对于所有式（11.26）成立。 最后一步证明式（11.26）有唯一解。该唯一解也是式（11.25）和式（11.24）的解，故式（11.24）和式（11.25）的解是相同的。 Rayleigh-Ritz法就是通过最小化积分来逼近解，它并不是对上所有函数最小化，而是对包含某些基函数的线性组合的比较少的函数集最小化。这些基函数是线性无关的，且满足 ， 对方程（11.24）的解的近似解是通过找到常数使最小化来得到的。 由方程（11.25），有 （11.26） 且知道要达到最小，将视为的函数，则必有 ， （11.27） 对式（11.27）进行微分，得到 代入方程（11.28）得到 ， （11.28） 方程（11.29）所描述的法方程产生一个线性方程组Ac=b,变量为，其中A为对称矩阵，其元素为 定义b为 基函数最简单的选择是分段线性多项式。第一步是选择满足 的点来对[0,1]进行划分。设（），定义基函数为 （11.29） 其中，。（见图11.4） 图11.4 由于函数是分段线性的，所以导数不