数值分析习题(含答案)

合集下载

数值分析试题与答案

一、单项选择题（每小题3分,共15分）1. 和分别作为π(de)近似数具有（）和（）位有效数字. A ．4和3 B ．3和2 C ．3和4 D ．4和42. 已知求积公式()()211211()(2)636f x dx f Af f ≈++⎰,则A ＝（）A ． 16B ．13C ．12D ．233. 通过点()()0011,,,x y x y (de)拉格朗日插值基函数()()01,l x l x 满足（）A ．()00l x ＝0,()110l x =B ．()00l x ＝0,()111l x =C ．()00l x ＝1,()111l x = D ．()00l x ＝1,()111l x =4. 设求方程()0f x =(de)根(de)牛顿法收敛,则它具有（）敛速.A ．超线性B ．平方C ．线性D ．三次5. 用列主元消元法解线性方程组1231231220223332x x x x x x x x ++=⎧⎪++=⎨⎪--=⎩ 作第一次消元后得到(de)第3个方程（）.A ．232x x -+=B ．232 1.5 3.5x x -+=C ．2323x x -+=D ．230.5 1.5x x -=-二、填空题（每小题3分,共15分）1. 设TX )4,3,2(-=, 则=1||||X ,2||||X = .2. 一阶均差()01,f x x =3. 已知3n =时,科茨系数()()()33301213,88C C C ===,那么()33C = 4. 因为方程()420x f x x =-+=在区间[]1,2上满足 ,所以()0f x =在区间内有根.5. 取步长0.1h =,用欧拉法解初值问题()211y y yx y ⎧'=+⎪⎨⎪=⎩(de)计算公式 .0,1,2分人三、计算题（每题15分,共60分）1. 已知函数211y x =+(de)一组数据：求分段线性插值函数,并计算()1.5f (de)近似值.1. 解 []0,1x ∈,()1010.510.50110x x L x x --=⨯+⨯=---[]1,2x ∈,()210.50.20.30.81221x x L x x --=⨯+⨯=-+--所以分段线性插值函数为()[][]10.50,10.80.31,2x x L x x x ⎧-∈⎪=⎨-∈⎪⎩ ()1.50.80.3 1.50.35L =-⨯=2. 已知线性方程组1231231231027.21028.35 4.2x x x x x x x x x --=⎧⎪-+-=⎨⎪--+=⎩（1）写出雅可比迭代公式、高斯－塞德尔迭代公式；（2）对于初始值()()00,0,0X =,应用雅可比迭代公式、高斯－塞德尔迭代公式分别计算()1X （保留小数点后五位数字）.1.解原方程组同解变形为1232133120.10.20.720.10.20.830.20.20.84x x x x x x x x x =++⎧⎪=-+⎨⎪=++⎩雅可比迭代公式为()()()()()()()()()1123121313120.10.20.720.10.20.830.20.20.84m m m m m m m m m x x x x x x x x x +++⎧=++⎪⎪=-+⎨⎪=++⎪⎩(0,1...)m =高斯－塞德尔迭代法公式()()()()()()()()()1123112131113120.10.20.720.10.20.830.20.20.84m m m m m m m m m x x x x x x x x x ++++++⎧=++⎪⎪=-+⎨⎪=++⎪⎩ (0,1...)m =用雅可比迭代公式得()()10.72000,0.83000,0.84000X =用高斯－塞德尔迭代公式得()()10.72000,0.90200,1.16440X =3. 用牛顿法求方程3310x x --=在[]1,2之间(de)近似根（1）请指出为什么初值应取2 （2）请用牛顿法求出近似根,精确到. 3. 解()331f x x x =--,()130f =-<,()210f =>()233f x x '=-,()12f x x ''=,()2240f =>,故取2x =作初始值4. 写出梯形公式和辛卜生公式,并用来分别计算积分111dxx+⎰.四、证明题（本题10分）确定下列求积公式中(de)待定系数,并证明确定后(de)求积公式具有3次代数精确度()()()()1010hhf x dx A f h A f A f h --=-++⎰证明：求积公式中含有三个待定系数,即101,,A A A -,将()21,,f x x x =分别代入求一、填空（共20分,每题2分）1. 设2.3149541...x *=,取5位有效数字,则所得(de)近似值x= .2.设一阶差商()()()21122114,321f x f x f x x x x --===---,()()()322332615,422f x f x f x x x x --===--则二阶差商 ()123,,______f x x x =3. 设(2,3,1)TX =--, 则2||||X = ,=∞||||X .4．求方程 21.250x x --= (de)近似根,用迭代公式 1.25x x =+,取初始值 01x =, 那么 1______x =。

数值分析课后习题答案

0 1
0 10 1 1 0 0 0 1
0 0 12 1 1 2 0 0 0

1 2
0 0 0 1 1 0
1 2

1 2

1 2
1
0 0 0 1 0

1 2

1 2

0
1 2

1 2
0
0
0
341 1 1
2-5.对矩阵A进行LDLT分解和GGT分解，并求解方程组
Ax=b，其中
16 4 8
1
A 4 5 4 , b 2
8 4 22
3
解
16 A 4
4 5
84
44 11
2-3(1).对矩阵A进行LU分解,并求解方程组Ax=b,其中
2 1 1 A1 3 2
4 ,b6
1 2 2
5
解
2 A 1
1 3
1 2

2 11
22
1
5 2
1

3 21来自，所以 A12
1
2 1 1

5 3
2-2(1).用列主元Gauss消元法解方程组
3 2 6x1 4 10 7 0x2 7 5 1 5x3 6
解
3 2 6 4 10 7 0 7 10 7 0 7

r1r2
消元

10 7 0 7 3 2 6 4 0 0.1 6 6.1
r=0.5101-n/3.162…<0.5101-n/3<0.01% 因此只需n=5.即取101/2=3.1623

数值分析习题含答案

2
x1 )
f (x0)
(x
x 0 )( x x0 x1
x1 )
f ' ( x0 )
(x ( x1
x0)
2 2
x0 )
f ( x1 )
R ( x)
其中 R(x) 由以下计算得到：构造辅助函数：
(t ) f (t ) N 2 (t ) (t (x x0 ) (t x0 ) ( x
2 2
x1 ) x1 )
f [ 2 ,2 ] =-2089 ，
0 1 2 7
0 1 7
f (x)
M ,
x
[ a , b ] ，证明：在任意相邻两节点间
R1 ( x )
1 8
Mh
2
。
x xi x xi M
1
f ( ) R1 i ( x ) 2 M 8 h 2,
h ,
2
x
8 ,n
[ xi , xi
1
]
R1 ( x )
max R1 i ( x )
1 2
s
2
[( x
xi
1
))( x
x
i
1 2
)( x
x i )]
e
4
h
3
[ s( s
1)( s
1)] 24
3 9
e h
4
3
10
6
3!
8
h
1 . 317
则用二次插值的步长应：
h
0 .6585
10
2
2-6 对区间 [a,b] 作步长为 h 的剖分，且做线性插值，其误差限为证明：区间上的误差限：误差限： 2-7 设 f ( x ) 解：自变量 1 2

