【正文】 ]={{2400,2600},{1200,1450},{150,200},{,},{,}}。//double bo[nn][2]={{2400,2800},{1200,1600},{150,200},{,},{,}}。double bo[nn][2]={{2400,2800},{1200,1600},{150,200},{,},{,}}。//double x[nn]={, , ,0}。double x[nn]={2600,1437,175,}。//原始數(shù)據(jù)double phamax,lacmin,lacmax。int subject() { int kk=0,i。 double ddd,llmin,llmax。 for (i=0。inn。i++){ if (x[i]bo[i][0]) {kk=1。 return kk。} if (x[i]bo[i][1]) {kk=1。 return kk。} } if (x[4]x[3]/3.) {kk=1。 return kk。} if (L/(2*x[1])) {kk=1。 return kk。} phamax=asin(L/(2*x[1]))。 if (phamax) {kk=1。 return kk。} ddd=x[0]*x[0]powl(x[1]*sin(x[3])x[2],2)。 if (ddd=0) {kk=1。 return kk。} llmin=sqrt(ddd)x[1]*cos(x[3])。 if (llmin=lmin) {kk=1。 return kk。} ddd=x[1]*sin(x[3]+phamax)x[2]。 llmax=x[1]*cos(x[3]+phamax)+sqrt(x[0]*x[0]ddd*ddd)。 lacmin=llmin*llminx[2]*x[2]。 lacmax=llmax*llmaxx[2]*x[2]。 ddd=x[0]*x[0]+x[1]*x[1]2*x[0]*x[1]*cos()。 if (lacminddd) {kk=1。 return kk。} ddd=x[0]*x[0]+x[1]*x[1]+2*x[0]*x[1]*cos()。 if (lacmaxddd) {kk=1。 return kk。} // coutphamax/pi dddd llmin llmaxendl。 // system(pause)。 return kk。}double function() { double ff,i1,pha,phaa,phak。 ff=0。 for (i1=1。i1=10。i1++) { pha=i1*phamax/10。 phaa=pha+x[3]。 phak=pha+x[4]。 ff=ff+sin(phaa)/(2*(sin(phak)*(sin(pha)sin(phaa))+cos(phak)*(cos(pha)cos(phaa))))。 } ff=ff/10。 coutffendl。 system(pause)。 return ff。} main() { int i,j,k,k1,nh,nl,nc,count=0。 double a,df。 char y。 const char convex=39。y39。 double xx[n2][nn],xo[nn],f[n2]。/* coutIf feasile set is not convex, plex method is used. Please input yendl。 ciny。 y=tolower(y)。 if (subject()) {coutX0 is not feasile 。 return 0。} coutsetiosflags(ios::left)。*/ //將X0放到復(fù)合形頂點的數(shù)組中l(wèi)oop1: for (i=0。inn。i++) xx[0][i]=x[i]。 // 用隨機函數(shù)產(chǎn)生復(fù)合形其余頂點 for (i=1。in2。i++) { for (j=0。jnn。j++) { xx[i][j]=bo[j][0]+rand()/32767.*(bo[j][1]bo[j][0])。 x[j]=xx[i][j]。 } if (!subject()) continue。 //向前面各點的形心收縮 for (j=0。jnn。j++) { xo[j]=。 for (k=0。ki。k++) xo[j]+=xx[k][j]。 xo[j]=xo[j]/i。 x[j]=xo[j]。 } if (subject()) { if (y!=convex) {coutFeasile is not convex. endl。 return 0。} } do{ for (j=0。jnn。j++) { xx[i][j]=xo[j]+*(xx[i][j]xo[j])。 x[j]=xx[i][j]。 } } while (subject())。 } //計算復(fù)合形個頂點的函數(shù)值loop2: for (i=0。in2。i++) { for (j=0。jnn。j++) x[j]=xx[i][j]。 f[i]=function()。 } //求最好點和最壞點 loop3: nh=0。 nl=0。 for (i=0。in2。i++) { if (f[i]f[nh]) nh=i。 if (f[i]f[nl]) nl=i。 } //打印每次復(fù)合形最好點和最壞點的函數(shù)值 coutfmax=setw(15)f[nh]fmin=setw(15)f[nl]。 //驗證收斂準(zhǔn)則 df=。 for (i=0。in2。i++) df+=powl((f[i]f[nl]),2)。 df=sqrt(df/(n21))。 if (fabs(f[nl])) df/=fabs(f[nl])。 count+=1。 if (dfee || count==5000*nn) { coutendlendlendlThe last resultsendlendl。 for (i=0。inn。i++) { if (i9) coutX[(i+1)]=setw(5)xx[nl][i] 。 else coutX[(i+1)]=setw(4)xx[nl][i] 。 if ((i+1)%4==0) coutendl。 } coutendlendl。 coutf(X*)=f[nl]endlendl。 coutcount=countendlendl。 return 0。 } //計算除最壞點以外復(fù)合形其余定點的形心并驗證形心是否是可行點 for (i=0。inn。i++) { xo[i]=。 for (j=0。jn2。j++) xo[i]+=xx[j][i]。 xo[i]=(xo[i]xx[nh][i])/(n21)。 x[i]=xo[i]。 } if (subject()) { if (y!=convex) {coutFeasile is not convex. endl。 return 0。} // 修改各變量的上下限 for (i=0。inn。i++) { x[i]=xx[nl][i]。 bo[i][0]=xx[nl][i]。 if (fabs(xx[nh][i]xx[nl][i])2) bo[i][1]=xx[nh][i]。 else bo[i][1]=bo[i][0]+2。 } goto loop1。 }loop4: a=。 //求反射點loop5: do { for (i=0。inn。i++) x[i]=xo[i]+a*(xo[i]xx[nh][i])。 a/=2.。 } while (subject())。 if (function()f[nh]) goto loop7。 if (aeb) goto loop5。 k1+=1。 if (k1=2) goto loop6。 //k1=2說明上次是用次壞點代替最壞點 //找次壞點,用次壞點代替最壞點 if (nh=0) nc=1。 else nc=0。 for (i=0。in2。i++) if ((i!=nh) amp。amp。 (f[i]f[nc])) nc=i。 for (i=0。inn。i++) xx[nh][i]=xx[nc][i]。 f[nh]=f[nc]。 goto loop4。 //向最好點縮邊loop6: for (i=0。in2。i++) { for (j=0。jnn。j++) { xx[i][j]=xx[nl][j]+*(xx[i][j]xx[nl][j])。 x[j]=xx[i][j]。 } if (!subject()) continue。 if (y!=convex) {coutFeasile is not convex. endl。 return 0。} // 修改各變量的上下限 for (j=0。jnn。j++) { x[j]=xx[nl][j]。 bo[j][0]=xx[nl][j]。 if (fabs(xx[nh][j]xx[nl][j])2) bo[j][1]=xx[nh][j]。 else bo[j][1]=bo[j][0]+2。} goto loop1。} k1=0。 goto loop2。 //求反射點代替最壞點loop7: k1=0。 f[nh]=function()。 for (i=0。inn。i++) xx[nh][i]=x[i]。 goto loop3。}致 謝從2014年11月開題到即將來到的畢業(yè)答辯，一共經(jīng)歷了將近半年的時間。畢業(yè)設(shè)計已經(jīng)基本完成了，作為一個本科畢業(yè)生，之前從來沒有自己獨立做過設(shè)計，由于經(jīng)驗的匱乏，難免有許多考慮不周全的地方，如果沒有導(dǎo)師的指引和公司單位同事的幫助，想要完成這個設(shè)計是很困難的。如今我發(fā)現(xiàn)，自身的成長離不開大家的支持和幫助，我很慶幸自己生活在這樣的空間。我要向那些幫助過我、鼓勵過我的人們致以我最誠摯的謝意。首先感謝樓紀(jì)國老師，這次畢業(yè)設(shè)計中，作為我的指導(dǎo)老師，雖然平日里工作繁多，但在我做畢業(yè)設(shè)計的每個階段，從選題到查閱資料，還有論文提綱的確定以及中期論文的修改、后期論文格式調(diào)整等各個環(huán)節(jié)中都給予了我悉心的指導(dǎo)。他嚴(yán)謹(jǐn)?shù)淖鍪聭B(tài)度，淵博的知識，還有平易近人的人格魅力將會對我影響深遠(yuǎn)。其次感謝下羅中華老師，他雖然不是我畢業(yè)設(shè)計的指導(dǎo)老師，但是在我畢業(yè)設(shè)計后期遇到的關(guān)于抓斗優(yōu)化的部分都是他指導(dǎo)我完成的。他指出了我之前優(yōu)化過程中很多的錯誤，并幫助我加以改善，最終我才能得到一個較理想的優(yōu)化結(jié)果，對我畢設(shè)的完成起了很關(guān)鍵的作用。最后我要感謝下我實習(xí)的公司，我的畢設(shè)題目就是來自于公司，他們?yōu)槲耶呍O(shè)的完成提供了很多的資料。編寫畢業(yè)論文是一個學(xué)習(xí)的過程，論文的完成，奠定了我的能力，它會幫助我在今后的工作與學(xué)習(xí)當(dāng)中更加有自信。最后，感謝這四年中陪我一起走過的各位老師和同學(xué)們，你們的信任和慰勉是我前進(jìn)的動力!也要感謝我的朋友們，感謝你們在我失意時給我鼓勵，在失落時給我支持，感謝你們和我一路走來，在此過程中我倍感溫暖!為了那些支持和愛護我的人，我會更加努力，好好工作!好好生活!在此我向我的指導(dǎo)老師，和小組的同學(xué)們致以真摯的謝意，謝謝你們的大力支持和幫助。