计算方法讲义：七 数值积分
计算方法数值积分
C1(2)
(1)1 2
4
21!1!0t(t2)dt6
类似可得，n=3时有四个Cotes系数
C 0 (3 ) 8 1 , C 1 (3 ) 8 3 , C 2 (3 ) 8 3 , C 3 (3 ) 8 1
n=4时，有五个Cotes系数
C 0 ( 4 ) 9 7 , 0 C 1 ( 4 ) 9 3 , 0 2 C 2 ( 4 ) 1 9 , 2 0 C 3 ( 4 ) 9 3 , 0 2 C 4 ( 4 ) 9 70
(k=0,1,2,…,4N)，h
ba N
4、复合Simpson公式算法
(1) 输入a,b,N (2) hba,sf(a),xa
2N
(3) 当 i=1,2, …,N时 做循环
① x=x+h ② s=s+4f(x) ③ x=x+h ④ s=s+2f(x)
(4) s h(s f (b)) 3
f()f(1 )f(2 ) f(n) N
1、复合梯形公式的余项
所以 IT Nb 1 a 2 h2f() , (a,b)
由 f (x) 在[a,b]上连续可知，f (x) 在[a,b]
上有界，于是存在常数M2，使
公式(5.6)称为等距节点内插求积公式。
求Ak
Ak
ablk(x)dx
b( n xxj a j0xk xj
)dx
jk
在等距节点前提下，做变换 t x a ，由 axb，可得 0tn 而x-xj=(t-j)h (j=0,1,2,…,n) ，xk-xhj=(k-j)h (j,k=0,1,2,…,n且j≠k)。 于是(5.5)式即为
A kk h !(( n 1 )n k k )! 0 n j n 0(tj)d t(b a )n(k !1 () n n kk)! 0 n j n 0(tj)dt
《数值积分方法》课件
数值积分的分类
按方法分类
可分为直接法和间接法。直接法如蒙特卡洛方法，间 接法如梯形法则、辛普森法则等。
按精确度分类
可分为低阶和高阶方法。低阶方法如梯形法则，高阶 方法如复合梯形法则、复合辛普森法则等。
按使用范围分类
可分为有限区间上的数值积分和无限区间上的数值积 分。
02
直接法
矩形法
总结词：简单直观
在金融建模中的应用
期权定价模型
数值积分方法可以用于求解期权定价模型，从而为金融衍生品定价提供依据。例如，二叉 树模型和蒙特卡洛模拟等。
利率衍生品定价
在利率衍生品定价中，数值积分方法可以用于求解利率期限结构模型，例如LIBOR市场模 型等。
风险管理
通过数值积分方法，可以对金融风险进行量化评估和管理。例如，计算VaR（风险价值） 和CVaR（条件风险价值）等指标，以评估投资组合的风险暴露程度。
自适应插值控制法
总结词
自适应插值控制法是一种通过插值技术来提 高数值积分精度的控制方法。
详细描述
在数值积分过程中，自适应插值控制法利用 插值技术对积分函数进行逼近，以提高数值 积分的精度。这种方法能够根据积分区间和 积分函数的特性，自动选择合适的插值方法 ，以获得更高的积分精度。同时，自适应插 值控制法还能够有效地处理复杂积分函数和
80%
算法设计与实现
数值积分方法的设计与实现是计 算数学的重要研究内容，推动了 科学计算的发展。
数值积分的概念
定义
数值积分是对函数在某个区间 上的定积分进行数值逼近的方 法。
思想
通过选取适当的积分点和权函 数，将定积分的计算转化为数 值逼近问题。
近似公式
常用的数值积分公式有梯形公 式、辛普森公式、复合梯形公 式、复合辛普森公式等。
数值积分课件 (《计算方法》)讲解共83页
36、殚精竭 虑将神 灵揉进 宪法， 总体上 来说， 法律就 会更好 。—— 马克·吐 温 37、纲纪废弃之日，便是暴政兴起之 时。— —威·皮 物特
38、若是没有公众舆论的支持，法律 是丝毫 没有力 量的。 ——菲 力普斯 39、一个判例造出另一个判例，它们 迅速累 聚，进 而变成 法律。 ——朱 尼厄斯
23、一切节省，归根到底都归结为时间的节省。——马克思 24、意志命运往往背道而驰，决心到最后会全部推倒。——莎士比亚
25、学习是劳动，是充满思想的劳动。——乌申斯基
谢谢！
40、人类法律，事物有规律，这是不 容忽视 的。— —爱献 生
21、要知道对好事的称颂过于夸大，也会招来人们的反感轻蔑和嫉妒。——培根 22、业精于勤，荒于嬉；行成于思，毁于随。——韩愈
第7章 数值积分
第七章 高斯数值积分法对于等参单元推导载荷列阵和刚度矩阵时，需计算如下形式的积分：其中被积函数一般比较复杂，甚至得不到显式。
因此，通常采用数值积分代替函数积分，即在单元内部选取某些点，先计算被积函数在这些点的函数值，然后用这些系数（称为加权系数，简称权）乘上这些函数值，再求总和作为近似积分值。
在有限元法中通常采用精度较高的高斯数值求积分法。
首先介绍一维高斯求积公式式中，()k f ξ是被积函数f 在积分点k ξ处的函数值；k w 是加权系数；n 是所选积分点的数目。
例如取一个积分点01=ξ（此时即1=n ），该点的函数值为1f （如图4.9 (a))，并取加权系数21=w ，则积分这是一种最简单的计算方法，只有当函数()ξf f =是一条直线时，即()ξf f =线之下是一个梯形才是精确的，若()ξf f =是任意曲线，则此计算结果是相当粗糙的。
为了改善精度，在11+≤≤-ξ范围内，取两个对称点1ξ，2ξ其函数值分别为()1ξf 和()2ξf 如图7.1（b ），但是横坐标1ξ、2ξ以及相应的权1w 和1w 需要确定。
为此设()ξf 为三次式，即则而由高斯求积公式于是由式（c ）和（d ）两式得即为了在3210,,,c c c c 取任意数值时式（d ）都是精确的，因此上式两边对应的系数必须相等，则有因此解得实根值得说明的是，上面确定的两个积分点的高斯求积公式（d ）对于被积函数是四次以下（不包括四次）的多项式是完全精确的，否则是近似的表达式。
另外，如图7.1（b ）所示，用两个矩形面积来表示函数()ξf 在区间[—1，十1]与轴ξ所围的面积，这就是式（d ）的几何意义。
图7.1 被积函数f 在积分点处的数值以相同的方法可以处理由3个函数值所组成的近似积分，如图7.1（c ）。
对不同的积分点数可确定相应的积分点坐标和加权系数，由此构成高斯积分表，见表7.1。
下面讨论二维、三维的高斯求积公式，对于二重积分可先对ξ积分，而把η视为常量，此时引入一维的高斯求积公式，则有再对η积分有将式（e ）代入式（f ），则可得二维的高斯求积公式用相同的方法可以导得三维的高斯求积公式在实际计算中，为了保证计算精度，并且不过分增加计算工作量，高斯积分中的积分点数n 通常可根据等参单元的节点数来选取，对于讨论的平面8节点等参单元和空间20节点等参单元都可以取3=n 。
数值积分PPT讲稿
),
0,
只要
~ f ( xk ) fk
(k 0,1,, n)
构造求积公式，原则上是一个确定参数
xk
和 A的k 代数问题.
b
n
f ( x)dx
a
Ak f ( xk ),
k 0
11
例 求a,b,c的值使下列求积公式的代数精度 达到最高。
1
f ( x)dx a f (1) bf (0) cf (1) 1
12
3. 插值型的求积公式
设给定一组节点
a x0 x1 x2 xn b,
x
初等函数表示的原函数； （2）当 f是(x由) 测量或数值计算给出的一张数据表.
这时，牛顿-莱布尼茨公式也不能直接运用. 因此有必要研究积分的数值计算问题. 由积分中值定理知，在积分区间 [a内, b存]在一点ξ，
成立
b
a f (x)dx (b a) f ( ),
3
就是说，底为 b 而a高为 f的(矩)形面积恰等于所求
曲边梯形的面积 I(图4-1).
图4-1
4
问题在于点ξ的具体位置一般是不知道的，因而难以
准确算出 f (的)值. 将 f 称(为) 区间
上[a的, b平]均高度.
这样，只要对平均高度 f提(供) 一种算法，相应地便
获得一种数值求积方法.
用两端点“高度“ f (a与) f的(b算) 术平均作为平均高度
j 0
15
注意到 lk (x j ) kj, 上式右端实际上等于 Ak , 因而
b
Ak a lk (x)dx
成立. 这样，有下面定理.
定理1
求积公式至少有 次代n数精度的
充分必要条件是，它是插值型的.
数值计算方法教案数值积分(有添加哦
数值积分教案教学目标：1. 理解数值积分的概念和意义；2. 掌握数值积分的基本方法和原理；3. 能够运用数值积分解决实际问题。
教学内容：1. 数值积分的概念和意义；2. 数值积分的基本方法：梯形法、辛普森法、高斯法等；3. 数值积分的原理：数值积分近似解的误差估计；4. 数值积分的应用：解决实际问题，如物理、工程等领域中的积分计算。
教学方法：1. 讲授法：讲解数值积分的概念、方法和应用；2. 案例分析法：分析实际问题，引导学生运用数值积分解决；3. 练习法：让学生通过练习题巩固所学知识。
教学准备：1. 教案、PPT、教学视频等教学资源；2. 计算器、电脑等教学工具。
教学过程：一、导入（5分钟）1. 引入数值积分的重要性，例如在物理、工程等领域中的应用；2. 引导学生思考如何利用数值方法近似计算积分值。
二、数值积分的概念和意义（10分钟）1. 讲解数值积分的定义；2. 解释数值积分的意义和作用；3. 举例说明数值积分在实际问题中的应用。
三、数值积分的基本方法（10分钟）1. 介绍梯形法、辛普森法和高斯法等基本方法；2. 讲解各种方法的原理和步骤；3. 通过实例演示数值积分的计算过程。
四、数值积分的原理（10分钟）1. 介绍数值积分近似解的误差估计；2. 解释误差估计的原理和意义；3. 引导学生思考如何选择合适的数值积分方法以减小误差。
五、数值积分的应用（10分钟）1. 分析实际问题，引导学生运用数值积分解决；2. 让学生通过练习题巩固所学知识；3. 引导学生思考数值积分在实际工程中的应用和限制。
教学评价：1. 课堂问答：检查学生对数值积分的概念和方法的理解；2. 练习题：评估学生对数值积分的应用能力；3. 课后作业：巩固学生对数值积分的掌握程度。
数值积分教案数值积分(有添加哦)六、梯形法的改进与应用（10分钟）1. 分析梯形法的局限性，如计算量大、精度低等问题；2. 介绍梯形法的改进方法，如自适应梯形法、辛普森法与梯形法的组合等；3. 通过实例讲解改进方法的原理和应用。
数值积分
b f(x)d x b L (x)d x b a [ f (a) 4 f ( a b ) f (b)] S a a 2 6 2
称为抛物线求积公式或Simpson公式. 几何意义:用抛物线围成的曲边梯形的面积代替围成的 曲边梯形面积
Newton-Cotes公式
将区间[a,b]n等分，其分点为xi=a+ih ， i=0,1,2,,n ， h=(b-a)/n，以这n+1个等距分点 为插值基点,作n次值多项式
n
称为牛顿—柯特斯（Newton-Cotes）系数.
当n=1时， Newton-Cotes公式(2)为梯形求积公 式
b
a
ba f ( x)dx [ f (a ) f (b )] T 2
当n=２时， Ne
b
a
ba ab f ( x)dx f (a ) 4 f 2 f (b ) S 6
xk 1 f(x)d x xk 1 xk [ f ( x ) f ( x )] x k k 1 2 k
则
a f(x)d x
b
n1 k 0
xk 1 f(x)d x n1 xk 1 xk [ f ( x ) f ( x )] x k k 1 k 0 2 k
Ln ( x)
b a
n i0
f ( xi )li ( x)
n b a i i0
f ( x)dx 求积系数
( l ( x)dx) f ( xi )
b Ai a li ( x)dx, (i 0,1,2,, n)
Newton-Cotes系数
作变量替换x=a+th，于是
计算方法 7 数值积分的基本方法
计算方法（2016/2017 第一学期） 西南科技大学 制造科学与工程学院
9
插值型求积公式
b
b
bn
f ( x )dx P( x )dx
a
a
a
f ( xk )lk ( x
n
f ( xk ) a lk ( x )dx
f ( xk )Ak
k 0
定义（代数精度） 设求积公式（1）对于一切次 数小于等于m的多项式
f ( x) 1, x, x 2 , , xm
或 f ( x) a0 a1 x a2 x2 am xm
是准确的，而对于次数为m+1的多项式是不准确 的，则称该求积公式具有m次代数精度（简称代 数由精定度义）可知，若求积公式(1)的代数精度为n，则 求积系数 A应k 满足线性方程组：
计算方法（2016/2017 第一学期） 西南科技大学 制造科学与工程学院
12
插值型求积公式
A0 A1 An b a
b2 a2
A0 x0 A0 x0n
A1 x1 A1 x1n
An xn
An xnn b
2
n1 n
20
本科毕设
2016年题目： ➢ 内啮合齿轮泵的设计 ➢ 一种移动式升降平台的设计 ➢ 木板贴纸生产线规划及上料装置设计 ➢ 通用小型斗式提升机设计 ➢ 码垛机结构设计 ➢ 齿轮泵的有限元力学分析
办公室地址：科技园（东9）A栋5-56
计算方法（2016/2017 第一学期） 西南科技大学 制造科学与工程学院
R( f=) 0,求积公式(1)能成为准确的等式。由于闭区间
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第七章 数值积分如果函数f(x)在区间[a,b]上连续，且原函数为F(x)，则可用牛顿―莱布尼兹公式：)()()(a F b F dx x f b a-=⎰来求得定积分。
然而很多函数无法用牛顿―莱布尼兹公式求积分。
一个简单被积函数，例如，其不定积分可能很复杂，见下面的MA TLAB 实例： >> syms a b c x>> int(sqrt(a+b*x+c*x*x),x)ans=1/4*(2*c*x+b)/c*(a+b*x+c*x^2)^(1/2)+1/2/c^(1/2)*log((1/2*b+c*x )/c^(1/2)+(a+b*x+c*x^2)^(1/2))*a-1/8/c^(3/2)*log((1/2*b+c*x)/c^(1/2)+(a+b*x+c*x^2)^(1/2))*b^2所以有必要研究简单、高效的计算定积分的方法（即数值积分方法）。
数值积分的基本思想是构造一个简单函数P n (x )来近似代替被积分函数f (x )，然后通过求⎰ba n dx x P )(得⎰ba dx x f )(的近似值。
7.1 插值型求积公式设⎰=ba dx x f I )(*，插值型求积公式就是构造插值多项式P n (x )，使⎰=≈ba n dx x P I I )(*。
构造以a ，b 为结点的线性插值多项式)()()(1b f ab ax a f b a b x x P --+--=，[])()()(21)()()(1b f a f a b dx b f a b a x a f b a b x dx x P T ba ba +-=⎥⎦⎤⎢⎣⎡--+--==⎰⎰称为梯形公式。
以a , 2ba c +=，b 为三个插值节点，构造二次插值多项式)())(())(( )())(())(()())(())(()(2b f c b a b c x a x c f b c a c b x a x a f b a c a b x c x x P ----+----+----=，则可以推出)()()()(2102b f c f a f dx x P S baλλλ++===⎰，)(61))(())((0a b dx b a c a b x c x ba-=----=⎰λ，)(64))(())((1a b dx b c a c b x a x ba-=----=⎰λ，)(61))(())((2a b dx c b a b c x a x b a -=----=⎰λ。
由此得公式：[])()(4)(6b fc f a f ab S ++-=，称为辛卜生（Sinpson ）求积公式。
根据经典拉格朗日插值公式)()()(0k nk k n x f x l x P ∑==，代入求定积分则有)()()()(0k nk b a k k nk ba k x f dx x l dx x f x l I ⋅==∑⎰∑⎰==，引入记号dx x l ba k k )(⎰=λ，)(0k nk k x f I ∑==λ，λk 为求积系数，x k 为求积节点。
注意：一积分结果为函数值的一个代数和，二是ab dx x l nk ba k -=∑⎰=0)(。
如果积分区间比较大，直接使用上述求积公式精度难以保证。
可对f (x )用分段抛物插值。
通常采取的办法是复化求积方法： （1）等分求积区间，比如取步长nab h -=，分[a, b]为n 等分，分点为kh x x k +=0，k = 0, 1, 2,…, n 。
（2）在区间 [x k , x k+1]上使用以上求积公式求得I k 。
（3）取和值∑-==10n k k I I 作为整个区间上的积分值。
将梯型公式和辛卜生公式应用于各子区间[]1,+k k x x ()1,,0-=n k 上得到子区间的定积分，再将子区间的定积分加起来得到整个区间的定积分近似值，相关公式称为复化梯型公式和复化辛卜生公式。
相对于复化梯型公式，复化辛卜生公式是一种精度较高的求积公式。
例如对于复化梯型公式，令[])()(21++=k k k x f x f hI ，则∑-==10n k k n I T [])()(2110+-=+=∑k k n k x f x f h ⎥⎦⎤⎢⎣⎡++=∑-=)()(2)(211b f x f a f h n k kx k 0 1/8 1/4 3/8 1/2 5/8 3/4 7/8 1 f (x k ) 4 3.93846 3.7647 3.50685 3.2 2.8764 2.46 2.26549 2 例 利用数据表计算积分 dx xI ⎰+=12*14( 1415926.3|arctg 410*===πx I )。
解：取n = 8用复化梯形公式：()13899.31872432852212832412812)0(21818=⎥⎦⎤+⎪⎭⎫⎝⎛+⎪⎭⎫ ⎝⎛+⎪⎭⎫ ⎝⎛+⎢⎣⎡⎪⎭⎫ ⎝⎛+⎪⎭⎫ ⎝⎛+⎪⎭⎫ ⎝⎛+⎪⎭⎫ ⎝⎛+⨯=f f f f f f f f f T 7.2 变步长梯形方法使用复化求积公式须给出合适的步长，步长太大精度难保证，步长太小会增加计算量，事先给出一个合适的步长是十分困难的。
递推公式避免了老节点的重复计算，使计算量减少了一半。
变步长积分法思想是将区间逐次对分，比较前后两次计算结果，若满足精度要求就停止，否则再次对分，直到到达精度要求为止。
设将区间[a, b] n 等分，共有n+1个分点，按复化梯形公式计算T n ，需要计算n+1个f (x )的值。
T 2n 的全部分点中有n+1个是原有的点。
小区间[x k , x k +1]经过二分增加分点21+k x 后，用复化梯形公式得积分为：[])()(2)(4121++++=k k k k x f x f x f hI ，因此有[][]∑∑∑∑-=+-=+-=+++-=+=++=++=110101112)(221)(2)()(4)()(2)(4212121n k k n n k k n k k k k k k n k n x f h T x f h x f x f h x f x f x f hT例 计算。
解：根据梯形公式和复化梯形公式2/)]1()0([1f f T +=，∑-=++=15.02)(221n k k n n x f h T T ，于是有n 1 2 4 8 16 32 T n 0.9397 0.9445 0.9456 0.9459 0.9461 0.94617.3 求积公式的误差P n (x )是f (x )的n 次插值多项式，当)(x f 本身就是次数不超过n 的多项式时)()(x P x f n ≡，求积公式)()()(*0*k nk k ba n ba x f dx x P dx x f I ∑⎰⎰====λ是精确的。
由于a b nk k-=∑=0λ，若f (x k )的舍入误差小于ε ，则)()()()()(*0**a b x f x f x f x f I I k k nk k k nk k k nk k-<-≤-=-∑∑∑===ελλλ。
所以舍入误差对数值积分的影响不大。
应用插值多项式余项定理)()!1()()(1)1(x n f x R n n +++=ωξ，对于插值多项式次数为1的情况有：))((2)()(b x a x f x R --''=ξ，可以证明梯形公式的截断误差为：3*)(12)())((2))((a b f dx b x a x x f T I ba-''=--''=-⎰ξξ-（注意：这里用到了积分中值定理：设)(x f 在区间[a,b]上连续，)(x g 在[a,b]区间上可积且不变号，则在[a,b]区间上至少有一个ξ满足⎰⎰=babadx x g f dx x g x f )()()()(ξ）。
将[a, b]区间n 等分，取nab h -=考虑复化梯形积分公式的误差，这个误差是n 个等分区间段上得误差之和，即12)()(-12)(-)(12-12)(-2310313*h f a b h f n f h h f T I n i i n i i n ξξξξ''-=''=''=''=-∑∑-=-=（注意：这里用到了介值定理，即对于连续函数)(x f 和自然数n ，存在ξ使。
7.4 收敛条件及收敛加速梯形法简单，但精度低，收敛的速度慢。
如何提高收敛速度？设I 是精确积分值，根据复化梯形公式的余项表达式可知：),(12)()(-2*b a h f a b T I n ∈''-=-ηη，，),(12)2/)(()(-*2*2*b a h f a b T I n ∈''-=-ηη，。
假定)()(*''≈''ηηf f , 则有41*2*=--n n T I T I 。
整理得：)(3122*n n n T T T I -=-。
可见只要二分前后T n 与T 2n 相当接近，就可以保证T 2n 的误差很小。
T 2n 的误差大致等于)(312n n T T -，用误差值作为T 2n 的补偿，可期望所得到的)(3122*n n n T T T I -+=，可能是更好的结果。
也可以这样考虑，将所有T n 看做构成一个函数T ，变量是h 2。
当h 趋近0时，T(h 2)接近I *，即T(0)=I *，连接点),2n n T h （和),121++n n T h （得一直线，其方程为：1221221221++++--+--=n nn n n n n n T h h h x T h h h x y ，延伸该直线与Y 轴相交（见下图），3-4T ~11221221221n n n n n n n n n n T T T h h h T h h h +++++=--+-=。
这就是一种迭代的加速。
图 一种迭代加速7.5 高斯型求积公式在插值型求积公式中，插值节点是事先固定的，有时还进一步限定是等距的。
是否可以在[a, b]上自由选择节点的位置，使精度提高？称：)()(0i ni k ba x f A dx x f ∑⎰=≈为一般求积公式，这里A k 为不依赖f (x )的常数。
若对任意不高于m 次的多项式精确成立，而对于x m +1不能精确成立，就说具有m 次代数精确度。
下面讨论最高代数精确度的求积公式（叫做高斯型求积公式）。
例 求形如 ⎰-+≈111100)()()(x f A x f A dx x f 的两点求积公式。
本题的解法很多，结果也不一定相同。