最速下降法和牛顿法求最小值点的算法及结果.docVIP

最速下降法计算：f=(x-1) +5(x-5) 程序代码如下： #includestdio.h #includemath.h double f1(double x,double y) { double r; r=(pow(x-1,2)+5*pow((y-5),2)); return r; } //最速下降法求最优解 void main() { double h=3,x0=3,x1,y0=6,y1,s,r0,r1; double e0=0.000001,e1=0.000001; int k=0; s=sqrt(pow(2*x0-2,2)+pow(10*y0-50,2)); printf(%d x1=%f x2=%f s=%f\n,k,x0,y0,s); while(se1) { x1=x0; y1=y0; r0=f1(x0,y0); h=3; //一维搜索，成功失败法 while(fabs(h)e0) { r1=f1((x1-h*2*(x1-1)),(y1-h*10*(y1-5))); if(r1r0) { x0=x1-h*2*(x1-1); y0=y1-h*10*(y1-5); r0=r1; h=h+h; } else if(fabs(h)e0) h=(-1)*h/4; } s=sqrt(pow(2*x0-2,2)+pow(10*y0-50,2)); k++; printf(%d x1=%f x2=%f s=%f\n,k,x0,y0,s); } printf(x1=%f x2=%f,x0,y0); } 初始值设为x1=3,x2=6时，运行结果如下图1-1： 图1-1 初始值设为x1=30000,x2=60000时，运行结果如下图1-2： 图1-2 牛顿法计算：f=(x-1) +5(x-5) 程序代码如下： #includestdio.h #includemath.h double f1(double x,double y) { return (pow(x-1,2)+5*pow(y-5,2)); } //牛顿法求最优解 void main() { double h=3,x0=3,x1,y0=6,y1,s,r0,r1; double e0=0.000001,e1=0.000001; int k=0; s=sqrt(pow(2*x0-2,2)+pow(10*y0-50,2)); printf(%d x=%f y=%f s=%f\n,k,x0,y0,s); while(se1) { x1=x0; y1=y0; r0=f1(x0,y0); h=3; //一维搜索 while(fabs(h)e0) { r1=f1((x1-h*(x1-1)),(y1-h*(y1-5))); if(r1r0) { x0=x1-h*(x1-1); y0=y1-h*(y1-5); r0=r1; h=h+h; } else if(fabs(h)e0) h=(-1)*h/4; } s=sqrt(pow(2*x0-2,2)+pow(10*y0-50,2)); k++; printf(%d x=%f y=%f s=%f\n,k,x0,y0,s); } printf(x=%f y=%f,x0,y0); } 初始值设为x1=3,x2=6时，运行结果如下图2-1： 图2-1 初始值设为x1=30000,x2=60000时，运行结果如下图2-2： 图2-2 最速下降法求：f=(x-1) +(x-1)+5(x-5)+5(x-5) 程序代码如下： #includestdio.h #includemath.h double f1(double x,double x1,double y,double y1) { double r; r=(pow(x-1,2)+pow(x1-1,2)+5*pow(y-5,2)+5*pow(y1-5,2)); return r; } //最速下降法求最优解 void main() { //x1,x2,x3,x4分别为多项式函数中四个分量的初始值 //x5,x6,x7,x8分别为对应的下一个值 double d=3,x1=3,x5,x2=4,x6,x3=8,x7,x4=6,x8,s,r0,r1; double e0=0.000001,e1=0.000001; int k=0; s=sqrt(pow(2*