基于Matlab的数值积分公式问题.doc

数值分析
学号:130080402015
学生所在学院:测试与光电工程学院
学生姓名:张翀
任课教师:郑华盛
教师所在学院:数信学院
基于Matlab的数值积分公式问题
张翀,测试与光电工程学院测试计量技术及仪器,130080402015
摘要:在求一些函数的定积分时,由于原函数十分复杂难以求出或用初等函数表达,导致积分很难精确求出,只能设法求其近似值,因此能够直接借助牛顿----莱布尼兹公式计算定积分的情形是不多的。数值积分就是解决此类问题的一种行之有效的方法。积分的数值计算是数值分析的一个重要分支;因此,探讨近似计算的数值积分方法是有着明显的实际意义的。
本文介绍了数值积分法的几种计算公式,如矩形求积公式、梯形求积公式和辛普森求积公式及相应的MATLAB命令,并给出了用 MATLAB编程求数值积分的实例。
关键词: MATLAB;数值积分;矩形求积公式;梯形求积公式;辛普森求积公式
目录
1引言..............................................................1
2数值积分算法介绍..................................................1
............................................1
................................................2
-科特斯求积公式.....................................5
..................................................7
3关于河流横断面积的数值积分问题....................................8
4问题的求解过程....................................................9
5基于MATLAB编程的各种求积公式对问题的求解.........................9
6总结.............................................................13
参考文献..........................................................14
附录...............................................................15
1 引言
实际问题当中常常需要计算积分。有些数值方法,如微分方程和积分方程的求解,也都和积分计算相联系。
在一元微积分学中,对于积分
,
只要找到被积函数f(x)原函数为F( x) ,求f(x)在该区间上的定积分便可用牛顿- 莱布尼兹公式求解,即
。
用牛顿- 莱布尼兹公式计算定积分的方法在理论上和解决实际问题中起到了很大的作用,但它并不能解决定积分计算的所有问题。在工程技术领域常遇到十分复杂的情况而无法用牛顿- :
某些被积函数f(x),其原函数无法用初等函数表示,如
, 等。
函数f(x)结构复杂,求其原函数非常困难。
函数f(x)的结构虽然简单且其原函数存在,但其原函数的结构相对复杂。
函数f(x)没有具体的表达式,只有一些由试验测试数据形成的表格或图形。
而在这些情况下,可采用“数值积分”的方法求出定积分(近似值) 。
2 数值积分算法介绍

大多数实际问题的积分是需要用数值积分方法求出近似结果的。数值积分原则上可以用于计算各种被积函数的定积分,无论被积函数是解析形式还是数表形式,其基本原理都是用多项式函数近似代替被积函数,用多项式的积分结果近似代替对被积函数的积分。由于所选多项式形式的不同,可以有许多种数值积分方法。而利用插值多项式来构造数值求积公式是最常用的一种方法。
对于积分,用一个容易积分的函数去代替被积函数,这样的自然以多项式为最佳,因为多项式能很好的逼近任何连续函数,而且容易求出其原函数。

在积分区间上取有限个点,作的次插值多项式,其中,为次插值基函数。用近似代替被积函数,
则得
()
若记
()
则得数值求积公式
()
其中称为求积系数,称为求积节点。则称该求积公式为插值型求积公式。
为了便于计算与应用,常将积分区间的等分点作为求积节点,这样构造出来的