机械优化设计随机法程序设计3页.doc

#include stdio.h#include stdlib.h#include math.hdouble objf(double x)double ff;ff=pow(x0-4),2)+pow(x1-5),2);return(ff);/*约束条件子程序*/double strain(double x)int i;double g3;g0=x0-1;g1=x1-1;g2=8-x0-x1;for (i=0;i3;i+)if(g0)return(0);return(1);double rn(double *r)int m;double s,u,v,p;s=65536.0;u=2053.0;v=13849.0;*r=u*(*r)+v;m=(int)(*r/s);*r=*r-m*s;p=*r/s;return(p);double randopti(double x,double x