欢迎来到天天文库
浏览记录
ID:21822086
大小:293.16 KB
页数:7页
时间:2018-10-24
《机械优化-复合型法c语言程序》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、XxL,ffL一次坏点xG代替最坏点XH(一)源程序如下:/*输入值选择n=2,k=3,本程序可以处理n为2或3,k为3或4的情况*/#include^include#include#defineEOle-5double**apply(int,int);doublef(double*);double*g(double*);booljudge(double*);intmain(){intn,k;inti,j,kl;int1;
2、doubletemporary;doublerestrain;doublereflect;合形法收敛控制精度*//*申请矩阵空间*//*目标函数*//*约束函数*//*可行点的判断*//*收敛条件*//*反射系数*/srand((unsigned)time(NULL));primf(”请输入目标函数的维数n:”);scanf(u%d",&n);printf(n请输入复合形的顶点数k:n);scanf("%d",&k);double林x=apply(k,n);double*y=(double*)calloc(
3、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));/*输入已知数裾*//*存放复合形顶点
4、*//*存放目标函数值*//*存放约束函数值*//*存放设计变量的下限*//*存放设计变量的上限*//*存放可行点中心*//*存放最坏点的反射点*/printf(”本程序中的所有输入,两个数之间用空格隔开,然后按enter键时不要长时间的按,否则,可能会出错n);printf(n请输入选定的策一个可行点*1(包含%(1个数):n,n);for(i=0;i5、nf(n%lf’,a+i);printf("请输入初选变量的上限b(包含%(1个数):",n);for(i=0;i6、=(',,*(a+n-1));for(i=0;i7、*(x+i))){for(j=l;j8、rary=*(*(x+i)+k1);*(*(x+i)+kl)=*(*(x+j)+kl);*(*(x+j)+kl)=temporary;}for(i=0;i
5、nf(n%lf’,a+i);printf("请输入初选变量的上限b(包含%(1个数):",n);for(i=0;i6、=(',,*(a+n-1));for(i=0;i7、*(x+i))){for(j=l;j8、rary=*(*(x+i)+k1);*(*(x+i)+kl)=*(*(x+j)+kl);*(*(x+j)+kl)=temporary;}for(i=0;i
6、=(',,*(a+n-1));for(i=0;i7、*(x+i))){for(j=l;j8、rary=*(*(x+i)+k1);*(*(x+i)+kl)=*(*(x+j)+kl);*(*(x+j)+kl)=temporary;}for(i=0;i
7、*(x+i))){for(j=l;j8、rary=*(*(x+i)+k1);*(*(x+i)+kl)=*(*(x+j)+kl);*(*(x+j)+kl)=temporary;}for(i=0;i
8、rary=*(*(x+i)+k1);*(*(x+i)+kl)=*(*(x+j)+kl);*(*(x+j)+kl)=temporary;}for(i=0;i
此文档下载收益归作者所有