数值计算方法上机实验.docVIP

PAGE 1 PAGE 16 重 庆 交 通 大 学 学 生 实 验 报 告 实验课程名称 数值计算方法 开课实验室 理学院实验室 学 院 土木建筑学院 专 业 班 级 09 级 测绘工程 专业 2 班 学 生 姓 名 吴必坤 学 号 开 课 时 间 2010 至 2011 学年第 二 学期 评分细则 评 分 考勤（10分） 排版格式（10分） 实验程序代码的正确性、完整性（40分） 实验结果及分析（30分） 实验创新（10分） 总 成 绩 教师签名 张文忠 目 录 TOC \o 1-3 \h \z \u 实验一 误差——定积分 2 实验二 非线性方程求根 4 （一） 牛顿迭代法 4 （二） 牛顿下山法 5 实验三 解线性代数方程组——列主元消元法 7 实验四 解线性代数方程组——追赶法 9 实验五 插值与拟合——三次样条差值 11 实验一 误差——定积分 实验内容： 用两种方法求解以下积分： , 初值 2. In=1-In-1*n 2 实验过程： 1利用C语言编程求解，程序如下： #includestdio.h #includemath.h void main() {double x1;int i,n; x1=100; for(i=50;i=0;i--) {x1=(1-x1)/i; printf(i=%d,xi=%lf\n,i,x1);} } 1.2.命令窗口中程序结果如下： 2.1利用C语言编程求解，程序如下： #includestdio.h #includemath.h void main() {double x1=1-exp(-1); int i,n,x0; scanf(%d,n); x0=1-exp(-1); printf(i=0,x0= 0.6320205671\n); for(i=1;i=n;i++) {x1=1-i*x1; printf(i=%d,xi=%12.9lf\n,i,x1);} } 1.2.命令窗口中程序结果如下： 实验结果分析： 从程序运行的结果可以看出，当输入n=14时，即运算到n=14时，结果与书上的结果并不一致，然而从n=1到n=10的运算结果是完全一致的，从n=11以后的运算结果可能是由于不同时期的计算机的处理数据的能力的不同或精度不同导致的。 实验二 非线性方程求根 实验内容： 用牛顿迭代法和牛顿下山法求解方程f（x）=，取近似值。 1．牛顿迭代法： ，直到。 2．牛顿下山法： ，（λ=1至合适），且直到| 实验过程： 1．用牛顿迭代法求解 1.1利用C语言编程求解，程序如下： #includestdio.h #includemath.h void main() {double x=-0.99,y; int n=0; printf(%d ,%lf\n,n,x); do {y=x; x=x-(x*x*x/3-x)/(x*x-1); n++; printf(n= %d ,x= %lf\n,n,x); }while(fabs(y-x)1e-5); } 1.2.命令窗口结果截屏如下： 2．用牛顿下山法求解 2.1利用C语言编程求解，程序如下： #includestdio.h #includemath.h double f(double x) { double t; t=x*x*x/3-x; return fabs(t); } void main() {double x=-0.99,y; int k=0; printf(%d ,%lf\n,k,x); do {y=x; int n=0; do{ x=y; x=x-pow(0.5,n)*(x*x*x/3-x)/(x*x-1); n++; }while(f(x)=f(y)); k++; printf(k=%d, n= %d, x= %lf\n,k,n,x); }while(fabs(x-y)1e-5); } 2.2.命令窗口结果截屏如下： 实验结果分析： 从运算截图中我们可以知道，一方面运用牛顿迭代法和牛顿下山迭代法的程序运算的结果是一致的，若出现小小的差异的话，可能是由于设置的变量类型不同导致，另一方面运用牛顿下山迭代法比运用一