计算方法与实习上机报告

计算方法与实习

——上机报告

学院：电子工程学院

2015.1.4

学号：13020130016 姓名： 刘 波

计算方法与实习上机报告

习题一：

1 舍入误差及稳定性 一、实验目的

（1）通过上机编程，复习巩固以前所学程序设计语言及上机操作指令； （2）通过上机计算，了解舍入误差所引起的数值不稳定性 二、实验内容

1、用两种不同的顺序计算?n?2，分析其误差的变化

n?1100002、已知连分数f?b0?a1，利用下面的算法计算f：

b1?a2/?b2?a3/(...?an/bn)?ai?1 (i?n?1,n?2,...,0) f?d0 di?1写一程序，读入n,b0,b1,...,bn,a1,...,an,计算并打印f dn?bn,di?bi?3、给出一个有效的算法和一个无效的算法计算积分

xnyn??dx (n?0,1,...,10)

04x?1N1?311?14、设SN??2，已知其精确值为????

2?2NN?1?j?2j?1（1）编制按从大到小的顺序计算SN的程序

1（2）编制按从小到大的顺序计算SN的程序

（3）按两种顺序分别计算S1000,S10000,S30000,并指出有效位数 三、实验步骤、程序设计、实验结果及分析

1、用两种不同的顺序计算?n，分析其误差的变化

n?110000?2（1）实验步骤：

分别从1~10000和从10000~1两种顺序进行计算，应包含的头文件有stdio.h和math.h （2）程序设计： a.顺序计算

#include #include void main() {

double sum=0; int n=1; while(1) {

sum=sum+(1/pow(n,2)); if(n00==0)

printf(\ if(n>=10000) break;

n++; }

printf(\ }

b.逆序计算

#include #include void main() {

double sum=0; int n=10000;

while(n!=0) {

sum=sum+(1/pow(n,2)); if(n 0==0)

printf(\ if(n<1) break;

n--; }

printf(\ }

（3）实验结果及分析： 程序运行结果： a.顺序计算