单纯性法C语言程序.docVIP

/* HELLO.C -- Hello, world */ #includestdio.h #includemath.h #define m 4 #define n 6 /*定义未知量的个数*/ float M=1000000.0; float A[m][n]; /*用于记录方程组的数目和系数;*/ float C[n]; /*用于存储目标函数中各个变量的系数*/ float b[m]; /*用于存储常约束条件中的常数*/ float CB[m]; /*用于存储基变量的系数*/ float seta[m]; /*存放出基与入基的变化情况*/ float delta[n]; /*存储检验数矩阵*/ float x[n]; int num[m]; /*用于存放出基与进基变量的情况*/ float ZB=0; /*记录目标函数值*/ void input(); void print(); int danchunxing1(); int danchunxing2(int a); void danchunxing3(int a,int b); /* 检查校验数delta 如果所有校验数delta=0，当前解为最优解 否则找出最大的校验数delta及所在的列号，即确定入基变量 */ int danchunxing1(){ int i=0; /*定义循环变量i*/ int k=0; /*定义最大校验数下表变量k*/ int flag=1; /*定义检查校验数标志flag，初始为1*/ float max=0; /*从第1列开始检查校验数，如果校验数delta存在大于0的，flag=0，并跳出循环*/ for(i=0;in;i++) if(delta[i]0){flag=0;break;} /*flag=1当前解为最优解，返回主函数-1*/ if(flag==1) return -1; /*从第1列的校验数开始检查，寻找最大的校验数delta并返回校验数最大的列号*/ for(i=0;in;i++){ if(delta[i]max) {max=delta[i];k=i;} } return k; } /* 确定入基变量 input：int a入基变量对应的列号 output：如果所有校验数delta=0，当前解为最优解；否则找出最小的seta及所在的行号 */ int danchunxing2(int a){ int i,k,j; int flag=0; float min; k=a; /*针对入基变量对应的列从第1行开始检查A[i][k]，如果所有的A[i][k]=0则该线性规划无最优解*/ for(i=0;im;i++) if(A[i][k]=0) flag=1; else {flag=0;break;} if(flag==1) {printf(\n This linear programming no optimal solution! \n); return -1;} /*针对入基变量对应的列从第1行开始计算seta*/ for(i=0;im;i++){ if(A[i][k]0) seta[i]=b[i]/A[i][k]; else seta[i]=M; } /*寻在最小的seta，确定入基变量*/ min=M; for(i=0;im;i++){ if(min=seta[i]) {min=seta[i];j=i;} } num[j]=k+1; CB[j]=C[k]; return j; } /* 迭代计算，重新获得单纯性表? */ void danchunxing3(int p,int q) { int i,j,c,l; float temp1,temp2,temp3; c=p;/*主元素所在行号*/ l=q;/*主元素所在列号*/ temp1=A[c][l]; /*记录主元素*/ /*主元素所在行的b的计算*/ b[c]=b[c]/temp1