坐标轮换法(麦华锋 东莞理工).docVIP

#includemath.h #includestdio.h #includestdlib.h int main(void) { double f(double x1,double x2); double fun(double x,double a0,double h0,double *a,double *b,int N); double fun2(double x,double a,double b,int N); double h0=1,p,a,b,y,x1[2]={0},x2[2]={0}; int i=0; printf(input:x2[1],x2[2]:\n); scanf(%f%f,x2[1],x1[2]); do { x1[1]=x2[1]; x1[2]=x2[2]; fun(x1[1],x1[2],h0,a,b,1); x2[2]=fun2(x1[1],a,b,1); fun(x2[2],x1[1],h0,a,b,2); x2[1] =fun2(x2[2],a,b,2); i=i+1; printf((%d) [%f,%f],i,x2[1],x2[2]); p=sqrt((x2[1]-x1[1])*(x2[1]-x1[1])+ (x2[2]-x1[2])*(x2[2]-x1[2])); printf(p=%f\n,p); }while(sqrt((x2[1]-x1[1])*(x2[1]-x1[1])+ (x2[2]-x1[2])*(x2[2]-x1[2]))=0.0001); y=f(x2[1],x2[2]); printf(x=[%f,%f],f(x)=%f\n\n\n,x2[1],x2[2],y); system(pause) ; return 0; } double f(double x1,double x2) { double y; y=2*x1*x1+3.1*x2*x2+2*x1*x2-2*x1+6*x2-1; return y; } void fun(double x,double a0,double h0,double *a,double *b,int N) { double y1,y2,y3,a2,a1,a3,h; a1=a0; h=h0; switch(N) { case 1:y1=f(x,a1); a2=a1+h; y2=f(x,a2);break; case 2:y1=f(a1,x); a2=a1+h; y2=f(a2,x);break; } if(y2y1) { h=-h;a3=a1;y3=y1; a1=a2;y1=y2; a2=a3;y2=y3; } a3=a2+h; switch(N) { case 1:y3=f(x,a3); break; case 2:y1=f(a3,x); break; } if(y3y2) { h=2*h; a1=a3;y1=y3; a3=a2;y3=y2; a2=a1;y2=y1; } else { if(h0) { *a=a3; *b=a1; } else { *a=a1; *b=a3; } printf(N=%d a=%f b=%f\n,N,a,b); } } double fun2(double x,double a,double b,int N) { double a1,a2,y1,y2,result; a1=a+0.382*(b-a); a2=a+0.618*(b-a); switch(N) { case 1:y1=f(x,a1); y2=f(x,a2);break; case 2:y1=f(a1,x); y2=f(a2,x);break; } while(fabs(b-a)0.0001) { if(y1=y2) { a=a1; a1=a2; y1=y2; a2=a+0.618*(b-a); switch(N) { case 1:y2=f(x,a2); break; case 2:y2=f(a2,x); break; } } else { b=a2; a2=a1; y2=y1; a1=a+0.382*(b-a); switch(N) { case 1:y1=f(x,a1); break; case 2:y1=f(a1,x); break; } } return 0; }