数值分析上机实习题及答案.docxVIP

PAGE \* MERGEFORMAT1 [数值分析] 2017-2018 第二学期上机实习题  1：编程计算， 其中c= 4.4942?10307， 给出并观察计算结果，若有问题，分析之。 解：matlab编程如下： 图一：编程 图二：运行结果 Matlab中，format long g 对双精度，显示15位定点或浮点格式，由上图可知，当输入较小的n值（n分别取10,100,1000,100000,100000000）的时候，结果后面的指数中总是含有e-308，这和题目中的C值很相似，我认为是由于分母中的C值相对于n值过大，出现了“大数吃小数”的现象，这是不符合算法原则的。 2：利用牛顿法求方程于区间的根，考虑不同初值下牛顿法的收敛情况。 解：牛顿法公式为： x 利用matlab编程 当输入初值=3的时候并不能收敛。 输入2.8依旧不能收敛。 考虑利用plot函数绘制函数f=x-lnx-2的图像，发现其在区间[2,4]内并无根。 3：给出 在xk=0+0.25k处的值yk, k=0,1,2,3,4. 请给出由节点 xk 确定的三次样条插值函数S(x), 使其满足条件：S/(0)=0, S/(1)=-0.074, 分析逼近效果如何？ 解：根据李庆扬《数值分析》第五版Page43的公式可以计算得到三次样条插值。