常州大学数值分析07-08试卷B及参考答案.docVIP

江苏工业学院2007~2008学年第 2 学期硕士生考试试题解答 一、（10分）举例说明如何在数值计算过程中防止相近数相减及避免“大数吃小数”。 答：1）防止相近数相减举例：当充分大时，即时，计算 会出现相近数相减， 可以用下述数学上等价的表达式 来计算，以避免相近数相减。 ………… 5 分 避免“大数吃小数”举例： 设为区间[0, 0.5]上的随机数，在字长为5的计算机上计算 时，如果采用上述给定的顺序计算，则会出现大数吃小数的现象；要避免大数吃小数，这里可以采用表达式：规定的顺序来计算即可。 ………… 5 分 注：学生的举例只要符合要求均可以算对。 二、(15分) （1）叙述Newton插值方法的方法思想； (2) 设, 试求的三次Newton插值多项式； (3) 利用上述插值公式近似计算. 解： (1) 牛顿插值方法是通过构造如下形式的多项式 其中通过Newton差商公式得到，且仅与有关，由此可以保证在增加节点时， 原先的计算量能够被充分利用。 ………… 6 分 (2) 根据列表函数可得差商表如下： 0 0 0 0 1.0000 1.0000 0 0 3.0000 2.0000 0.5000 0 5.0000 2.0000 0 -0.1667 的三次Newton插值多项式为: ………… 6 分 (3) ………… 3 分 三、（15分） （1）简要叙述求非线性方程根的迭代法的方法思想。 （2）选用适当的迭代方法求方程在附近的一个根, 精度 为。 解：（1）求非线性方程根的迭代法的方法思想： 将方程改写成 由给定的初始近似解，给出如下迭代公式 ，， 如果上述迭代序列收敛，即 则为方程的根。 ………… 7 分 （2）将方程改写成 由此可得到相应的迭代公式 。 由于上述迭代公式的迭代函数在处的导数的绝对值小于1，因此上述迭代公式在附近具有局部收敛性。 ………… 4 分 通过上述迭代可以得到上述方程根的数值近似计算结果如下： 0 2.5000 1 2.5198 2 2.5313 3 2.5379 4 2.5417 5 2.5439 ………… 4 分 四、（10分）叙述确定二次函数拟合下述列表函数的步骤 …… …… 解：第一步，构造函数 ………… 5 分 第二步，令 上述方程是关于所求参数的线性方程组。 第三步，解上述方程组可得所求参数，由此可得到用二次函数拟合上述列表函数的最小二乘解。 ………… 5 分 五、（15分）叙述复化Simpson积分公式计算的方法思想，并用复化Simpson公式计算积分 , 其中。 解：方法思想：由截断误差可知, 当区间长度b－a较大时, 梯形求积公式的误差较大. 为此，利用积分关于区间具有可加性, 将[a，b]区间上的积分, 分成若干小区间上的积分, 以此来减少积分区间长度引起的误差. 这就引入了复合求积公式. 具体如下： 设分点将区间[a，b]分成n等分，则 将每个小区间上的积分都用Simpson公式给出，则得计算定积分的复化Simpson公式如下： ………… 7 分 利用上述公式，取时，可得如下计算结果： s = 0.7023 ………… 8 分 六、（15分）说明求解微分方程初值问题的Runge-kutta方法的方法思想，并选用适当的方法解微分方程初值问题 的数值解, 取步长. 解：对于方程 龙格-库塔（Runge－Kutta）方法利用微分中值: 其中. 由于未知, 从而未知, 故我们不能直接应用上述公式解微分方程. 在实际计算时, 我们只好寻找的适当近似, 并由此得到相应的算法. 例如: 若用近似就可以得到Euler公式; 更一般地, 我们可以考虑用区间上某几个点处的值的加权平均作为的近似, 并满足一定误差要求，由此给出相应的公式: 这就是Runge-Kutta(龙格-库塔)算法的思想. ………… 7 分 上述方程利用Runge-lutta方法计算得到的数值解为： x=0.2时, k1= 0.0000 k2= 0.2002 k3= 0.1982 k4= 0.3935 y( 1)= 1.0397 x=0.4时, k1= 0.3934 k2= 0.5811 k3= 0.5756 k4= 0.