数值分析模拟试题.docVIP

2011数值分析模拟试题 （教师：南华大学王礼广） 数值分析模拟试题一 一、(15分)已知函数且的离散值如下表： 0 1 2 0．32 0．34 0．36 0．314567 0．333487 0．352274 试用拉格朗日二次插值公式求（保留小数点后6位）并估计截断误差。 （提示：拉格朗日n次插值公式：， 误差公式：） 解： 方法一：用拉格朗日插值公式计算： 计算完毕。 方法二：用牛顿插值公式计算： 计算完毕。 （15分）求在[-1，1]上次数不超过3的最佳一致逼近多项式。 （提示：第一类切比雪夫多项式）。 解：设最佳一致逼近多项式为，则有 ，（10分） ，计算完毕。 三、（15分）已知实验数据如下： Xi 19 25 31 38 44 Yi 19．0 32．3 49．0 73．3 97．8 求形如的经验公式，并估计均方误差。 （拟合公式提示： ，， ，） 令 ，， 解 ，， ， ， ， ，，计算完毕。 四、（15分）用四点高斯-勒让德公式求定积分（要求小数点后4位）。 （提示：四点公式系数：） 解： 令，得 用四点高斯-勒让德公式求得 ，计算完毕。 （注：精确值为）。 （15分）对方程构造一种高斯-塞德尔迭代公式， 说明其收敛的理由； 以为初值，迭代到。 (提示：高斯-塞德尔迭代： ，，) 解： ,,, ,,, D=diag(diag(A)),L=-tril(A,-1),U=-triu(A,1),B=inv(D-L)*U, f=inv(D-L)*b. , λ1=0，λ2=-0.1125 + 0.1654i, λ3=-0.1125 - 0.1654i， 谱半径，所以迭代 收敛, （5分）。取 ,(k=0,1,2,3,4,5)得下表 k Xk ||Xk-Xk-1||∞ ||Ax-b||2 0 (0 0 0) 1.65 13.0384 1 (1.6 1.65 1.075) 2.5373 4.8747 2 (0.7250 0.8937 1.0231) 1.1577 1.5678 3 (1.0379 0.9979 0.9918) 0.3312 0.1877 4 (1.0025 1.0047 1.0009) 0.0372 0.0292 5 (0.9979 0.9990 1.0001) 0.0073 0.0123 6 (1.0004 1.0000 0.9999) 0.0026 0.0019 实际取k=5即可. 计算完毕。 注：精确解为 ， （25分）对于初值问题 证明用梯形公式是绝对稳定的； 取步长，用（1）的格式计算其数值解（小数点后保留5位）。 （提示：梯形公式：） 解： ，解得 对任意的，，故有，误差是缩小的，所以格式绝对稳定。（5分）取，则， ，由计算得： ， ， ， ， ， 计算完毕。 数值分析模拟试题二 一、(20分)按照Hermit插值的思想，用构造基函数的方法求一个次数不超过2的插值多项式使,并计算的值 (计算保留4位有效数字). （提示：Hermit插值公式：, 其中, ,） 解:设基函数为,则有 , 由得 由得, 由得, 所以 ,，计算完毕。 二、(20分)求在上的二次最佳平方逼近多项式及其平方误差(计算保留4位有效数字).(提示:对于函数在上的最佳平方逼近的法方程为,其中记号,平方误差为) 解：,,,,;,,解得,所以 平方误差=3.0228e-005;计算完毕。 三、(10分)选取适当的步长，用计算定积分的复化梯形公式 (余项)或二分递推梯形公式(,事后误差估计)来计算定积分，使其积分的误差不超过0.01。 解: 由于, 设,得,取有, ,，计算完毕。 四、(10分)对于方程，因右向量的误差成为新的方程时，试用其误差估计公式，关于矩阵和向量的最大模，求出解的相对误差（计算保留4位有效数字）。 解：因为. , ，计算完毕。 五、(20分)对于非线性方程组，试用牛顿法构造一个求解迭代格式，并对于初值，计算(计算保留小数点后4位,提示:方程组。 解：牛顿法迭代公式为:, 其中为的雅可比式) 令,则 ,, 所以 , 计算完毕。 六、(20分)对于一阶初值问题 请用向前欧拉公式给出其求解格式； 请讨论自变量步长的稳定区间； 对于计算到的数值解(计算保留小数点后4位)。 （提示：一阶问题的向前欧拉公式为：） 解：, (1), (2)由(1),,要使格式稳定,必须,即,解得为步长的稳定区间. (3), 已知初值,由(1)及得 计算完毕。 数值分析模拟试题三