C语言优化算法：直线摸索与共轭梯度法应用

C-一维直线搜索2008年12月05日 星期五 10:04type *fl_s(type *p,type *q)              /*一维直线搜索（不精确）*/
{
type r1=0.2,r2=0.6,r3=2,mp2,mp;
type t=1,g1,g2,g3=0,fs1,fs2,ra=0,rb=5000;
type a[num],a1[num],b[num],g[num];
int i,sig=1,k=0;
for(i=0;i<num;i++){a[i]=*p;b[i]=*q;p++;q++;} 
for(i=0,mp2=0;i<num;i++){mp2=mp2+(b[i]*b[i]);}
mp=sqrt(mp2);
for(i=0;i<num;i++)b[i]=b[i]/mp;
fs1=f_c(&a); 
p=f_d(&a); 
for(i=0;i<num;i++) {g[i]=*p; p++; }

for(i=0,g1=0;i<num;i++) {g1=b[i]*g[i]+g1; }
while(sig)
{k++;
   for(i=0;i<num;i++) a1[i]=a[i]+t*b[i];
   fs2=f_c(&a1); 
   p=f_d(&a1); 
   for(i=0;i<num;i++){g[i]=*p; p++;}
   for(i=0,g2=0;i<num;i++) {g2=b[i]*g[i]+g2; }
   if(g2==g3)sig=0;
   if(g2<(r2*g1)){ra=t; t=(((ra+rb)/r3)>2*ra)?(2*ra):((ra+rb)/r3);g3=g2;}
   else if(fs1-fs2<-r1*t*g1){rb=t;t=(ra+rb)/r3;}
        else sig=0;
   }
return(&a1);

 

C-共轭梯度法(Fletcher-Reeves)2008年12月05日 星期五 10:06Fletcher-reeves()                              /*共轭梯度法*/
{
FILE *fp;
char filename[40]={"E:\date.txt"};
type f,f0,*q,rg,rx=0,rm=0,rf,h1=0.000001,h2=0.000001,h3=0.00001,e=0.001;
type x[num],p0[num],p[num],g[num],g0[num],m[num],rg0,rpg,rp,*tt,mp,mp2;
int i,k,count=0;
type x0[num]={3,9,1,9,6};            /*设置初始点*/

fp=fopen(filename,"a");
printf("\n%s has created",filename);
getch();
for(i=0;i<num;i++)x[i]=x0[i];
fprintf(fp,"\n方程维数 %d",num);             
fprintf(fp,"\n初值 x0 ");                    /*写入文件*/
for(i=0;i<num;i++)fprintf(fp," %f ",x0[i]);
f=f_c(&x); 
q=f_d(&x); 
for(i=0;i<num;i++){g[i]=*q;q++;}
for(i=0;i<num;i++){x0[i]=x[i];p0[i]=-g[i];}
while(1)
{count++;
for(i=0;i<num;i++){x0[i]=x[i];g0[i]=g[i];}   
f0=f;
tt=fl_s(&x0,&p0);
for(i=0;i<num;i++,tt++)x[i]=*tt;
f=f_c(&x); 
q=f_d(&x); 
for(i=0;i<num;i++){g[i]=*q; q++; }
printf("\nmain-x");