無約束最優(yōu)化方法可變單純形法(simplex)Nelder-MeadWord版

1、無約束最優(yōu)化方法可變單純形法（simplex）Nelder-Mead 可愛的饅頭本程序是用編寫的，從編寫的算例來看，應該是沒有問題的。所采用的原理和步驟是參考華南理工大學出版社蔣金山等編寫的最優(yōu)化計算方法第章第三節(jié)可變單純形法。歡迎各位批評指正。#include<stdio.h>#include<math.h>#include<iostream.h>int i,j;double d3100=0,1,0,0,0,0,1,0,0,0,0,1,f100;/d為單純形的頂點，本算例中未知數(shù)個數(shù)為3，則頂點個數(shù)為4double g,h,l,q,s=1,t=2,u=0.

2、5,v=0.0001,y=0;/s為反射系數(shù)，t為擴展系數(shù)，u為壓縮系數(shù)，v為允許誤差int o,F,r,D,e,lj=0,N=4;/N為頂點的個數(shù)，o為最大值點的位置，F(xiàn)為最小值點的位置，r為次大值點的位置void function1(int e)/求解函數(shù)fe     fe=(d0e-3)*(d0e-3)+2*(d1e+2)*(d1e+2)+(d2e-4)*(d2e-4);/函數(shù)為f=(x1-3)2+2(x2+2)2+(x3-4)2,求其最小值 void function2()   while(+lj)<100)/

3、最大迭代次數(shù)   for(i=0,g=fi;i<N-1;i+)/求最大值點   if(g<fi+1)     g=fi+1;o=i+1;    else            if(i=0)    o=i;  for(i=0,h=fi;i<N-1;i+)/求最小值點

4、60;  if(h>fi+1)     h=fi+1;F=i+1;    else   if(i=0)   F=i;     for(i=0,l=fi;i<N-1;i+)/求次大值點         if(i=o&&i=0)     

5、;l=fi+1;   r=i+1;   continue;    if(l<fi+1&&(i+1)!=o)     l=fi+1;r=i+1;    else   if(i=0)   r=i;     for(i=0;i<N-1;i+)/求除最大值點，其余點的形心 

6、60; diN=0;  for(j=0;j<N;j+)              if(j!=o)                   diN=diN+dij;        

7、;       diN=diN/(N-1);/平均值  function1(N);    for(i=0;i<N-1;i+)   diN+1=diN+(diN-dio)*s;/反射  function1(N+1); if(fN+1<fF)   for(i=0;i<N-1;i+)     diN+2=diN+(diN+1-diN)*

8、t;/擴展       function1(N+2);        if(fN+2<fN+1)              for(i=0;i<N-1;i+)         dio=diN+2;    &#

9、160;    fo=fN+2;   y=0;            for(i=0;i<N;i+)          y=y+(fi-fN)*(fi-fN);         y=sqrt(y/(N*1.0);   

10、60;  if(y<v)       for(i=0;i<N-1;i+)         printf("d%d%d=%f",i,F,diF);     printf("n");          printf("f%

11、d=%f",F,fF);    break;      else       function2();    break;                      &#

12、160;   else              for(i=0;i<N-1;i+)          dio=diN+1;         fo=fN+1;   y=0;      

13、      for(i=0;i<N;i+)          y=y+(fi-fN)*(fi-fN);         y=sqrt(y/(N*1.0);      if(y<v)       for(i=0;i<N-1;i+) &

14、#160;       printf("d%d%d=%f",i,F,diF);     printf("n");        printf("f%d=%f",F,fF);    break;      else  &#

15、160;    function2();    break;          if(fN+1>=fF&&fN+1<=fr)          for(i=0;i<N-1;i+)       dio=diN+1;    

16、60; fo=fN+1;   y=0;         for(i=0;i<N;i+)       y=y+(fi-fN)*(fi-fN);      y=sqrt(y/(N*1.0);   if(y<v)       for(i=0;i<N-1;i+)    

17、     printf("d%d%d=%f",i,F,diF);     printf("n");        printf("f%d=%f",F,fF);    break;      else      

18、 function2();    break;     if(fN+1>fr)   if(fN+1>fo)/壓縮     D=o;    else     D=N+1;          for(i=0;i<N-1;i+)

19、     diN+3=diN+(diD-diN)*u;    function1(N+3);  if(fN+3<=fD)              for(i=0;i<N-1;i+)          dio=diN+3;   

20、      fo=fN+3;   y=0;   for(i=0;i<N;i+)          y=y+(fi-fN)*(fi-fN);         y=sqrt(y/(N*1.0);      if(y<v)   

21、60;   for(i=0;i<N-1;i+)         printf("d%d%d=%f",i,F,diF);     printf("n");        printf("f%d=%f",F,fF);    break; 

22、     else       function2();    break;       else     for(i=0;i<N-1;i+)       for(j=0;j<N;j+)    &#

23、160;    dij=dij+0.5*(diF-dij);/收縮          for(i=0;i<N;i+)       function1(i);      y=0;            for(i=0;i<N;i+)