常微分方程边值问题的数值解法重点.docVIP

常微分方程论文 题目： 常微分方程边值问题的数值解法 组 长： 数学132文洲 组 员： 数学131王琦 数学132姚瑶 信息132郭斌 院 （系）： 理学院 指导教师： 岳宗敏 时 间： 2015年6月9日 常微分方程边值问题的数值解法 摘要：作为一类定解问题，补充条件由以自变量取某些值时，未知函数及其导数的值而定，称其为边值条件。许多物理和数学问题都归结为边值问题。本文介绍边值问题的待定常数法和格林函数。 关键词：边值问题 待定常数法 格林函数 Abstract: as a kind of definite solution problems, the supplementary conditions by took the certain values in the independent variable, the value of the unknown function and its derivative, referred to as boundary value conditions. Many physical and mathematical problems boil down to boundary value problems. In this paper, the boundary value problem of the method of undetermined constants and greens function. Keywords: boundary value problem Method of undetermined constants Greens function 11.1 引言 在很多实际问题中都会遇到求解常微分方程边值问题. 考虑如下形式的二阶常微分方程 , , (11.1.1) 在如下三种边界条件下的定解问题: 第一种边界条件: , (11.1.2) 第二种边界条件: , (11.1.2) 第三种边界条件: , (11.1.13) 其中. 常微分方程边值问题有很多不同解法, 这里只介绍打靶方法和有限差分方法. 11.2 打靶算法 将边值问题转化成初值问题来求解，即根据边界条件(11.1.2)，也就是说，反复是调整初始时刻的斜率值，使得初值问题的积分曲线能“命中”。 例11.1 利用打靶法求解两点边值问题 解：(1) 为应用打靶方法，需要假定初值来先求解初值问题，取初始参数 ， 解 令，，则上述初值问题变为 ， 由标准的R-K 方法，取，可得。 （2）再令， 解 ， 取。 仍由标准的R-K 方法，可得。 （3）令 再解， 取，得。 （4）令，得 （5）类似有，， 因为， 可以认为即为所求。 11.3 有限差分方法 原理：微分方程数值解的方法。基本思想是把连续的定解区域用有限个离散点构成的网格来代替， 这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似， 积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组 ?， 解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。, , (11.3.1) , . (11.3.2) 不失一般性, 将求解区间等份, 分点为, ,…, , 步长, 这里 称为结点. 在内部结点上, 记 , , , , 并分别用一阶和二阶中心差分来代替一阶导数和二阶导数, , . 将上面各式代入方程(11.3.1), 略去项, 得到关系式: , 将上式改写成 , , (11.3.3) 在上式中 , , ,, (11.3.4) 利用边界条件 , 可以得到如下形式的方程组: (11.3.5) 其中 (11.3.6) (11.3.7)