非线性方程求解实验报告.docxVIP

课程实验报告课程名称计算方法班级实验日期姓名学号实验成绩实验名称实验三 非线性方程求解实验目的及要求1、掌握二分法与牛顿法的基本原理及应用2、用二分法计算方程在内的根的近似值3、用迭代计算方程在内的根的近似值。实验环境Microsoft Visual C++ 6.0实验内容二分法：假定在上连续，且在内仅有一实根取区间中点，若，则恰为其根，否则，根据是否成立，可判断出根所属的新的有根子区间或，为节省内存，仍称其为。运算重复进行，直到满足精度要求为止，即。式中为新的有根子区间的端点。迭代法：为初始猜测，则由递推关系产生逼近解的迭代序列，这个递推公式就是Newton法。当距较近时，很快收敛于。但当选择不当时，会导致发散。故我们事先规定迭代的最多次数。若超过这个次数，还不收敛，则停止迭代另选初值。算法描述及实验步骤1、实验1实验步骤：（1）程序设计（2）计算实例（3）撰写实验报告2、实验2实验步骤：（1）程序设计（2）计算实例（3）撰写实验报告调试过程及实验结果二分法求根实验结果如图：迭代法的求根实验结果如图：总结通过本次实验，使我掌握了使用二分法求非线性方程的根，加深了对二分法的了解。通过本次实验，使我掌握了使用迭代法求非线性方程的根的方法，能够使用使用迭代法正确的求出求非线性方程的根。附录二分法程序清单：#includestdio.h#includemath.h#includeconio.hdouble f(double x){double y = 0;y = sin(x) - x*x/2;return y;}void main(){double a,b,y1,y2,temp;int e1,e2=1,n,i=0;printf(请输入根的区间：);scanf(%lf,a);scanf(%lf,b);y1=f(a);y2=f(b);if((y1*y2)0){printf(要求小数点后的位数：);scanf(%d,e1);while(e10){e2=e2*10;e1--;}n=(int)(log((b-a)*e2)/log(2));printf(对方程sinx-x2/2=0开始进行二分：\n);while(n0){printf(第%d次二分结果:,i);temp=f((a+b)/2);if((f(a)*temp)0){b=(a+b)/2;printf(a=%f,b=%f\n,a,b);}if((f(b)*temp)0){a=(a+b)/2;printf(a=%f,b=%f\n,a,b);}i++;n--;}}printf(方程的解是：%6.5f\n,(a+b)/2);}迭代法的求根程序清单：#includestdio.h#includemath.h#includeconio.hdouble f(double x){double y = 0;y = x*x*x-x-1;return y;}double fdao(double x){double a=1.0/3;return pow(x+1,a);}void main(){double x1,x2,a,b,e2=1;int i = 0,e1;printf(请输入根的区间：);scanf(%lf,a);scanf(%lf,b);if((f(a)*f(b))0){printf(要求小数点后的位数：);scanf(%d,e1);while(e10){e2=e2/10;e1--;}e2=e2/2;printf(e2=%f,e2);x1=a;printf(第%d次迭代后结果是：%20.19f\n,i,x1);x2=f(x1);while(fabs(x2-x1)e2){i++;printf(第%d次迭代后结果是：%20.19f\n,i,x2);x1=x2;x2=fdao(x1);}printf(共进行%d次迭代过程\n,i);printf(结果是：%20.19f\n,x2);}elseprintf(无法进行迭代!);}