您的位置:百味书屋 > 试题答案 > 数值分析习题集及答案 正文

数值分析习题集及答案

2016-11-05 10:51:35 来源网站:百味书屋

篇一:数值分析习题集及答案[1]

数值分析习题集

(适合课程《数值方法A》和《数值方法B》)

长沙理工大学

第一章 绪 论

1. 设x>0,x的相对误差为δ,求lnx的误差.

2. 设x的相对误差为2%,求x的相对误差.

3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出

它们是几位有效数字: 4. 利用公式(3.3)求下列各近似值的误差限:

*****x1?1.1021,x2?0.031,x3?385.6,x4?56.430,x5?7?1.0.

n

************

(i)x1?x2?x4,(ii)x1x2x3,(iii)x2/x4,其中x1,x2,x3,x4

均为第3题所给的数.

5. 计算球体积要使相对误差限为1%,问度量半径R时允许的相对误差限是多少? 6. 设Y0?28,按递推公式

( n=1,2,…)

Y计算到Y100.

27.982(五位有效数字),试问计算100将有多大误差?

Yn?Yn?12

7. 求方程x?56x?1?0的两个根,使它至少具有四位有效数字

27.982).

8. 当N充分大时,怎样求

?

??

N

1dx1?x2?

2

9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝? 10. 设

误差增加,而相对误差却减小. 11. 序列

S?

12gt2假定g是准确的,而对t的测量有±0.1秒的误差,证明当t增加时S的绝对

{yn}满足递推关系yn?10yn?1?1(n=1,2,…),

若y0??1.41(三位有效数字),y10时误差有多大?这个计算过程稳定吗?

计算到

6

12.

计算f?1),

?1.4,利用下列等式计算,哪一个得到的结果最好?

3

?? 13.

f(x)?ln(x,求f(30)的值.若开平方用六位函数表,问求对数时误差有多大?若

改用另一等价公式

计算,求对数时误差有多大?

ln(x???ln(x

14. 试用消元法解方程组

?

x1?1010x2?1010;x1?x2?2.

假定只用三位数计算,问结果是否可靠?

1?absinc,0?c?22,且测量a ,b ,c 的误差分别为15. 已知三角形面积其中c为弧度,

?a,?b,?c.证明面积的误差?s满足

s?

?s?a?b?c

???.sabc

第二章 插值法

1. 根据(2.2)定义的范德蒙行列式,令

1

Vn(x)?Vn(x0,x1,?,xn?1,x)?

11

证明Vn(x)是n次多项式,它的根是x0,?,xn?1,且

x0?xn?1x

2

x0

???

nx0

?x2

?xn

2nxn?xn?1?1

Vn(x)?Vn?1(x0,x1,?,xn?1)(x?x0)?(x?xn?1).

2. 当x= 1 , -1 , 2 时, f(x)= 0 , -3 , 4 ,求f(x)的二次插值多项式.

3.

4. 给出cos x,0°≤x ≤90°的函数表,步长h =1′=(1/60)°,若函数表具有5位有效数字,

研究用线性插值求cos x 近似值时的总误差界.

maxl2(x)x?x?khx0?x?x3k05. 设,k=0,1,2,3,求.

xj

6. 设

为互异节点(j=0,1,…,n),求证:

i) ii) 7. 设

?xl(x)?x

kjjj?0n

n

k

(k?0,1,?,n);

?(x

j?0

j

?x)klj(x)???k?1,2,?,n).

2

f(x)?C?a,b?且f(a)?f(b)?0,求证maxa?x?b

x

?6

1

