坐标轮换法C语言相关程序

坐标轮换法C语言相关程序坐标轮换法C语言相关程序坐标轮换法C语言相关程序坐标轮换法C语言相关程序编制仅供参考审核批准生效日期地址：电话：传真:邮编:坐标轮换法C语言相关程序#include<>#include<>floatf(floatx[]){floatz;z=4*(x[0]-5)*(x[0]-5)+(x[1]-6)*(x[1]-6);returnz;}floaty,y1,y2,y3;floatx1[2],x2[2],x3[2],s[2],a[2],b[2];inti,n=2;huangjin()/*黄金分割法确定每轮最优值*/{floate=,h0=,r=,h=h0;/*外推法确定区间*/{y1=f(x1);for(i=0;i<n;i++)s[i]=s[i]*h;for(i=0;i<n;i++)x2[i]=x1[i]+s[i];y2=f(x2);if(y2>y1){for(i=0;i<n;i++){s[i]=-s[i];x3[i]=x1[i];y3=y1;}for(i=0;i<n;i++){x1[i]=x2[i];y1=y2;x2[i]=x3[i];y2=y3;}}for(i=0;i<n;i++){x3[i]=x2[i]+s[i];y3=f(x3);}while(y3<y2){for(i=0;i<n;i++)s[i]=s[i]*2;for(i=0;i<n;i++){x1[i]=x2[i];y1=y2;x2[i]=x3[i];y2=y3;}for(i=0;i<n;i++){x3[i]=x2[i]+s[i];y3=f(x3);}}for(i=0;i<n;i++){a[i]=x1[i];b[i]=x3[i];}}for(i=0;i<n;i++){x1[i]=b[i]-r*(b[i]-a[i]);x2[i]=a[i]+r*(b[i]-a[i]);}y1=f(x1);y2=f(x2);if(y1>=y2){for(i=0;i<n;i++)a[i]=x1[i];for(i=0;i<n;i++)x1[i]=x2[i];y1=y2;for(i=0;i<n;i++)x2[i]=a[i]+r*(b[i]-a[i]);y2=f(x2);}else{for(i=0;i<n;i++)b[i]=x2[i];for(i=0;i<n;i++)x2[i]=x1[i];y2=y1;for(i=0;i<n;i++)x1[i]=b[i]-r*(b[i]-a[i]);y1=f(x1);}while(abs(y2-y1)>=e){if(y1>=y2){for(i=0;i<n;i++)a[i]=x1[i];for(i=0;i<n;i++)x1[i]=x2[i];y1=y2;for(i=0;i<n;i++)x2[i]=a[i]+r*(b[i]-a[i]);y2=f(x2);}else{for(i=0;i<n;i++)b[i]=x2[i];for(i=0;i<n;i++)x2[i]=x1[i];y2=y1;for(i=0;i<n;i++)x1[i]=b[i]-r*(b[i]-a[i]);y1=f(x1);}}for(i=0;i<n;i++)x1[i]=*(a[i]+b[i]);}main(){floatx0[2],e0=,sum;x0[0]=;x0[1]=;s[0]=0;s[1]=1;do{if(s[1]<1){s[0]=1;s[1]=0;}else{s[0]=0;s[1]=1;}for(i=0;i<n;i++)x1[i]=x0[i];huangjin();sum=0;for(i=0;i<n;i++)sum+=(x1[i]-x0[i])*(x1[i]-x0[i]);if(sqrt(sum)<e0)break;for(i=0;i<n;i++)x0[i]=x1[i];}while(1);printf("thebestres