微积分的数值计算方法

微积分的数值计算方法

7.1 微积分计算存在的问题/数值积分的基本概念 微分计算问题

求函数的导数(微分)，原则上没有问题。当然，这是指所求函数为连续形式且导数存在的情形。但如果函数一表格形式给出，要求函数在某点的导数值；或者是希望某点的导数值只用其附近离散点上的函数值近似地表示，这就是新问题了，它称为微分的数值计算，或称为数值微分。 2．定积分计算问题

计算函数在上的定积分

当被积函数的原函数能用有限形式给出时，可用积分基本公式来计算：

然而，问题在于：① 的原函数或者很难找到，或者根本不存在；②可能给出一个函数表；③仅仅知道是某个无穷级数的和或某个微分方程的解等等。这就迫使人们不得不寻求定积分的近似计算，也称数值积分。 3．数值积分的基本形式

数值积分的基本做法是构造形式如下的近似公式

(7.1.1)

或记成 (7.1.2)

f ],[b a dx x f I b

a ⎰=)(f )(x F )()()(a F

b F dx x f I b

a -==⎰f f f ∑⎰=≈n

k k

k

b

a

x f A dx x f 0

)()(∑⎰=+=n

k n k k b

a f R x f A dx x f 0

][)()(

和 分别成为上的的数值求积

公式及其余项(截断误差)，和分别称为求积节点和求积系数(求积系数与被积函数无关)。

这种求积公式的特点是把求积过(极限过程)程转化为乘法与加法的代数运算。构造这种求积公式需要做的工作是：确定节点及系数，估计余项以及讨论

的算法设计及其数值稳定性。

4．插值型求积公式

如何构造求积公式呢？基本的技术是用被积函数的Lagrange 插值多项式近似代替，也即对上指定的个节点

及相应的函数值

，作

代入(7.1.2)式等号左边有

或写成形如(7.1.2)式的一般形式：

∑==n

k k k x f A I 0

*

)(][f R n ],[b a f k x k A ),,1,0(n k Λ=k x k A ),,1,0(n k Λ=][f R n *I f )(x L n f ],[b a 1+n b

x n ≤<⋯⋯<<≤10x x a )(,),(),(10n x f x f x f Λ)()()!

1(1)()()()()()1()1(0x f n x f x l x R x L x f n n k n

k k n n ++=++=+=∑ωξ⎰⎰⎰+=b

a n

b a n b a dx x R dx x L dx x f )()()(⎰∑⎰++=++

=b a

n n k n

k b

a k dx x x f n x f dx x l )())(()!

1(1)(])([)1()

1(0

ωξ∑⎰=+=n

k n k k b

a f R x f A dx x f 0

][)()(

(7.1.4)

其中 (7.1.5)

(7.1.6)

称(7.1.4)为插值型求积公式。下面将要介绍的几种实用求积公式无不都是插值型公式的某种具体形式。

由上述定义，用余项公式可以衡量数值求积公式的精确程度。不过，由余项公式(7.1.6)可见，如果为次数的多项式，则中有,故，从而公式(7.1.4)成为等式 ，这就是说，当被积函数为次数的多项式时，其相应的插值型求积公式不是近似公式，而是准确公式。据此，人们引入了另一个衡量数值求积公式近似程度好坏的“代数精度” 概念。

5. 代数精度

定义7.1.1 如果数值求积公式，当是次多项式(可表示为)时，均有

，也即有，便称求积公式

至少具有次代数精度；如果当是次多

⎰=b

a k k dx x l A )(),,1,0(n k Λ=⎰+++=

b a

n n n dx x x f n f R )())(()!

1(1][)1()

1(ωξf n ≤][f R n 0)1(=+n f 0][=f R n ∑⎰==n

k k k b

a x f A dx x f 0

)()(f n

≤∑⎰=≈n

k k k b

a x f A dx x f 0

)()()

(x f m ~0m x x x f ,,,1)(Λ=∑⎰==n

k k

k

b

a

x f A dx x f 0

)()(0][=f R n

∑=n

k k

k

x f A 0

)(m )(x f 1+m

项式(可表示为)时，只能是，

也即有 ，便称求积公式具有 次代数精

度。

定理 7.1.1 有个节点的插值型求积公式(7.1.4)，至少具有次代数精度。

例7.1.1 确定节点和系数，使得下列形式的求积公式

具有3次代数精度。

解 利用代数精度定义 取 令 计算得 即得至少具有3次代数精度的求积公式为

验证上式两端对，有 左边=2/5 ，右边=1/3

两者不相等，故可知所得公式正好为3次代数精度。

1

)(+=m x x f ∑⎰=≈n

k k k b

a x f A dx x f 0

)()(0][≠f R n ∑=n

k k k x f A 0

)(m 1+n n 321,,x x x A )]()()([)(3211

1x f x f x f A dx x f ++=⎰-32,,,1)(x x x x f =⎪⎪⎩⎪⎪⎨⎧++=++=++=++=⎰⎰⎰⎰----11

3

332313112

3222121

13211

1)()()

()111(1x x x A dx x x x x A dx x x x x A xdx A dx 21

,0,21,323

21==-==x x x A )]2

1

()0()21([32)(1

1f f f dx x f ++-

≈⎰-4)(x x f =