数值分析 数值积分与数值微分.pptVIP

数值分析 §3 龙贝格算法 §3 龙贝格算法 §3 龙贝格算法 §3 龙贝格算法 §3 龙贝格算法 §3 龙贝格算法 例： ? + ? 1 0 1 1 0 0 ) ( ) ( ) ( x f A x f A dx x f x Step 1：构造正交多项式?2 设 c bx x x a x x x + + = + = = 2 2 1 0 ) ( , ) ( , 1 ) ( j j j ? ? 5 3 - = a 0 ) ( 1 0 = + ? dx a x x 0 ) , ( 1 0 = j j ? ? = + + - ? = = + + ? = 1 0 2 1 1 0 2 2 0 0 ) )( 5 3 ( 0 ) , ( 0 ) ( 0 ) , ( dx c bx x x x dx c bx x x j j j j 21 5 9 10 = - = c b 即： Step 2：求?2 = 0 的 2 个根，即为 Gauss 点 x0 ，x1 Step 3：代入 f (x) = 1, x 以求解 A0 ，A1 解线性方程组，简单。 结果与前一方法相同： ? 利用此公式计算 的值 二、Gauss-Legendre求积公式 §4 高斯型求积公式 1.Legendre多项式 其具有前面提到的正交的性质，即对于任意次数不超过n的多项式q(x)，成立 因此，多项式Pn+1(x)的零点就是相应的Gauss-Legendre公式的高斯点。 Legendre多项式的前几项如下： §4 高斯型求积公式 Legendre多项式的性质 性质1： Legendre多项式的首项系数为 性质2：当n为奇数时，Pn(x)为奇函数，当n为偶数时， Pn(x)为偶函数。 性质3：对一切次数不高于n-1次的多项式q(x),有 §4 高斯型求积公式 §4 高斯型求积公式 2.Gauss-Legendre求积公式 （1）当n=0时，P1(x)=x，其零点为x0=0，易得A0=2 的高斯－勒让德求积公式是2f(0) 其代数精确度为1。 （2）当n=1时， ，其零点为 ，此时A0=1，A1=1， 的高斯－勒让德求积公式是 其代数精确度为1。 以此类推，可得到n=2,3…时的高斯－勒让德求积公式。课本中表4.5列出了高斯－勒让德求积公式的节点和系数。 §4 高斯型求积公式 三、高斯公式的余项 §5 数值微分 导数是用极限来定义的，如果一个函数是初等函数，我们可以用求导法则来求得其导函数或在某点的导数值。如果我们只知道一个函数在若干已知离散点上的函数值，则我们必须利用数值方法来求解函数的导数，这就是数值微分。 一、插值型求导公式 已知f(x)在n+1个点上的函数值，我们可以构造插值多项式Pn(x)来近似代替原函数f(x)，从而可以利用P′n(x)来代替f ′(x)，这就称为插值型求导公式。 §5 数值微分 插值型求导公式的截断误差： 因为 所以 数值微分的截断误差为 由于我们只讨论在节点上的数值微分，而 ＝0 于是 一般地 Company name LOGO 主讲 侯晓慧 第四章 数值积分与数值微分 其中F(x)是f(x)的原函数之一，可用不定积分求得。而在实际问题中，大量函数的原函数不容易或根本无法求出。例如，概率统计中常用的概率积分 ,及积分 等根本无法用初等函数来表示其原函数，因而也无法精确计算其定积分，只能运用数值积分。 计算定积分可用牛顿－莱布尼兹公式计算 §1 数值积分和代数精确度 一、数值积分公式 在区间[a, b]上的定积分 其某个数值积分公式就是在区间[a, b]内取n+1个点x0,x1,…,xn。利用被积函数f(x)在这n+1个点的函数值的某一种线性组合来近似作为待定求定积分的值，即 右端公式称为左边定积分的某个数值积分公式。其中xk为积分节点，Ak为相应的求积系数。因此，一个数值积分公式关键在于积分节点xk的选取和积分系数Ak的决定，其中Ak与被积函数f(x)无关。 §1 数值积分和代数精确度 二、积分公式的代数精确度 定义：若积分 的数值积分公式 对于任意一个次数不高于m次的多项式都精确成立，而存在一个m+1次多项式使之不精确成立，则称该数值积分公式的代数精确度为m。 对于代数精确度为m的求积公式，若f(x)是不超过m次的代数多项式，求积