当前位置:文档之家› 计算方法-论文

计算方法-论文

浅论拉格朗日与牛顿插值法一、课程简介计算方法是一种以计算机为工具,研究和解决有精确解而计算公式无法用手工完成和理论上有解而没有计算公式的数学问题的数值近似解的方法。

在实际中,数学与科学技术一向有着密切关系并相互影响,科学技术各领域的问题通过建立数学模型和数学产生密切的联系,并以各种形式应用于科学与工程领域。

而所建立的这些数学模型,在许多情况下,要获得精确解是十分困难的,甚至是不可能的,这就使得研究各种数学问题的近似解变的非常重要了,计算方法就是这样一门课程,一门专门用来研究各种数学问题的近似解的一门课程。

计算方法的一般步骤四:实际问题抽象出实际问题的物理模型,再有物理模型具体出数学模型,根据相关的数值方法利用计算机计算出结果。

从一般的过程可以看出,计算方法应该具有数学类课程的抽象性和严谨性的理论特性和实验课程的实用性和实验性的技术特征等。

随着计算机的飞速发展,数值计算方法已深入到计算物理、计算力学、计算化学、计算生物学、计算机经济学等各个领域,并且在航天航空、地质勘探、桥梁设计、天气预报和字形字样设计等实际问题领域得到广泛的应用。

二、主要内容《计算方法》这门课程可以分为三大块:数值逼近,数值代数,常微分方程。

1.数值逼近模块这模块的知识点主要分布在第一章到第三章。

第一章:数值计算中的误差。

主要的知识点是绝对误差和绝对误差限、相对误差和相对误差限、有效数字等概念的引入和计算绝对误差和绝对误差限、相对误差和相对误差限及有效数字的方法。

第二章:插值法。

在这一章中,主要的就是拉格朗日插值法与牛顿插值法的讲述。

拉格朗日插值法中核心就是去求插值结点的插值基函数,牛顿插值法中核心就是计算插值结点的差商,还有就是截断误差的说明。

第三章:曲线拟合的最小二乘法。

重点是最小二乘法的法则和法方程组列写,如何利用法方程组去求一个多项式各项的系数。

最小二乘法是与插值方法是有区别的,它不要求过所有的结点,只要靠近这些点,尽可能的表现出这些点的趋势就行了。

2.数值代数模块这一部分内容主要在第四章至第七章。

第四章:数值积分。

主要说的是插值型的数值积分的公式和积分系数。

刚开始讲了牛顿-柯特斯插值求积公式,包括梯形公式、Simpson公式、Cotes公式-系数、代数精度和截断误差。

然后就是复合的牛顿-柯特斯求积公式,包括复合的梯形公式、复合的Simpson公式、各个复合公式的收敛阶和它们各自的截断误差。

最后讲的是龙贝格算法的计算思想和公式的讲述。

第五章:非线性方程的数值解法。

在这一章中主要就是向我们介绍了四种非线性方程求根的迭代法,即为二分法、牛顿切线法、牛顿下山法和正割法。

牛顿切线法、牛顿下山法和正割法种方法的迭代公式是怎样的,各自的收敛阶,及它们相互之间的比较。

第六章:方程组的数值解法。

本章的内容讲的都是求解方程组的值,可以分为两类:一类是求解方程组的精确值的方法,即高斯列主元消去法、LU分解法和高斯消去法;另一类是求解方程组的近似解的方法,即Jacobi迭代法、S-R迭代法和SOR迭代法。

用迭代法求解方程组要判断所用的方法是否收敛,引入了矩阵的范数,迭代法迭代矩阵谱半径的求解,条件数及病态方程等知识。

3.常微分方程这个是在第七章:常微分方程的数值解法。

在这一章中讲的就是欧拉方法的介绍,由初值,利用欧拉方法去计算微分方程的值。

主要的内容就是欧拉公式、向后欧拉公式和改进的欧拉公式。

三、重点与难点1.数值逼近这一部分的重点与难点就是两种插值方法(即拉格朗日插值法和牛顿插值法)和插值条件。

在拉格朗日插值法中要知道如何去求每个插值结点的基函数,计算基函数是拉格朗日插值法的核心部分,并且要理解基函数的定义和插值余项。

在牛顿插值法中,要知道怎么去求差商,求差商是牛顿插值法的核心。

在这块的知识点中主要是要掌握好这两种插值法,利用它们去解决实际中的一些问题,知道它们的优缺点,根据实际的问题去选择用哪种方法解决实际中的问题。

2.数值代数这部分主要的是求积分近似值;求解非线性方程的解主要的三种方法(牛顿切线法、牛顿下山法和正割法);求方程组解的五中方法(高斯消去法、LU分解法、雅可比迭代法、高斯—塞德尔迭代法和SOR迭代法)以及这几种方法的收敛性是怎么样的,如何判断用这几种方法解方程组的根就是收敛的。

3.常微分方程在最后一章中主要的是掌握欧拉公式和改进的欧拉公式,学会怎么用欧拉公式和改进的欧拉公式来常微分方程的值。

