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

成 绩 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方法从而通过对参数的选择实现对稳定性和收敛性的控制。 下面分别是h=0.05，h=0.01时的结果输出： （h=0.05时的结果，与h同阶） （h=0.01时的结果） 从输出结果可以看出，该算法的全局误差都为阶。与理论结果一致。随着h的减小，误差也相应减小。当然这是建立在误差阶为的基础上的。 3.4.2 下面分别是改进的Euler法在h=0.1，0.05.0.01时的结果截图： （h=0.1） （h=0.05） （h=0.01，其中e-006表示） 经分析，改进的Euler法的全局误差为阶。而上面的结果则很好的符合了这个理论。而h越小，则误差越小。且其误差比Euler法要好很多。这个也和理论相符。 3.4.3 下面分别是Adams外插法在h=0.1，0.05，0.01时的输出结果： （h=0.1） （h=0.05） （h=0.01） 根据题目要求，全局误差要为3阶，应该选取k=2.而当h=0.1及0.01时，出现部分结果的误差为两阶。与理论结果不太符合。而h越小，则误差越小这个理论还是符合的。说明这个数值格式不是很尽人意。而由于u