实验一非线方程组求解实验报告.docVIP

计算方法实验报告1．【实验题目2．【实验目的3．【实验内容 ，针对该方程，可以构造多种迭 代法，如： 取初始值 ，取 ，分别用以上迭代格式作实验，记录各算法的迭代过程 4. 【实验 ，分别用 方案1 用二分法求解； 方案2 用牛顿法求解； 5. 【算法描述】 二分法算法步骤 计算有根区间的端点a,b及预先给定的精度e。 计算中点(a+b)/2。 若f(x)f(a)0,则xb,转向4)；否则，xa，转向4). 若b-ae，则输出满足精度的根x,结束；否则转向2)。 牛顿法迭代法的计算步骤 给出初始近根及精度e。 计算。 若|-|e，则转向4)；否则转向2)。 输出满足精度的根，结束。 【源程序（带注释）】 二分法 #includestdio.h #includemath.h #includeconio.h #includewindows.h float f(float x) { float a; a=x*x*x-x-1; return a; } /*求函数值，如果求其它函数，只需改成其它函数即可*/ main() { float a,b,e,x; /* a,b 分别表示有根区间的左、右端点, e 是精度要求,x 区间中点值*/ system(CLS);//清屏 printf(对代数方程x^3-x-1=0，分别用\n方案1 用二分法求解\n); printf( \n please input data a =); scanf(%f,a); printf( \n please input data b=); scanf(%f,b); if(f(a)*f(b)0) { while(f(x)!=0) { x=(a+b)/2; if(f(x)*f(a)0) { b=x; if(fabs(b-a)0.000001) break; else continue; } else { a=x; if(fabs(b-a)0.000001)break; else continue; } } printf(\n); x=(b+a)/2; printf(the root of f(x)=0 is x=%f\n,x); } else printf(\ not root! afresh input\n); /*表示 [a,b] 区间无根，重新选择有根区间*/ getch(); return(x); } 牛顿法 #includestdio.h #includemath.h #includeconio.h #includewindows.h #define maxrept 1000 /*最大迭代次数*/ float f(float x) { float a; a=x*x*x-x-1; return a; /*函数f(x) */ } float df(float x) { return(1+exp(-x)); /* 函数f(x)的导数 ) (x f ′ */ } float iterate(float x) { float x1; x1=x-f(x)/df(x); /* 牛顿迭代函数 iterate(x)=x-f(x) / ) (x f ′ */ return(x1); } main() { float x0,x1,d; int k=0; // clrscr(); system(CLS); printf(对代数方程x^3-x-1=0，分别用\n方案2 用牛顿法求解\n); printf(\n please input x0=); /* 输入迭代初值 x0 */ scanf(%f,x0); printf(\n k xk\n); printf(\ %d %f\n,k,x0); do { k++; x1=iterate(x0); printf( %d %f\n,k,x1); d=fabs(x1-x0); x0=x1; } while((d=0.000001)(kmaxrept)); if(kmaxrept) printf(the root of f(x)=0 is x=%f, k=%d\n,