数值分析习题与答案试卷及答案.docxVIP

数值分析习题与答案试卷及答案

考试时间：______分钟总分：______分姓名：______

一、

1.设计算法计算某根f(x)=0的近似值，若该算法具有线性收敛速度，且初始近似值x0与真根x*的距离为0.1，要使近似值达到10^(-4)的精度，大约需要迭代多少次？

2.解释什么是舍入误差和截断误差，并举例说明它们在数值计算中可能产生的影响。

3.给定方程x^3-x-1=0，试用二分法求其在区间[1,2]内的根，要求误差不超过10^(-3)。请写出关键步骤和计算结果。

4.写出牛顿迭代法求解f(x)=0的公式，并分析其收敛性的条件。

二、

1.已知矩阵A=[[3,2],[1,2]]，试用高斯消元法将A化为行阶梯形矩阵。

2.什么是LU分解？若矩阵A可以进行LU分解为A=LU，其中L是单位下三角矩阵，U是上三角矩阵，请写出解线性方程组Ax=b的LU分解方法。

3.写出雅可比迭代法和高斯-赛德尔迭代法求解线性方程组Ax=b的迭代公式，并简述判断这两种迭代法收敛的基本条件。

4.对于线性方程组Ax=b，若雅可比迭代法收敛，则高斯-赛德尔迭代法是否一定收敛？请说明理由。

三、

1.给定数据点(x0,y0),(x1,y1),...,(xn,yn)，写出拉格朗日插值多项式L(x)的定义式。

2.什么是插值误差？请给出拉格朗日插值误差的估计公式，并说明其中各符号的含义。

3.什么是分段线性插值？与拉格朗日插值相比，它有哪些优缺点？

4.给定数据点(0,1),(1,3),(2,2),(3,4)，求通过这四个点的三次样条插值函数S(x)的自然边界条件，并写出边界条件公式。

四、

1.写出用梯形法则计算定积分∫[a,b]f(x)dx的公式。

2.什么是辛普森法则？它与梯形法则相比，在精度上有何优势？请说明其公式。

3.什么是复化梯形法则？请写出计算∫[a,b]f(x)dx的复化梯形法则公式（n等分）。

4.试用辛普森法则计算∫[0,1]e^xdx，要求n=4。请写出计算过程和结果。

五、

1.写出欧拉法求解初值问题y=f(x,y),y(x0)=y0的公式。

2.什么是改进欧拉法（Heun法）？请写出其计算公式，并说明它与欧拉法相比有何改进。

3.什么是龙格-库塔法（RK4）？请写出用RK4方法求解初值问题y=f(x,y),y(x0)=y0的公式。

4.比较欧拉法、改进欧拉法和RK4法在求解初值问题时，在相同步长下的收敛速度和精度。

试卷答案

一、

1.因为算法具有线性收敛速度，所以误差满足|x_{k+1}-x*|≈c|x_k-x*|，其中c为收敛常数(0c1)。设初始误差为|x0-x*|=ε0=0.1，目标误差为εT=10^(-4)。迭代k次后的误差为|x_k-x*|≈c^k*ε0。要求|x_k-x*|εT，即c^k*ε0εT。取对数得k*ln(c)+ln(ε0)ln(εT)。由于c1，ln(c)0，所以k(ln(εT)-ln(ε0))/ln(c)。取c=1/2作为保守估计（实际c值未知，收敛速度取决于算法和问题），k(ln(10^(-4))-ln(0.1))/ln(0.5)=(ln(10^(-4)/0.1))/ln(0.5)=(ln(10^(-3)))/ln(0.5)≈6.906/-0.693≈9.97。因此，大约需要迭代10次。

2.舍入误差是在数值计算中由于计算机表示有限位数而丢失的精度。例如，将1/3存储为有限小数。截断误差是由于用近似公式（如有限项级数、数值积分公式）代替精确公式而引起的误差。例如，用sin(x)≈x替代sin(x)的精确值。两者都可能导致计算结果与真值存在偏差。

3.f(1)=1-1-1=-10，f(2)=8-2-1=50。根在[1,2]内。中点x1=(1+2)/2=1.5。f(1.5)=3.375-1.5-1=0.8750。根在[1,1.5]内。x2=(1+1.5)/2=1.25。f(1.25)=1.953125-1.25-1=-0.2968750。根在[1.25,1.5]