f(x)?(b?a)2maxf?(x8a?x?b

x

8. 在?4?x?4上给出f(x)?e的等距节点函数表,若用二次插值求e的近似值,要使截

断误差不超过10,问使用函数表的步长h应取多少? 9. 若yn?2,求?yn及?yn.

10. 如果f(x)是m次多项式,记?f(x)?f(x?h)?f(x),证明f(x)的k阶差分

n

4

4

?kf(x)(0?k?m)是m?k次多项式,并且?m?lf(x)?0(l为正整数).

11. 证明?(fkgk)?fk?gk?gk?1?fk. 12. 证明k?0

n?1n?1

n?1

?f?g

k

k

?fngn?f0g0??gk?1?fk.

k?0

13. 证明

??

j?0

2

yj??yn??y0.

n?1n

14. 若f(x)?a0?a1x???an?1x?anx有n个不同实根x1,x2,?,xn,证明

?f?(x)

j?1

j

n

xkj

?

?

0,0?k?n?2;

?1an,k?n?1.

15. 证明n阶均差有下列性质: i)

若F(x)?cf(x),则

F?x0,x1,?,xn??cf?x0,x1,?,xn?;

Fx,x,?,xn??f?x0,x1,?,xn??g?x0,x1,?,xn?.

ii) 若F(x)?f(x)?g(x),则?01

74f?20,21,?,27?f?20,21,?,28?f(x)?x?x?3x?1????. 16. ,求及

17. 证明两点三次埃尔米特插值余项是

并由此求出分段三次埃尔米特插值的误差限.

R3(x)?f(4)(?)(x?xk)2(x?xk?1)2/4!,??(xk,xk?1)

18. 求一个次数不高于4次的多项式P(x),使它满足P(0)?P(?k?1)并由此求出分段三次

埃尔米特插值的误差限. 19. 试求出一个最高次数不高于4次的函数多项式P(x),以便使它能够满足以下边界条件

P(0)?P?(0)?0,P(1)?P?(1)?1,P(2)?1.

20. 设

f(x)?C?a,b?,把?a,b?分为n等分,试构造一个台阶形的零次分段插值函数?n(x)

并证明当n??时,?n(x)在?a,b?上一致收敛到f(x).

2

f(x)?1/(1?x),在?5?x?5上取n?10,按等距节点求分段线性插值函数Ih(x),21. 设

计算各节点间中点处的Ih(x)与f(x)的值,并估计误差.

a,b?上的分段线性插值函数Ih(x),并估计误差.

22. 求f(x)?x在?

24

23. 求f(x)?x在?a,b?上的分段埃尔米特插值,并估计误差.

试求三次样条插值并满足条件 i) ii)

2

f(x)?C?a,b?,S(x)是三次样条函数,证明 25. 若

S?(0.25)?1.0000,S?(0.53)?0.6868; S?(0.25)?S?(0.53)?0.

i)

??f?(x)?dx???S?(x)?dx???f?(x)?S?(x)?dx?2?

a

a

a

b

2

b

2

b

2

b

a

S?(x)?f?(x)?S?(x)?dx

;

ii) 若f(xi)?S(xi)(i?0,1,?,n),式中xi为插值节点,且a?x0?x1???xn?b,则

?

b

a

S?(x)?f?(x)?S?(x)?dx?S?(b)?f?(b)?S?(b)??S?(a)?f?(a)?S?(a)?

.

26. 编出计算三次样条函数S(x)系数及其在插值节点中点的值的程序框图(S(x)可用(8.7)

式的表达式).

第三章 函数逼近与计算

1. (a)利用区间变换推出区间为?

a,b?的伯恩斯坦多项式.

?上求1次和三次伯恩斯坦多项式并画出图形,并与相应的(b)对f(x)?sinx在?

马克劳林级数部分和误差做比较. 2. 求证:

0,?/2

(a)当m?f(x)?M时,m?Bn(f,x)?M. (b)当f(x)?x时,Bn(f,x)?x.

0,2??的最佳一致逼近多项式.

3. 在次数不超过6的多项式中,求f(x)?sin4x在?

a,b?上连续,求f(x)的零次最佳一致逼近多项式.

4. 假设f(x)在?

5. 选取常数a,使0?x?1

maxx3?ax

达到极小,又问这个解是否唯一?

0,?/2?上的最佳一次逼近多项式,并估计误差.

6. 求f(x)?sinx在?

0,17. 求f(x)?e在??上的最佳一次逼近多项式.

