西安交通大学计算方法b上机报告(完整题目和程序).docVIP

交通大学 计算方法B 上机编程报告 学号: XXX 姓名: XXX 专业: 工程热物理 班级: 硕 XXXX date 2015/12/10 Word资料 注：本上机报告使用的程序语言均为 Matlab语言，为本人独立完成! 1.计算以下和式:，要求:_1 1.计算以下和式: ，要求: n 016n 8n 1 8n 4 8n 5 8n 6 若保留11个有效数字，给出计算结果，并评价计算的算法; 若要保留30个有效数字，则又将如何进行计算。 解题思想 (1)先根据精度要求估计所需累加的项数 n,使用后验误差估计方法，条 件为:1 4 2 1n o16 件为: 1 4 2 1 n o16n 8n 1 8n 4 8n 5 1 8n 6 10-m (m为有效数 为了避免有效数字字位数) 为了避免有效数字 (2)在该题中S的和式存在两个相近的数相减的问题, 损失，在计算中改变了运算顺序，分别将正数和负数分别相加，然后 再将其和相加。 为了避免大数吃小数的问题，本题先计算出保留目标有效数字所需要 的迭代次数，然后采用倒序相加的方法提高计算精度。 算法实现的结构 an 1 4 16n 8n 1 2. for n=0,1,2,- …,i If an 10 m end 1. S1=s2=0 ； 2 1_ 1 (8n 4 8n 5 8n 6〉 for n=i,i-1,i-2,…,0a仁 4心6八 n