单纯形法C语言程序代码.pdfVIP

//单纯形法程序（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;

第1页共5页

floatmin;

k=a;

for(i=0;im;i++)

if(A[i][k]=0)

flag=1;//加上分号

else{flag=0;break;}

if(flag==1)if(flag==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++)

第2页共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()//用于输出，每