(偏微分方程数值解实验报告.docVIP

成 绩 09信计2011-2012（二） 偏微分方程数值解法实验 实验题目　 Euler法与Adams内外插法的偏微分数值解法 实验题目 2012年5月2日 学生姓名　 夏秋涵、梁化义、武良鹏、赵露青、唐琦 所在班级　 09信计3 指导教师　 杨 扬 徐州工程学院数学与物理科学学院 目??? 录..........................................3 1.1 Euler数值积分法................................3 1.2 Adams方法 二、背景..........................................3 三、实验..........................................4 3.1 预备知识........................................4 3.2 实验题目........................................4 3.3 实验过程分析....................................4 3.4 实验结果与误差分析..............................6 四、总结…………………………………………………10 五、参考书文献....................................10 一、摘要 1.1 Euler数值积分法： 分析数学的重要分支之一。包括一个自变量和它的未知函数以及未知函数的微商的等式叫做常...当时的许多实际问题只能用数值方法求近似解，欧拉折线法便是这方面工作的开端。19世纪初期和中期是数学发展史上的一个转折时期，分析基础的重建、复变函数，等等 ... 1.2 Adams方法?延迟积分微分方程广泛应用于物理学、生物学、生态学及控制科学等科学领域,这类方程由于通常很难获得理论解的解析式,因此研究这类方程的数值方法是十分有意义的。为求解这些方程,学者们构造了许多的数值计算方法,如:Runge-Kutta方法、线性多步法、Rosonbrock方法等。 Adams方法是线性多步方法的一种特殊形式,由于Adams方法具有良好的稳定性和方便的变阶变步长方法,因而在求解刚性问题时Adams方法显现出其突出的优势。然而,在使用Adams方法时不可避免的要遇到一个问题:步长和阶的选择。选择适当的阶和步长是非常重要的,这将直接影响Adams方法的精度和效率。因此Adams方法的步长和阶的选择也成为了许多研究的主题。 本文首先研究了求解一类Volterra离散-分布型延迟积分微分方程的扩展Adams方法。我们扩展了一般隐式Adams方法求解延迟问题,采用同阶的求积公式离散方程中分布型延迟函数,并用Newton迭代方法执行函数迭代。 其次本文还构造了求解刚性延迟微分方程的变阶变步长Adams方法。利用Nordsieck方法实现变阶变步长策略。再引入小参数改造Adams方法从而通过对参数的选择实现对稳定性和收敛性的控制。 用Euler法和改进的Euler法求解，其中步长h=0.1,0.05,0.01 用三阶Adams外插法及内插法求解，步长h=0.1,0.05,0.01 3.3实验过程分析： 3.3.1 本次实验代码用c++语言编写，编译器为VC++6.0.Euler法最简单的数值积分法， 数值格式也很简单。在将数学语言转化为计算机语言时，可以直接转化，无需作其他计算或转换。在程序代码中，其数值格式用C++代码表示为 u[j]=u[j-1]+h*(-5)*u[j-1]; 在编写程序时，先输入步长h的值，针对h将区间(0,1)等分成N=1/h份。然后计算u在这些节点上的近似值。经计算可以得出微分方程满足边值条件的解为。计算真解在节点上的值。输出时，只是输出及数值解误差-.对其优良性进行分析。按照书本上的分析，Euler法的局部截断误差为阶，故其全局误差为阶。因此可以预期它的逼近效果不会很理想。 3.3.2 进行改进的Euler法格式实验时，由于右端项只与u有关，t没有显式出现，故进行迭代时，可以将右端的移项到左边，合并