数值分析题库答案(含详细解题步骤)

第 1 页/共 22 页1. 正方形的边长大约为100cm ，应怎样测量才干使面积误差不超过1cm 22. 已测得某场地长l 的值为110=*l m ，宽d 的值为80=*d m ，已知 2.0≤-*l l m,1.0≤-*d d m, 试求面积ld s =的绝对误差限与相对误差限.3.为使π的相对误差小于0.001%,至少应取几位有效数字？4.设x的相对误差界为δ，求n x的相对误差界.5.设有3个近似数a=2.31，b=1.93，c=2.24，它们都有3位有效数字，试计算p=a+bc的误差界和相对误差界，并问p的计算结果能有几位有效数字？第 3 页/共 22 页6. 已知333487.034.0sin ,314567.032.0sin ==，请用线性插值计算3367.0sin 的值，并预计截断误差.7. 已知sin0.32=0.314567, sin0.34=0.333487, sin0.36= 0.352274,用抛物插值计算sin0.3367的值, 并预计误差.8. 已知16243sin ,sin πππ===请用抛物插值求sin50的值,并预计误差9. . .6,8,7,4,1)(,5,4,3,2,1求四次牛顿插值多项式时设当==i i x f x第 5 页/共 22 页10. 已知4)2(,3)1(,0)1(=-=-=f f f ，求函数)(x f 过这3点的2次牛顿插值多项式.11. 设x x f =)(，并已知483240.1)2.2(,449138.1)1.2(,414214.1)0.2(===f f f ，试用二次牛顿插值多项式计算(2.15)f 的近似值，并研究其误差12. 设],[)(b a x f 在上有四阶延续导数，试求满意条件)2,1,0()()(==i x f x P i i 及)()(11x f x P '='的插值多项式及其余项表达式.13. 给定3201219(),,1,,44f x x x x x ====试求()f x 在1944⎡⎤⎢⎥⎣⎦，上的三次埃尔米特插值多项式()P x ,使它满意11()()(0,1,2),()(),i i P x f x i P x f x ''===并写出余项第 7 页/共 22 页表达式.14. 设],1,0[,23)(2∈++=x x x x f 试求)(x f 在]1,0[上关于,,1{,1)(x span x =Φ=ρ}2x 的最佳平方逼近多项式15.已知实验数据如下：用最小二乘法求形如y=a+bx2的拟合曲线，并计算均方误差.16.已知数据表如下第 9 页/共 22 页x i 1 2 3 4 5 y iωi4 4.56 8 8.5 2 1 3 1 1试用最小二乘法求多项式曲线与此数据组拟合17. .1)(},1{span ,1]41[)(的最佳平方逼近多项式中的关于上的在在求==Φ=x x x x f ρ18. 决定求积公式⎰++≈10110)1()(32)0()(f A x f f A dx x f 中的待定参数110,,A x A , 使其代数精度尽量高，并指出所决定的求积公式的代数精度.19. 用复化辛普森公式计算积分⎰=10dx e I x , 问区间[0,1]应分多少等分才干使截断误差不超过?10215-⨯第 11 页/共 22 页20. 利用下表中给出的数据，分离用复化梯形公式和复化辛甫生公式计算定积分dx x I ln 21⎰=的近似值（要求结果保留到小数点后六位）21. 用复化梯形公式和复化辛甫生公式计算积分⎰=6.28.1)(dx x f I ,函数)(x f 在某些节点上的值如下图：（本题共14分）22. 决定公式⎰+≈101100)()()(x f A x f A dx x f x 的系数1010,,,x x A A ,使其具有最高代数精度23. 决定求积公式⎰++≈1110)1()(32)0()(f A x f f A dx x f 中的待定参数110,,A x A ,使其代数精度尽量高，并指出所决定的求积公式的代数精度第 13 页/共 22 页24.用LU 分解法求解以下方程组 (10分)123123142521831520x x x ⎛⎫⎛⎫⎛⎫ ⎪⎪ ⎪= ⎪⎪ ⎪ ⎪⎪ ⎪⎝⎭⎝⎭⎝⎭25.用LU 分解法求解以下方程组⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛8892121514131615141321x x x26. 用LU 分解法求解以下方程组⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫⎝⎛542631531321321x x x27. 设方程组b Ax =，其中⎪⎪⎪⎭⎫⎝⎛-=220122101A ，Tb ⎪⎭⎫ ⎝⎛-=32,31,21, 已知它有解Tx ⎪⎭⎫⎝⎛-=0,31,21，若右端有小扰动61021-∞⨯=bδ，试预计由此引起的解的相对误差.第 15 页/共 22 页28. 设方程组b Ax =，其中212 1.0001A -⎛⎫= ⎪-⎝⎭，11.0001b -⎛⎫= ⎪⎝⎭,当右端向量b 有误差00.0001δ⎛⎫= ⎪⎝⎭b 时，试预计由此引起的解的相对误差(用∞范数计算)29. 给定b Ax =，其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=111a a a a a a A 证实：(1) 当121<<-a 时，A 对称正定，从而GS 法收敛. (2) 惟独当2121<<-a 时，J 法收敛.30. 对于线性方程组⎪⎩⎪⎨⎧-=+-=-+=+1242043 16343232121x x x x x x x ，列出求解此方程组的Jacobi 迭代格式，并判断是否收敛。

《数值分析》所有参考答案

解：
等价三角方程组
，，
11.设计算机具有4位字长。分别用Gauss消去法和列主元Gauss消去法解下列方程组，并比较所得的结果。
解：Gauss消去法
回代
列主元Gauss消去
15.用列主元三角分解法求解方程组。其中
A= ,
解：
等价三角方程组
回代得
，，，
16.已知，求，，。
解：
, ,
17.设。证明
，(II)
，
当时
当时
迭代格式(II)对任意均收敛
3) ，
构造迭代格式 (III)
，
当时
当时
迭代格式(III)对任意均收敛
4)
取格式(III)
，，，
4.用简单迭代格式求方程的所有实根，精确至有3位有效数。
解：
当时，，
1 2
当时
，
，
，，
1)
迭代格式，
,
当时，，
任取迭代格式收敛于
是中的一种向量范数。
解：
当时存在使得
，
,
所给为上的一个范数
18.设。证明
(1) ;
(2) ;
(3) 。
解：(1)
(2)
(3)
19.设
A=
求，，及，。
解： ,
Newton迭代格式
,
20.设为上任意两种矩阵(算子)范数，证明存在常数
，使得
对一切均成立。
解：由向量范数的等价性知道存在正常数使得
,
=0.187622
[23.015625 , 23.015625+0.187622]

