工程数学模型及数值方法例题 修改页码.doc

《工程数学模型及数值方法》解答 一、 计算的数值稳定性 1. 对于一个采用迭代计算的问题 假设初值，如果采用新的近似初值进行计算，分析这种计算方法是否稳定？并说明如果按照新的初值计算到，产生的误差有多大？ 解：设以初值递推得到的近似计算值为，设以初值递推得到的近似计算值为，则。从而得到递推计算的第步的误差是初始误差的倍，，尽管初始误差不大，但是随着的增大， 就可以非常大。从而这种算法是一种不稳定的算法。 如果按照新的初值计算到，产生的误差为： 2. 在实际计算中，人们除了分析算法是否稳定之外，还要求在计算中尽量避免误差危害，防止有效数字的损失.通过研究和实际计算，前人总结出了一些行之有效的基本的避免误差危害的原则，阐述有哪些基本原则？ 答：避免相近的数相减；避免量级相差太大的两数相除；避免大数和小数相加减；简化计算步骤。 思考： 1. 对于一个采用迭代计算的问题 假设初值，如果采用新的近似初值进行计算，分析这种计算方法是否稳定？并说明如果按照新的初值计算到，产生的误差有多大？ 2. 对于积分，容易导出计算公式；若给出的一个近似值，则可以逐步递推计算出的近似值。分析这种计算方法的稳定性？（书30页例2.3.1） 二、 插值和曲线拟合问题 假设表示时间，表示不同时间点的某支股票价格.现在连续观测了2012个时间点股票价格数据，下表仅仅给出了部分数据. 1.0 2.0 3.0 4.0 1.9 2.7 4.8 5.3 为了分析股票价格随时间变化的趋势，可以考虑采用插值或者曲线拟合的方法来进行分析. 利用上述已知数据做Lagrange插值，请写出插值函数的表达式； 利用上述数据进行最小二乘曲线拟合，试求出曲线拟合的表达式； 如果上述数据表中把2012个时间点的数据全部给出，是否可以利用2012个数据进进行插值计算？简要说明原因； 如果上述数据表中把2012个时间点的数据全部给出，是否可以利用2012个数据进行曲线拟合？简要说明原因； 对于给定的数据做插值和进行曲线拟合，两种方法有什么差异？ 对于表中的实验数据，利用数值微分求函数在节点处的一阶导数和二阶导数。（书89-93页） 解 1. 取， 则Lagrange插值函数的表达式为： 2. 在t-y坐标系内画出其散点图，这些点落在一直线附近，则选择线性函数作拟合曲线，即，这里。由于数据表中没有给出权重，所以我们取，则 于是有 ，解得 故所求拟合曲线为 3. 不能。因为对于任意的插值节点，当时，插值多项式不一定收敛到，即对于节点个数多的时候容易出现数值振荡（龙格现象）会发生；如果采用样条插值，计算量又太大 4. 可以。因为这些数据组成的点在坐标系中可能落在一直线或曲线附近，所以可以选择合适的直线或曲线进行拟合。 5. 对于给定的数据做插值时，所求曲线必须通过给定的个点；但对于给定的数据进行曲线拟合时，所求曲线可以不通过给定的个点，而只要求在给定点 上的误差（也称为残差） 按某种标准最小。 6. 用两点公式，取，则 若取，则 用三点公式，取，则 ， 三、 采用数值积分方法计算积分 1. 复化梯形公式（）； 2. 两点高斯-勒让德（Gauss-Legendre）公式； 3. 数值积分的基本思想是什么？简要说明可以考虑什么方法计算广义积分？ 解 1. 解 将区间分为2等分，先计算出各个节点处的函数值. 令，则有，， 则? 2. 先做变量代换，这时原积分变成了， 则用两点高斯勒让德公式可得 原式 3. 数值积分的基本思想是依据积分中值定理，将积分求值问题归结为平均高度 的函数值的计算。即， 在区间[a,b]上适当选取某些节点，，然后用加权平均得到平均高度的近似值，这样构造出的求积公式（也称为求积格式）具有形式 ，式中 称为求积节点，?称为求积系数，亦称伴随节点的权。权仅仅与节点的选取有关，而不依赖于被积函数 f (x)的具体形式。 计算广义积分的方法有高斯—拉盖尔求积公式方法（书87页）。 思考：如何用复化梯形公式（）和 两点高斯-勒让德（Gauss-Legendre）公式计算积分？ 四、 方程组的求解问题 1. 采用直接法求一个线性代数方程组的近似解时，有一种重要方法称为矩阵的LU分解，说明这种分解的含义并解释这种分解在方程组求解中起什么作用？ 2. 对于方程组 写出求解方程组的Jacobi迭代和Gauss-Seidel迭代法的分量形式的迭代公式； 判定上述Jacobi迭代法和Gauss-Seidel迭代法的收敛性； 如果取初值，利用（2）中收敛的迭代法迭代计算两次； 对于一个非线性方程组的求解，该如何求解，请解释它们的大致含义？ 解 1. 矩阵的LU 分解定理 设为n 阶矩阵， 如果的顺序主子式，(