x

?1,1?上与零偏差最小?r是否唯一?

8. 如何选取r,使p(x)?x?r在?

2

0,19. 设f(x)?x?3x?1,在??上求三次最佳逼近多项式.

4

3

***

T(x)?T(2x?1),x?0,1??T(x),T(x),T(x),T3(x). nn01210. 令,求

11. 试证12. 在?

?T

*

n

(x)?

是在?0,1

?上带权

??

的正交多项式.

?1,1?上利用插值极小化求1f(x)?tg?1x的三次近似最佳逼近多项式.

?

x?1,1?上的插值极小化近似最佳逼近多项式为Ln(x),若f?Ln

13. 设f(x)?e在?

有界,

证明对任何n?1,存在常数?n、?n,使

?nTn?1(x)?f(x)?Ln(x)??nTn?1(x)(?1?x?1).

112331541655

?(x)?1?x?x?x?x?x?1,1??28243843840,试将?(x)降低到3次多14. 设在上

项式并估计误差. 15. 在?

?1,1?上利用幂级数项数求f(x)?sinx的3次逼近多项式,使误差不超过0.005.

?a,a?上的连续奇(偶)函数,证明不管n是奇数或偶数,f(x)的最佳逼近多项式

16. f(x)是?

Fn*(x)?Hn也是奇(偶)函数.

?ax?b?sinx?dx为最小.并与1题及6题的一次逼近多项式误差作比较.

17. 求a、b使?

1

g(x)?C?a,b?,定义 f(x)18. 、

2

?

2

(a)(f,g)??f?(x)g?(x)dx;(b)(f,g)??f?(x)g?(x)dx?f(a)g(a);

a

a

bb

问它们是否构成内积?

1

x6

?01?x19. 用许瓦兹不等式(4.5)估计的上界,并用积分中值定理估计同一积分的上下界,

并比较其结果.

20. 选择a,使下列积分取得最小值:21. 设空间

?

1

?1

(x?ax2)2dx,?x?ax2dx

?1

1

.

???span?1,x?,?2?span?x100,x101?

,分别在?1、?2上求出一个元素,使得其为

x2?C?0,1?的最佳平方逼近,并比较其结果.

?1?span?1,x2,x4?f(x)?x?1,1??22. 在上,求在上的最佳平方逼近.

sin(n?1)arccosxun(x)?

23.

是第二类切比雪夫多项式,证明它有递推关系

un?1?x??2xun?x??un?1?x?.

24. 将

近多项式并画出误差图形,再计算均方误差.

f(x)?sin

1

x??1,1?2在上按勒让德多项式及切比雪夫多项式展开,求三次最佳平方逼

?1,1?上展成切比雪夫级数.

25. 把f(x)?arccosx在?

26.

y?a?bx.

2

27.

用最小二乘拟合求.

29. 编出用正交多项式做最小二乘拟合的程序框图. 30. 编出改进FFT算法的程序框图. 31. 现给出一张记录?

xk???4,3,2,1,0,1,2,3?,试用改进FFT算法求出序列?xk?的离散频谱

?Ck?(k?0,1,?,7).

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

篇二:数值分析习题解答1

第一章 引论(习题)

2.证明:x的相对误差约等于x的相对误差的1/2.

证明记 f(x)?

Er(f)?x ,则 ?x?x*

x(x?x)*x?x*x

?x?x*1??Er(x).□ *x2x?xx

3.设实数a的t位?进制浮点机器数表示为fl(a). 试证明fl(a?b)?(a?b)/(1??),|?|?

其中的记号*表示+、-、?、/ 中一种运算.

证明: 令: ??11?t?, 2(a?b)?fl(a?b) fl(a?b)

c?1可估计: |fl(a?b)|??

故: |?|? (c为a?b阶码), 1c?tc?111?t???? 22

于是:fl(a?b)?(a?b)(1??). □

4.改变下列表达式使计算结果比较精确:

(1)

(2)