数值分析练习题加答案(一)

数值分析期末考试一、设80~=x ，若要确保其近似数的相对误差限为0.1%，则它的近似数x 至少取几位有效数字？（4分）解：设x 有n 位有效数字。

因为98180648=<<=，所以可得x 的第一位有效数字为8（1分）又因为21101011000110821--⨯=<⨯⨯≤n ε，令321=⇒-=-n n ，可知x 至少具有3位有效数字（3分）。

二、求矩阵A 的条件数1)(A Cond （4分）。

其中⎥⎦⎤⎢⎣⎡=4231A 解：⎥⎦⎤⎢⎣⎡--=-5.05.1121A （1分） 1A =7（1分） 2711=-A （1分）249)(1=A Cond （1分）三、用列主元Gauss 消元法法求解以下方程组（6分）942822032321321321=++-=++--=+-x x x x x x x x x解：→⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---→⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----→⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----5.245.2405.35.230914220321821191429142821120321 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---→⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---8175835005,245.24091425.33.2305.245.2409142(4分) 等价三角方程组为：⎪⎪⎩⎪⎪⎨⎧-=-=+-=++,8175835,5.245.24,942332321x x x x x x （1分）回代得1,3,5123==-=x x x （1分）四、设.0,2,3,1,103)(3210234=-===-+-=x x x x x x x x f 1）求以3210,,,x x x x 为节3次Lagrange 多项式；（6分） 2）求以3210,,,x x x x 为节3次Newton 多项式；（6分）3）给出以上插值多项式的插值余项的表达式（3分）解：由0,2,3,13210=-===x x x x 可得10)(,34)(,1)(,11)(3210-==-=-=x f x f x f x f即得： +------+------=))()(())()(()())()(())()(()()(312101320130201032103x x x x x x x x x x x x x f x x x x x x x x x x x x x f x L=------+------))()(())()(()())()(())()(()(23130321033212023102x x x x x x x x x x x x x f x x x x x x x x x x x x x f+-+--+-⨯-+-+--+-⨯-)03)(23)(13()0)(2)(1()1()01)(21)(31()0)(2)(3(11x x x x x x326610.)20)(30)(10()2)(3)(1()10()02)(32)(12()0)(3)(1(34x x x x x x x x x -+--=+--+--⨯-+---------⨯2）计算差商表如下：i x )(i x f 一阶差商二阶差商三阶差商1 -11 3 -1 5 -2 34 -7 4 0-10-225-1则=+-----+-+-=)2)(3)(1()3)(1(4)1(511)(3x x x x x x x N326610x x x -+--3）)2)(3)(1())()()((!4)()(3210)4(3+--=----=x x x x x x x x x x x x f x R ξ五、给定方程组b Ax =，其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=100131w w w w A 。

数值分析习题（含答案）

数值分析习题（含答案）第一章绪论姓名学号班级习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。

1 若误差限为5105.0-?,那么近似数0.003400有几位有效数字？（有效数字的计算）解：2*103400.0-?=x ，325*10211021---?=?≤-x x 故具有3位有效数字。

2 14159.3=π具有4位有效数字的近似值是多少？（有效数字的计算）解：10314159.0?= π，欲使其近似值*π具有4位有效数字，必需41*1021-?≤-ππ，3*310211021--?+≤≤?-πππ，即14209.314109.3*≤≤π即取（3.14109 , 3.14209）之间的任意数，都具有4位有效数字。

3 已知2031.1=a ，978.0=b 是经过四舍五入后得到的近似值，问b a +，b a ?有几位有效数字？（有效数字的计算）解：3*1021-?≤-aa ，2*1021-?≤-b b ，而1811.2=+b a ，1766.1=?b a 2123****102110211021)()(---?≤?+?≤-+-≤+-+b b a a b a b a故b a +至少具有2位有效数字。

2123*****10210065.01022031.1102978.0)()(---?≤=?+?≤-+-≤-b b a a a b ba ab 故b a ?至少具有2位有效数字。

4 设0>x ，x 的相对误差为δ，求x ln 的误差和相对误差？（误差的计算）解：已知δ=-**xx x ，则误差为δ=-=-***ln ln xx x x x则相对误差为******ln ln 1ln ln ln xxx x xxx x δ=-=-5测得某圆柱体高度h 的值为cm h 20*=，底面半径r 的值为cm r 5*=，已知cm h h 2.0||*≤-，cm r r 1.0||*≤-，求圆柱体体积h r v2π=的绝对误差限与相对误差限。

数值分析练习题附答案

1
2-3 对矩阵 A 进行 LDLT 分解和 GGT 分解，求解方程组 Ax=b,其中
16 4 8
1
A=( 4 5 −4) ， b＝(2)
8 −4 22
3
解：（注：课本 P26 P27 根平方法）
设 L=(l i j )，D=diag(di)，对 k=1,2,…,n,
其中��=��－∑��=−11 ��2��
��31=(��31 − ∑0��=1 ��3��1�� )/ ��1=186=12
��32=(��32
−
∑1��=1
��3��2��
�� )/
6.6667
，得 ��3 ＝ 1.78570
−1 209
��4
0
��4
0.47847
(
56
−1
780 (��5) 209)
(200)
(��5) ( 53.718 )
1 −1
4
1 −4
15
��1
25
��2
6.6667再由1源自− 15561
− 56
209

x (k1) 1

1 5
(12

2 x2( k )

x (k) 3
)

2 5
x (k) 2

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

第一章绪论姓名学号班级习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。

1 若误差限为5105.0-⨯,那么近似数0.003400有几位有效数字？（有效数字的计算）解：2*103400.0-⨯=x ，325*10211021---⨯=⨯≤-x x 故具有3位有效数字。

2 14159.3=π具有4位有效数字的近似值是多少？（有效数字的计算）解：10314159.0⨯= π，欲使其近似值*π具有4位有效数字，必需41*1021-⨯≤-ππ，3*310211021--⨯+≤≤⨯-πππ，即14209.314109.3*≤≤π即取（3.14109 , 3.14209）之间的任意数，都具有4位有效数字。

3 已知2031.1=a ，978.0=b 是经过四舍五入后得到的近似值，问b a +，b a ⨯有几位有效数字？（有效数字的计算）解：3*1021-⨯≤-aa ，2*1021-⨯≤-b b ，而1811.2=+b a ，1766.1=⨯b a 2123****102110211021)()(---⨯≤⨯+⨯≤-+-≤+-+b b a a b a b a故b a +至少具有2位有效数字。

2123*****10210065.01022031.1102978.0)()(---⨯≤=⨯+⨯≤-+-≤-b b a a a b b a ab 故b a ⨯至少具有2位有效数字。

4 设0>x ，x 的相对误差为δ，求x ln 的误差和相对误差？（误差的计算）解：已知δ=-**xx x ，则误差为 δ=-=-***ln ln xx x x x则相对误差为******ln ln 1ln ln ln xxx x xxx x δ=-=-5测得某圆柱体高度h 的值为cm h 20*=，底面半径r 的值为cm r 5*=，已知cm h h 2.0||*≤-，cm r r 1.0||*≤-，求圆柱体体积h r v2π=的绝对误差限与相对误差限。

（误差限的计算）解：*2******2),(),(h h r r r h r r h v r h v -+-≤-ππ绝对误差限为πππ252.051.02052)5,20(),(2=⨯⋅+⨯⋅⋅⋅≤-v r h v相对误差限为%420120525)5,20()5,20(),(2==⋅⋅≤-ππv v r h v 6 设x 的相对误差为%a ,求nx y =的相对误差。

（函数误差的计算）解：%**a x x x =-，)%(******na xx x nxx x yy y nnn =-≤-=-7计算球的体积，为了使体积的相对误差限为%1，问度量半径r 时允许的相对误差限为多大？（函数误差的计算）解：球体积为 334)(r r v ⋅⋅=π，3**34)(r r v ⋅⋅=π 欲使%13344)()()(**3**2***=-=⋅⋅-⋅⋅=-r r r r r r r r v r v r v ππ，必须%31**=-r r r 。

8 设⎰-=11dx e x eI x n n ，求证：（1）)2,1,0(11 =-=-n nI I n n（2）利用（1）中的公式正向递推计算时误差逐步增大；反向递推计算时误差逐步减小。

（计算方法的比较选择）解：11111011011111][-------=-=-==⎰⎰⎰n x n xn x nx n n nI dx e x nedx e xn ex e de x eI 111101)1(----=-==⎰e e e dx e e I x如果初始误差为*000I I -=ε，若是向前递推，有0221*11*!)1()1()1()1()1(εεεεn n n n nI nI I I nn n n n n n n -==--=-=---=-=----可见，初始误差0ε的绝对值被逐步地扩大了。

如果是向后递推n n I nn I 111-=-，其误差为 n n n I I εεεε!)1(211)1(11)1111()1111(221*110-==⋅-=-=---= 可见，初始误差n ε的绝对值被逐步减少了。

第二章插值法姓名学号班级习题主要考察点：拉格朗日插值法的构造，均差的计算，牛顿插值和埃尔米特插值构造，插值余项的计算和应用。

1 已知1)2(,1)1(,2)1(===-f f f ，求)(x f 的拉氏插值多项式。

（拉格朗日插值）解法一（待定系数法）：设c bx ax x L ++=2)(，由插值条件，有⎪⎩⎪⎨⎧=++=++=+-12412c b a c b a c b a 解得：3/4,2/1,6/1=-==c b a 。

故 342161)(2+-=x x x L 。

解法二（基函数法）：由插值条件，有1)12)(12()1)(1(1)21)((11()2)(1(2)21)(11()2)(1()(⋅-+-++⋅-+-++⋅------=x x x x x x x L)1)(1(31)2)(1(21)2)(1(31-++-+---=x x x x x x 3421612+-=x x 2 已知9,4,10===x x x y ，用线性插值求7的近似值。

（拉格朗日线性插值）解：由插值节点与被插函数，可知，240==y ，391==y ，其线性插值函数为565134942949)(+=⋅--+⋅--=x x x x L 7的近似值为6.25135657)7(≈=+=L 。

3 若),...1,0(n j x j =为互异节点，且有)())(())(()())(())(()(11101110n j j j j j j j n j j j x x x x x x x x x x x x x x x x x x x x x l ----------=+-+-试证明),...1,0()(0n k x x l xnj k jk j =≡∑=。

（拉格朗日插值基函数的性质）解：考虑辅助函数∑=-=nj k jk j x x l xx F 0)()(，其中，n k ≤≤0，),(∞-∞∈x 。

)(x F 是次数不超过n 的多项式，在节点i x x =（n i ≤≤0）处，有)()()(0=-=-=-=∑=ki k i k i i i k i nj k i i j k j i x x x x l x x x l x x F 这表明，)(x F 有n+1个互异实根。

故0)(≡x F ，从而∑=≡nj k jkj x x l x)(对于任意的n k ≤≤0均成立。

4 已知352274.036.0sin ,333487.034.0sin ,314567.032.0sin ===，用抛物线插值计算3367.0sin 的值并估计截断误差。

（拉格朗日二次插值）解：由插值条件，其抛物线插值函数为314567.0)36.032.0)(34.032.0()36.0)(34.0()(⋅----=x x x L333487.0)36.034.0)(32.034.0()36.0)(32.0(⋅----+x x352274.0)34.036.0)(32.036.0()34.0)(32.0(⋅----+x x将3367.0=x 代入，计算可得：3304.0)3367.0(≈L 。

其余项为：)36.0)(34.0)(32.0(!3sin )(----=x x x x r ξ其中，36.032.0<<ξ )36.0)(34.0)(32.0(61)(---≤x x x x r 故误差的上界为：71014.2)36.03367.0)(34.03367.0)(32.03367.0(61)3367.0(-⨯≤---≤r 。

5 用余弦函数x cos 在00=x ，41π=x ，22π=x 三个节点处的值，写出二次拉格朗日插值多项式, 并近似计算6cos π及其绝对误差与相对误差，且与误差余项估计值比较。

（拉格朗日二次插值）解：由插值条件，二次拉格朗日插值多项式为0)4/2/)(02/()4/)(0(21)2/4/)(04/()2/)(0(1)2/0)(4/0()2/)(4/()(⋅----+⋅----+⋅----=ππππππππππππx x x x x x x L22)2/(28)2/)(4/(8πππππ----=x x x x8508.09242)2/6/(6/28)2/6/)(4/6/(8)6(22≈+=----=ππππππππππL 绝对误差为：0153.01828439924223)6(6cos≈--=+-=-ππL 相对误差为：0179.028428439)6()6(6cos≈+--=-ππL L余项为：)2/)(4/(!3sin )(ππξ--=x x x x r ，其中，2/0πξ<< 其余项的上界为：)2/)(4/(61)(ππ--≤x x x x r 0239.06)26)(46(661)6(43≈=--≤πππππππr 比较可知，实际计算所得的绝对误差较余项公式所估计出的值要小一些。

6 已知函数值212)6(,82)4(,46)3(,10)1(,6)0(=====f f f f f ，求函数的四阶均差]6,4,3,1,0[f 和二阶均差]3,1,4[f 。

（均差的计算）从表中可查得：]6,4,3,1,0[=f 。

故6]3,1,4[=f 。

其实，根据均差的对称性，6]4,3,1[]3,1,4[==f f ，该值在第一个表中就可以查到。

7 设)())(()(10n x x x x x x x f ---= 求][1,0p x x x f 之值，其中1+≤n p ，而节点)1,1,0(+=n i x i 互异。

（均差的计算）解：由均差可以表示成为函数值的线性组合，有∑=-+-------=pi p i p i i i i i i i i p x x x x x x x x x x x x x f x x x f 0111101,0))(())(())(()(][而0)(=i x f p i ≤≤0，故0][1,0=p x x x f 。

8 如下函数值表建立不超过三次的牛顿插值多项式。

（牛顿插值多项式的构造）解：先构造均差表故 )2)(1(4)1(381)(----++=x x x x x x x N 。

9求一个次数小于等于三次多项式)(x p ，满足如下插值条件：2)1(=p ，4)2(=p ，3)2(='p ，12)3(=p 。

（插值多项式的构造）解法一（待定系数法）：设d cx bx ax x p +++=23)(，则c bx ax x p ++='23)(2，由插值条件，有⎪⎪⎩⎪⎪⎨⎧=+++=++=+++=+++123927341242482d c b a c b a d c b a d c b a 解得：6,15,9,2-==-==d c b a 。