上海大学数值分析历届考题

-数值分析历届考题03-04学年秋季学期一 简答题每题5分1. 数值计算中要注意哪些问题。答：第一、两个相近的数应防止相减。第二、绝对值很小的数应防止作除数。第三、注意选取适当的算法减少运算次数。第四、两个绝对值相差很大的数运算时，注意机器零的问题。第五、注意算法的收敛性和稳定性。2. 用迭代法求解非线性方程时，迭代收敛的条件是什么，可以用什么方法来确定初值。答：对于非线性方程其迭代格式为，如果满足：（1） 当时，；（2） 在上连续，且对任意的都有。则有结论：对任意给定的，由迭代格式，k=0,1,2,产生的序列收敛于，即迭代收敛。可以用二分法来确定初值。3. 用消元法求解线性方程组时，为什么要选主元。答：因为用简单高斯消元法求得的近似解与准确解相差甚远，其主要原因是绝对值很小的数作除数，导致了误差的快速增长。为了防止这种情况的发生，我们可以通过行交换，在需要消元的列中，取绝对值最大者作为主对角线元素即主元，计算效果将得到改善。4. 矩阵的条件数是什么，它对求解线性方程组有什么影响。答：对于n阶可逆方阵A，正实数|称为A的条件数，记为cond(A)。条件数对于线性方程组A*=b的影响如下：，其中为A准确时b产生的误差；，其中为b准确时A产生的误差。5. 把以下二阶常微分方程的初值问题化为一阶常微分方程组，并写出求解该方程的改进Euler方法。答：令则，其中。所以用改进的Euler方法表示为：，。二 20分给出数据表*012f(*)212f(*)-1求一个满足插值条件的三次插值多项式，并写出余项公式。解：先求出满足函数值插值条件，i=0，1，2的二次插值多项式。i*f(*)一阶差商二阶差商102211-132211由牛顿插值公式：令，其中A是待定常数，则，由条件，代入可得：；所以。其插值余项为，其中。三 20分给出数据表*0.10.20.40.5y10.80240.61740.53023用最小二乘法求拟合曲线保存3位小数。解：对于曲线，令，得。把*，y的数据转换为t，z的数据取3位有效数字：t=1/*2.002.505.0010.0z=1/y1.891.621.251.00对于，其法方程组为：；其中：，数据代入后得法方程组为；解得。所以拟合曲线为。四 15分确定以下求积公式的系数，使公式成为Guass型求积公式。解：通过待定系数法：当时，有1当时，有2当时，有3由此得到一个关于未知数，的线性方程组：；解得。五 20分证明：对任意参数t以下求解常微分方程初值问题的算法，其局部截断误差都是c：。证：令，则1对作泰勒展开得：。代入到1式中：由于在的条件下。即对任意参数t，上述求解微分方程初值问题的算法其局部截断误差都是。六 16分证明：以下求解常微分方程初值问题的数值方法，其局部截断误差为。证：在的条件下将上述两式代入中，可得：由于在的条件下。所以上述求解微分方程初值问题的算法其局部截断误差都是。05-06学年秋季学期一 简答题每题4分，共20分1. 设*=0.06020，y=0.0418是按四舍五入得到的近似值，则*+y，*y的绝对误差限，相对误差限，有效数字各是多少。答：，；，所以*+y三位有效，；，所以*/y三位有效，2. 同03-04学年秋季学期第一题33. 在解线性方程组时，原始数据的误差对解的影响如何；对病态方程组可以采用什么方法处理。答：原始数据的误差对于线性方程组A*=b的影响如下：，其中为A准确时b产生的误差；，其中为b准确时A产生的误差；其中cond(A)=|为条件数。对于病态方程组，可以使用迭代改善的方法处理。4. 给出三个等距节点，及其相应的函数值，试导出二阶数值导数的计算公式。答：以这三个点为节点的根本插值多项式为：，；求二阶导得：，；设，i=0，1，2。则。5. 用数值方法求解常微分方程时，怎样选择适宜的步长。答：先选取一个步长h，计算和，如果，则将步长逐次减半，直到为止。如果对于初始步长h，就有，则尝试将步长逐次加倍，知道满足的最大步长。二 16分给出数据表*123f(*)2412f(*)3求一个3次插值多项式；并证明其余项公式为解：先求出满足函数值插值条件，i=0，1，2的二次插值多项式。i*f(*)一阶差商二阶差商1122242331283由牛顿插值公式：令，其中A是待定常数，则，由条件，代入可得：；所以。由插值条件可知，是R(*)的二重零点，和是R(*)的单重零点，所以，其中K(*)是待定函数。令，当的4阶导数连续时，反复用罗尔定理，可得，所以。三 16分给出一组数据*1.001.251.501.752.00Y8.467.456.535.795.10用最小二乘法求拟合曲线。解：对于曲线，两边取对数得：令，则可得到：把*，y的数据转换为t，z的数据取3位有效数字：t=1/*0.5000.5710.6670.8001.00z=lny1.631.761.882.012.14对于，其法方程组为：；其中：，数据代入后得法方程组为；解得。所以拟合曲线为。四 16分用龙贝格方法求以下积分，要求5位有效数字。解：；。五 16分对于非线性方程f(*)=0，求证：改进的牛顿迭代格式：，k=0，1，在单根附近是至少三阶收敛的。并判别该方法对重根是几阶收敛。解：1在单根的情况下，设是的单重根。，所以是的二重零点，即该迭代格式是三阶收敛的。2在重根的情况下，设是的m重根。m1则，且，同理：这时：由于m为大于1的整数，所以显然，所以在重根情况下题设迭代法线性收敛。一阶收敛06-07学年冬季学期一、 简答题(每题4分，共20分)1. 设*=-0.0307，y=1.230是按四舍五入得到的近似值，则*-y，*/y的绝对误差限，相对误差限，有效数字各是多少。答：，；，所以*-y三位有效，；，所以*/y三位有效，2. 插值型数值积分方法的根本原理是什么，其截断误差是什么。答：根本原理：，其中是的n次插值多项式。截断误差：3. 写出求解非线性方程组，i=1,2,n一般迭代法的迭代格式和收敛条件。答：一般迭代法的格式：，i=1,2,n，其中：是的等价方程。当，时收敛。4. 同03-04秋季学期第一题45. 把以下二阶常微分方程的初值问题化为一阶常微分方程组的初值问题，并写出数值求解的欧拉格式。答：令则，其中。所以用欧拉形式表示为：，i=0,1,2,n-1。二、 16分给出数据表*012f(*)129f(*)3用3次插值多项式求f(1.5)的近似值，并估计误差：解：先求出满足函数值插值条件，i=0，1，2的二次插值多项式。i*f(*)一阶差商二阶差商101212132973由牛顿插值公式：令，其中A是待定常数，则，由条件，代入可得：；所以。三、 16分给出一组数据*1.001.251.501.752.00y5.105.796.537.458.46用最小二乘法求拟合曲线。解：对于曲线，两边取对数得：令，则可得到：把*，y的数据转换为t，z的数据取3位有效数字：*1.001.251.501.752.00z=lny1.631.761.882.012.14对于，其法方程组为：；其中：，数据代入后得法方程组为；解得。所以拟合曲线为。四、 16分用任意一种方法求以下积分，要求5位有效数字。解：用龙贝格方法求解题目中的积分：；。五、 16分同05-06学年秋第五题六、 16分求参数a，b，c使以下求解常微分方程初值问题的数值方法的局部截断误差到达。解：令则1对作泰勒展开得：代入到1式中可以得到：由于在的条件下假设要使，必须满足：，解得。. z.