机械优化-复合型法C语言程序.docVIP

源程序如下： /*输入值选择n=2,k=3,本程序可以处理n为2或3，k为3或4的情况*/ #include stdio.h #include stdlib.h #include time.h #include math.h #define E0 1e-5 /*复合形法收敛控制精度*/ double **apply(int,int); /*申请矩阵空间*/ double f(double *); /*目标函数*/ double *g(double *); /*约束函数*/ bool judge(double *); /*可行点的判断*/ int main() { int n,k; int i,j,k1; int l; double temporary; double restrain; /*收敛条件*/ double reflect; /*反射系数*/ srand((unsigned)time(NULL)); printf(请输入目标函数的维数n:); /*输入已知数据*/ scanf(%d,n); printf(请输入复合形的顶点数k:); scanf(%d,k); double **x=apply(k,n); /*存放复合形顶点*/ double *y=(double *)calloc(k,sizeof(double)); /*存放目标函数值*/ double *p=(double *)calloc(3,sizeof(double)); /*存放约束函数值*/ double *a=(double *)calloc(n,sizeof(double)); /*存放设计变量的下限*/ double *b=(double *)calloc(n,sizeof(double)); /*存放设计变量的上限*/ double *x_c=(double *)calloc(n,sizeof(double)); /*存放可行点中心*/ double *x_r=(double *)calloc(n,sizeof(double)); /*存放最坏点的反射点*/ printf(本程序中的所有输入，两个数之间用空格隔开，然后按enter键时不要长时间的按，否则，可能会出错\n); printf(请输入选定的第一个可行点x1(包含%d个数):,n); for(i=0;in;i++) scanf(%lf,*x+i); printf(请输入初选变量的下限a(包含%d个数):,n); for(i=0;in;i++) scanf(%lf,a+i); printf(请输入初选变量的上限b(包含%d个数):,n); for(i=0;in;i++) scanf(%lf,b+i); printf(输出输入结果为:\nn=%d,k=%d,x1=(,n,k); /*输出已知数据*/ for(i=0;in-1;i++) printf(%.5lf ,*(*x+i)); printf(%.5lf)\na=(,*(*x+n-1)); for(i=0;in-1;i++) printf(%f ,*(a+i)); printf(%.5lf),b=(,*(a+n-1)); for(i=0;in-1;i++) printf(%f ,*(b+i)); printf(%.5lf)\n,*(b+n-1)); L1: for(i=1;ik;i++) /*随机得到其余(k-1)个可行点*/ for(j=0;jn;j++) *(*(x+i)+j)=*(a+j)+(double)(rand()%10000)/10000*(*(b+j)-*(a+j)); l=1; for(i=1;ik;i++)