) 中图学技术大荸学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 第二章数值微分和数值积分
数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 第二章 数值微分和数值积分
) 中图学技术大荸学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 数值微分 1.函数fx)以离散点列给出时,而要求我们给出导数值, 2.函数f(x)过于复杂 这两种情况都要求我们用数值的方法求函数的导数值 微积分中,关于导数的定义如下: f(x)=lm f(x+h)-f(x f(x)-f(x-h) f(x+h)-f(x h→>0 h h→0 h h→>0 2h 自然,而又简单的方法就是,取极限的近似值,即差商
数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 数值微分 h f x h f x h h f x f x h h f x h f x f x h h h 2 ( ) ( ) lim ( ) ( ) lim ( ) ( ) '( ) lim 0 0 0 + − − = − − = + − = → → → 1. 函数f(x)以离散点列给出时，而要求我们给出导数值， 2. 函数f(x)过于复杂 这两种情况都要求我们用数值的方法求函数的导数值 微积分中，关于导数的定义如下： 自然，而又简单的方法就是，取极限的近似值，即差商
) 中图学技术大荸学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 向前差商 f(xo +h)-f(ro) 由 Taylor展开 f(x+h)=f(x)+h(x0)+f"(5x≤5≤x0+h 2 因此,有误差 R(x)=f(x0) f(ro +h)fo ∫"'()=O(h
数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS h f x h f x f x ( ) ( ) '( ) 0 0 0 + − 由Taylor展开 f x x h h f x + h = f x + hf x + 0 0 + 2 0 0 0 ''( ), 2! ( ) ( ) '( ) 因此，有误差 ''( ) ( ) 2! ( ) ( ) ( ) '( ) 0 0 0 f O h h h f x h f x R x f x = − = + − = − 向前差商
) 中图学技术大荸学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 向后差商 f(o) f(o)-f(xo -h h 由 Taylor展开 h f(x。-h)=f(x)-b(x)+f"(5,x0≤5≤x+h 因此,有误差 R(x)=f(ro)f(xo)-f(xo-h)h =f"()=O(h h
数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS h f x f x h f x ( ) ( ) '( ) 0 0 0 − − 由Taylor展开 f x x h h f x − h = f x − hf x + 0 0 + 2 0 0 0 ''( ), 2! ( ) ( ) '( ) 因此，有误差 ''( ) ( ) 2! ( ) ( ) ( ) '( ) 0 0 0 f O h h h f x f x h R x f x = = − − = − 向后差商
中图苔技术大荸数学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 中心差商 f(x0+h)-f(x0-h) 2h 由 Taylor展开 h f(xo+h)=f(o)+hf(xo)+af(xo)+af (s,xo sssxo +h 2 f(x。-h)=f(x0)-h(x0)+f"(x0)-f"(52),x-h≤52≤x 因此,有误差 R(x)=f(x0) (xo+h)-f(o-h) 2h I八"(1)+f"(2)=f"(5)=O(h2)
数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS h f x h f x h f x 2 ( ) ( ) '( ) 0 0 0 + − − 由Taylor展开 2 3 0 0 0 0 1 0 1 0 2 3 0 0 0 0 2 0 2 0 ( ) ( ) '( ) ''( ) '''( ), 2! 3! ( ) ( ) '( ) ''( ) '''( ), 2! 3! h h f x h f x hf x f x f x x h h h f x h f x hf x f x f x h x + = + + + + − = − + − − 因此，有误差 '''( ) ( ) 6 [ '''( ) '''( )] 12 2 ( ) ( ) ( ) '( ) 2 2 1 2 2 0 0 0 f O h h f f h h f x h f x h R x f x = + = = + − − = − 中心差商
) 中图学技术大荸学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 由误差表达式,h越小,误差越小,但同时舍入误差增大,所以,有个最佳步长 我们可以用事后误差估计的方法来确定 设D(h),D(h/2)分别为步长为hh/2的差商公式。则 h D(h)-D()<E 时的步长h2就是合适的步长 f(x)-D(h)=O(h) f(x)-d(h) O(h) 2 f(x)-D(h/2)=O(h/2)f(x)-D(h/2)O(h/2) f"(x)-D(h)=2f(x)-2D(h/2) f"(x)-D(h/2)=D(h)-D(h/2)
数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 由误差表达式，h越小，误差越小，但同时舍入误差增大，所以，有个最佳步长 我们可以用事后误差估计的方法来确定 设D(h),D(h/2)分别为步长为h,h/2的差商公式。则 − ) 2 ( ) ( h D h D 时的步长h/2就是合适的步长 '( ) ( ) ( ) '( ) ( / 2) ( / 2) f x D h O h f x D h O h − = − = '( ) ( ) ( ) 2 '( ) ( / 2) ( / 2) f x D h O h f x D h O h − = − f x D h f x D h '( ) ( ) 2 '( ) 2 ( / 2) − = − f x D h D h D h '( ) ( / 2) ( ) ( / 2) − = −
) 中图学技术大荸学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 例: f(ax=exp(x) f(1.15)|R(x) f'(1.15)R(x) 0.10 3.16300.00480053.1590-00008 0.09 3.1622-0.00400.043.1588-0.0006 0.08 3.1613|-0.00310.033.1583-0.0001 0m7|816070051002 3.1575-0.0007 3.1600 0.00180013.1550-0.0032
数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS f(x)=exp(x) h f’(1.15) R(x) h f’(1.15) R(x) 0.10 3.1630 -0.0048 0.05 3.1590 -0.0008 0.09 3.1622 -0.0040 0.04 3.1588 -0.0006 0.08 3.1613 -0.0031 0.03 3.1583 -0.0001 0.07 3.1607 -0.0025 0.02 3.1575 -0.0007 0.06 3.1600 -0.0018 0.01 3.1550 -0.0032 例：
) 中图学技术大荸学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 插值型数值微分 插值是建立逼近函数的手段,用以研究原函数的性质。因此,可以用插值 函数的导数近似为原函数的导数 f6)(x)≈Ln6(x) 误差 (n+1) R,(x) (5 O,(x)=f(x)-L,(x) n+ (n+1) R,(x) (2 0(x dx(n+1)
数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 插值是建立逼近函数的手段，用以研究原函数的性质。因此，可以用插值 函数的导数近似为原函数的导数 ( ) ( ) ( ) ( ) f x L x k n k 误差 ( ) ( ) ( ) ( 1)! ( ) ( ) ( 1) x f x L x n f R x n n n n = − + = + + = + ( ) ( 1)! ( ) ( ) ( 1) ( ) x n f dx d R x n n k k k n 插值型数值微分
) 中图学技术大荸学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 例: 给定点列{x,f(x)=0且x2-x1=x1-x0=h,求 ∫(x2),f(x)(x) 解: X-Xr-x (x) f(x0)+ (x-x0)(x-x2) f(x1)+ x-x0(x-x1) f(x2) 2h h ch x-Xtx-x x-xo+x-x xxo+x-x x1)+ 2h2 2h
数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 给定点列 2 0 ( , ( )) i i i= x f x 且 x2 − x1 = x1 − x0 = h ，求 '( ), '( ), '( ) 2 1 0 f x f x f x 解： ( ) 2 ( )( ) ( ) ( )( ) ( ) 2 ( )( ) ( ) 2 2 0 1 2 1 0 2 2 0 1 2 2 f x h x x x x f x h x x x x f x h x x x x L x − − + − − − + − − = ( ) 2 ( ) ( ) ( ) ( ) 2 ( ) ' ( ) 2 2 0 1 2 1 0 2 2 0 1 2 2 f x h x x x x f x h x x x x f x h x x x x L x − + − + − − + − + − + − = 例：
) 中图学技术大荸学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS f(x)≈L2(x0) 2h (-3f(x)+4f(x1)-f(x2)+,f"(5) f(x)≈L2(x1)=1(-f(x)+f(x2)-f"() f(x)≈L2(x2)=(f(x)-4f(x)+3f(x2)+,f"( f"(x)≈L2(x)=2(f(x)-2f(x)+f(x2)+[-"(5)+() f"(x)≈L"(x2)=(f(x)-2f(x)+f(x2),f((5) f"(x2)≈L"(x2)=2(f(x)-2f(x)+f(x2)+[f"(5)-f((52) Taylor)展开分析,可以知道,它们都是O(h2)称为三点公式
数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS ( ) 2 0 2 0 0 1 2 1 '( ) ' ( ) 3 ( ) 4 ( ) ( ) '''( ) 2 3 h f x L x f x f x f x f h = − + − + ( ) 2 1 2 1 0 2 1 '( ) ' ( ) ( ) ( ) '''( ) 2 6 h f x L x f x f x f h = − + − ( ) 2 2 2 2 0 1 2 1 '( ) ' ( ) ( ) 4 ( ) 3 ( ) '''( ) 2 3 h f x L x f x f x f x f h = − + + Taylor展开分析，可以知道，它们都是 ( ) 2 O h 称为三点公式 ( ) 2 (4) 0 2 0 0 1 2 1 2 2 1 ''( ) '' ( ) ( ) 2 ( ) ( ) [ '''( ) ( )] 6 h f x L x f x f x f x hf f h = − + + − + ( ) 2 (4) 1 1 2 0 1 2 2 1 ''( ) '' ( ) ( ) 2 ( ) ( ) ( ) 12 h f x L x f x f x f x f h = − + − ( ) 2 (4) 2 2 2 0 1 2 1 2 2 1 ''( ) '' ( ) ( ) 2 ( ) ( ) [ '''( ) ( )] 6 h f x L x f x f x f x hf f h = − + + −