数值计算方法教案

《计算方法》教案

课程名称：计算方法

适用专业：医学信息技术

适用年级：二年级_______

任课教师：张利萍

编写时间：2011年8月

新疆医科大学工程学院张利萍

教案目录

《计算方法》教学大纲.................................4

一、课程的性质与任务........................................................4

二、课程的教学内容、基本要求及学时分配.....................................4

三、课程改革与特色..........................................................5

四、推荐教材及参考书........................................................5

《计算方法》教学日历.................错误!未定义书签。

第一章绪论...........................................6

第1讲绪论有效数字........................................................6

第2讲误差...............................................................

第二章线性方程组的直接法............................14

第3讲直接法、高斯消去法..................................................14

第4讲高斯列主元消去法....................................................22

第5讲平方根法、追赶法....................................................29

第三章插值法与最小二乘法...........................31

第6讲机械求积、插值型求积公式...........................................32

第7讲牛顿柯特斯公式、复化求积公式.......................................37

第8讲高斯公式、数值微分..................................................42

第9讲

第10讲

第12讲

第四章数值积分与数值微分...........................48

第11讲欧拉公式、改进的欧拉公式..........................................48

第12讲龙格库塔方法、亚当姆斯方法........................................52

第13讲收敛性与稳定性、方程组与高阶方程..................................56

第14讲

第15讲

第五章微分常微分方程的差分方法.....................59

第16讲迭代收敛性与迭代加速...............................................60

第17讲牛顿法、弦截法.....................................................64

第18讲

第19讲

第20讲

第六章线性方程组的迭代法...........................67

第21讲迭代公式的建立...................................................68

2

第22讲

第23讲

第24讲向量范数、迭代收敛性71

第25讲

3

《计算方法》教学大纲

课程名称：计算方法/ComputerNumericalAnalysisB

学时/学分：54/4

先修课程：高等数学、线性代数、高级语言程序设计（如：Matlab语言）

适用专业：计算机科学与技术、信息管理与信息系统

开课学院（部）、系（教研室）：医学工程技术学院、医学信息技术专业

一、课程的性质与任务

计算方法是一门专业必修课。当前，由于科学技术的快速发展和计算机的广泛应用，

学习和掌握计算机上常用的数值计算方法及有关的基础理论知识,并能用某种高级语言（如

Matlab语言）将这些常用算法编程实现，这对于计算机专业的学生来说是非常重要的。

本课程着重介绍进行科学建设所必须掌握的一些最基本、最常用的算法，向高等院校

有关专业的学生普及计算方法的知识.

二、课程的教学内容、基本要求及学时分配

（一）教学内容

1.引论

数值分析的研究对象、误差及有关概念、数值计算中应注意的一些原则。

2.线性代数方程组的数值解法

Gauss消去法、Gauss消去法的矩阵形式、主元消去法、三角分解法、迭代法、迭代法

的收敛条件及误差估计。

3.插值方法

Lagrange插值、Newton插值、分段插值、Hermite插值、三次样条插值、数据拟合的

最小二乘法。

4.数值积分与微分

机械求积、Newton-Cotes求积公式、复化求积、Romberg求积算法、Gauss求积公式、

数值微分。

5.常微分方程初值问题的数值解法

Euler方法及其改进、龙格-库塔（Runge-Kulta）方法、线性多步法、收敛性与稳定性、

一阶方程组与高阶方程。

6.方程求根的数值方法

二分法、迭代法、迭代过程的加速、Newton迭代法、Newton迭代法的几种变形。

（二）基本要求

1.了解数值分析的研究对象、掌握误差及有关概念。

2.正确理解使用数值方法求方程的解的基本思想、数学原理、算法设计。

3.了解插值是数值逼近的重要方法之一，正确理解每一种算法的基本思想、计算公式、

算法设计、程序框图设计和源程序。

4.掌握数值积分的数学原理和程序设计方法。

5.能够使用数值方法解决一价常微分方程的初值问题。

6.理解和掌握使用数值方法对线性方程组求解的算法设计。

（三）学时分配

本课程的理论教学时数为54学时分配如下表:

