第5章 数值积分与数值微分方法
关于定积分计算,已经有较多方法,如公式法、分步积分法等,但实际问题中,经常出现不能用通常这些积分方法计算的定积分问题。怎样把这些通常方法失效的定积分在一定精度下快速计算出来,特别是通过计算机编程计算出来就是本章研究的内容。
此外,怎样根据函数在若干个点处的函数值去求该函数的导数近似值也是本章介绍的内容。
本章涉及的方法有Newton-Cotes求积公式、Gauss求积公式、复化求积公式、Romberg求积公式和数值微分。b5E2RGbCAP 5.1 引 例
人造地球卫星轨道可视为平面上的椭圆。我国的第一颗人造地球卫星近地点距离地球表面439km,远地点距地球表面2384km,地球半径为6371km,求该卫星的轨道长度。p1EanqFDPw 本问题可用椭圆参数方程
来描述人造地球卫星的轨道,式中a, b分别为椭圆的长短轴,该轨道的长度L就是如下参数方程弧长积分
但这个积分是椭圆积分,不能用解读方法计算。 5.2问题的描述与基本概念
要想用计算机来计算定积分是如下和式的极限
,应对其做离散化处理。注意到
要离散化,做
1) 去掉极限号
2) 将
取为具体的
值
用待定系数
代替
且任取
3) 为减少离散化带来的误差,将于是就得到
定义5.1 若存在实数
都有
<5.1)
则称式(5.1>为一个数值求积公式。式中求积节点,而称
称为求积系数,称为
<5.2)
为求积余项或求积公式<5.1)的截断误差。
从定义可以看到,数值求积公式依赖于求积节点个数n、求积节点和求积系数,这三个量有一个发生变化,则产生不同的求积公式。DXDiTa9E3d 定义5.2 若求积公式
对所有不超过m次的多项式某一个m+1次多项式代数精度为m。
有
有求积余项
,而对
,则称该求积公式的
一般,一个求积公式的代数精度越大,则该求积公式越好。 确定代数精度的方法
依次取
代入公式
并验证若第一个使为m-1。
是否成立。
不成立的k值为m,则对应的代数精度
例 5.1确定求积公式
的代数精度。 解 取
代入求积公式有
易验证
,但
本题求积公式代数精度为3。 例5.2确定下面求积公式
,故
的参数A,B,C,使它具有尽可能高的代数精度,并指出相应的代数精度。
解本题要先求出具体的求积公式,然后再判断所求公式的代数精度。
公式有3个待定参数,h不是求积公式的参数,故利用3个条件得到的3个等式关系就可以解决求出具体求积公式的问题。RTCrpUDGiT 依次取
代入求积公式并取等号,有
数值积分与数值微分课件
![](/skin/haowen/images/icon_star.png)
![](/skin/haowen/images/icon_star.png)
![](/skin/haowen/images/icon_star.png)
![](/skin/haowen/images/icon_star.png)