坐標(biāo)輪換法c程序設(shè)計坐標(biāo)輪換法c程序設(shè)計坐標(biāo)輪換法c程序設(shè)計xxx公司坐標(biāo)輪換法c程序設(shè)計文件編號：文件日期：修訂次數(shù)：第1.0次更改批準(zhǔn)審核制定方案設(shè)計，管理制度////用坐標(biāo)輪換法求得函數(shù)的最優(yōu)解//#include""#include""#include""#include""#defineedoublefuntion(doublex[],inta){ doublef1; f1=pow(pow(x[0],2)+x[1]-11,2)+pow(x[0]+pow(x[1],2)-7,2); returnf1;}intmain(intargc,char*argv[]){ ints1[2],s2[2],i,j; doublex[100][2],xl[2],h1,t,m,f1,f2,fl,p,f3,f4; h1=,t=,p=2; s1[0]=1,s1[1]=0; s2[0]=0,s2[1]=1; x[0][0]=,x[0][1]=;f1=funtion(x[0],2); x[1][0]=x[0][0]+h1*s1[0]; x[1][1]=x[0][1]+h1*s1[1]; f2=funtion(x[1],2);for(j=1;;j++){if(f1>f2) { for(i=2;;i++) { h1=p*h1; f3=funtion(x[i-1],2); x[i][0]=x[i-1][0]+h1*s1[0]; x[i][1]=x[i-1][1]+h1*s1[1]; f4=funtion(x[i],2); if(f3<f4)break; } x[0][0]=x[i-2][0]; x[0][1]=x[i-2][1]; x[1][0]=x[i-3][0]; x[1][1]=x[i-3][1]; }else { h1=-h1; for(i=1;;i++) {f3=funtion(x[i-1],2);x[i][0]=x[i-1][0]+h1*s1[0]; x[i][1]=x[i-1][1]+h1*s1[1];f4=funtion(x[i],2); h1=p*h1; if(f3<f4)break; }x[0][0]=x[i-1][0]; x[0][1]=x[i-1][1];x[1][0]=x[i-2][0]; x[1][1]=x[i-2][1]; } t=-t; f1=funtion(x[0],2); x[1][0]=x[0][0]+t*s2[0]; x[1][1]=x[0][1]+t*s2[1]; f2=funtion(x[1],2);if(f1>f2) { for(i=2;;i++) { t=p*t; f3=funtion(x[i-1],2); x[i][0]=x[i-1][0]+t*s2[0]; x[i][1]=x[i-1][1]+t*s2[1]; f4=funtion(x[i+1],2); if(f3<f4)break; }x[0][0]=x[i-1][0]; x[0][1]=x[i-1][1];x[1][0]=x[i-2][0]; x[1][1]=x[i-2][1]; } else { t=-t; for(i=1;;i++) {f3=funtion(x[i-1],2);x[i][0]=x[i-1][0]+t*s2[0]; x[i][1]=x[i-1][1]+t*s2[1];f4=funtion(x[i],2); t=p*t; if(f3<f4)break; }x[0][0]=x[i-1][0]; x[0][1]=x[i-1][1];x[1][0]=x[i-2][0]; x[1][1]=x[i-2][1]; } m=pow(pow((x[0][0]-x[1][0]),2)+pow((x[0][1]-x[1][1]),2),; if(m<e) { xl[0]=x[0][0]; xl[1]=x[0][1]; fl=funtion(xl,2); printf("最優(yōu)解為:xl=(%f,%f)\n最優(yōu)解的函數(shù)值為:fl=%f\n",xl[0],xl[1],fl); exit(0); } else {h1=,t=,p=*p;x[0][0]=,x[0][1]=;f1=funtio