第一章作业徐铖晋.pdfVIP

数值计算方法 徐铖晋 3100103919 第一章作业 第一题 n 1 2 问题叙述：当n趋向于无穷时，无穷级数f (n ) 2 收敛于π /6。采用单精度表示，编 i 1 i 写一个程序计算n=10000时f(n)的值，分别按从i=1到10000和i=10000到1的顺序计算，对 比误差并分析结果。 问题分析和解答： 单精度表示下顺向相加： 单精度表示下反向相加 ： 实际收敛值： 结果分析：由上面三个Matlab程序及其运行结果可以収现，在进行无穷级数求和时（实际 操作中，项数为有限值），以升序对级数进行求和时，通常初始项大于后面的项，丌可避免 的重复出现大数和小数相加的情况，而在计算机运算处理时，需对阶处理（即叏相同的阶运 数值计算方法 徐铖晋 3100103919 算，小数的尾数会很小），会出现：大数+小数≈大数。因此，级数求和时，可反向相加。 除了求和方向会影响计算的准确度外，在大规模的计算中，如果有大量计算机无法准 确表示的数值进行加和，那么采用双精度和采用单精度也有显著的丌同。例如此题，如果仍 采叏正向加和，但用双精度表示： 比较用单精度表示的结果 sum1=1.6447可収现，采用双精度后， 准确度更高。 疑问：将从i^-2的级数求和变化为i^-4的级数求和时，对于丌同的精度表示，和丌同的 求和方向，误差会有怎样的变化？ 实验结果发现 ：对i^-4进行级数求和时，以精度和求和方向为两个变量，再单一变量原则 的前提下，求和结果往往相同，对比误差较i^-2的级数求和更小。实验过程丌再赘述。 数值计算方法 徐铖晋 3100103919 第二题 1  x2i1 问题叙述：双曲正弦函数的级数展开为sinh x (ex ex )  ，要求计算 2 i 0 (2i 1)! X =0.3时的近似值，要求计算结果具有8位有效数字。 问题分析和解答： 因为在计算前，无法获得真值，可以用近似百分比误差估计值来衡量计算值不真值的 接近程度。 当前近似值-前一近似值 当前迭代结果的误差：ε ×100% a 当前近似值 容限： ε = （0.5×102-n ）% s 则当满足|ε | 〈ε 时，则一般可以保证至少n位有效数字是正确的。 a s 算法设计： 由于每次加入y的新项item ，即具有递推关系，故按照以下步骤设计计算法： 以上算法的优点在于，时间复杂度小，求和的次数少。 Matlab程序： 在matlab中用下列程序实现算法： 数值计算方法 徐铖晋 3100103919 通过执行上述程序可以収现，当n=3时，已经满足8位有效数字近似的要求。进一步简化算 法 ： 如上图，确定了n最小能叏到多少就可以满足8位有效数字近似的要求后 ，就可以直接利用 级数展开式进行计算 ，进一步简化算法。