《典型例题与习题4》由会员分享，可在线阅读，更多相关《典型例题与习题4（16页珍藏版）》请在金锄头文库上搜索。
1、数值求积公式及代数精度数值求积公式及代数精度 数值求导方法与截断误差数值求导方法与截断误差 一阶常微分方程数值法一阶常微分方程数值法 局部截断误差与精度局部截断误差与精度 数值分析典型例题典型例题 IV 插值型求积公式插值型求积公式: 求积系数求积系数 ), 2 , 1 , 0(,)(njdxxlAbajj )()(0fRxfAdxxfnjjjba 求积余项求积余项 bannbandxxnfdxxLxffR)()!1()()()(1)1( 等距结点插值型求积公式称为等距结点插值型求积公式称为Newton-Cotes公式公式, 偶数阶偶数阶Newton-Cotes公式至少有公式至少有(n+1)阶
2、代数精度阶代数精度 2/16 求积结点求积结点 bxxxan 101.梯形公式梯形公式 )(12)()()(2)(3 fabbfafabdxxfba )(2)()(2)(11 njbajhafbfafhdxxf复合梯形求积公式复合梯形求积公式 令令h=(b-a)/n ),(, )(12)(23 bafnabfR 求积余项求积余项 3/16 2. 辛卜生公式辛卜生公式 )()2(4)(6)(bfbafafabdxxfba )(290)()4( 55 fabfR 求积余项求积余项: 两点高斯型数值求积公式两点高斯型数值求积公式 )()()(31 3111ffdxxf)()()(2322322aba
3、bfababfabdxxfba4/16 练习练习: 复合辛卜生公式求积余项复合辛卜生公式求积余项? )()()()(hOhafhafaf一阶向前差商一阶向前差商 )()()()(hOhhafafaf一阶向后差商一阶向后差商 )()()()()(2 22hOhhafafhafaf 二阶中心差商二阶中心差商 )()()()(2 2hOhhafhafaf一阶中心差商一阶中心差商 5/16 外推算法外推算法 )()()(hxfhxfhhG214 22 1hhxfhG )()(142411 mmmmmhGhG hG)()( )()()()()(12m mhOhGxf 16/4/)()2/(4 22 1h
4、hxfhG 3)()2/(4hGhG 4/3)(4 2hxf 练习练习：二阶中心差商的外推公式？：二阶中心差商的外推公式？ 6/16 000 )(),(yxyxxyxfy1. Euler方法方法 ), 2 , 1 , 0( , ),(),(1100 nyxhfyyhxxxyynnnnnn常微分方程初值问题常微分方程初值问题 ),(),(2111 nnnnnnyxfyxfhyy2. 梯形公式梯形公式: 7/16 ),( nnnnyxhfyy1),(),(1112nnnnnnyxfyxfhyy预测预测- -校正公式校正公式 局部截断误差局部截断误差 设设 yn= y(xn), 称称Rn+1=y(x
5、n+1) - yn+1为局部截断误差为局部截断误差 常表示为常表示为: O(hp+1)， p 称为单步法的精度阶数称为单步法的精度阶数 又称为修正的又称为修正的Euler公式公式 yn+1= yn+ 0.5h k1+ k2 k1=f(xn,yn), k2=f(xn+h, yn+hk1) 8/16 Ex1.推导左矩形求积公式推导左矩形求积公式 2)(2)()()()(abfafabdxxfba 2)(2)()()()(abfbfabdxxfba 3)(24)()2()()(abfbafabdxxfba 令令 uadxxfuF)()(F(u)= F(a) + (u-a)F(a) +0.5(u-a)
6、2F ”( ) )()(),()(, 0)( fFafaFaF 2)(2)()()()(abfafabdxxfba 练习练习: 9/16 Ex3. 求复合中矩形公式求复合中矩形公式 的求积误差的求积误差? 10)5 . 0()(njbahjafhdxxfEx2.复合左矩形求积公式的求积误差复合左矩形求积公式的求积误差 njjnjbafhjhafhdxxf1210)(2)()( 设被积函数在积分区间上的一阶导数连续设被积函数在积分区间上的一阶导数连续,由连续函数由连续函数 介值定理介值定理 ),()()(11baffnnjj ),()(2)()(2)()(2122 bafnabfnabxfRnj
7、j 10/16 Ex4.利用复合梯形公式计算积分利用复合梯形公式计算积分 10sindxxxI使其截断误差不超过使其截断误差不超过 0.510-3,应算多少次函数值？应算多少次函数值？ 提示提示： 10)cos(sin)(dtxtxxxf练习练习: : 给定积分给定积分 当要求误差小于当要求误差小于10-3时用复合梯形公式和时用复合梯形公式和Simpson公公 式计算时式计算时, 需要计算多少次函数值？需要计算多少次函数值？ dxxex 31sin11/16 Ex5. 验证，复合梯形公式与复合验证，复合梯形公式与复合Simpson 公式之公式之 间有如下关系间有如下关系 43122mmmTTS
8、 )(2)()(211 mjmjhafbfafhTmabh )5 . 0(2112 mjmmhjhafhTT )2(23143112 mjmmmhjhafhTTT )(4)(2)()(31121122 mjjmjjmxfxfbfafhmabhm22 mjjmjjmSxfxfxfh221211222)()(4)(3 12/16 Ex6. 定积分定积分 的计算问题可化为初值问题的计算问题可化为初值问题 y= f (t) , y(a)=0 试证明用试证明用Euler公式计算结果为公式计算结果为 其中其中, h = (b a )/N, tn= a + n h ( n = 0,1,2, , N) bad
9、xxf)(10Nnnhtfby)()(Ex7. 试证明试证明4阶阶Range-Kutta公式解公式解a, b内初值问题内初值问题 y= f (x) , y(a)=0 结果有结果有: 其中其中, h = (b a )/N, xn= a + n h ( n = 0,1,2, , N) 1012/1)()(4)(6)(Nnnnnxfxfxfhby13/16 14/16 Ex 8 将积分上限函数将积分上限函数 xdttxxf 022)exp()exp()(转化为常微分方程初值问题。并推导用二阶和四阶龙转化为常微分方程初值问题。并推导用二阶和四阶龙 格格- -库塔方法求解的计算公式库塔方法求解的计算公式
10、 Ex9. 初值问题初值问题 有解有解y(x)=0.5a x2 + b x 。若取若取 xn = nh，yn为欧为欧 拉方法得到的数值解，试证明拉方法得到的数值解，试证明 y(xn) yn = 0.5 a h xn 0)0(ybaxy若取若取 xn = nh，yn为用梯形公式计算所得的数值解，为用梯形公式计算所得的数值解， 记记y(xn)为初值问题的在为初值问题的在x=xn处的解析解。试证明处的解析解。试证明: y(xn) = yn Ex10. 初值问题初值问题 0)0(ybaxy15/16 Ex11. 推导用于求解初值问题的推导用于求解初值问题的Adamas显格式显格式: 16/16 ),(),(321112nnnnnnyxfyxfhyy 23nnyy),(5),(16),(23121122nnnnnnyxfyxfyxfh