(3) 11?x?,1?2x1?xx?1?xx?1,x对|x|??1; 对x??1; 1?cosx,x

2对x?0,|x|??1. 解 (1) 2x

(2) (1?x)(1?2x). . x

(x?x?x?x)

1?cosxsin2xsinx?? (3) . □ xx(1?cosx)1?cosx

6.设a?0.937关于精确数x有3位有效数字,估计a的相对误差. 对于f(x)??x,估计f(a)对于f(x)的误差和相对误差.

解 a的相对误差:由于|E(x)|?x?a?1x?a, ?10?3. Er(x)?2x

11Er(x)?10?2??10?2.(Th1) 2?918

f(a)对于f(x)的误差和相对误差. |E(f)|?|?x??a|=a?x

?x??a??3???102

2?0.25=10?3 |Er(f)|?10

?3?a?4?10?3.□

9.序列{yn}满足递推关系:yn?1?100.01yn?yn?1. 取y0?1,y1?0.01及y0?1?10?5,

定的. y1?0.01,试分别计算y5,从而说明该递推公式对于计算是不稳

解递推关系: yn?1?100.01yn?yn?1

(1) 取初值 y0?1, y1?0.01 计算

可得: y2?100.01?10

?6?2?1?1.0001?1?10?4 ?8y3?10 , y4?10

(2) 取初值 0?1?10

记: ?n?yn?n, ?5 , y5?10?2?10 , … , 1?10,

序列 ??n? ,满足递推关系,且?0??10?5 , ?1?0

?n?1?100.01?n??n?1, 于是: ?2?10?5,

?3?100.01?10?5, ?4?(100.01)2?10?5?10?5,

?5?(100.01)?103?5?200.02?10?5,

n?2 可见随着 ?n 的主项 (100.01)?10?5 的增长,说明该递推关系式是不稳定的.

篇三:《数值分析》习题1

习题1

1. 以下各表示的近似数,问具有几位有效数字?并将它舍入成有效数。

(1)x1=451.023, x1=451.01;

*(2)x2=-0.045 113, x2=-0.045 18; *

*(3)x3=23.421 3,x3=23.460 4;

(4)x*

4=1, x4=0.333 3; 3

*(5)x5=23.496, x5=23.494;

*(6)x6=96×105,x6=96.1×105;

*(7)x7=0.000 96, x7=0.96×10?3;

*(8)x8=-8 700, x8=-8 700.3。

解:(1) x1?451.023 x1?451.01

*x1*1?x1?0.013??10?1,x1具有4位有效数字。x1?451.0 2

*(2) x2??0.045 113 x2??0.045 18

11*?10?4?x2?x2?0.045 18?0.045113=0.000 067??10?3 22

x2具有2位有效数字,x2??0.045

*(3)x3?23.4213 x3?23.4604 *x31?x3?23.4213?23.4604?23.4604?23.4213?0.0391??10?1 2x3 具有3位有效数字,x3?23.4 (不能写为23.5)

(4) x4?*1 ,x4?0.3333 3

*x4?x4?0.000033???10?4 ,x4具有4位有效数字,x4?0.3333

*?23.496,x5?23.494 (5) x512

1*x5?x5?23.496?23.494?0.002??10?2 2

x5 具有4位有效数字, x5?23.50 (不能写为23.49)

*?96?105?0.96?107 x6?96.1?105?0.961?107 (6) x6

1*x6?x6?0.001?10?7??10?2?10?7 2

x6具有2位有效数字,x6?0.96?107?96?105

*?0.00096 x7?0.96?10?3 (7) x7

**x7?0.96?10?3 x7?x7?0 x7精确

*??8700 x8??870.03 (8) x8

*x81?x8?0.3??100x8具有4位有效数字,x8??8700精确 2

2.以下各数均为有效数字:

(1) 0.1062 + 0.947;

(3)2.747?6.83; (2)23.46―12.753; (4)1.473 / 0.064 。

问经过上述运算后,准确结果所在的最小区间分别是什么?

