计算方法上机作业集合.docx

第一次第二次上机作业上机作业：1.在Matlab上执行： 5.1-5-0.1和 1.5-1-0.5给出执行结果，并简要分析一下产生现象的原因。解：执行结果如下：在Matlab中，小数值很难用二进制进行描述。由于计算精度的影响，相近两数相减会出现误差。2.（课本181页第一题）解：（1）n=0时，积分得=ln6-ln5，编写如下图代码从以上代码显示的结果可以看出，的近似值为0.012712966517465（2）=dx,可得dxdxdx,得,则有,取,以此逆序估算。代码段及结果如下图所示结果是从逆序输出至，所以得到的近似值为0.088392216030227。（3）从估计的过程更为可靠。首先根据积分得表达式是可知，被积函数随着n的增大，其所围面积应当是逐步减小的，即积分值应是随着n的递增二单调减小的，（1）中输出的值不满足这一条件，（2）满足。设表示的近似值，-=(根据递推公式可以导出此式),可以看出，随着n的增大，误差也在增大，所以顺序估计时，算法不稳定性逐渐增大，逆序估计情况则刚好相反，误差不断减小，算法逐渐趋于稳定。2.（课本181页第二题）（1）上机代码如图所示求得近似根为0.09058（2）上机代码如图所示得近似根为0.09064；（3）牛顿法上机代码如下计算所得近似解为0.09091第三次上机作业上机作业181页第四题线性方程组为=顺序消元法A=[1.1348,3.832