求约束最优化算法

dyjun888666

哪位能给我一个求解约束最优化问题的算法，如梯度投影法，乘子法等的程序，我已经有罚函数法了，需要另外一种算法进行比较，谢谢了。也可发到我的邮箱dyjun888666@163.com

风花雪月

乘子法，网上收集到的一个c程序

1. #include  <stdio.h>
   
2. #include  <math.h>
   
3. #define  L  10
   
4. 
   
5. void FX(int m,int n,float x[L][L],float z,float s[L][L],float x0[L][L])
   
6. {int i,j;
   
7. for(i=0;i<m;i++)        {for(j=0;j<n;j++)
   
8.                 { x0[i][j]=x[i][j]+z*s[i][j];}}}       
   
9. void grad(int m,float input[L],float x[L][L],float g0[L][L])
   
10. {int i,j;
    
11.   for(i=0;i<m;i++)
    
12.   { for(j=0;j<m;j++)        { if(j!=i)
    
13.                     g0[i][0]=2*input[i]*x[i][0]+input[2]*x[j][0]+input[i+3];}}}
    
14. void qiudao(float input[L],float s[L][L],float x[L][L],float a,float b,float *z)
    
15. {a=input[0]*s[0][0]*s[0][0]+input[1]*s[1][0]*s[1][0]+input[2]*s[0][0]*s[1][0];
    
16. b=2*input[0]*x[0][0]*s[0][0]+2*input[1]*x[1][0]*s[1][0]+input[2]*(x[0][0]*s[1][0]
    
17.          +x[1][0]*s[0][0])+input[3]*s[0][0]+input[4]*s[1][0];
    
18.         a=2*a;*z=-1*b/a;}
    
19. void gongetidu(float input[L],float x[L][L])
    
20. {int i,j,k=0,stop=0;float s[L][L]={0};float g0[L][L]={0},g1[L][L]={0},a=0,b=0,z=0,w=0;
    
21.     while(stop==0)
    
22.         {        grad(2,input,x,g0);
    
23.                 if((g0[0][0]*g0[0][0]+g0[1][0]*g0[1][0])<0.000001)
    
24.                 {stop=1;break;}
    
25.                 for(i=0;i<2;i++)
    
26.                 {for(j=0;j<1;j++) {s[i][j]=-1*g0[i][j];}}
    
27.                 qiudao(input,s,x,a,b,&z);
    
28.            FX(2,1,x,z,s,x);
    
29.        grad(2,input,x,g1);
    
30.           w=(g1[0][0]*g1[0][0]+g1[1][0]*g1[1][0])/(g0[0][0]*g0[0][0]+g0[1][0]*g0[1][0]);
    
31.          for(i=0;i<2;i++)
    
32.                  {for(j=0;j<1;j++){s[i][j]=-1*g1[i][j]+w*s[i][j];}}
    
33.                 qiudao(input,s,x,a,b,&z);
    
34.                 FX(2,1,x,z,s,x);
    
35.                 k++;
    
36.         }
    
37. }
    
38. void main()
    
39. {float f[L]={0},x[L][L]={0};
    
40. float r=0.25,a=5,u=0,c=2,e=0.000001,h,h1,h2;
    
41. int i,j,k=0,stop=1;
    
42. while(stop==1)
    
43. { f[0]=0.5+c/2;f[1]=0.166667+c/2;f[2]=c;f[3]=u-c;
    
44.    f[4]=u-c;f[5]=(u*u+u*c)/(2*c);h1=x[0][0]+x[1][0]-1;
    
45.    gongetidu(f,x);
    
46.    if(fabs(x[0][0]+x[1][0]-1)<e)stop=0;          
    
47.    h2=x[0][0]+x[1][0]-1;
    
48.    h=h1/h2;
    
49.    if(h>r) c=a*c;
    
50.    u=u+c*h2;
    
51.    k++;}
    
52. k--;
    
53. printf("k=%d\n",k);
    
54. printf("x*=");
    
55. {for(i=0;i<2;i++)
    
56.     {for(j=0;j<1;j++)
    
57.         printf("%f  ",x[i][j]);}
    
58. printf(" \n");}
    
59. 
    
60. }

复制代码

donglei06

楼主能不能把求解约束最优化问题的算法罚函数法程序发我邮箱啊,谢谢!
急用!我的邮箱是:donglei06@163.com
能运行出结果的,小弟将不胜感激!!

风花雪月

原帖由 donglei06 于 2007-5-16 21:38 发表

                               
登录/注册后可看大图

楼主能不能把求解约束最优化问题的算法罚函数法程序发我邮箱啊,谢谢!
急用!我的邮箱是:donglei06@163.com
能运行出结果的,小弟将不胜感激!!

搜索本版