单纯型法的c语言实现！

风花雪月 · 发表于 2005-9-22 09:17

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

/*本程序是单纯形法，参考书目:河海大学<吴凤平>运筹学方法与应用*/ #include<stdio.h> #include<math.h> #define X 5 #define Y 7 void xi_max(int *m2,int *mn1,float *c,int *is,int *ir,int *j0,float (*a)[X][Y]) { int j; *c=0; for(j=1;j<=*is;j++) { if((*a)[*ir][j]-*c>0) { *c=(*a)[*ir][j]; *j0=j; } } } /***************** 参数说明 **********************/ /* m_约束方程个数(基变量个数),n_非基变量个数 */ /* m2-m+2整个变量,(*a)[X][Y]存放初始数据 */ /* (*k)[]存放基变量脚标,(*x)存放基变量最优值 */ /********* <a href="http://happyyangxu.home.sunbo.net*******/" target="_blank" >http://happyyangxu.home.sunbo.net*******/</A> int xi_sm(int m,int n,int m2,int mn1,int l1,float (*a)[X][Y], int(*k)[],float(*x)[]) { int i,m1,mn,j0,i0,j; float c,g; m1=m+1; mn=m+n; for(i=1;i<=m;i++) (*k)=n+i; leap1: if(l1-1==0) xi_max(&m2,&mn1,&c,&mn,&m1,&j0,a); else { leap2: if(l1-50==0) xi_max(&m2,&mn,&c,&n,&m2,&j0,a); else xi_max(&m2,&mn1,&c,&mn,&m2,&j0,a); } c-=1e-8; if((c<=0)&&(l1==1)&&((*a)[m1][mn1]-1e-8>0)) { printf("\n\t*********Not min&&No solution**********\n"); return(-1); } if((c<0)&&(l1==1)) { l1=50; goto leap2; } if(c<=0) { for(i=1;i<=m;i++) (*x)=(*a)[mn1]; printf("\n\t********Optimal solution********\n"); for(i=1;i<=m;i++) printf("\ni=%d,x=%f\n",(*k),(*x)); printf("\nF=%f\n",(*a)[m2][mn1]); return 0; } c=1e8; for(i=1;i<=m;i++) { if((*a)[j0]>1e-8) { g=(*a)[mn1]/(*a)[j0]; if(g-c<0); { c=g; i0=i; } } } if(c==1e+8) { printf("\n\t*******Lp no solution********\n"); return -3; } (*k)[i0]=j0; for(j=1;j<=mn1;j++) { if((j==j0)||(l1==50)&&(n<j)&&(j<mn1)) continue; g=(*a)[i0][j]/(*a)[i0][j0]; (*a)[i0][j]=g; for(i=1;i<=m2;i++) { if((i==i0)||(!(l1==1)&(i==m1))) continue; (*a)[j]=(*a)[j]-(*a)[j0]*g; } } for(i=1;i<=m2;i++) (*a)[j0]=0; (*a)[i0][j0]=1; goto leap1; } main() { float a[5][7]={{0,0,0, 0, 0,0,0}, {0,3,-4,3, 1,0,12}, {0,3,0, 6, 0,1,12}, {0,0,0, 0, 0,0,0}, {0,69,0,144,0,0,300}}; int k[3]; float x[3]; clrscr(); xi_sm(2,3,4,6,0,&a,&k,&x); }

aya331 · 发表于 2005-10-10 22:11

版主可以讲讲没部迭代的原理吗。[em01]

风花雪月 · 发表于 2005-10-11 11:34

<DIV class=quote>以下是引用aya331在2005-10-10 22:11:22的发言： 版主可以讲讲没部迭代的原理吗。[em01]</DIV>
没部迭代是什么意思？没太看懂，请原谅

aya331 · 发表于 2005-10-11 22:51

这个程序应该是把单纯型表中两阶法计算过程用c实现的程序。我的意思是可否解释下程序每部运算的过程，谢谢了。

johhan · 发表于 2005-10-16 23:31

谢谢分享好资源

sprine · 发表于 2006-4-1 16:29

编译出错了

binry · 发表于 2006-4-6 16:24

xiexie[em01][em01]

账号		自动登录	找回密码
密码			我要加入

[C/C++] 单纯型法的c语言实现！

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（aya331）版主可以讲讲没部迭代的原理吗。[em0...