速学环节

课程辐、、、学时讲课

引论4

线性代数方程组的数值解法6

插值方法12

数值积分与微分10

常微分方程初值问题的数值解法10

方程求根的数值方法10

总复习2

合计54

（四）课程内容的重点、难点

重点：Lagrange插值、Newton插值、分段插值、Heimite插值、三次样条插值、机械

求积、Newlon-Coies求积公式、复化求积、Romberg求积算法。

难点：Gauss消去法、Gauss消去法的矩阵形式、主元消去法、三角分解法、迭代法、

迭代法的收敛条件及误差估计。

三、课程改革与特色

本课程是一门重要的专业基础课。数值计算方法既是一门古老的学科，又是一门新兴

的学科。电子计算机的产生和发展极大地促进了数值计算方法的发展。只有把数值计算方

法和程序设计紧密结合起来，把算法变为计算机能直接执行的程序，才能真正使计算机帮

助人们解决各种复杂的计算任务。

本课程试图将数值计算方法和程序设计方法学融为一体，这也是一种尝试。

四、推荐教材及参考书

推荐教材：《计算机数值方法》（第三版），主编：施吉林、刘淑珍、陈桂芝，出版社:

高等教育出版社，出版时间：2005年3月

参考书：

《数值计算方法和算法》，主编：张韵华、奚梅成、陈效群，出版社：科学出版社，出

版时间：2002年3月

《NumericalAnalysis》，主编：RichardL.Burden,出版社：高等教育出版社影印，出

版或修订时间：2003

《数值分析》，主编：金聪、、熊盛武，出版社：武汉理工大学出版社，出版时间：2003

年8月

5

第一章绪论

一、教学目标及基本要求

通过对本章的学习，使学生对了解涉及工程和科学实验中常见的数学问题，其中包括

线性方程组、函数插值、离散数据的拟合、微积分、微分方程等，这些问题是其他数学问

题的基础。

二、教学内容及学时分配

本章主要介绍数值分析的研究对象及误差的概念。具体内容如下：

第1-2学时讲授内容：计算方法的研究内容、对象与特点；误差的基本概念。

三、教学重点难点

1.教学重点：误差、误差种类：误差分析：误差与有效数字的关系.

2.教学难点：误差分析、误差与有效数字的关系。

四、教学中应注意的问题

多媒体课堂教学为主。适当提问，加深学生对概念的理解。

第1讲绪论

基本求解步骤

编程上机

计算结果

数学模型是通过科学实验或者观察分析一系列数据后，用数学作为工具近似地描述客观事

物的一种数学表达式。在数学模型中，往往包含了若干参量，这些物理参数通常由实验仪

器测得，根据仪器的精密程度，物理参数的确定也会产生一定的误差。

在建立了数学模型之后，并不能立刻用计算机直接求解，还必须寻找用计算机计算这些数

学模型的数值方法，即将数学模型中的连续变量离散化，转化成一系列相应的算法步骤，

编制出正确的计算程序，再上机计算得出满意的数值结果。

算法：从给定的已知量出发，经过有限次四则运算及规定的运算顺序，最后求出未知量的

数值解，这样构成的完整计算步骤称为算法。

计算多项式p(x)=31+4x2-2x+6的值。

算法1：由X计算出x',3后再进行计算。

需乘法5次，加法3次。

6

〃(x)=x[x(3x+4)—2J+6

需乘法3次，加法3次。

一般地，计算n次多项式的值

nnx

巴(%)=anx+an_xx~+…++4

P_,、«i〃(〃+i)

如若按4"1有14次乘法运算，计算K(x)共需"2++〃=1—次乘法和〃次加

法运算。

采用：秦九韶算法(1247)有递推公式：

%)=工(旧.《(『+%)+噎+.+4)+%从内往外一层一层计算，社巳表示第k层

以=(...(atlx+a,^)x+...+an_k.x)x+an_k

