funct.c

运筹学非线性规划的c语言编写的外点惩罚函数法子程序

#include "stdio.h"
#include "stdlib.h"
#include "math.h"
const int kkg=2;
const int qkg=1;
double r0;
double f(double x[])
{double ff;
 ff=pow((x[0]-4),2)+pow((x[0]-5),2);
 return(ff);
}
void strain(double x[],double g[])
{double dlt=0.001;
 g[0]=x[0]-1-dlt;
 g[0]=x[1]-4-dlt;
 g[2]=7-x[0]-x[1];
}
double objf(double p[])
{int i;
 double ff,sg,*g;
 g=(double *)malloc((kkg+qkg)*sizeof(double));
 sg=0;
 strain(p,g);
 for(i=0;i<kkg;i++)
  if(*(g+i)<0)
   sg=sg+pow((*(g+i)),2);
  if(qkg!=0)
   for(i=kkg;i<kkg+qkg;i++)
     sg=sg+pow((*(g+i)),2);
 free(g);
 ff=f(p)+r0*sg;
 return(ff);
}