数值分析数值积分与数值微分应用.ppt

* Gauss型数值求积公式 正交多项式及其零点 数值微分方法 数值积分与数值微分应用 《数值分析》 22 ? ? ? ? 高斯型数值求积公式 考虑两点插值型求积公式 为使代数精度尽可能高,取 f(x)=1, x, x2, x3 (1) (2) (3) (4) (4)-(2)×x02 x12= x02 (3)-(1) ×x02 x02=1/3 2/18 显然,A0=1,A1=1.代数精度为3的数值求积公式为 对于[a, b]区间上的定积分,构造变换 t∈[-1, 1] 3/18 定义 如果求积结点x0, x1,······,xn,使插值型求积公式 的代数精度为2n+1,则称该求积公式为Gauss型求积公式. 称这些求积结点为Gauss点. 定理7.2 如果多项式wn+1(x)=(x – x0) (x – x1)···(x – xn) 与任意的不超过n次的多项式P(x) 正交，即 则, wn+1(x)的所有零点x0, x1 ,······, xn 是Gauss点 证明: 设f (x)是任意(2n+1)次多项式 , 由多项式除法 4/18 构造插值型求积公式,有 由于 Q(xk)= Q(xk)+wk+1(xk)P(xk) = f (xk) 其中,P(x) ,Q(x) 均为n次多项式. 两端积分，得 其中, 所以 插值结点为wn+1(x)的零点 5/18 例 证明多项式 是[–1,1]上正交多项式. 证: 显然 两点Gauss公式 得Gauss点 插值公式: 6/18 三点Gauss数值求积公式 Legendre多项式递推式 7/18 例.用两点Gauss公式计算 解:作变换 x = 0.5( t + 1 ), 则 取 = 0.9460411 MATLAB函数: sinint(1)=0.946083070367183 8/18 例.测得一个运动物体的距离D(t)数据如下 t 8.0 9.0 10.0 11.0 12.0 D(t) 17.45 21.46 25.75 30.30 35.08 用数值微分求速率v(10) 9/18 Tylor展开 方法 一阶向前差商 一阶向后差商 10/18 二阶中心差商 一阶中心差商 11/18 隐式方法: 设 xk= a + k h ,( k =0,1,···,n) 值 令 mk = f ’(xk) ( k=0,1,···,n) k=1,2,···,n-1 12/18 Lagrange插值函数方法 xk = x0 + kh (k=0，1，2) 13/18 *