单纯形法程序Word文档格式.docx
《单纯形法程序Word文档格式.docx》由会员分享,可在线阅读,更多相关《单纯形法程序Word文档格式.docx(14页珍藏版)》请在冰豆网上搜索。
m;
i++)
for(j=0;
j<
j++)
b[i][j]=0.0;
for(i=1;
i++)b[i][i]=1.0;
printf("
*********************************\n"
);
{for(j=0;
n;
printf("
%6.1f"
a[i][j]);
\n"
}
i++)cb[i]=a[0][jj[i]];
for(j=1;
{y[j]=0.0;
for(i=1;
i++)y[j]+=cb[i]*b[i][j];
}
{
for(r=1;
r<
r++)
if(jj[r]==j){cj[j]=0.0;
break;
cj[j]=a[0][j];
i++)cj[j]=cj[j]-y[i]*a[i][j];
maxcj=cj[1];
maxjj=1;
for(j=2;
if(maxcj<
cj[j])
{maxcj=cj[j];
maxjj=j;
k=maxjj;
j++)printf("
cj[j]);
k=%d\n"
k);
while(maxcj>
0)
{count++;
{tb[i]=0.0;
r++)tb[i]+=b[i][r]*a[r][0];
{tp[i]=0.0;
r++)tp[i]+=b[i][r]*a[r][k];
/*for(i=1;
i++)printf("
tb[i]);
tp[i]);
tb--tp\n"
*/
theta=1000.0;
for(j=1;
{if(tp[j]>
if(theta>
tb[j]/tp[j])
{theta=tb[j]/tp[j];
mini=j;
/*printf("
LLL=%dTheta=%f\n"
mini,theta);
l=mini;
lk=tp[l];
jj[l]=k;
for(i=0;
e[i][j]=0.0;
i++)e[i][i]=1.0;
i++)e[i][l]=(-1)*tp[i]/lk;
e[l][l]=1.0/lk;
{for(j=1;
e[i][j]);
/*for(j=1;
tp[j]);
*/
{t[i][j]=0.0;
r++)t[i][j]+=e[i][r]*b[r][j];
{b[i][j]=t[i][j];
b[i][j]);
x%d\n"
jj[i]);
maxcj=cj[1];
for(j=2;
k=maxjj;
z=0.0;
{cb[i]=a[0][jj[i]];
x%d=%6.1f\n"
jj[i],tb[i]);
z+=tb[i]*cb[i];
Max_z=%5.1fTHE_count=%d\n"
z,count);
}_
/*p37--li1-13*/
math.h"
inti,j,r,k,l,jj[4],m,n,maxjj,mini,count=0;
floata[4][6]={{0,700,1200,0,0,0},{360,9,4,1,0,0,},{200,4,5,0,1,0},{300,3,10,0,0,1}};
m=4;
n=6;
jj[1]=3;
jj[2]=4;
jj[3]=5;
0.0001)
k=%dMAXcj=%f\n"
k,maxcj);
}/*while--end*/
/*p47--6(5)*/
inti,j,r,k,l,jj[4],m,n,maxjj,mini,count=0,ok=0;
floatb[4][4],e[4][4],t[4][4],y[4],maxcj,tb[4],tp[4],cb[4],cj[8],theta,lk,z;
floata[4][8]={{0,-1,-1,3,0,0,-1000,-1000},{11,1,-2,1,1,0,0,0},{3,2,1,-4,0,-1,1,0},{1,1,0,-2,0,0,0,1}};
n=8;
jj[2]=6;
jj[3]=7;
0.001&
&
count<
5)
/*for(i=1;
i++)if(tp[i]<
=0)ok=-1;
if(ok==-1)
{printf("
***z\030\354***\n"
return;
}*/
0.01)
/*getchar();
/*p47---6(3)*/
floata[4][7]={{0,-4,-1,0,0,-1000,-1000},{3,3,1,0,0,1,0},{6,4,3,-1,0,0,1},{4,1,2,0,1,0,0}};
n=7;
jj[1]=5;
jj[3]=4;
0.001)
%6.1f