数值分析期末试题

数值分析期末试题
数值分析期末试题

第一套

一、(8分)用列主元素消去法解下列方程组:

???

??=++-=+--=+-11

2123454

321321321x x x x x x x x x

二、(10分)依据下列数据构造插值多项式:y(0)=1,y(1)= —2,y '(0)=1, y '(1)=—4

三、(12分)分别用梯形公式和辛普生公式构造 复化的梯形公式、复化的辛普生公式并利用复化的梯形公式、复化的辛普生公式计算下列积分:

?9

1

dx

x n=4

四、(10分)证明对任意参数t ,下列龙格-库塔方法是二阶的。

五、(14分)用牛顿法构造求c 公式,并利用牛顿法求115。保留有效数字五位。

六、(10分)方程组AX=B 其中A=?????

????

?10101a a a a 试就AX=B 建立雅可比迭代法和高斯-赛德尔迭代法,并讨论a 取何值时

迭代收斂。

七、(10分)试确定常数A,B,C,a,使得数值积分公式?-++-≈2

2

)(}0{)()(a Cf Bf a Af dx x f 有尽可能多的

代数精确度。并求该公式的代数精确度。 八、{6分}

证明:

A ≤ 其中A 为矩阵,V 为向量.

第二套

一、(8分)用列主元素消去法解下列方程组:

???

??=++=+-=+3

2221

43321

32132x x x x x x x x

二、(12分)依据下列数据构造插值多项式:y(0)=y '(0)=0,

y(1)=y '(1)= 1,y(2)=1

三、(14分)分别用梯形公式和辛普生公式构造 复化的梯形公式、复化的辛普生公式,并利用

?

??????

-+-+=++==++=+1

3121231)1(,)1((),(),()(2

hk t y h t x f k thk y th x f k y x f k k k h y y n n n n n n n n

复化的梯形公式、

复化的辛普生公式及其下表计算下列积分:

?2

/s i n πx d x

四、(12分)证明下列龙格-库塔方法是三阶的。

五、(10分)试确定常数A,B,C 使得数值积分公式

?++≈2

)

2()1()0()(Cf Bf Af dx x f

共 2 页 第 2 页

有尽可能多的代数精确度。并求该公式的代数精确度。

六、(14分)用牛顿法构造求c 1

公式,验证其收敛性。并求1/ e(保留4位有效数字)。 七、{10分}证明:设非负函数N(x )=x 为R n 上任意向量范数,则N(x )是x 分量x 1,x 2,…x n 的

连续函数.

参考答案

一、解:(8分)

???

??=++=+-=+3

2221

43321

32132x x x x x x x x

增广矩阵:

??????????→??????????-→??????????-12003/13/4102/312/112/102/3014302/312/11321221111431 (4分) 解得:x 1=2/3, x 2=-1/3 x 3=1./2 (8分) 二、解:(12分)

注:直接待定系数简单,或者用牛顿茶商

设 P(x)=φ0(x)y(0)+φ1(x) y(1)+φ2(x)y(2)+ψ0(x) y’(0)+ψ1(x) y’(1) (4分)

解得:

1(x)=x 2(x-2)2 φ2(x)=(1/12)x 2(x-1)2 ψ1(x)=-x 2(x-1)(x-2) (4分) P(x)= φ1(x) y(1)+φ2(x)y(2)+ψ1(x) y’(1)= φ1(x) +φ2(x)+ψ1(x)

??????

?++=++==++=+)

3/2,3/2()3/,3/(),()3(423121131hk y h x f k hk y h x f k y x f k k k h y y n n n n n n n n

= x 2(x-2)2+(1/12)x 2(x-1)2 +x 2(x-1)(x-2) (4分)

三、解:(14分) 推证复化的梯形公式 (3分)

推证复化的辛普生公式 (3分)

利用复化的梯形公式

?2

/0

sin πxdx

=0.96593 利用复化的辛普生公式

?2

/0

sin πxdx

=1.000003

四、(12分)证明:

k 3=f(x n ,y n )+2h/3f’(x n ,y n )+(2h/3)2f’’(x n ,y n )/2+0(h 2) (4分) y n+1=y n +h/4(3 k 3+k 1)= y n + h f(x n ,y n )+h 2f’(x n ,y n )/2+h 3/6f’’(x n ,y n ) +0(h 3) (8分) y n+1*= y n + h y n ’ +h 2y n ’’/2+h 3/6 y n ’’’ +0(h 3) y n+1 -y n+1*=0(h 3)

则该公式是三阶的 (12分)

五、解:(10分) 将1,x,x 2代入原式得A+B+C=2 B+2C=2 B+4C=8/3

解得:A=1/3, B=4/3 C =1/3

?

++≈

2

)2(31

}1{34)0(31)(f f f dx x f (8分)

代数精确度为2 (10分)。

六、证明:(14分)1/x-c=0

X k+1=x k -)()

(k k x f x f '=x k (2-cx k ) X k+1-1/c=-c(x k -1/c)2

设r k =1-cx k r k+1=r k 2 反复递推 r k =02

r k

(8分)

若选初值0

七、{10分}证明:设x =

∑=n

i i

i e

x 1

y =

∑=n

i i

i e

y 1

(4分)

)

(0)()()(1

1

∑∑=∞

==→-≤--≤-=-n

i i i n i i i e c y

x c e y x y x y x y N x N

?

?????

?

++=++==++=+)

3/2,3/2()3/,3/(),()3(423121131hk y h x f k hk y h x f k y x f k k k h y y n n n n n n n n

..(10分)

第三套

一、 (10分)利用列主元素消去法解方程:

??????????=?????????????????????---453311294642321x x x

二、 (15分)证明下面龙格-库塔方法是三阶的:

)

43

,43()

2,2()

,()

432(9

3213211h y h x f k h

y h x f k y x f k k k k h

y y n n n n n n n n ++=++==+++=+

三、 (10分)求3次插值多项式使:P(0)=3, P(1)=5,4)0(='P ,6)1(='P ,

四、 (20分)确定下面公式中的a,b ,使其代数精确度尽量高,并指出其代数精确度的次

数:

)]()([)()]()([2)(2b f a f a b a b f a f a

b dx x f b

a

'-'-++-≈

?

五、(20分)分别利用梯形公式和Simpson 公式推导复化的梯形公式和Simpson 公式,并分别

利用复化的梯形公式和Simpson 公式计算积分

?9

1

dx

x (n=8)

六、(15分)用二分法求方程f(x)=x 3+4x 2-10在区间[1,1.5]上的根。(1)要得到具有3位有效数

的近似根,须作几次二分;(2)用二分法求具有3位有效数的近似根。 七、(10分)设?是n n R ?中的任意范数,n n R A ?∈,则有A A ≤)(ρ

参考答案

五、(10分)利用列主元素消去法解方程:

解:

???????????????

?

-→?????

?????---4112

5

)45(021521

0529445331129)4(642 (5分)

x 1=139/20, x 2=5/2, x 3=-3/20 (10分)

六、(15分)证明下面龙格-库塔方法是三阶的:

证:)(61

)(21)()()(321ξy h x y h x y h x y x y n n n n '''+''+

'+=+(5分)

))

(,(21

)21(),(21),(22ξξy f h y x f h y x f k n n n n ''+'+=(9分) ))

(,(21

)43(),(43),(23ξξy f h y x f h y x f k n n n n ''+'+=(13分) ∴y(x n+1)- y n+1=o(h 3

) (15分)

七、(10分)求3次插值多项式使:P(0)=3, P(1)=5,4)0(='P ,6)1(='P ,

解:设)()()()()(2211

21103x p x p x p x p x p ??φφ'+'++= (2分) 0)1(,0)0,0)1(,1)0(1(111='='==φφφφ

0)1(,0)0(,1)1(,0)0(22

22='='==φφφφ 0)1(,1)0(,0)1(,0)0(11

11='='==???? 1)1(,0)0(,0)1(,0)0(22

22='='==???? (6分) =∴)(3x p 3+4x-2x 2+6x 2(x-1) (10分)

八、(20分)确定下面公式中的a,b ,使其代数精确度尽量高,并指出其代数精确度的次数:

)]()([)()]()([2)(2b f a f a b a b f a f a

b dx x f b

a

'-'-++-≈

?

解:将1,x ,x 2,,x 3

代入)]()([)()]()([2)(2b f a f a b a b f a f a

b dx x f b

a

'-'-++-≈

?

(4分)

得]

22[)(][2)(31222

33b a a b a b a a b a b --++-=-(10分)

]33[)(][2)(412223344b a a b a b a a b a b --++-=-

a=b=1/2(15分)

将1,x ,x 2,,x 3,x 4,x 5代入公式的两端,可得该公式具有4次代数精确度。(20分)

五、(20分)分别利用梯形公式和Simpson 公式推导复化的梯形公式和Simpson 公式,并分别利用复化的梯

形公式和Simpson 公式计算积分

?9

1

dx

x (n=8)

证: 利用梯形公式推导复化的梯形公式(5分)

Simpson 公式推导复化Simpson 公式(10分)

解:利用复化的梯形公式

?9

1

dx

x (n=8) =17.22774 (15分)

Simpson 公式计算积分

?9

1

dx

x (n=8)=17.32222 (20分)

六、(15分)用二分法求方程f(x)=x 3+4x 2-10在区间[1,1.5]上的根。(1)要得到具有3位有效数的近似根,

须作几次;(2)用二分法求具有3位有效数的近似根。 解:须作3次(5分)

将[1,1.5] [1,1.25], [1.25,1.5] f(1)<0, f(1.25) <0, (8分) 将[1.25,1.5] 二分为[1.25,1.375],[1.375,1.5] f(1.375) >0, (10分) 将[1.25,1.375]二分为[1.25,1.3125],[1.3125,1.375] f(1.3125) <0(12分)

∴[1.3125,1.375]的 中点为方程f(x)=x 3+4x 2-10的近似根(15分)

七、设

?是n n R ?中的任意范数,n n R A ?∈,则有A A ≤)(ρ

证: 设λ是的任意特征值,x 为相应的向量, (2分)

则x Ax λ=, x A Ax x x ≤==λλ (8分)∴A A ≤)(ρ(10分)

数值分析学期期末考试试题与答案(A)

期末考试试卷(A 卷) 2007学年第二学期 考试科目: 数值分析 考试时间:120 分钟 学号 姓名 年级专业 一、判断题(每小题2分,共10分) 1. 用计算机求 1000 1000 1 1 n n =∑时,应按照n 从小到大的顺序相加。 ( ) 2. 为了减少误差,进行计算。 ( ) 3. 用数值微分公式中求导数值时,步长越小计算就越精确。 ( ) 4. 采用龙格-库塔法求解常微分方程的初值问题时,公式阶数越高,数值解越精确。( ) 5. 用迭代法解线性方程组时,迭代能否收敛与初始向量的选择、系数矩阵及其演变方式有 关,与常数项无关。 ( ) 二、填空题(每空2分,共36分) 1. 已知数a 的有效数为0.01,则它的绝对误差限为________,相对误差限为_________. 2. 设1010021,5,1301A x -????????=-=-????????-???? 则1A =_____,2x =______,Ax ∞ =_____. 3. 已知5 3 ()245,f x x x x =+-则[1,1,0]f -= ,[3,2,1,1,2,3]f ---= . 4. 为使求积公式 1 1231 ()()(0)33 f x dx A f A f A f -≈- ++? 的代数精度尽量高,应使1A = ,2A = ,3A = ,此时公式具有 次的代数精度。 5. n 阶方阵A 的谱半径()A ρ与它的任意一种范数A 的关系是 . 6. 用迭代法解线性方程组AX B =时,使迭代公式(1) ()(0,1,2,)k k X MX N k +=+=产 生的向量序列{ }() k X 收敛的充分必要条件是 . 7. 使用消元法解线性方程组AX B =时,系数矩阵A 可以分解为下三角矩阵L 和上三角矩

数值分析典型习题

特别声明:考试时需带计 算器作辅助计算 1.2015x *=是经四舍五入得到的近似值,则其相对误差* r e ≤-31 104 ?. 2. 01(),(), ,()n l x l x l x 是以01,, ,n x x x 为节点的拉格朗日插值基函数,则 3.设(0)1(1)3(2)4(3)2f =,f =,f =,f =,[0123]f =,,,1 3 - . 4. 利用Simpson 公式求?2 1 2dx x = 7.3 5. 设求积公式1 0()d (),(1)n k k k f x x A f x n ≈≥∑?=是Gauss 型求积公式,则3 n k k k A x == ∑1 .4 6. 数值微分公式(2)(2) ()i i i f x h f x h f x h +≈ --'的截断误差为 2().O h 7. 设1101A ?? = ??? ,则A 的谱半径()A ρ= 1 ,A 的条件数1cond ()A = 4. 8. 用牛顿下山法求解方程3 03 x x -=根的迭代公式是 2 13 3(1),3n n n n x x x x x λ+-=-- 下山条件是 1()().n n f x f x +< 9.对任意初始向量(0)x 及任意向量f ,线性方程组的迭代公式(1)()(0,1,2,)k k k +=+=x Bx f ,迭代序列()k x 收敛于方程组的精确解x *的充分必要条件是()1.ρ

11:数值分析试题2009~2010

中国石油大学(北京)2009--2010学年第一学期 研究生期末考试试题A (闭卷考试) 课程名称:数值分析 注:计算题取小数点后四位 一、填空题(共30分,每空3分) 1、 已知x =0.004532是由准确数a 经四舍五入得到的近似值,则x 的绝对误差 界为_______________。 2、数值微分公式()() '()i i i f x h f x f x h +-≈ 的截断误差为 。 3、已知向量T x =,求Householder 变换阵H ,使(2,0)T Hx =-。 H = 。 4、利用三点高斯求积公式 1 1 ()0.5556(0.7746)0.8889(0)0.5556(0.7746) f x d x f f f -≈-++? 导出求积分 4 0()f x dx ?的三点高斯求积公式 。 5、4 2 ()523,[0.1,0.2,0.3,0.4,0.5]_____.f x x x f =+-=若则 6、以n + 1个互异节点x k ( k =0,1,…,n ),(n >1)为插值节点的 Lagrange 插值基函数为l k (x)( k =0,1,…,n ),则 (0)(1)__________.n k k k l x =+=∑ 7、已知3()P x 是用极小化插值法得到的cos x 在[0,4]上的三次插值多项式,则3()P x 的 截断误差上界为3()cos ()R x x P x =-≤_________. 8、已知向量(3,2,5)T x =-,求Gauss 变换阵L ,使(3,0,0)T Lx =。L =_________. 9、设3 2 ()(7)f x x =-, 给出求方程()0f x =根的二阶收敛的迭代格式_________。

数值分析典型例题

第一章典型例题 例3 ln2=0.…,精确到10-3的近似值是多少 解 精确到10-3=,即绝对误差限是=, 故至少要保留小数点后三位才可以。ln2 第二章典型例题 例1 用顺序消去法解线性方程组 ??? ??1 -=4+2+4=+2+31 -=4++2321 321321x x x x x x x x x 解 顺序消元 ?? ?? ??????---???→???????????---????→???????????--=-?+-?+-?+1717005.555.00141 25.025.105.555.001412142141231412]b A [)3()2/1()2/3(231312r r r r r r M 于是有同解方程组 ?? ? ??-==--=++17175.555.0142332321x x x x x x 回代得解 x 3=-1, x 2=1,x 1=1,原线性方程组的解为X =(1,1,-1)T 例2 取初始向量X (0)=(0,0,0)T ,用雅可比迭代法求解线性方程组 ??? ??5 =+2+23=++1=2-2+321 321321x x x x x x x x x 解 建立迭代格式 ???????+--=+--=++-=+++5223122) (2)(1)1(3 ) (3)(1)1(2 ) (3)(2)1(1k k k k k k k k k x x x x x x x x x (k =1,2,3,…)

第1次迭代,k =0 X (0)=0,得到X (1)=(1,3,5)T 第2次迭代,k =1 ???????-=+?-?-=-=+--==+?+?-=3 532123 351515232)2(3) 2(2)2(1x x x X (2)=(5,-3,-3)T 第3次迭代,k =2 ???????=+-?-?-==+---==+-?+-?-=1 5)3(2521 3)3(511)3(2)3(2)2(3) 3(2)3(1x x x X (3)=(1,1,1)T 第4次迭代,k =3 ???????=+?-?-==+--==+?+?-=1 512121 311111212)2(3) 2(2)2(1x x x X (4)=(1,1,1)T 例4 证明例2的线性方程组,雅可比迭代法收敛,而高斯-赛德尔迭代法发散。 证明 例2中线性方程组的系数矩阵为 A =?? ?? ? ?????-122111221 于是 D =?? ?? ??????100010001 D -1=D ??????????=022001000L ~ ????? ?????-=000100220U ~ 雅可比迭代矩阵为

数值分析期末考试复习题及其答案.doc

数值分析期末考试复习题及其答案 1. 已知325413.0,325413* 2* 1==X X 都有6位有效数字,求绝对误差限。(4分) 解: 由已知可知,n=6 5.01021 ,0,6,10325413.0016*1=?= =-=?=ε绝对误差限n k k X 2分 620* 21021,6,0,10325413.0-?=-=-=?=ε绝对误差限n k k X 2分 2. 已知?????=001A 220 - ???? ?440求21,,A A A ∞ (6分) 解: {},88,4,1max 1==A 1分 {},66,6,1max ==∞A 1分 () A A A T max 2λ= 1分 ?????=001A A T 420 ?? ?? ? -420?????001 220 - ?????440=?????001 080 ???? ?3200 2分 {}3232,8,1max )(max ==A A T λ 1分 24322==A 3. 设3 2 )()(a x x f -= (6分) ① 写出f(x)=0解的Newton 迭代格式 ② 当a 为何值时,)(1k k x x ?=+ (k=0,1……)产生的序列{}k x 收敛于2 解: ①Newton 迭代格式为: x a x x x a x a x x a x x x f x f x x k k k k k k k k k k 665)(665)(6)()(')(2 2 32 1 += +=---=-=+? 3分

②时迭代收敛即当222,112 10)2(',665)('2<<-<-=-=a a x a x ?? 3分 4. 给定线性方程组Ax=b ,其中:? ??=1 3A ??? 22,??????-=13b 用迭代公式)()()()1(k k k Ax b x x -+=+α(k=0,1……)求解Ax=b ,问取什么实数α,可使迭代收 敛 (8分) 解: 所给迭代公式的迭代矩阵为?? ? --? ??--=-=ααααα21231A I B 2分 其特征方程为 0) 21(2)31(=----= -αλα ααλλB I 2分 即,解得αλαλ41,121-=-= 2分 要使其满足题意,须使1)(

数值分析典型习题资料

数值分析典型习题

特别声明:考试时需带计 算器作辅助计算 1.2015x *=是经四舍五入得到的近似值,则其相对误差* r e ≤ -31 104 ?. 2. 01(),(),,()n l x l x l x L 是以01,,,n x x x L 为节点的拉格朗日插值基函数,则 3.设(0)1(1)3(2)4(3)2f =,f =,f =,f =,[0123]f =,,,1 3 - . 4. 利用Simpson 公式求?2 1 2dx x = 7.3 5. 设求积公式1 0()d (),(1)n k k k f x x A f x n ≈≥∑?=是Gauss 型求积公式,则3 n k k k A x == ∑1 .4 6. 数值微分公式(2)(2) ()i i i f x h f x h f x h +≈ --'的截断误差为 2().O h 7. 设1101A ?? = ??? ,则A 的谱半径()A ρ= 1 ,A 的条件数1cond ()A = 4. 8. 用牛顿下山法求解方程3 03 x x -=根的迭代公式是 2 13 3(1),3n n n n x x x x x λ+-=-- 下山条件是 1()().n n f x f x +< 9.对任意初始向量(0)x 及任意向量f ,线性方程组的迭代公式(1)()(0,1,2,)k k k +=+=L x Bx f ,迭代序列()k x 收敛于方程组的精确解x *的充分必要条件是()1.ρ

数值分析典型例题

第一章典型例题 例3…,精确到10-3的近似值是多少? 解 精确到10-3=,即绝对误差限是?=, 故至少要保留小数点后三位才 可以。ln2? 第二章典型例题 例1 用顺序消去法解线性方程组 解 顺序消元 于是有同解方程组 回代得解 x 3=-1, x 2=1,x 1=1,原线性方程组的解为X =(1,1,-1)T 例2 取初始向量X (0)=(0,0,0)T ,用雅可比迭代法求解线性方程组 解 建立迭代格式 ??? ????+--=+--=++-=+++5223122)(2)(1)1(3) (3)(1)1(2 )(3)(2)1(1k k k k k k k k k x x x x x x x x x (k =1,2,3,…) 第1次迭代,k =0 X (0)=0,得到X (1)=(1,3,5)T 第2次迭代,k =1 X (2)=(5,-3,-3)T 第3次迭代,k =2 X (3)=(1,1,1)T 第4次迭代,k =3

X (4)=(1,1,1)T 例4 证明例2的线性方程组,雅可比迭代法收敛,而高斯-赛德尔迭 代法发散。 证明 例2中线性方程组的系数矩阵为 A =?? ?? ? ?????-122111221 于是 D =?? ?? ??????100010001 D -1 =D ?? ?? ? ?????=022001000L ~ ?? ?? ? ?????-=000100220U ~ 雅可比迭代矩阵为 B 0=?? ?? ? ?????--=??????????-??????????-=+--022101220022101220100010001)U ~L ~(D 1 得到矩阵B 0的特征根03,2,1=λ,根据迭代基本定理4,雅可比迭代法收敛。 高斯-赛德尔迭代矩阵为 G =-U ~ )L ~D (1-+ =-?? ?? ??????----=??????????-??????????---=??????????-??????????-2003202200001002201200110010001002201220110011 解得特征根为?1=0,?2,3=2。由迭代基本定理4知,高斯-赛德尔迭代发散。 例5 填空选择题: 1. 用高斯列主元消去法解线性方程组 作第1次消元后的第2,3个方程分别为 。

数值分析期末复习资料

数值分析期末复习资料

数值分析期末复习 题型:一、填空 二、判断 三、解答(计算) 四、证明 第一章 误差与有效数字 一、 有效数字 1、 定义:若近似值x*的误差限是某一位的半个单位,该位到x*的第一位非零数字共有n 位,就说 x*有n 位有效数字。 2、 两点理解: (1) 四舍五入的一定是有效数字 (2) 绝对误差不会超过末位数字的半个单位eg. 3、 定理1(P6):若x*具有n 位有效数字,则其相对误差限为 4、 考点: (1)计算有效数字位数:一个根据定义理解,一个根据定理1(P7例题3) 二、 避免误差危害原则 1、 原则: (1) 避免大数吃小数(方法:从小到大相加;利用韦达定理:x1*x2= c / a ) (2) 避免相近数相减(方法:有理化)eg. 或 (3) 减少运算次数(方法:秦九韶算法)eg.P20习题14 三、 数值运算的误差估计 1、 公式: (1) 一元函数:|ε*( f (x *))| ≈ | f ’(x *)|·|ε*(x )|或其变形公式求相对误差(两边同时 除以f (x *)) eg.P19习题1、2、5 (2) 多元函数(P8)eg. P8例4,P19习题4 *(1) 11 102n r a ε--≤?;x εx εx εx ++=-+();1ln ln ln ??? ? ??+=-+x εx εx x cos 1-2sin 22x =

第二章 插值法 一、 插值条件 1、 定义:在区间[a,b]上,给定n+1个点,a ≤x 0<x 1<…<x n ≤b 的函数值 yi=f(xi),求次数不超过n 的多项式P(x),使 2、 定理:满足插值条件、n+1个点、点互异、多项式次数≤n 的P(x)存在且唯一 二、 拉格朗日插值及其余项 1、 n 次插值基函数表达式(P26(2.8)) 2、 插值多项式表达式(P26(2.9)) 3、 插值余项(P26(2.12)):用于误差估计 4、 插值基函数性质(P27(2.17及2.18))eg.P28例1 三、 差商(均差)及牛顿插值多项式 1、 差商性质(P30): (1) 可表示为函数值的线性组合 (2) 差商的对称性:差商与节点的排列次序无关 (3) 均差与导数的关系(P31(3.5)) 2、 均差表计算及牛顿插值多项式 四、埃尔米特插值(书P36) 两种解法: (1) 用定义做:设P 3(x)=ax 3+bx 2+cx+d ,将已知条件代入求解(4个条件:节点函数值、导数值相 等各2个) (2) 牛顿法(借助差商):重节点eg.P49习题14 五、三次样条插值定义 n i y x P i i n ,,2,1,0)( ==

数值分析典型例题

数值分析典型例题 例1 对下列各数写出具有5位有效数字的近似值。236.478, 0.00234711, 9.000024, 9.0000343 10?. 解:按照定义,以上各数具有5位有效数字的近似值分别为:236.478, 0.0023471, 9.0000, 9.0000310?。 注意: *x =9.000024的5位有效数字是9.0000而不是9,因为9 是1位有效数字。 例2 指出下列各数具有几位有效数字。2.0004, -0.00200, -9000, 9310?, 23 10-?。 解:按照定义,以上各数的有效数字位数分别为5, 3, 4,1,1 例3 已测得某物体行程* s 的近似值s=800m ,所需时间* s 的近似值为t=35s ,若已知m s s s t t 5.0||,05.0||**≤-≤-,试求平均速度v 的绝对误差和相对误差限。 解:因为t s v /=,所以)()(1)()()(2t e t s s e t t e t v s e s v v e -=??+??≈ 从 而 05.00469.035 800 5.0351|)(||||)(|1|)(|22≤≈+?≤+≤t e t s s e t v e 同样v v e v e r )()(≈)()()()(t e s e t e v t t v s e v s s v r r r -=??+??= 所以00205.035 05 .08005.0|)(||)(||)(|≈+≤+≤t e s e v e r r r 因此绝对误差限和相对误差限分别为0.05和0.00205。 例4试建立积分20,,1,05 =+=n dx x x I n n 的递推关系,并研究它的误差 传递。 解:151 --= n n I n I ……………………………………………..…...(1) 5ln 6ln 0-=I ,计算出0I 后可通过(1)依次递推计算出1I ,…,20I 。 但是计算0I 时有误差0e ,由此计算出的1I ,…,20I 也有误差,由(1)可 知近似值之间的递推关系为 151 --= n n I n I ……………………………………………….…..(2) (1)-(2)可得 01)5(5e e e n n n -=-=-,由0I 计算n I 时误差被放大了n 5倍。所以(1)不稳 定。 (1) 可以改写为 n I I n n 51 511+ -=- ……………………………………… (3) 如果能先求出20I ,则依次可以求出19I ,…,0I ,计算20I 时有误差,这样根据(3)计算19I ,…,0I 就有误差,误差传播为 n n n e e ?? ? ??-=-511 ,误差依次减少。 例5 用二分法求解方程012)(23=+--=x x x x f 在区间[0,1]内的1个实根,要求有3为有效数字。 解:因为0)1()0(

数值计算方法期末考精彩试题

1. 已知函数 21 1y x = +的一组数据: 求分段 线性插值函数,并计算 () 1.5f 的近似值. 计算题1.答案 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.35 L =-?=% 4. 写出梯形公式和辛卜生公式,并用来分别计算积分1 01 1dx x +?. 计算题4.答案 4 解 梯形公式 ()()()2b a b a f x dx f a f b -≈ ?+???? 应用梯形公式得 1 01111 []0.75121011dx x ≈+=+++? 辛卜生公式为

确定下列求积公式中的待定系数,并证明确定后的求积公式具有3次代数精确度 ()()()() 1010h h f x dx A f h A f A f h --=-++? 证明题答案

故 ( )()()()40333h h h h f x dx f h f f h -= -++? 具有三次代数精确度。 1.设 3 2 01219 (), , 1, 44f x x x x x ==== (1)试求()f x 在 19,44???? ??上的三次Hermite 插值多项式()x H 使满足''11()(), 0,1,2,... ()()j j H x f x j H x f x === () x H 以升幂形式给出。 (2)写出余项()()()R x f x H x =-的表达式 计算题1.答案 1、(1) ()32142632331 22545045025x x x x H =- ++- (2) ()522191919()(1)(),()(,) 4!164444R x x x x x ξξξ-=---=∈ 3.试确定常数A ,B ,C 和 a ,使得数值积分公式 有尽可能高的代数精度。试问所得的数值积分公式代数精度是多少?它是否为Gauss 型的? 计算题3.答案

数值分析期末试卷

数值分析2006 — 2007学年第学期考试 课程名称:计算方法 A 卷 考试方式:开卷[] 闭卷[V ] 半开卷[] IV 类 充要条件是a 满足 二、(18分)已知函数表如下 1?设 f(0) = 0, f (1) =16 , f( 2) =46,则 f [0,1]= ,f[0,1,2]二 2 ?设 AJ <2 -3 -1 ,则X ,A := A 1 1 j — 3 ?计算积分 xdx ,取4位有效数字。用梯形公式求得的近似值为 "0.5 (辛普森)公式求得的近似值为 ,用 Spsn 4?设f (x )二xe x -3,求方程f (x ) =0近似根的牛顿迭代公式是 ,它的收 敛阶是 5 ?要使求积公式 1 1 [f (x)dx 拓一(0) + A , f (x 1)具有2次代数精度,则 捲= _________________ , 0 4 6 ?求解线性方程组 x 1 ax 2 = 4 , 12_3 (其中a 为实数)的高斯一赛德尔迭代格式收敛的 10 11 12 13 In x 2.3026 2.3979 2.4849 2.5649

三、(20分)构造如下插值型求积公式,确定其中的待定系数,使其代数精度尽可能高, 并指出所得公式的代数精度。 2 f (x)dx : A o f (0) A f (1) A2f(2) o

X 2 4 6 8 y 2 11 28 40 五、(14分)为求方程X ’ -X 2 -1 =0在X o =1.5附近的一个根,将方程改写为下列等价 形式,并建立相应的迭代公式: 试问上述两种迭代公式在 x 0 =1.5附近都收敛吗?为什么?说明理由。 (1)X =1 ?丄,迭代公式 X 1 X k 1 = 1 - X k (2) X 2二1 ,迭代公式 X —1 2 (X k ); X k 1

数值分析期末试题

数值分析期末试题 一、填空题(20102=?分) (1)设??? ? ? ??? ??---=28 3 012 251A ,则=∞ A ______13_______。 (2)对于方程组?? ?=-=-3 4101522121x x x x ,Jacobi 迭代法的迭代矩阵是=J B ?? ? ? ??05.25.20。 (3)3*x 的相对误差约是*x 的相对误差的 3 1倍。 (4)求方程)(x f x =根的牛顿迭代公式是) ('1)(1n n n n n x f x f x x x +-- =+。 (5)设1)(3 -+=x x x f ,则差商=]3,2,1,0[f 1 。 (6)设n n ?矩阵G 的特征值是n λλλ,,,21 ,则矩阵G 的谱半径=)(G ρi n i λ≤≤1max 。 (7)已知?? ? ? ??=1021 A ,则条件数=∞ )(A Cond 9 (8)为了提高数值计算精度,当正数x 充分大时,应将)1ln(2 -- x x 改写为 )1ln(2 ++ -x x 。 (9)n 个求积节点的插值型求积公式的代数精确度至少为1-n 次。 (10)拟合三点))(,(11x f x ,))(,(22x f x ,))(,(33x f x 的水平直线是)(3 1 3 1 ∑== i i x f y 。 二、(10分)证明:方程组? ?? ??=-+=++=+-1 211 2321321321x x x x x x x x x 使用Jacobi 迭代法求解不收敛性。 证明:Jacobi 迭代法的迭代矩阵为 ???? ? ?????---=05 .05 .01015.05.00J B J B 的特征多项式为

2012研究生数值分析课期末考试复习题及答案

一、填空 1. 设 2.3149541...x * =,取5位有效数字,则所得的近似值x= 2.3150 . 2.设一阶差商 ()()()21122114 ,321f x f x f x x x x --= = =---, ()()()322332 615 ,422f x f x f x x x x --= = =-- 则二阶差商 ()123,,______ f x x x =11/6 3. 设(2,3,1)T X =--, 则2||||X = 14 ,=∞||||X 3 。p49 4. 4.求方程 2 1.250x x --= 的近似根,用迭代公式 1.25x x =+,取初始值 01 x =, 那么 1______x =。 1.5 5.解初始值问题 00 '(,)()y f x y y x y =?? =?近似解的梯形公式是 1______k y +≈。 ()()[]11,,2 ++++k k k k k y x f y x f h y 6、 1151A ??= ? -??,则A 的谱半径 = 6 。 7、设 2()35, , 0,1,2,... , k f x x x kh k =+== ,则 []12,,n n n f x x x ++= —————— ————3 和 []123,,,n n n n f x x x x +++= _______________0_____ 。 8、 若线性代数方程组AX=b 的系数矩阵A 为严格对角占优阵,则雅可比迭代和高斯-塞德尔迭代都 收敛 。 9、解常微分方程初值问题的欧拉(Euler )方法的局部截断误差为_______O(h ) ___。

数值分析期末试题

一、(8分)用列主元素消去法解下列方程组: ??? ??=++-=+--=+-11 2123454 321321321x x x x x x x x x 二、(10分)依据下列数据构造插值多项式:y(0)=1,y(1)= —2,y '(0)=1, y '(1)=—4 三、(12分)分别用梯形公式和辛普生公式构造 复化的梯形公式、复化的辛普生公式并利用复化的梯形公式、复化的辛普生公式计算下列积分: ? 9 1dx x n=4 四、(10分)证明对任意参数t ,下列龙格-库塔方法是二阶的。 五、(14分)用牛顿法构造求c 公式,并利用牛顿法求115。保留有效数字五位。 六、(10分)方程组AX=B 其中A=????????? ?10101a a a a 试就AX=B 建立雅可比迭代法和高斯-赛德尔迭代法,并讨论a 取何值时 迭代收斂。 七、(10分)试确定常数A,B,C,a,使得数值积分公式?-++-≈2 2 ) (}0{)()(a Cf Bf a Af dx x f 有尽可能多的 代数精确度。并求该公式的代数精确度。 八、{6分} 证明: A ≤ 其中A 为矩阵,V 为向量. 第二套 一、(8分)用列主元素消去法解下列方程组: ??? ??=++=+-=+3 2221 43321 32132x x x x x x x x 二、(12分)依据下列数据构造插值多项式:y(0)=y '(0)=0, y(1)=y '(1)= 1,y(2)=1 三、(14分)分别用梯形公式和辛普生公式构造 复化的梯形公式、复化的辛普生公式,并利用复化的梯形公式、 复化的辛普生公式及其下表计算下列积分: ?2 /0 sin πxdx ????? ? ? -+-+=++==++=+1 3121231)1(,)1(() ,(),()(2 hk t y h t x f k thk y th x f k y x f k k k h y y n n n n n n n n

数值分析教案

数值分析教案 土建学院 工程力学系 2014年2月 一、课程基本信息 1、课程英文名称:Numerical Analysis

2、课程类别:专业基础课程 3、课程学时:总学时32 4、学分:2 5、先修课程:《高等数学》、《线性代数》、《C 语言》 6、适用专业:工程力学 二、课程的目的与任务: 数值分析是工程力学专业的重要理论基础课程,是现代数学的一个重要分支。其主要任务是介绍进行科学计算的理论方法,即在计算机上对来自科学研究和工程实际中的数学问题进行数值计算和分析的理论和方法。通过本课程的学习,不仅使学生初步掌握数值分析的基本理论知识,而且使学生具备一定的科学计算的能力、分析问题和解决问题的能力,为学习后继课程以及将来从事科学计算、计算机应用和科学研究等工作奠定必要的数学基础。 三、课程的基本要求: 1.掌握数值分析的常用的基本的数值计算方法 2.掌握数值分析的基本理论、分析方法和原理 3.能利用计算机解决科学和工程中的某些数值计算应用问题,增强学生综合运用知识的能力 4.了解科学计算的发展方向和应用前景 四、教学内容、要求及学时分配: (一) 理论教学: 引论(2学时) 第一讲(1-2节) 1.教学内容: 数值分析(计算方法)这门课程的形成背景及主要研究内容、研究方法、主要特点;算法的有关概念及要求;误差的来源、意义、及其有关概念。数值计算中应注意的一些问题。 2.重点难点: 算法设计及其表达法;误差的基本概念。数值计算中应注意的一些问题。3.教学目标: 了解数值分析的基本概念;掌握误差的基本概念:误差、相对误差、误差限、相对误差限、有效数字;理解有效数字与误差的关系。学会选用相对较好的数值计算方法。

数值分析试题及答案

数值分析试题 一、 填空题(2 0×2′) 1. ?? ????-=? ?????-=32,1223X A 设x =是精确值x *=的近似值,则x 有 2 位有效数字。 2. 若f (x )=x 7-x 3+1,则f [20,21,22,23,24,25,26,27]= 1 , f [20,21,22,23,24,25,26,27,28]= 0 。 3. 设,‖A ‖∞=___5 ____,‖X ‖∞=__ 3_____, ‖AX ‖∞≤_15_ __。 4. 非线性方程f (x )=0的迭代函数x =?(x )在有解区间满足 |?’(x )| <1 ,则使用该迭代 函数的迭代解法一定是局部收敛的。 5. 区间[a ,b ]上的三次样条插值函数S (x )在[a ,b ]上具有直到 2 阶的连续导数。 6. 当插值节点为等距分布时,若所求节点靠近首节点,应该选用等距节点下牛顿差 商公式的 前插公式 ,若所求节点靠近尾节点,应该选用等距节点下牛顿差商公式的 后插公式 ;如果要估计结果的舍入误差,应该选用插值公式中的 拉格朗日插值公式 。 7. 拉格朗日插值公式中f (x i )的系数a i (x )的特点是:=∑=n i i x a 0)( 1 ;所以 当系数a i (x )满足 a i (x )>1 ,计算时不会放大f (x i )的误差。 8. 要使 20的近似值的相对误差小于%,至少要取 4 位有效数字。 9. 对任意初始向量X (0)及任意向量g ,线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…) 收敛于方程组的精确解x *的充分必要条件是 ?(B)<1 。 10. 由下列数据所确定的插值多项式的次数最高是 5 。 11. 牛顿下山法的下山条件为 |f(xn+1)|<|f(xn)| 。

数值分析典型例题

数值分析典型例题 Revised as of 23 November 2020

第一章典型例题 例3 ln2=0.…,精确到10-3的近似值是多少 解 精确到10-3=,即绝对误差限是=, 故至少要保留小数点后 三位才可以。ln2 第二章典型例题 例1 用顺序消去法解线性方程组 ??? ??1 -=4+2+4=+2+31-=4++2321 321321x x x x x x x x x 解 顺序消元 ?? ?? ??????---???→???????????---????→???????????--=-?+-?+-?+1717005.555.00141 25.025.105.555.001412142141231412]b A [)3()2/1()2/3(231312r r r r r r 于是有同解方程组 ?? ? ??-==--=++17175.555.01 42332321x x x x x x 回代得解 x 3=-1, x 2=1,x 1=1,原线性方程组的解为X =(1,1,-1)T 例2 取初始向量X (0)=(0,0,0)T ,用雅可比迭代法求解线性方程组 ??? ??5 =+2+23=++1=2-2+321 321321x x x x x x x x x 解 建立迭代格式 ???????+--=+--=++-=+++5223122) (2)(1)1(3 ) (3)(1)1(2 ) (3)(2)1(1k k k k k k k k k x x x x x x x x x (k =1,2,3,…)

第1次迭代,k =0 X (0)=0,得到X (1)=(1,3,5)T 第2次迭代,k =1 ???????-=+?-?-=-=+--==+?+?-=3 532123 351515232)2(3) 2(2)2(1x x x X (2)=(5,-3,-3)T 第3次迭代,k =2 ???????=+-?-?-==+---==+-?+-?-=1 5)3(2521 3)3(511)3(2)3(2)2(3) 3(2) 3(1x x x X (3)=(1,1,1)T 第4次迭代,k =3 ???????=+?-?-==+--==+?+?-=1 512121 311111212)2(3)2(2) 2(1x x x X (4)=(1,1,1)T 例4 证明例2的线性方程组,雅可比迭代法收敛,而高斯- 赛德尔迭代法发散。 证明 例2中线性方程组的系数矩阵为 A =?? ?? ? ?????-122111221 于是 D =??????????100010001 D -1=D ??????????=022001000L ~ ?? ?? ? ?????-=000100220U ~ 雅可比迭代矩阵为

数值分析期末复习题

一、填空题 1.设真值x=983350,则其近似值y=98000的有效数字的位数 ,绝对误差为 , 相对误差为 。 2.x=0.1062,y=0.947,计算x+y 其有效数字的位数为 。 3.对f(x)=x 3 +x+1,差商f[0,1,2,3]= ;f[0,1,2,3,4]= 。 4.设f(x)可微,求方程x=f(x)根的牛顿迭代法格式是 。 5.设方程x=?(x)有根x * ,且设?(x)在含x * 的区间(a,b)内可导,设x 0∈(a,b)则迭代格式x k+1=?(x k )收敛的充要条件为 。 6.求解线性方程组Ax=b 的迭代格式x (k+1)=Jx (k)+f 收敛的充要条件为 。 7.??? ? ??=01100 1001001....A ,||A||∝= ,cond(A)∝= 。 8.n 次Legendre 多项式的最高次项系数为 。 9.中矩形公式:)()2( )(a b b a f dx x f b a -+=?的代数精度为 。 10.求积公式:)1(2 1)0()(10 f f dx x f '+ ≈?的代数精度为 。 11.在区间[1,2]上满足插值条件? ??==3)2(1 )1(P P 的一次多项式P(x)= 。 12.设∑ == n k k k n x f A f I 0 )()(是函数f(x)在区间[a,b]上的插值型型求积公式,则 ∑=n k k A = 。 13.梯形公式和改进的Euler 公式都是 阶精度的。 二、计算题 1.利用矩阵的高斯消元法,解方程组??? ??=++=++=++20 53182521432321 321321x x x x x x x x x 2.设有函数值表 试求各阶差商,并写出Newton 插值多项式。

数值分析试题及答案汇总

数值分析试题及答案汇 总 TTA standardization office【TTA 5AB- TTAK 08- TTA 2C】

数值分析试题 一、填空题(2 0×2′) 1. ?? ????-=? ?????-=32,1223X A 设x =是精确值x *=的近似值,则x 有 2 位有效数字。 2. 若f (x )=x 7-x 3+1,则f [20,21,22,23,24,25,26,27]= 1 , f [20,21,22,23,24,25,26,27,28]= 0 。 3. 设,‖A ‖∞=___5 ____,‖X ‖∞=__ 3_____, ‖AX ‖∞≤_15_ __。 4. 非线性方程f (x )=0的迭代函数x =(x )在有解区间满足 |’(x )| <1 ,则使用该迭代函数 的迭代解法一定是局部收敛的。 5. 区间[a ,b ]上的三次样条插值函数S (x )在[a ,b ]上具有直到 2 阶的连续导数。 6. 当插值节点为等距分布时,若所求节点靠近首节点,应该选用等距节点下牛顿差 商公式的 前插公式 ,若所求节点靠近尾节点,应该选用等距节点下牛顿差商公式的 后插公式 ;如果要估计结果的舍入误差,应该选用插值公式中的 拉格朗日插值公式 。 7. 拉格朗日插值公式中f (x i )的系数a i (x )的特点是:=∑=n i i x a 0)( 1 ;所以当系数 a i (x )满足 a i (x )>1 ,计算时不会放大f (x i )的误差。 8. 要使 20的近似值的相对误差小于%,至少要取 4 位有效数字。 9. 对任意初始向量X (0)及任意向量g ,线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…) 收敛于方程组的精确解x *的充分必要条件是 (B)<1 。 10. 由下列数据所确定的插值多项式的次数最高是 5 。

相关文档
最新文档