计算方法复习题.pdfVIP

计算法复习题

软13计算法复习题

1、对下的计算式做适当的等价变换，以避免两个相近的数相减时的精度损失。

（1）)ln()1ln(xx-+，其中x较

（2）xx-+12，其中x较

222、已知函数程0)ln(3)(=--=xxxf有正根，请完成以下的作：

（1）分析并选定个含有这正根的区间[a0,b0]，以便于分法求解；

（2）验证在[a0,b0]上分法求根的可性，并计算逐步缩的区间[a1,b1]和[a2,b2]；

（3）若考虑简单迭代法求此根，试构造个在[a0,b0]上能保证收敛的迭代式)(1kkxx?=+。解：（1）把程的根看成

y=3-x和y=ln(x)的交点，经分析可取含根区间[1.0,3.0]（2）经验算可得f(1.0)*f(3.0)0，另f’(x)在[1.0,3.0]上不变号，f(x)单

调，分法可（3）迭代式)ln(31kkxx-=+从迭代收敛定理两作完整讨论，知迭代式能保证收敛

3、Doolittle分解法求解线性程组

=?564221231112321xxx（要求写明求解过程）。解：（1）先对系数矩阵A作LU分解得

A=LU=??

5/32/32/511

215/32/112/11

（2）由LY=B解出Y=(4,4,3/5)T，由UX=Y解出X=(1,1,1)T

4、关于某函数y=f(x)，已知如下表所的批数据

（1）由上表中的数据构建差商表，并求出各阶差商；（2）分别点、三点顿插值法计算f(0.75)的近似值；

（3）若bx

aey=来拟合这批数据，试求出系数a和b（提：两边取然对数得lny=lna+bx，

令u=lny，问题转化为求拟合直线u=lna+bx）；

（4）分别复化梯形积分和复化普森积分计算

20

)(dxxf的近似值。

=5.75ln5.7555ba得a=1，b=1，故有x

ex

f=)(（4）分别复化梯形积分公式和复化普森积分公式计算

5、若Jacobi迭代法求解线性程组??

=++-=-+=+

5zyxzyxzyx：

（1）能否从系数矩阵判定Jacobi迭代求解是收敛的？请说明原因；（2）写出经过等价变换得到的Jacobi迭代格式fBXX

kk+=+1；

（3）求出迭代矩阵B的范数∞

B

和列范数1B，并说明B能否保证收敛。

6、规范化幂法求矩阵?

=1403A的按模最特征值，使误差不超过1

105.0-?。初始向量取为V(0)=(1,1)T。（另：若给出规范化幂法迭代计算的向量序列，你是否掌握根据向量序列的收敛情

况计算按模最特征值和特征向量的法。）7、改进欧拉法求初值问题==0

.1)0.0(/yxy

dxdy在区间[0.0,1.0]上的解，取步长h=0.2。计算结果保留到

数点后3位。

8、）对于函数)1()(xxxxf-+=，按下两种法计算)1000

(f的近似值，分别讨论两个结果的绝对误差限和有效数字的位数，并说明产差别的原因。（特别注意：计算过程按四位舍

法进。例如2

103162.01000?≈，2

103164.01001?≈）

（1）直接按表达式计算；（2）按等价变换式)1/()(xxxxf+

+=计算。

8、答题要点

精确值f(1000)=0.1580743(102)

（1）f1(1000)≈1000*(0.3164-0.3162)*102=0.2*102，与精确值较得绝对误差限ε1=0.5*101

，得有效数字位数为1位；

（2）f2(1000)≈1000/(0.3164*102+0.3162*102)≈0.1581*102

，与精确值较得绝对误差限为

ε2=0.5*10*10-2

,得有效数字的位数为4位。

原因在于直接按表达式计算时两个相近的数相减导致有效数字位数减少误差增

9、已知函数程052)(3

=-+=xxxf在区间[2，3]上有根（令a0=2，b0=3）：（1）验证在此区间上分法求根的可性，并计算逐步缩的区

间[a1,b1]和[a2,b2]；（2）若简单迭代法求此根，试分析并构造个在[a