[匕=%TX+4T

Vo=4

需乘法n次，加法n次，存储单元n+3个。

对算法所要考虑的问题,包括如下：

•计算速度

例如，求解一个20阶线性方程组，用消元法需3000次乘法运算；而用克莱姆法则要进行

9.7X1O20次运算,如用每秒1亿次乘法运算的计算机要30万年。

7

•存储量

大型问题必要考虑计算机的数据存贮。

•数值稳定性

在大量计算中,舍入误差是积累还是能控制,这与算法有关。

实际算法往往表现为某种无穷递推过程

算法的精度控制

方程根的二分法求解

/(幻在［〃,切上单调连续，f(a)f(b)<0,根据连续函数性质，/(处在［©力内一定有实的零点,

即方程〃幻=0在［。,川内一定有唯一实根。解实根为/

若/(/)=0,则%为所求根

否则若f(a)J(Xo)<。，则根在区间［。，须)］，取q=x0

若/S)/(%)<0,则根在区间岛，勿，取4=Xo，b［=b

[a,b]n[«1，/?!]z>...n[ak,bk]Tt...

每一区间为前一区间的一半，有根区间［4,4］长度%一见=-(b-a)

2

,一(4

§1.2预备知识和误差

(1)误差的来源

实际问题"建立数学模型”研究计算方法》编程上机计算解结果。

模型误差：在建立数学模型过程中,不可能将所有因素均考虑,必然要进行必要的简化,这

就带来了与实际问题的误差。

测量误差：测量已知参数时，数据带来的误差。

截断误差：模型的准确解与某种数值方法的准确解之间的误差称为截断误差或方法误差。

舍入误差：计算机的字长是有限的,每一步运算均需四舍五入，由此产出的误差称舍入误

差。如：n.1/3,……取小数点8位、16位。

［截断误差的实例］

21一+

己知e"=1+x+—X4--"3+・..+------X十

2!7i!

求e-的近似值，并估计误差。

解：利用展开式的前三项，取n=2,

6T«14-(-1)4-1(-1)2=0.5

由Qy如公式:

/(x)=f(x0)+/'(x0)(x-x(l)+

+k(i。-即a-"

8

"+l

=Ov®<1

5+1)!

\R\=Ie-1-O.S|^—<1.7*IO-1

213!截断误差为：0.17

［舍入误差的实例］

1.492x1.066=1.590472,设在一台虚构的4位数字的计算机上计算

1.492x1.066«1.590,舍入误差为0.000472。

数值计算方法主要讨论截断误差和舍入误差的影响，不讨论模型误差和测量误差。

三、误差的基本概念

(1)误差与误差限

误差不可避免，设以工代表数K*的近似值，称《二”一/是近似值大的绝对误差。简称误

差。误差是有量纲的，可正可负。

误差通常是无法计算的，但可以估计出它的一个上界。即

卜一‘称£是近似值X的误差限，或称精度，即

**

X-8<x<x+8

O

(2)相对误差与相对误差限

e_x*—x

绝对误差并不能完全反应精度，称♦-X为近似值x的相对误差，记作。相对

误差是个相对数,是无量纲的,也可正可负。

相对误差的估计图",「，称£，为相对误差限，即

(3)有效数字

定义：如果近似值X的误差限是3(某一数位的半个单位)，则称X准确到小数点后n

位,并从第一个非零的数字到这一位的所有数字均为有效数字。

如：n=3.1415926535,

3.14有三位有效数字,误差限e=0.005;

3.1416有五位有效数字,误差限为0.00005o

(4)有效数字与误差限的关系：

x有n位有效数字，标准形式为.x=±KTx0.生生…%其中a(i=l,2,…)是0~9之间的

整数，且qW0,如果误差|x-x|<^xl(F"JV/V〃，称x为/的具有1位有效数值的

近似值.

(5)有效数字与相对误差的关系：

9

标准形式为x=±UTxO.q/…耳，则:

M

a)若寸有〃位有效数字，J_xio'-

kI2q

1.?10止"

品甯WxlO1

若邑？!<一!一xio〜

b)以12(4+1)，则x"有〃位有效数字

,n

证:lx-x*-----------x10l-nxx*|4―!—x10~x(a.+l)x10-'=-x10*”

112(4+1)2(q+l)'2

例，已知乃=3.14159265..,试问其近似值内=3.1,x2=3.14,x3=3.1415,A:3=3.1416

各有几位有效数字？并给出它们的误差限和相对误差限。

e,=|^-^|«0.04<^10十分位以前都是有效数字，有两位有效数字

1-2",

e\r-<—xlO=-xl0

2x36

♦=归一々核0.002<-xl02有三位有效数字

年一天|<^—XlO1-3=-xl0-2

2x36

3

|^--x3|«0.00009<^xl0-,有四位有效数字

<—!—xio1-4=-xio-3

，「

22x36

/=年一%|。0.00001<^xW4,有五位有效数字

^,叱小叱

例：为使二*的相对误差小于0.001%,至少应取几位有效数字?

解:

£「M」一XlO-1<0.001%

2al

〃>6—k)g6,即〃之6,取〃=6,则"*=3.14159

