第二次求根方法牛顿法等实验.pptxVIP

会计学1第二次求根方法牛顿法等实验 第一题：问题分析： 先利用作图工具绘出方程草图第1页/共20页 方程e^(-x*x)=cos(x)+1在定义域内的值在3附近。利用牛顿法求时F(x)=(cos(x)-exp(-x*x)+1);F’(x)=(2*x*( exp(-x*x))-sin(x)) ;根据题意代入x=1或x=0均得不合题意。代入x=2时x*=3.131081 符合题意算法设计 如下第2页/共20页 #include #include float F( float b) { float a; a=(cos(b)-exp(-b*b)+1); return(a);} float dao( float b) {double a; a=(2*b*( exp(-b*b))-sin(b)); return(a);} void main() { double m,b[30];int i; m=0; for(i=1;i30;i++) {b[0]=2; b[i]=(b[i-1]-F(b[i-1])/dao(b[i-1])); m=b[i]; printf( %f\n ,m); } } 第3页/共20页 测试结果第4页/共20页 方程e^(-x*x)=cos(x)+1在定义域内的值在3附近。利用割线法求时F(x)=(cos(x)-exp(-x*x)+1) F(x)’=(F(x)-F(x-1))/(x-(x-1)) 算法设计 如下第5页/共20页 #include #include float F( float b) { float a; a=(cos(b)-exp(-b*b)+1); return(a);} float dao( float b,float c) {double a; a=(F(b)-F(c))/(b-c); return(a);} void main() {double m,b[30];int i; m=0; for(i=2;i30;i++) {b[0]=0;b[1]=1; b[i]=(b[i-1]-F(b[i-1])/dao(b[i-1],b[i-2])); m=b[i]; printf( %f\n ,m); } } 第6页/共20页 测试结果 第7页/共20页 第二题：问题分析： 先利用作图工具绘出方程草图第8页/共20页 牛顿法分析：如函数的图像所示：根在1到2 之间。使用牛顿法，取x=1作为初值。收敛至x(k)与x(k+1)的差别满足精度要求即可求出它的根。Newton法是收敛速度较快的方法。算法设计如下第9页/共20页 #include #include float F( float b) { float a; a=(b*(b*(b+2)+10)-20); return(a);} float dao( float b) {double a; a=(b*(3*b+4)+10); return(a);} void main() {double m,b[30];int i; m=0; for(i=1;i20;i++) {b[0]=1; b[i]=(b[i-1]-F(b[i-1])/dao(b[i-1])); m=b[i]; printf( %f\n ,m); } } 第10页/共20页 测试结果第11页/共20页 割线法#include #include float F( float b) { float a; a=(b*(b*(b+2)+10)-20); return(a);} float dao( float b,float