单纯形法C语言程序代码

第1页共5页//单纯形法程序（C语言）#includestdio.h#includemath.h#definem4/*定义约束条件方程组的个数*/#definen6/*定义未知量的个数*/floatM=1000000.0;floatA[m][n];/*用于记录方程组的数目和系数;*/floatC[n];/*用于存储目标函数中各个变量的系数*/floatb[m];/*用于存储常约束条件中的常数*/floatCB[m];/*用于存储基变量的系数*/floatseta[m];/*存放出基与入基的变化情况*/floatdelta[n];/*存储检验数矩阵*/floatx[n];/*存储决策变量*/intnum[m];/*用于存放出基与进基变量的情况*/floatZB=0;/*记录目标函数值*/voidinput();voidprint();intdanchunxing1();intdanchunxing2(inta);voiddanchunxing3(inta,intb);intdanchunxing1()//把void改为int{inti,k=0;intflag=0;floatmax=0;//把min改为max，因为这里是求目标函数最大值for(i=0;in;i++)if(delta[i]=0)flag=1;else{flag=0;break;}if(flag==1)return-1;for(i=0;in;i++){if(maxdelta[i])//把min改为max{max=delta[i];k=i;}}returnk;//找到检验数最小的一列i}intdanchunxing2(inta){inti,k,j;intflag=0;第2页共5页floatmin;k=a;for(i=0;im;i++)if(A[i][k]=0)flag=1;//加上分号else{flag=0;break;}if(flag==1){printf(\n该线性规划无最优解!\n);return-1;}for(i=0;im;i++){if(A[i][k]0)seta[i]=b[i]/A[i][k];elseseta[i]=M;}min=M;for(i=0;im;i++){if(min=seta[i]){min=seta[i];j=i;}}num[j]=k+1;CB[j]=C[k];returnj;}voiddanchunxing3(intp,intq){inti,j,c,l;floattemp1,temp2,temp3;c=p;/*行号*/l=q;/*列号*/temp1=A[c][l];b[c]=b[c]/temp1;for(j=0;jn;j++)A[c][j]=A[c][j]/temp1;for(i=0;im;i++){if(i!=c)if(A[i][l]!=0){temp2=A[i][l];b[i]=b[i]-b[c]*temp2;for(j=0;jn;j++)第3页共5页A[i][j]=A[i][j]-A[c][j]*temp2;}}temp3=delta[l];for(i=0;in;i++)delta[i]=delta[i]-A[c][i]*temp3;}voidprint()//用于输出，每个数长度至少为8位且保留2位小数{inti,j=0;printf(\n--------------------------------------------------------------------------\n);for(i=0;im;i++){printf(%8.2f\tX(%d)%8.2f,CB[i],num[i],b[i]);for(j=0;jn;j++)printf(%8.2f,A[i][j]);printf(\n);}printf(\n--------------------------------------------------------------------------\n);printf(\t\t\t);for(i=0;in;i++)printf(%8.2f,delta[i]);printf(\n--------------------------------------------------------------------------\n);}voidinput()//用于输入线性规划问题，分为约束方程组的系数矩阵A[i][j]、初始基变量的数字代码num矩阵、方程组右端的常数矩阵b[i]、目标函数的系矩阵C[i]{//加上花括号inti,j;/*循环变量*/intk;printf(请输入方程组的系数矩阵A(%d行%d列):\n,m,n);for(i=0;im;i++)for(j=0;jn;j++)scanf(%f,&A[i][j]);printf(\n请输入初始基变量的数字代码num矩阵:\n);for(i=0;im;i++)scanf(%d,&num[i]);printf(\n请输入方程组右边的值矩阵b:\n);for(i=0;im;i++)scanf(%f,&b[i]);printf(\n请输入目标函数各个变量的系数所构成的系数阵C:\n);for(i=0;in;i++)scanf(%f,&C[i]);第4页共5页for(i=0;in;i++)delta[i]=C[i];for(i=0;im;i++){k=num[i]-1;CB[i]=C[k];}}main()//程序运行的起点{inti,j=0;intp,q,temp;input();printf(\n--------------------------------------------------------------------------\n);printf(\tCB\tXB\tb\t);for(i=0;in;i++)printf(X(%d)\t,i+1);for(i=0;in;i++)x[i]=0;printf(\n);while(1){q=danchunxing1();if(q==-1){print();printf(\n所得解已经是最优解!\n);printf(\n最优解为:\n);for(j=0;jm;j++){temp=num[j]-1;x[temp]=b[j];}for(i=0;in;i++){printf(x%d=%.2f,i+1,x[i]);//输出最优解ZB=ZB+x[i]*C[i];//计算最优解对应的目标函数值}printf(ZB=%.2f,ZB);//输出目标函数最大值ZB，保留两位小数break;}print();p=danchunxing2(q);第5页共5页printf(\np=%d,q=%d,p,q);if(q==-1)break;danchunxing3(p,q);}}