10

§1.3数值计算的若干原则

1,避免两相近数相减

当x较大时，计算工T-6，可先转化为'-6=-^J—尸

VX+I4-VX

/(x)=G&x=2得导数值f⑵X,2+%二J2i,精确值尸Q)=0.353553

2%

人k八［组，”\V2+^-V2^A1.4491-1.3784n_„.n

令h=0.1得/(2)=-----------------------«------------------------«0.35350

2h0.2

人」AAAA,田V2+A-V2^h1.4142-1.4142八

令h=0.0001得f(2)*-----------------------a------------------------=0

2h0.0002

计算"c°s*,x=l,分子出现相近数相减，可转换为

sinx

1—cosxsinxh、、a

—；----=--------,再计算

sinx1-cosx

2.避免绝对值太小的数做除数

分母接近零的数会产生溢出错误,因而产生大的误差,此时可以用数学公式化简后再做.

V=,___],_==ViooT+Viooo

Viooi-Viooo

3.要防止大数“吃掉”小数

计算机在进行算术计算时，首先要把参加运算的数对阶，即把两数都写成绝对值小于1,而

阶码相同的数。如：%=1。9+1必须改写成：x=0.1x1()1°+0.0000000001x1010如果

计算机只能表示8位小数，则算出xnO.lxlOio,大数吃掉了小数。这种情形是要尽量避

免的。

4.简化计算步骤，提高计算效率

简化计算步躲是提高程序执行速度的关键，它不仅可以节省时间，还能减少舍入误差。

例4：设A、B、C、D分别是10x20、20x50、50x1、1x100的矩阵，试按不同的算法求

矩阵乘积E=ABCD.

解：由矩阵乘法的结合律，可有如下算法

1.E=((AB)C)D.计算量N=11500flop

2.E=A(B(CD)).计算量N=125000flop

3.E=(A(BC))D.计算量N=2200flop

5.要使用数值稳定的算法

我们已经知道，所谓算法的稳定性，是指误差的传播可以得到控制，在用计算机解决实际

问题时，运算次数成千上万。如果误差的传播得不到控制，那么误差的累积会使问题的解

答成为荒谬的，尤其是某些病态问题(如病态方程组)，舍入误差对其计算结果往往有非常

严重的影响。因此，在选择计算方案时，要特别谨慎。

考察方程组

11

11

3~6

13

x解为X]=1,电=1,X=1

242n3

JX347

34560

四舍五入系数后，解为M=1.09,%=0-484，%3=149

尽管系数变动不大，但求出得解却变动很大，这类问题称为病态的。

例：蝴蝶效应（气象学家洛伦兹，1963）

——南美洲亚马孙河流域热带雨林中的一只蝴蝶翅膀一拍，偶尔扇动几下翅膀，可能在两

周后引起美国德克萨斯引起一场龙卷风？！

12

13

第二章插值法

一、教学目标及基本要求

通过对本章的学习，使学生掌握插值法计算常见的数学问题。

二、教学内容及学时分配

本章主要介绍数值分析的插值法。具体内容如下：

第3-4学时讲授内容：问题的提法、拉格朗日插值公式。第5-6学时讲授内容：插值

