单纯形法C程序源代码_精品文档Word文件下载.doc

单纯形法C程序源代码_精品文档Word文件下载.doc

《单纯形法C程序源代码_精品文档Word文件下载.doc》由会员分享，可在线阅读，更多相关《单纯形法C程序源代码_精品文档Word文件下载.doc（7页珍藏版）》请在冰豆网上搜索。

/*存放出基与入基的变化情况*/

floatdelta[n];

/*存储检验数矩阵*/

floatx[n];

intnum[m];

/*用于存放出基与进基变量的情况*/

floatZB=0;

/*记录目标函数值*/

voidinput（）;

voidprint（）;

intdanchunxing1（）;

intdanchunxing2（inta）;

voiddanchunxing3（inta,intb）;

intdanchunxing1（）

{

inti,k=0;

intflag=0;

floatmin=0;

for（i=0;

i<

n;

i++）

if（delta[i]>

=0）

flag=1;

else{flag=0;

break;

}

if（flag==1）

return-1;

if（min>

delta[i]）

{min=delta[i];

k=i;

returnk;

intdanchunxing2（inta）

inti,k,j;

floatmin;

k=a;

m;

if（A[i][k]<

{printf（"

\n该线性规划无最优解!

\n"

）;

return-1;

if（A[i][k]>

0）

seta[i]=b[i]/A[i][k];

elseseta[i]=M;

min=M;

=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;

j<

j++）

A[c][j]=A[c][j]/temp1;

if（i!

=c）

if（A[i][l]!

{

temp2=A[i][l];

b[i]=b[i]-b[c]*temp2;

for（j=0;

A[i][j]=A[i][j]-A[c][j]*temp2;

}

temp3=delta[l];

delta[i]=delta[i]-A[c][i]*temp3;

voidprint（）

inti,j=0;

printf（"

\n--------------------------------------------------------------------------\n"

printf（"

%8.2f\tX（%d）%8.2f"

CB[i],num[i],b[i]）;

%8.2f"

A[i][j]）;

\t\t\t"

%8.2f"

delta[i]）;

voidinput（）

inti,j;

/*循环变量*/

intk;

请输入方程组的系数矩阵A（%d行%d列）:

m,n）;

scanf（"

%f"

&

A[i][j]）;

\n请输入初始基变量的数字代码num矩阵:

%d"

num[i]）;

\n请输入方程组右边的值矩阵b:

b[i]）;

\n请输入目标函数各个变量的系数所构成的系数阵C:

C[i]）;

delta[i]=C[i];

k=num[i]-1;

CB[i]=C[k];

voidmain（）

intp,q,temp;

input（）;

\tCB\tXB\tb\t"

X（%d）\t"

i+1）;

x[i]=0;

while

（1）

q=danchunxing1（）;

if（q==-1）

print（）;

\n所得解已经是最优解!

\n最优解为:

temp=num[j]-1;

x[temp]=b[j];

for（i=0;

x%d=%.2f"

i+1,x[i]）;

ZB=ZB-x[i]*C[i];

ZB=%.2f"

ZB）;

break;

p=danchunxing2（q）;

\np=%d,q=%d"

p,q）;

if（q==-1）break;

danchunxing3（p,q）;

运行结果如下：

请输入方程组的系数矩阵A（3行5列）:

12100

40010

04001

请输入初始基变量的数字代码num矩阵:

345

请输入方程组右边的值矩阵b:

81612

请输入目标函数各个变量的系数所构成的系数阵C:

-2-3000

--------------------------------------------------------------------------

CB 

XB 

b 

X

（1） 

X

（2） 

X（3） 

X（4） 

X（5）

0.00 

8.00 

1.00 

2.00 

0.00

16.00 

4.00 

X（5） 

12.00 

1.00

-2.00 

-3.00 

p=2,q=1

-0.50

3.00 

0.25

0.75

p=0,q=0

-4.00 

2.00

-