计算机实习任务报告.docVIP

计算机实习 C语言编程计算实习报告 学号： 姓名： 班级： 专业： 2011年12月23 一、实习目的 通过16学时的计算机实习（C语言编程计算），掌握一种技能。 二、实习课题 1、用牛顿迭代法求在附近的实根，取四位有效数字。 2、用高斯消去法求解下述线性方程组： 三、算法分析 ㈠、牛顿迭代法 设r是f(x) = 0的根，选取x0作为r初始近似值，过（x0,f(x0)）y = f(x)的切线L，L的方程为y = f(x0)+f(x0)(x-x0)，求出L与x轴交点的横坐标 x1 = x0-f(x0)/f(x0)，称x1为r的一次近似值。过点（x1,f(x1)）曲线y = f(x)的切线，并求该切线与x轴交点的横坐标 x2 = x1-f(x1)/f(x1)，称x2为r的二次近似值。重复以上过程，得r的近似值序列，其中x(n+1)=x(n)－f(x(n))/f(x(n))，称为r的n+1次近似值，上式称为牛顿迭代公式。 消元过程如下： 2.将每行的数字减去第一行的第一个数的某个倍数。 重复上述两个步骤直至增广矩阵变成一个行梯阵式，再用代入法就可以求得这个方程组的解。 在程序中用一维数组分别存储系数矩阵，常数项矩阵。 高斯消元法的算法复杂度是O(n3)；这就是说，如果系数矩阵的是n × n，那么高斯消元法所需要的计算量大约与n3成比例。对于较大的线性方程组，一般不会直接使用高斯消去法，而采用相应的迭代加速法。 四、程序流程图 牛顿迭代法： 高斯消元法： 五、程序运行结果截图 牛顿迭代法结果： 高斯消去法结果： 附：（源程序） 牛顿迭代法源程序Newton.cpp： #includestdio.h #includemath.h　 double func (double x) { //函数f(x)=x^3-3*x-1 return x*x*x-3*x-1; } double func1(double x) { //导函数f(x)=3*x^2-3 return 3*x*x-3.0; } int Newton(double *x,double precision,int macyc) { double x1,x0; int k; x0=*x; for (k=0;kmacyc;k++) { if (func1(x0)==0.0 ) { printf (迭代过程中导数为0!\n); return 0; } x1 = x0 - func(x0)/func1(x0); //进行牛顿迭代计算 if (fabs(x1-x0)precision || fabs(func(x1))precision) //达到结束条件 { *x = x1; //返回结果 return 1; } else //未达到结束条件 { x0 = x1; //准备下一次迭代 } } printf (迭代次数超过预期!\n); //迭代次数达到，仍没有达到精度 return 0; } int main() { double x,precision; int maxcyc; printf(输入初始迭代值x0:); scanf(%lf,x); printf(输入最大迭代次数:); scanf(%d,maxcyc); printf(迭代要求的精度:); scanf(%lf,precision); if(Newton(x,precision,maxcyc)==1) //若函数返回值为1 printf(该值附近的根为:%.3lf\n,x);; else //若函数返回值为0 printf(迭代失败!\n); return 0; } 高斯消去法源程序GuassElimination.cpp： #include stdio.h void GuassCommon(double *pA, double *pB, double *pX, int n) { double akk, mik, temp; int k, i, j; for (k = 0; k n - 1; k++) //化成行阶梯形 { akk = pA[k * n + k]; if (akk 0.00001 akk -0.00001) { printf(\nError: a[k, k] = 0\n); break; } fo