余项、牛顿插值公式。第7-8学时讲授内容：曲线拟合。

三、教学重点难点

1.教学重点：插值方法的由来、拉格朗日插值公式、牛顿插值公式、曲线拟合。

2.教学难点：拉格朗口插值公式、牛顿插值公式。

四、教学中应注意的问题

多媒体课堂教学为主。适当提问，加深学生对概念的理解。

第2讲拉格朗日插值公式

众所周知，反映自然规律的数量关系的函数有三种表示方法：

A.解析表达式

f(x)=x3-2x-5

(开普勒(Kepler)方程)％=>一esiny.。

悬链线方程：y=4cos(x")。

B.图象法

C.表格法

14

Xy

0.924-0.008513725

0.928-0.003822324

09320.000343434

0.9360.005532443

0.9400.012976643

1、插值法对于一组离散点（％J（X））,（，=0,1,2,...,〃），选定一个便于计算的简单

函数P（M，如多项式函数，要求尸㈤满足「区）=/（茗），由此确定函数P（幻作

为*幻的近似函数，然后通过处理P（幻获得关于《幻的结果。这就是插值方法。

2、曲线拟合选定近似函数P㈤时，不要求近似函数P（刈必须满足

尸（七）=/（匕），而只要求在某种意义下（最小二乘法原理），使近似函数尸⑶在

这些点上的总偏差量最小，这类方法成为曲线拟合。

§1.1多项式插值问题的一般提法

1插值法的概念：

假设函数尸f（x）是［46］上的实值函数，的用,…,为是5］上加1个互异的

点，f（x）在这些点上的取值分别为必，儿…,以

求一个确定的函数尸（才），使之满足：

产（%）二%（/=0,1,2,­,n）（1）

称沏为,“.,心为插值节点，关系式⑴称为插值原则，函数PG）称为函数y¥（x）

的插值函数，区间［为3称为插值区间。

2泰勒插值：

人们熟悉的泰勒展开方法其实就是一种插值方法，泰勒多项式：

23=/*0）+/'（元0）*-/）+-。0）2+.•.+―/）"（1）

'2I!T%）r,v.，％

与刈在点与邻近会很好的逼近f（x）o

泰勒余项定理：

定理1假设《刈在含有点%的区间［a,b］内有直到n+1阶导数，则当

例时，对于式（1）给出的匕⑴，成立

/（幻一X。严

（〃+1）!

15

其中J介于与与X之间，因而J£[〃,/?]。

所谓泰勒插值指下述问题：

问题1求作n次多项式月⑴,使满足？，„)=带),％=0,1,2,…刀，端为

一组已给数据。

易看出，上述插值问题的解就是泰勒多项式(l)o

例1例题分析：

求作力幻=正在/=100的一次和二次泰勒多项式，利用它们计算斤的近似

值并估算误差。

解：

l/23/2-5/2

fix)=4x,f'(x)=^x~ff"(x)=^-X~,/"W=|x

248

/xJ=10,/'(xo)=l/2O,/"(/)=-1/4000,/优)=3/8000000

yw=正在/=IOO的一次泰勒多项式是

6(X)=/(%)+,尸(/Xx-X(J=5+0.05X

7=115时Vil?=/(115)«^(x)=10.75

根据定理1可估计误差

22

|/(x)-Pi(x)|="(x-A0)<，；(x-x0)<0.028125<0.05

误差小于十分位的一半，故十分位及前面的数字为有效数字，所以结果有三

位有效数字。

修正R(x)可进一步得到二次泰勒公式

鸟(幻=《。)+^^。一%)2

VH5=/(115)«^(x)=10.75-0.028125=10.721875

,-,|/"'(X0)|--Q

3

,(x)—P2(x)|=l2。_/)<]-(X-x0)<0.0006328125<0.005

误差小于百分位的一半，故百分位及前面的数字为有效数字，所以结果有四

位有效数字。

泰勒插值是一种有效的插值方法，对函数要求严格(要足够光滑，存在高阶

导数)，要计算函数的高阶导数，而高阶导数的计算对计算机来说就很困难；

另外，计算过程不能形成机械重复的过程，不利于计算机程序实现。

§1.2拉格朗日(Lagrange)插值

1多项式插值的存在惟一性：

多项式导数易于计算，函数表达式简单，计算机易于计算，故考虑用多项式

函数彳型插值函数来模拟实色函数。

从如下数据表着手,并假定七。。《二/4

X：XoX\X2...尤〃

y-yoyiy?...y〃

16

求〃次多项式々（］）=〃0+&11+...+々〃］：使得：

P（x）=yi（2=0,1,2,•••,n）。

根据插值条件，有：

P（M）=%+。用+…+4石=%

P（$）=%+4%+…+ax；=%

<n

P区）二旬+4升+…+=yn（i）

显然，这是一个关于。。，弓-一〃〃的〃丹元线性方程组，其系数矩阵的行

列式为

1玉）…玉）

匕（/，2，…，5）=：）7

1士…<

/•1rxV（x,x,•••,%„）=n（x；-x；）0

注意到插值节点必"=1，2,…，〃）两两相异，而"ft0MX"'

故方程组（1）有惟一解，4，•••"〃，于是满足插值条件的多项式存在且惟一。

定理由加1个不同插值节点％，*1，•••，工〃可以惟一确定一个n次多项式

匕（元）=%+4工+・・・+满足插值条件Pn（N）=yo

从理论上说，由方程组（1）可以求出〃。，4,…〃〃的惟一解，从而确定?（回。但

从数值计算上看，当〃较大时求解线性方程组的工作量较大且不便应用。

解方程组（1）需计算n+1个n阶行列式，每个n阶行列式为n!项之和，每

项乂是n个元素的乘积，需n-1次乘法，所以求解需要（〃+1）〃!（〃-1）次乘法，

当n较大时，计算量非常大。

为解决此问题，现已提出了不少构造2（幻的巧妙办法。

2Lagrange插值的基函数构造法

首先讨论炉1时的情形。

已知*0,%，为，y,求乙（%）=%+“I］使得4（/）=%；4%）=X

显然4（X）是过（％），%）和（％，M）两点的一条直线。

由点斜式容易求得

17

L1(x)=y0+-—―(x-x0)

1x0-xJyx,-XJ/=o

V-------Y-------)V-------Y-------)

4G)AG)

其中，4(x),(，=0,l)具有如下特点：

7o(xo)=l;/o(x)=O

4(/)=。；4a)=i

称其为线性插值基函数。。(“)可以通过函数4(%)，("二°」)组合得出，且组

合系数恰为所给数据y0,y.o

再讨论声2时的情形。

显然4(%)是过(/，％)、(用，％)、(%，j2)三点的一条抛物线。

y

°Xo.V1.X2X

仿照线性插值基函数的构造方法，令

/(幻二(xfXxr)

0(x0-x,)(x0-x2)

/(x)^U-xQ)(x-x2)

a-5)a-x2)

/式©=(…。)*7

(工2-%)(%一%)

其中，4G),(i二°,L2)具有如下特点：

小/)=l;/o(x,)=O;/o(x2)=0

</1(xo)=O;/l(x1)=l;/](x2)=0

/2(x0)=0"2(%)=0;/2(x2)=1

称其为抛物重插值基函数(如下面所示)0

18

于是,

。一%）（了一马）,

L(X)=

2（x0-x,）（x0-x2）0

*一%）（尢一9）

（石一工0）。-x2）

+（二?「）斗⑴%

（X,一演）（无，一%）r=0

最后讨论一版情形。

求乙（而使得L（M）=%（7=0,1,2,-,/?）o

令〃诙多项式插值基函数为：

〃（）

43=—X-X4.

4.（x）,（i=0,19・・・,〃）具有如下特点:

l,i=j

4（勺）=%.=<【。"打

于是，满足插值条件的〃次多项式可以直接写为：

i=0j^i（X,—XJ»=0

j=01J

我们称£〃（x）为Lagrange多项式，4（“）其Lagrange插值基函数。

19

■给定%=>+1,/=0,1,2,3,4,5.下面哪个是&(x)的图像?

3插值余项

如图所示，其截断误差尺5)=f(x)-£，(x),称为Lagrange插值多项式的余

项。

20

定理假设F5)在［a,b］上有连续的直到小1阶导数，且在不同插值节点

%,玉，・・・，%〃取值为/(%)=%,Ln(x)是经过插值样点(%,%),"=0,1,…㈤

的Lagrange插值多项式，若引进记号：

q+1(X)=(X)(X_%…(X_当)=—苍)

1=0

则当勿时，有如下的误差估计：

4。)-fM-W-9~~n。一七)

(〃+1)!1-0

=—儒多“⑴”)

证明：因为此(若)=/(七)一4(%)=°«=0,1，….)

于是可假定凡(X)具有如下形式：

n

RnM=-x0)(x-X,)•--U-x„)=k(x)U(x-xr.)

1=0

将X看作(a,b)上的一个固定点，作辅助函数

9(f)=/(/)—Ln(t)—k(x)(Z-x0)(Z—X))•••(/—xn)

=/(0-4(f)-Mx)加一七)

i=O

容易看出，。⑺有乂％不…，毛共加2个相异零点，且在［a,b］上存在加1阶导数。

根据罗尔，“⑺在。⑺的两个零点之间至少有一个零点，故。'⑺在［a,b］上至少

有加1个零点。如此类推,“川)⑺在(&b)上至少有1个零点1使得

小〃+1)

产⑹=f向皤)_£片©_攵⑶%而口n("%)|y

ati=o

=0

n

注意到4是〃次多项式，4"%)三°;口”刈的首项为f叫

)=5+i)!

故力(〃+~=。o由上述方程解得

（”+l）（g）

&（幻=

5+1）!

产）⑸〃

凡㈤=

于是

21

4例题

例1己知函数尸f(x)的观察数据为

Xk-2045

yk51-31

试构造F(x)的拉格朗日多项式〃(力，并计算人一1)。

解先构造基函数

x(x-4)(x-5)__x(x-4)(x-5)

(-2-0)(-2-4)(-2-5)-84

(x+2)(x-4)(x-5)_(x+2)(x-4)(x-5)

(0-(-2))(0-4)(0-5)~40-

.、(x+2)x(”5)x(x+2)(x-5)

iwz=--------=-------

2(4+2)(4-0)(4-5)24

。+2»(升—2)。-4)_(x+2)x(x—4)

J(5+2)(5-0)(5-4)35-

所求三次多项式为

3

Z3(%)-JO

-5J("4X"5)(x+2)a-4)Q-5)

84+40

(_3X，(x+2)(x-9(x+2)x(x-4)

~24+35-

上1

-15421

-155,24

〃———+1——

4214217

第3讲牛顿公式

§1.4差商与差分及其性质

1差商的概念：

称％一不为函数f(x)的一阶差商；

/[XpXj-Zlx^xJ

」一--

/[x0,Xpx2]=--

称马一%为函数f(x)的二阶差商;

22

rrrrri_/[X]，…,―/[%,.•・,七・/

JL人0,人]，…，人〃J—一

一般地，称为函数F（x）的〃阶

差商；

特别地，定义力%］二/（与）为函数/U）关于先的零阶差商。

由此可知，高阶差商总是由比它低一阶的的两个差商组合而成。

2差商性质

（a）性质1"阶差商可以表示成加1个函数值为'''.•"的线性组合，

即

/K，・・・代］二互——V——、-------「——;

,=。（X,.-%。）（西一百）…（％"—％_|'）（七一苦+1）…（％—X,）

该性质说明：4阶差商/［%，%，…，%］计算是由函数值代为），/'（幻，…人天）线

性组合而。

如：f[xQ,x[9x2]=/[xpx0,x2]=/[x2,xpx0].

九外,无|］二/（“）一/。°）=f（'°）I/（"）

%一元0%一%%7。

/知豆,引=&1见二色闻

%7。

ZU；）-/Ui）_/U,）2/（小）

超一玉X,-A-_/（x）,/（x,）

-------------------0-------0------

七一•%七一%X一%

fgf（M）［

与一小

/（见）।।/（W）

（$一x,）（x0-x2）（%一及）（％一与）（x2-%）（左-%）

（b）性质2（对称性）:差商与节点的顺序无关。即

/区,3］=小"（）］,

这一点可以从性质1看出。

3利用差商表计算差商

利用差商的递推定义,可以用递推来计算差商。

差商表：

23

一阶差商二阶差商三阶差商

,八阳）

八%）

为/（巧）小”X」

/[XpX2]/[”0,“1，/]

工2f（x2）

/[x2,x3]/[x19x2,x3]

/（/）/[x0,xnx2,x3]

如要计算四阶差商,应再增加一个节点,表中还要增加一行。

4差分的概念

定义设函数尸/U）在等距节点为=%+^a=°』，…，"）上的函数值『（为）=£,

其中，力为常数称作步长。称

▽工可/1

九一九

/力汨也2尸2，-2

分别为F（x）在％处以力为步长的一阶向前差分，一阶向后差分和一阶中心差分。

称符号/、▽、6分别为向前差分算子，向后差分算子和中心差分算子。

f+-C--

在节点等距情况I,差商%用差分表示，设步长力=匕+1-匕，有

Jyxi，演+i）--------------7

“川一天八

//、/（苞+1，巧+2）一/（如演+1）1/AA、1A2

f5,x/+1，z+2）=-------------------------------=T7T（一△”•）=R△M

xj+2-Xj2h2h

一般形式（数学归纳法可证）

f5,xM,•.，,+«）=M

§1.5牛顿插值公式

1.牛顿插值公式的构造

24

Lagrange插值虽然易算，但若要增加一个节点时，全部基函数1人6都

需重新算过。本节介绍另外一种方法-牛顿插值法，并用它解决上面所述问题。

由线性插值

N|（x）=y（）+^―^-（x-x0）,令。0==~—=+a1（x-x0）

二次插值能否写成

N2(x)=al)+ai(x-x^+a^x-x^x-x^

由条件N2ao)=%，(再)=y,N2(X2)=y2得

为一）‘。y-y。

推广得

N“(x)=4+4(x-x0)+4(x—MX%-%)

+…+Q〃(X—拓)…(X-X〃T),

其中，〃。.…，〃〃为待定系数。如何求"o，4尸・.，4??

/卬引/⑴一小。）

所以/（工）=/（/）+/〔九，/］（%_%）（0）

〃X,Xo，xJ=

/[X,XO]=/[XO,X1]+/[X,XO,X1](X-X1)⑴

又

/n1=

x-x2

/(x,x0,x1]=/[x0,x1,x,]+/[x,x0,x1,x2](x-x2)⑵

一般地，,%]=/阮知不…''/一/[%，再，…，品

f[x,XQ,X]xn_J=f[xQ,玉xn]+f[x,xQ,X[，…，xn](x-xn)(n)

将式(n)代入式(n-1),...，式⑵代入式(1),式(1)代入式(0),

25

如此可得:

/(x)=/(x0)+/rx0,x1](x-x0)

+/[x(pX],](x—*0)(2-X[)+•••

+/[x0,xI,-,xn](x-x0)(x-x1)-(x-xzi_1)

+/[^x0,x1,-,rj(x-x0)(x-x1)-(x-x„)

尤为注意的是：最后一项中，差商部分含有X,乃是余项部分，记作此(X)；而

前面小1项中,差商部分都不含有X,因而前面加1项是关于X的〃次多项式,

记作N”(x),这就是牛顿插值公式。

2算例

例1:当n=l嘲，

f(x)=/(x0)+f[xQ,xJ(x-工0)+/[x,x0,x1](x-x0)(x-芭)

其中，’

^i(^)=/(x0)+/[x0,x1](x-x0)

=典+生*1。)

/一再

O

这就是牛顿一次插值多项式，也就是点斜式直线方程。

当n=2时，

/(x)=/(x0)+/[x0,Xj](x-x0)+/[x0,x1,x2](x-x0)(x-x1)

+/[x,x0,x1,x2](x-x0)(x-x1