四、拉格朗日与牛顿插值法由于在生产和科研中出现的函数是多种多样的,所以常常会遇到这样的情况:在某个实际问题中,虽然可以断定所考虑的函数f(x)在区间[a,b]上存在且连续,但却难以找到它的解析表达式,只能通过实验和观测得到在有限个点的函数值(即一张函数表)。

显然,要利用这张函数表来分析函数f(x)的性态,甚至直接求出其他一些点上的函数值可能是非常地困难。

在有些情况下,虽然可以写出函数f(x)的解析表达式,但是由于结构相当复杂,使用起来很是不方便。

面对这些情况,总是希望根据所得的函数表(或结构复杂的解析表达式),构造某个简单函数P(x)作为f(x)的近似。

插值法就是为了解决此类问题的一种古老的确实目前常用的方法,它不仅直接广泛地应用于生产实际和科学研究中,而且也是进一步学习数值计算方法的基础。

拉格朗日插值法和牛顿插值法就是两种常用的简便的插值法。

在这里主要的就是说一说这两种插值法的理论和比较,它们是属于数值逼近模块的知识。

1. 拉格朗日插值法在求满足插值条件n 次插值多项式P n (x)之前,先考虑一个简单的插值问题:对结点x i (i=0,1,…,n)中任一点x k (0≤k ≤n),作一n 次多项式l k (x),使它在该点上取值为1,而在其余点x(i=0,1,…k-1,k+1…,n)上取值为零,即L 型插值多项式:2. 牛顿插值法由线性代数知,任何一个不高于n 次多项式,都可以表示成函数1,x-x 0,(x-x 0)(x-x 1),…,(x-x 0)(x-x 1)…(x-x n-1)的线性组合。

既可以把满足插值条件P(x i )=y i (i=0,1,…,n)的n 次插值多项式写成如下形式:a 0+ a 1(x-x 0)+ a 2(x-x 0)( x-x 1)+…+ a n ( x-x 0)( x-x 1)…(x-x n-1)其中,a k 为待定系数。

这种形式的插值多项式成为牛顿插值多项式,记为N n (x)。

对个n+1个互异节点:x 0, x 1,…,x n011011()()()()()()11()()()()k k k k n k k k k k i k i k n l x a x x x x x x x x l x a x x x x x x x x -+-+=----=⇒=---- 又:011,,)n x =两个互异的插值节点(x 01010110() , ()x x x x l x l x x x x x --==--插值基函数:10011()()()()()L x f x l x f x l x =+线性插值函数:0()()()()nn n kk k L x P x lx f x ===∑101(1)12()[,] ()[,]()()()()()()()(1)!n n n n n n n n f x C a b a x x x b P x x a b fR x f x P x f x L x x n ξω+++∈≤<<≤∀∈=-=-=+ 定理:若,且节点,则插值多项式对有:3. 两者的比较牛顿插值的误差不要求函数的高阶导数的存在,所以更具有一般性。

它对f(x)是由离散点给出的函数情形或f(x)的导数不存在的情形均适用。

拉格朗日插值法公式结构紧凑,在理论分析中方便,但是如果遇到结点的增减,所有的数据需要全部重算,没有承袭性。

而牛顿插值就避免了这一缺点,这样的话,在用计算机计算是就可以大量的节省乘除运算的次数,减少了计算的时间,所以可以说对于一些结构相当复杂的函数f(x),牛顿插值法比拉格朗日插值法要占有一定的优势。

五、心得与体会通过这学期对《计算方法》这门课程的学习对一些知识有了较清晰的认识和了解,印象也比较深刻。

计算方法中所提到的各种方法都有其自己所使用的范围和使用所需的注意事项,同样在相同问题的处理上,不同的处理方法的选择可能会造成两种截然不同的处理结果,也可能会造成误差的扩大化或者不是最有效的解决方法。

误差在数值计算中是不可避免的,误差的传播和积累直接影响到计算结果的精度。

在研究算法的同时,必须注重误差的分析,使建立起来的算法科学有效。

绝大多数情况下不存在绝对的严格和精确,在考虑数值算法时要能将误差限制在许可的范围之内,这种算法就是数值稳定的。

一般情况下除了选用较好的计算方法来防止误差传播和积累以外,还可以选用稳定性较好的计算公式、简化计算步骤和公式、合理安排运算顺序(避免大数“淹没”小数;多个数相加时,其绝对值小的先加;多个数相乘时,其有效位数多者先乘)、避免两相近数相减和绝对值太小的数作为除数等。

在这门课中,我觉得主要的是掌握解决实际问题的方法和思想,计算方法告诉了我们处理问题的方法,而我们在实际运用中就需要通过自己的判断和经验来选择实用性最好、使用最简单、误差最小的方法来解决实际问题。

10011012010122010100()()[,](,)(,)[,,][,,][,,][,,]n n n n f x f x f x x x x f x x f x x f x x x x x f x x f x x f x x x x --=--=--=- 称为n 阶差商称为1阶差商称为2阶差商)()()()( 10010---++-+=∴n n n x x x x a x x a a x N ⎪⎪⎩⎪⎪⎨⎧=--++-+==--+-+==-+===-)()()()()()())(()()()()()()()(1001021202202102101101000n n n n n n n n n n n x f x x x x a x x a a x N x f x x x x a x x a a x N x f x x a a x N x f a x N。

相关主题