《第五章数值微积分》由会员分享，可在线阅读，更多相关《第五章数值微积分（27页珍藏版）》请在金锄头文库上搜索。
1、通知通知数值分析上机考试时间改为数值分析上机考试时间改为2010年年1月月6日日(第第18周周三周周三)上上午午10:05-11:05共共60分钟分钟在图文信息中心在图文信息中心3号机房号机房引子引子定积分计算：用定积分计算：用原函数原函数( (不定积分不定积分) )无法找到原函数无法找到原函数F(x)怎么办？怎么办？第五章第五章数值微积分数值微积分 5.1数值积分公式数值积分公式5.2数值积分的余项数值积分的余项5.3复化求积法与步长的选取复化求积法与步长的选取5.4数值微分法数值微分法5.1 数值积分公式数值积分公式机械求积机械求积Newton-Cotes公式公式代数精度代数精度Gauss
2、求积公式求积公式1. 机械求积机械求积原理：定积分原理：定积分曲边梯形的面积曲边梯形的面积理论基础：积分中值定理理论基础：积分中值定理f( )=?矩形公式与梯形公式矩形公式与梯形公式左左矩形矩形右矩形右矩形中矩形中矩形梯形梯形机械求积一般公式机械求积一般公式 问题问题适当取适当取求积节点求积节点和和求积系数求积系数A0,An，计算函数值计算函数值f(x0),f(xn),近似近似解解误差误差T(f )-Q(f )2. Newton-Cotes公式公式 插值型求积公式插值型求积公式:P(x)是是f(x)的一个插值函数的一个插值函数(linear,Lagrange,Hermite,spline等等)
3、Newton-Cotes公式公式:采用等距节点采用等距节点Lagrange插值插值Cotes系数系数Ci(仅依赖于仅依赖于n, i)变量代换x=a+th低阶低阶Newton-Cotes公式公式梯形公式梯形公式(n=1)Simpson公式公式(n=2)Cotes公公式式 (n=4)数值稳定性数值稳定性:n8时，时，Cotes系数非负系数非负3代数精度代数精度 定义：定义：若机械求积公式若机械求积公式对所有对所有幂函数幂函数f(x)=1,x,x2xm准确，则称它具有准确，则称它具有m次次代数精度。代数精度。性质：若具有性质：若具有m次代数精度，次代数精度，则对所有次数则对所有次数不超过不超过m次的
4、多项式准确。次的多项式准确。代数精度：梯形公式代数精度：梯形公式(n=1)1次，次，Simpson公公式式(n=2)3次，次，Cotes公式公式 (n=4)5次。次。n为奇数时，为奇数时，n阶阶Newton-Cotes公式的代数精度为公式的代数精度为n；n为偶数时，为偶数时，n阶阶Newton-Cotes公式的代数精度为公式的代数精度为n+1。用代数精度构造插值公式用代数精度构造插值公式例题例题求求A1,A2及及x2，使求积公式使求积公式代数精度尽量高代数精度尽量高解：解：得得A1=1/4,A2=3/4，x2=2/34 Gauss求积公式求积公式 考虑将节点也视为待定参数，此时机械求积公式考虑
5、将节点也视为待定参数，此时机械求积公式的待定参数达的待定参数达2n+2个，从而可期望代数精度达到个，从而可期望代数精度达到2n+1，称此类高精度的求积公式为称此类高精度的求积公式为Gauss公式公式，而对应，而对应节点称为节点称为Gauss点点。一点一点Gauss(n=0)(中矩形公式中矩形公式)-1,1上的两点上的两点Gauss公式公式-1,1上的上的Gauss点点n n次次Legendre多项式多项式定理定理5.1-1,1上上n n阶阶Gauss点恰为点恰为n n次次Legendre多项式的根。多项式的根。nGauss点点求积系数求积系数代数精度代数精度002111，1325/9，8/9，
6、5/95一般区间一般区间a,b上的积分上的积分变换到变换到-1,1一般区间一般区间a,b上的两点上的两点Gauss公式公式5.2数值积分的余项数值积分的余项 引理引理5.1（积分中值定理积分中值定理）若）若f(x),g(x)均在均在a,b上连续且不变号上连续且不变号,则存在则存在a,b使使左矩形公式余项左矩形公式余项(证明证明:用用Taylor公式公式)中矩形公式余项中矩形公式余项(证明证明:用用Taylor公式公式)低阶情形低阶情形梯形公式余项梯形公式余项(证明证明:用积分中值定理用积分中值定理)Simpson公式余项公式余项(证明证明:用积分中值定理用积分中值定理+Hermite插值插值)
7、一般情况一般情况Newton-Cotes系列公式余项系列公式余项Gauss系列公式余项系列公式余项5.3 复化求积法与步长的选取复化求积法与步长的选取 复化求积原理复化求积原理定步长梯形法定步长梯形法2阶收敛性定步长定步长Simpson法法P119例例5.13(Simpson法法精度高精度高)4阶收敛性变步长梯形法变步长梯形法 递推关系递推关系逐级计算而在增加新节点时逐级计算而在增加新节点时,不浪费原先的不浪费原先的计算量计算量,并且可由并且可由|T2n(f) Tn(f)| 控制控制计算精度。计算精度。xi-1xixi-1/2Romberg公式公式 由|R2n(f) Rn(f)| 控制计算精度
8、 , j i= 1,2,由|Tii (f) Ti- 1 i-1 (f)| 控制计算精度. 自适应步长法自适应步长法 根据被积函数的陡缓自动选择局部步长根据被积函数的陡缓自动选择局部步长考虑某区间考虑某区间ak,bk,记记hk= bk ak,从从a,b开始按开始按 =|0.1(S=|0.1(S2 2-S-S1 1)|)| 检查精度检查精度 ，若满足精度则以，若满足精度则以S2为计算结果，否则分成两为计算结果，否则分成两个小区间各自重复逐步上述过程，每个小区间个小区间各自重复逐步上述过程，每个小区间精度用精度用 /2。这样重复下去，直至每个分段部。这样重复下去，直至每个分段部分达到相应精度（步长为
9、分达到相应精度（步长为h=(b-a)/2k时精度时精度 /2k）不同段的步长可能是不一样的，积分结果为每不同段的步长可能是不一样的，积分结果为每一小段积分的总和。一小段积分的总和。例例5 .15这里共使用了6个区间，调用函数13次，如果用等步长Simpson法达到该精度，需要调用函数17次。主要原因是自适应步长利用了函数的陡缓自动选择局部步长，变化快的地方细分，变化慢的地方粗分 5.4 数值微分法数值微分法 1差商法差商法向前差商向前差商公式公式向后差商向后差商公式公式中心差商中心差商公式公式差商公式比较差商公式比较中心差商精度比较高中心差商精度比较高插值型求导公式插值型求导公式 两点公式与三点公式两点公式与三点公式差商公式+余项