解:(1) x1=0.1062,x2=0.947,x1+x2=1.0532

11?4 e(x1)??10,e(x2)??10?3 22

e(x1?x2)?e(x1)?e(x2)?e(x1)?e(x2)?

=0.00055 11?10?4??10?3 22

**x1?x2?[1.0532?0.00055,1.0532+0.00055]=[1.05265,1.05375]

(2) x1=23.46, x2??12.753 x1?x2?10.707

11?2,??10??10?3 e(x1)e(x2)22

e(x1?x2)?e(x1)?e(x2)?e(x1)?e(x2)

11?2??10??10?3=0.0055 22

**x1?x2?[10.707?0.0055, 10.707+0.0055]=[10.7015,10.7125]

(3) x1?2.747 x2?6.83 x1x2?18.76201,

11?3e(x1)??10, e(x2)??10?2 22

e(x1x2)?x2e(x1)?x1e(x2)?x2e(x1)?x1e(x2)

111?6.83??10?3?2.747??10?2??10?2?(0.683+2.747)=0.01715 222

**x1x2?[18.76201?0.01715,18.76201?0.01715]?[18.74486,18.77916]

(4) x1?1.473 ,x2?0.064 ,x1x2?23.015625

x111xe(x2) e(x1)??10?3,e(x2)??10?3 e(1)?e(x1)?1

22x222x2

e(x12)?x1111.4731?3?3e(x1)?1e(x)? ??10???10222x20.06422x20.064

=0.187622

*x1*x2?[23.015625?0.187622, 23.015625+0.187622]

=[22.828003 , 23.203247]

3.对一元2次方程x2?40x?1?0,如果399?19.975具有5位有效数字,求其具有5位有效数字的根。

解:x2?40x?1?0

x2?40x?400?399

**x1?20? ,x2?20?399?1

20?

1记 x*? ,x?19.975e(x)??10?3 2

1?x20=20+19.975=39.975 x1?e(x1)?e(x2)??10?3 2

? x1具有5位有效数字。

111x2????0.0250156347? 20?x20?19.97539.975

e(x2)??e(x)

(20?x)2 ,

1?3?10e(x)1?6?6e(x2)?? ?0.313?10??10222(20?x)39.975

因而 x2具有5位有效数字。 x2?0.025016

也可根据 x1x2?1 得到 x2?11??0.0250156347? x139.975

1?6?10e(x)e(x)e(x2)??21 e(x2)?21? 2x139.975x1

4.若x1?0.937具有3位有效数字,问x1的相对误差限是多?设f(x)??x,求f(x1)的绝对误差限和相对误差限。

1 解:x1?0.937e(x1)??10?3 2

1?10?3e(x)er(x1)???0.534?10?3 x10.937

f(x)??x ,f?(x)??1 2?x

e(f)?f?(x)e(x)??11?e(x) , 2?x

e(f(x1))?11111?e(x1)????10?3?0.996?10?3 2?x12?0.9372

e(f)11e(x), er(f)????21?xf

er(f(x1))?11111?e(x1)????10?3 21?x121?0.9372

=0.00397?3.97?10?3

5.取2.01?1.42,2.00?1.41试按A?2.01?2.00和A?0.(2.01?2.00)两种算法求A的值,并分别求出两种算法所得A的近似值的绝对误差限和相对误差限,问两种结果各至少具有几位有效数字?

**?2.01 ,x1?1.42 ,x2 解:1) 记 x1=2.00 ,x2?1.41

则 e(x1)?11?10?2 ,e(x2)??10?2 22

A* ?2.01?2.00?1.42?1.41?0.01

A1?1.42?1.41?0.01

e(A1)?e(x1?x2)?e(x1)?e(x2)

11e(A1)?e(x1)?e(x2)?e(x1)?e(x2)??10?2??10?2?10?2 22


数值分析习题集及答案》出自:百味书屋
链接地址:http://www.850500.com/news/18681.html
转载请保留,谢谢!
查看更多相关内容>>数值分析习题集及答案
相关文章
推荐范文