西安交大计算方法上机报告.pdfVIP

计算方法 上机实习题目报告 班级：材料s3076 班 姓名：丁明帅 学号：3113305029 100000 1 6 1：计算 ,要求误差小于 ,给出实现算法。 S  2 10 k 1 k 【实现思路】 -6 设当k 值为i 时S-S 10 ,则只需要 i 100000 1 S  Si  2 i1 k 100000 100000 1  1 1          i1 k k  1 i1 k  1 k  1 1 6   10 i 100000 1 即可得到i 5 6 90909.0909 10  10 因此，只需取i=90910 即满足题目的近似要求。 【算法依据】 本题比较简单，主要通过for 循环语句对所需要计算的90910 个数带入公式 进行计算，然后运行程序即可得出S 的近似解。 【算法步骤】 1 ：给定初始值条件s, n 。 2 ：用for 循环对1/(s^2)从1 到90910 进行求和，并将每次的值加到初值为零的 sum 上。 3 ：输出sum 的结果。 【matlab 运行程序】 函数question1.m function question1 s=1;%起始值 n=90910;%项数 format long, sum=0;%累加的和 for i=1:n 3 / 46 sum = sum+1/(s^2); s=s+1; end display([原式求和结果为:,num2str(sum,16)]); 【运行结果】 question1 原式求和结果为:1.644923067018747 【结果分析】 本题通过理论近似计算算出当求和项大于90909 时，所计算出的S ’近似值即与精确值 6 S 之间的误差小于10 ，因而无需将求和计算到第 100000 项，这样就节省了计算机计算的 时间，同时也达到了计算精度。计算出的结果为S ’=1.644923067018747，由于误差项在10-7 ， 故取S ’=1.6449231。即S  S =1.6449231 【问题总结】 本题在刚开始分析的时候由于想到从1 到100000 计算S 的精确值的计算量对于电脑来 说并不是很大，因而曾先计算出S 的精确值，然后用while 语句对每一步计算过后值再进行 判断，如果近似值已经与精确值的差小于给定误差，则停止计算。那么计算出的近似值即是 所要求的解。但后来想到这样并没有减少计算量，反而比求解精确值的计算量还要大，然后 就将这样的求解模式否定了。另一个问题就是在求解过程中发现输出的结果的小数个数太少， 无法达到相应的精度，后来通过调整变量类型解决了这个问题。 2：编写实现对N 阶非奇矩