《数值分析典型例题与习题4.课件》由会员分享，可在线阅读，更多相关《数值分析典型例题与习题4.课件（16页珍藏版）》请在金锄头文库上搜索。
1、数值求积公式及代数精度数值求积公式及代数精度数值求导方法与截断误差数值求导方法与截断误差一阶常微分方程数值法一阶常微分方程数值法局部截断误差与精度局部截断误差与精度数值分析典型例题典型例题 IV插值型求积公式插值型求积公式:求积系数求积系数求积余项求积余项 等距结点插值型求积公式称为等距结点插值型求积公式称为Newton-Cotes公式公式,偶数阶偶数阶Newton-Cotes公式至少有公式至少有(n+1)阶代数精度阶代数精度2/16求积结点求积结点1.梯形公式梯形公式复合梯形求积公式复合梯形求积公式令令h=(b-a)/n求积余项求积余项3/162.辛卜生公式辛卜生公式求积余项求积余项:两点高
2、斯型数值求积公式两点高斯型数值求积公式4/16练习练习:复合辛卜生公式求积余项复合辛卜生公式求积余项?一阶一阶向前差商向前差商一阶一阶向后差商向后差商二阶中心差商二阶中心差商一阶中心差商一阶中心差商5/16外外推算法推算法练习练习：二阶中心差商的外推公式？：二阶中心差商的外推公式？6/161.Euler方法方法 常微分方程初值问题常微分方程初值问题2.梯形公式梯形公式:7/16预测预测- -校正公式校正公式局部截断误差局部截断误差设设yn= y(xn),称称Rn+1=y(xn+1)-yn+1为为局部截断误差局部截断误差常表示为常表示为:O(hp+1)， p 称为单步法的精度阶数称为单步法的精度
3、阶数又称为修正的又称为修正的Euler公式公式yn+1= yn+0.5hk1+k2 k1=f(xn,yn),k2=f(xn+h,yn+hk1)8/16Ex1.推导左矩形求积公式推导左矩形求积公式令令F(u)=F(a)+(u-a)F(a)+0.5(u-a)2F ”( )练习练习:9/16Ex3.求复合中矩形公式求复合中矩形公式的求积误差的求积误差?Ex2.复合左矩形求积公式的求积误差复合左矩形求积公式的求积误差设被积设被积函数在积分区间上的一阶导数连续函数在积分区间上的一阶导数连续,由连续函由连续函数介值定理数介值定理10/16Ex4.利用复合梯形公式计算积分利用复合梯形公式计算积分使其截断误差
4、不超过使其截断误差不超过0.510-3,应算多少次函数值应算多少次函数值？提示提示：练习练习: :给定积分给定积分当要求误差小于当要求误差小于10-3时用复合梯形公式和时用复合梯形公式和Simpson公公式计算时式计算时,需要计算多少次函数值？需要计算多少次函数值？11/16Ex5.验证，复合梯形公式与复合验证，复合梯形公式与复合Simpson 公式之公式之间有如下关系间有如下关系12/16Ex6.定积分定积分的计算问题可化为初值问题的计算问题可化为初值问题y= f (t),y(a)=0试证明用试证明用Euler公式计算结果为公式计算结果为其中其中,h = (b a )/N, tn= a +
5、n h ( n = 0,1,2,N)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)13/1614/16Ex8 将积分上限函数将积分上限函数转化为常微分方程初值问题。并推导用二阶和四阶龙转化为常微分方程初值问题。并推导用二阶和四阶龙格格- -库塔方法求解的计算公式库塔方法求解的计算公式Ex9.初值问题初值问题有解有解y(x)=0.5a x2 + b x。若取若取 xn = nh，yn为欧为欧拉方法得到的数值解，试证明拉方法得到的数值解，试证明y(xn)yn =0.5 a h xn若取若取 xn = nh，yn为用梯形公式计算所得的数值解，为用梯形公式计算所得的数值解，记记y(xn)为初值问题的在为初值问题的在x=xn处的解析解。试证明处的解析解。试证明:y(xn)=ynEx10.初值问题初值问题15/16Ex11.推导用于求解初值问题的推导用于求解初值问题的Adamas显格式显格式:16/16
