数值积分矩形公式的复化及误差分析
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
, （6）
4矩形公式的递推化
虽然复化求积方法对提高精度是行之有效的,但在使用求积公式之前必须给出合适的步长,步长太长,精度难以保证,步长太小,又会导致计算量的增加．而事先给出一个合适的步长往往是困难的,那到底怎样选取步长才是合适的呢
实际计算中常常采用变步长的方案,即在步长逐次分半（即步长二分）的过程中,反复利用复化求积公式进行计算,直至所求得的积分值满足精度要求为止．
由于 在 上连续, 均为 的内点,所以由中值定理知,存在一点 ,使得 ,所以有
, （5）
称为复化左矩形公式的误差估计式 ,下标 表示将区间 划分为 等分．
3．2复化中矩形公式的误差估计公式
类似于复化左矩形的误差公式,同样可得复化中矩形公式的误差估计公式
其中 ,由于 在 上连续, 均为 的内点,所以由中值定理知,存在一点 ,使得 ,所以有
4．1左矩形公式的递推公式及误差
变步长过程中左矩形的计算规律：
将求积区间 分成 等份,则一共有 个分点,按左矩形公式计算 ,需要提供 个函数值,如果将积分区间再分一次,则分点增至 个,将二分前后两个积分值联系起来加以考虑,注意到每个子区间 经过二分只增加了一个分点 ,利用复化的左矩形公式求得该子区间上的积分值为 ,其中代 表二分前的步长,将每个子区间上的积分值相加得
（8）
根据中矩形公式的误差公式 得,积分值的截断误差大致与 成正比,因此当步长二分后,误差将减至原有误差的 ,即有 ,移项整理得 ,同样,当二分前后的两个积分值 与 相差很近时,就可以保证计算结果误差很小,积分近似值 的误差大致等于 ,因此,如果用这个 误差值作为 的一种补偿,则可以得到
可能结果比较理想．
Keywords:the rectangle formula;the left rectangle formula; error analysis; complex formula; the formula of recurrence
1引言
以前我们在进行积分运算时,都是先对被积函数求出其原函数,然后代值进行计算,但不是每个被积函数都是能轻易找到其原函数的,有的甚至找不到它的原函数,这就要求我们找出另外一种方法来研究积分运算．
首先我们来定义即将用到的左矩形公式和中矩形公式：
对于积分 ,由积分中值定理知, 一点 ,使得 =(b-a) ( )
（7）
从而根据左矩形公式的误差公式 得,积分值的截断误差大致与h成正比,因此当步长二分后,误差将减至原有误差的 ,即有 ,移项整理得 ,由此可见只要二分前后的两个积分值 与 相当接近,就可以保证计算结果误差很小,积分近似值 的误差大致等于 ,因此,如果用这个误差值作为 的一种补偿,可以期望得到的
可能是更好的结果．
5矩圆公式
由右图可见,这样分割后,形成一些小网格,以上一些工作我们就是
通过计算这些小的矩形条的面积之和进而估计出曲线 在 上所围
的面积．那么除此之外还有无别的近似计算方法呢
首先,我们试想,如右图所示,把网格顶端的一些剩下的不全的网格
近似为底为 ,高为 的三角形,那么前面我们按照左矩形公式
算得的矩形条的面积就为 ,整理后为 ,那么
所谓复化[1]就是指将一个积分的积分区间 划分为 等分,在每一个小区间 上应用左、中矩形公式求出积分值 ,然后对 求和,近似估计出积分 的积分值的算法．
2．1复化左矩形公式
将积分区间 划分为 等分,步长 ,分点
对每一个小区间 采用左矩形公式有
（3）
称为复化左矩形求积公式,下标 表示将区间 划分为 等分．
．若用区间左端点a的函数值 (a)作为 ( )的近似值,则得到我们熟悉的左矩形公式：
,其积分余项 （1）
．若改用区间中点 的函数值 作为 的近似值,则得到中矩形公式：
,其积分余项 （2）
由于我们导出的左矩形公式和中矩形公式对积分值的近似估计误差很大,所以我们采用复化求积公式来近似估计积分的准确值．
2复化公式
（11）
这就是我们所得到的左矩圆公式．
同理按照中矩形公式得到的小矩形条的面积为 ,用它来近似积分的准
确值可得到
（12）
这就是我们所得到的中矩圆公式．
例应用复化矩形公式（3）和（4）计算以及递推公式（7）与（8）和矩圆公式（11）与（12）计算积分 的近似值,并与其准确值作相应的比较．
解：设 , ,分点个数 ;
当 取不同值时各种算法对积分的估计值与近似解的比较
所用
公式
精确解
通过上表容易看出，当步长 逐渐变小时，不论是复化公式还是递推公式，它们对准确值的逼近效果都显著提高，即 越小，逼近效果越好；另一方面容易看出，中矩形公式比左矩形对准确值的近似程度更高，当然其复化公式的近似程度也比左矩形复化公式的精确度高；还有我们最后推出的(7)式与(8)式，它比起各自的复化公式来，逼近效果也相对较好，同样地中矩形公式的复化公式比(7)式的逼近效果要好．由(11)的结果可知在除 之外，它的计算结果是比较理想的；明显的问题是(12)式的计算结果与准确值的差距特别大，因为对于复化的中矩形公式而言，精确度已经是比较好的了，那么，如果我们再去作(12)式那样的逼近，势必导致出现大的波动．这样的公式是不完美的，所以对这样的公式完全可以舍去．只要我们选取合适的步长，分别利用它们各自的递推公式算出的近似值比它们自己的复化公式精确度要高很多．
(9)
这就是我们所熟知的梯形公式,而梯形公式对准确值的逼近程度要优于左矩形公式和右矩形公式,所以这样的假设与估计是成立的,同理我们对上面算得的中矩形公式也可以加上这个小三角形的面积而得到与准确值更为近的值．
(10)
下面我们来讨论另一种情形:即把上述所描写的三角形换成半径为 (或 ),为了计算方便,可以直接看成半径为 )的圆,那么按照前面我们推导得到的左矩形公式,计算得到小矩形条的面积为 ,用它来近似积分的准确值可得到
4．2中矩形公式的递推公式及误差
同理对中矩形公式也一样,将求积区间 分成 等份,则一共有 个分点,按中矩形公式计算 ,需要提供 个函数值,如果将积分区间再分一次,则分点增至2 +1个,将二分前后两个积分值联系起来加以考虑,注意到每个子区间 经过二分只增加了一个分点 ,在上述二分后的子区间上利用复化的中矩形公式求得该子区间上的积分值为 ,同样 代表二分前的步长,将每个子区间上的积小区间 采用中矩形公式,且令 ,则有
（4）
称为复化中矩形求积公式,下标 表示将区间 划分为 等分．[3]
3复化公式的误差分析
3．1复化左矩形公式的误差估计公式
由(1)式对每个小区间有误差估计式
其中 介于 , 之间,将上式代入(3)中则有
从而复化左矩形公式的误差估计式为
[6]张韵华，奚梅成，陈效群.数值计算方法与算法[M].2版.北京:科学出版社，2000.
Numerical integration of rectangular complex formula and error analysis
Zhang Xiaoxia
(CollegeofMathematicsandInformationSciencesNorthwestNormalUniversityLanzhouGansu730070)
[2]王仁宏.数值逼近[M].北京:高等教育出版社，1999.
[3]李岳生，黄友谦.数值逼近[M].北京:人民教育出版社，1978.
[4]李晓红，堵秀风，张永胜，王延臣.计算方法[M].北京:北京航空航天大学出版社，2006.
[5]马东升，雷永军.数值计算方法[M].2版.武汉:机械工业出版社，2001.
Abstract:Firstly, wededuced in therectangleformula,the left rectangle formula,and then carry out restoration ofthem. Becausethe result is notsatisfied with our expectations, so wetry to find theirrecursive formulas andanalyze theirerrors. Finally, wegiveanexample ofthe some appropriateformulas for the change of error;
数值积分矩形公式的复化及误差分析
张晓霞
（西北师范大学数学与信息科学学院甘肃兰州730070）
摘要：先推导得出中矩形公式、左矩形公式,然后对其进行复化，但由于结果不理想,再对两个公式进行递推,求出它们的递推化公式以及对其误差进行分析,最后举例说明几种逼近公式误差的变化情况;
关键词：中矩形公式,左矩形公式,误差分析,复化公式,公式递推化
6结论：
通过本文的论述，得出复化求积公式比原近似公式的精确度高；同样地，复化中矩形公式的逼近效果比复化左矩形公式对准确值的逼近效果好;另外，通过公式的递推化之后，我们得出递推化的公式比复化求积公式的精确度高;理所当然，随着 的不断增大，误差逐渐减小，当 到一定程度大时，会无限接近准确值．
参考文献：
[1]李庆杨，王能超，易大义.数值分析[M].4版.武汉:华中科技大学出版社，2006.
4矩形公式的递推化
虽然复化求积方法对提高精度是行之有效的,但在使用求积公式之前必须给出合适的步长,步长太长,精度难以保证,步长太小,又会导致计算量的增加．而事先给出一个合适的步长往往是困难的,那到底怎样选取步长才是合适的呢
实际计算中常常采用变步长的方案,即在步长逐次分半（即步长二分）的过程中,反复利用复化求积公式进行计算,直至所求得的积分值满足精度要求为止．
由于 在 上连续, 均为 的内点,所以由中值定理知,存在一点 ,使得 ,所以有
, （5）
称为复化左矩形公式的误差估计式 ,下标 表示将区间 划分为 等分．
3．2复化中矩形公式的误差估计公式
类似于复化左矩形的误差公式,同样可得复化中矩形公式的误差估计公式
其中 ,由于 在 上连续, 均为 的内点,所以由中值定理知,存在一点 ,使得 ,所以有
4．1左矩形公式的递推公式及误差
变步长过程中左矩形的计算规律：
将求积区间 分成 等份,则一共有 个分点,按左矩形公式计算 ,需要提供 个函数值,如果将积分区间再分一次,则分点增至 个,将二分前后两个积分值联系起来加以考虑,注意到每个子区间 经过二分只增加了一个分点 ,利用复化的左矩形公式求得该子区间上的积分值为 ,其中代 表二分前的步长,将每个子区间上的积分值相加得
（8）
根据中矩形公式的误差公式 得,积分值的截断误差大致与 成正比,因此当步长二分后,误差将减至原有误差的 ,即有 ,移项整理得 ,同样,当二分前后的两个积分值 与 相差很近时,就可以保证计算结果误差很小,积分近似值 的误差大致等于 ,因此,如果用这个 误差值作为 的一种补偿,则可以得到
可能结果比较理想．
Keywords:the rectangle formula;the left rectangle formula; error analysis; complex formula; the formula of recurrence
1引言
以前我们在进行积分运算时,都是先对被积函数求出其原函数,然后代值进行计算,但不是每个被积函数都是能轻易找到其原函数的,有的甚至找不到它的原函数,这就要求我们找出另外一种方法来研究积分运算．
首先我们来定义即将用到的左矩形公式和中矩形公式：
对于积分 ,由积分中值定理知, 一点 ,使得 =(b-a) ( )
（7）
从而根据左矩形公式的误差公式 得,积分值的截断误差大致与h成正比,因此当步长二分后,误差将减至原有误差的 ,即有 ,移项整理得 ,由此可见只要二分前后的两个积分值 与 相当接近,就可以保证计算结果误差很小,积分近似值 的误差大致等于 ,因此,如果用这个误差值作为 的一种补偿,可以期望得到的
可能是更好的结果．
5矩圆公式
由右图可见,这样分割后,形成一些小网格,以上一些工作我们就是
通过计算这些小的矩形条的面积之和进而估计出曲线 在 上所围
的面积．那么除此之外还有无别的近似计算方法呢
首先,我们试想,如右图所示,把网格顶端的一些剩下的不全的网格
近似为底为 ,高为 的三角形,那么前面我们按照左矩形公式
算得的矩形条的面积就为 ,整理后为 ,那么
所谓复化[1]就是指将一个积分的积分区间 划分为 等分,在每一个小区间 上应用左、中矩形公式求出积分值 ,然后对 求和,近似估计出积分 的积分值的算法．
2．1复化左矩形公式
将积分区间 划分为 等分,步长 ,分点
对每一个小区间 采用左矩形公式有
（3）
称为复化左矩形求积公式,下标 表示将区间 划分为 等分．
．若用区间左端点a的函数值 (a)作为 ( )的近似值,则得到我们熟悉的左矩形公式：
,其积分余项 （1）
．若改用区间中点 的函数值 作为 的近似值,则得到中矩形公式：
,其积分余项 （2）
由于我们导出的左矩形公式和中矩形公式对积分值的近似估计误差很大,所以我们采用复化求积公式来近似估计积分的准确值．
2复化公式
（11）
这就是我们所得到的左矩圆公式．
同理按照中矩形公式得到的小矩形条的面积为 ,用它来近似积分的准
确值可得到
（12）
这就是我们所得到的中矩圆公式．
例应用复化矩形公式（3）和（4）计算以及递推公式（7）与（8）和矩圆公式（11）与（12）计算积分 的近似值,并与其准确值作相应的比较．
解：设 , ,分点个数 ;
当 取不同值时各种算法对积分的估计值与近似解的比较
所用
公式
精确解
通过上表容易看出，当步长 逐渐变小时，不论是复化公式还是递推公式，它们对准确值的逼近效果都显著提高，即 越小，逼近效果越好；另一方面容易看出，中矩形公式比左矩形对准确值的近似程度更高，当然其复化公式的近似程度也比左矩形复化公式的精确度高；还有我们最后推出的(7)式与(8)式，它比起各自的复化公式来，逼近效果也相对较好，同样地中矩形公式的复化公式比(7)式的逼近效果要好．由(11)的结果可知在除 之外，它的计算结果是比较理想的；明显的问题是(12)式的计算结果与准确值的差距特别大，因为对于复化的中矩形公式而言，精确度已经是比较好的了，那么，如果我们再去作(12)式那样的逼近，势必导致出现大的波动．这样的公式是不完美的，所以对这样的公式完全可以舍去．只要我们选取合适的步长，分别利用它们各自的递推公式算出的近似值比它们自己的复化公式精确度要高很多．
(9)
这就是我们所熟知的梯形公式,而梯形公式对准确值的逼近程度要优于左矩形公式和右矩形公式,所以这样的假设与估计是成立的,同理我们对上面算得的中矩形公式也可以加上这个小三角形的面积而得到与准确值更为近的值．
(10)
下面我们来讨论另一种情形:即把上述所描写的三角形换成半径为 (或 ),为了计算方便,可以直接看成半径为 )的圆,那么按照前面我们推导得到的左矩形公式,计算得到小矩形条的面积为 ,用它来近似积分的准确值可得到
4．2中矩形公式的递推公式及误差
同理对中矩形公式也一样,将求积区间 分成 等份,则一共有 个分点,按中矩形公式计算 ,需要提供 个函数值,如果将积分区间再分一次,则分点增至2 +1个,将二分前后两个积分值联系起来加以考虑,注意到每个子区间 经过二分只增加了一个分点 ,在上述二分后的子区间上利用复化的中矩形公式求得该子区间上的积分值为 ,同样 代表二分前的步长,将每个子区间上的积小区间 采用中矩形公式,且令 ,则有
（4）
称为复化中矩形求积公式,下标 表示将区间 划分为 等分．[3]
3复化公式的误差分析
3．1复化左矩形公式的误差估计公式
由(1)式对每个小区间有误差估计式
其中 介于 , 之间,将上式代入(3)中则有
从而复化左矩形公式的误差估计式为
[6]张韵华，奚梅成，陈效群.数值计算方法与算法[M].2版.北京:科学出版社，2000.
Numerical integration of rectangular complex formula and error analysis
Zhang Xiaoxia
(CollegeofMathematicsandInformationSciencesNorthwestNormalUniversityLanzhouGansu730070)
[2]王仁宏.数值逼近[M].北京:高等教育出版社，1999.
[3]李岳生，黄友谦.数值逼近[M].北京:人民教育出版社，1978.
[4]李晓红，堵秀风，张永胜，王延臣.计算方法[M].北京:北京航空航天大学出版社，2006.
[5]马东升，雷永军.数值计算方法[M].2版.武汉:机械工业出版社，2001.
Abstract:Firstly, wededuced in therectangleformula,the left rectangle formula,and then carry out restoration ofthem. Becausethe result is notsatisfied with our expectations, so wetry to find theirrecursive formulas andanalyze theirerrors. Finally, wegiveanexample ofthe some appropriateformulas for the change of error;
数值积分矩形公式的复化及误差分析
张晓霞
（西北师范大学数学与信息科学学院甘肃兰州730070）
摘要：先推导得出中矩形公式、左矩形公式,然后对其进行复化，但由于结果不理想,再对两个公式进行递推,求出它们的递推化公式以及对其误差进行分析,最后举例说明几种逼近公式误差的变化情况;
关键词：中矩形公式,左矩形公式,误差分析,复化公式,公式递推化
6结论：
通过本文的论述，得出复化求积公式比原近似公式的精确度高；同样地，复化中矩形公式的逼近效果比复化左矩形公式对准确值的逼近效果好;另外，通过公式的递推化之后，我们得出递推化的公式比复化求积公式的精确度高;理所当然，随着 的不断增大，误差逐渐减小，当 到一定程度大时，会无限接近准确值．
参考文献：
[1]李庆杨，王能超，易大义.数值分析[M].4版.武汉:华中科技大学出版社，2006.