機(jī)械優(yōu)化設(shè)計二維搜索中坐標(biāo)輪換法求解C++程序

二維搜索題目：y=1.9*xl*xl+l.l*x2*x2-1.2*xl*x2+ll*xl-13*x2+19

初始點(diǎn)(1,2),精度0.000001運(yùn)行結(jié)果:x(-1.241777,5.232949)y=-21.838148⑥E:\C語M二維搜T\erwei.exeinput:x2tl],x2(2]:13N=1a=7.292000b=8.708000N=2a=-0.854000b=-0.146000[-0.855358,6.455154]p=3.921791N=1a=4.601154b=5.309155N=2a—1.091358b—0.619358[-1.175135,5.442792]p=l.061666N=1a=5.206792b=5.678792N=2a=-l.411135b=-0.939135[-1.230134,5.268620]pM.182649N=1a=5.032620b=5.504620N=2a—1.466134b—0.994134[-1.239196,5.239353]pM.030638N=1a=5.003353b=5.475353N=2a=-l.475196b=-l.003196[-1.241777,5.233534]p=0.006366N=1『4.997534b=5.469534N=2a—1.477777b—1.005777[-1.241777,5.232949]p=0.000585N=1a=4.996949b=5.468949N=2a—1.477777b=-l.005777[-1.241777,5.232949]p=0.000000x=[-1.241777,5.232949]f(x)=-21.838148#iiiclude<math.h>#include<stdio.h>#include<stdlib.h>/*頭文件，第3個配合主程序最后一句，是讓程序畫面停留功能*/floatf(floatxl,floatx2)(floaty;y=L9*xl*xl+l.l*x2*x2-1.2*xl*x2+ll*xl-13*x2+19;returny; /*二元二次方程，以功能塊形式定義，等待調(diào)用*/}voidfun(floatx,floataO.floathO.float*a,float*b,intN)/*外推法尋找一維搜索區(qū)間，當(dāng)N=1時，x2被一個常數(shù)代替，成為對xl進(jìn)行一維搜索*/(floatyl,y2,y3,al,a2,a3,h;al=a0;h=hO;switch(N)(case1:yl=f(x,al);a2=al+h;y2=Rx,a2);break;case2:yl=Ral,x);a2=al+h;y2=Ra2,x);break;}if(y2>yl)(h=-h;a3=al;y3=yl;chf: al=a2;yl=y2; /*尋找高低高*/a2=a3;y2=y3;}a3=a2+h;switch(N)(case1:y3=f(x,a3);break;case2:y3=Ra3,x);break;}If(y3<y2)(h=2*h;gotochf;}else(lf(h<0)JI*a=a3;*b=al;}else(*a=al;*b=a3;} /*確定區(qū)間上下限a,b*/}printf(”N=%da=%fb=%f\n”,N,a,b);}floatfuii2(floatx^floata.floatbjntN)/*黃金分割法程序*/( floata1,a2,yl,y2,result;al=a+0.382*(b-a);a2=a+0.618*(b-a);switch(N) /*黃金分割法法尋找一維搜索區(qū)間，當(dāng)N=1時，x2被一個常數(shù)代替，該功能塊只對一個維度進(jìn)行分割求解*/(case1:yl=f(x,al);y2=f(x,a2);break;case2:yl=Ral,x);y2=f(a2,x);break;}wlule(fabs((b-a))>0.000001) /*把黃金分割的精度進(jìn)行設(shè)置*/if(yl>=y2)

a=al;al=a2;yl=y2;a2=a+0.618*(b-a);switch(N){case1:y2=Rx,a2);break;case2:y2=Ra2,x);break;))else)elseb=a2;a2=al;y2=yi;b=a2;a2=al;y2=yi;al=a+0.382*(b?a);switch(N){case1:yl=Rx,al);break;case2:yl=Ral,x);break;))}result=(a+b)/2;letuin(result);J*最后把某一方向上的“下一點(diǎn)坐標(biāo)”求出*J/*主函數(shù)*/voidmain()/*主函數(shù)*/voidmain()floataO,110=1,p,a.b,y,xl[2]={0},x2[2]={0};/*用aO把暫不需要做搜索的維度上的坐標(biāo)值以常數(shù)形式送入*/inti=0;pnntfC'input:x2[l],x2[2]:\ii");scanf(”％必f',&x2[l],&x2⑵)；do{xl[l]=x2[l];xl[2]=x2[2]; /*把前一點(diǎn)的xl,x2方向的坐標(biāo)值，分別傳遞給新點(diǎn)，也就是把新得到的點(diǎn)作為舊點(diǎn)進(jìn)行下一輪搜索*/fiin(xl[l],xl[2],h0,&a,&b,l); /*N=1,x2方向被以常數(shù)形式固定，送入xl方向的初值xll,按照初始步長為1進(jìn)行xl方向搜索*/x2[2]=fiiii2(xl[l],a,b,l);/*由給定前一點(diǎn)的xl方向的值，得到此方向下一個點(diǎn)的值*/fim(x2[2],xl[l],h0,&a,&b,2);x2[l]=fiin2(x2[2],a,b,2);i=i++;pnntf(”(％d)[%f,%f]”,i,x2[l],x2[2]);P=sqrt((x2[l]-xl[l])*(x2[l]-xl[l])+(x2[2]-xl[2])*(x2[2]-xl[2]));pnntf(”p=%f\n”,p);/*新舊兩點(diǎn)之間的距離p*/}wlule(sqn((x2[1]-x1[1])*(x2[1]-x1[1])+(X2⑵-xl[2])*(x2[2]-xl[2]))>=0.