算法及编程语言 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断

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

单纯型法的c语言实现！

/*本程序是单纯形法，参考书目:河海大学<吴凤平>运筹学方法与应用*/ #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)) { int j; *c=0; for(j=1;j<=*is;j++) { if((*a)[*ir]-*c>0) { *c=(*a)[*ir]; *j0=j; } } } /***************** 参数说明 **********************/ /* m_约束方程个数(基变量个数),n_非基变量个数 */ /* m2-m+2整个变量,(*a)存放初始数据 */ /* (*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), 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)-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); 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)); return 0; } c=1e8; for(i=1;i<=m;i++) { if((*a)>1e-8) { g=(*a)/(*a); if(g-c<0); { c=g; i0=i; } } } if(c==1e+8) { printf("\n\t*******Lp no solution********\n"); return -3; } (*k)=j0; for(j=1;j<=mn1;j++) { if((j==j0)||(l1==50)&&(n<j)&&(j<mn1)) continue; g=(*a)/(*a); (*a)=g; for(i=1;i<=m2;i++) { if((i==i0)||(!(l1==1)&(i==m1))) continue; (*a)=(*a)-(*a)*g; } } for(i=1;i<=m2;i++) (*a)=0; (*a)=1; goto leap1; } main() { float a={{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; float x; clrscr(); xi_sm(2,3,4,6,0,&a,&k,&x); }

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

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

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

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

<DIV class=quote>以下是引用aya331在2005-10-10 22:11:22的发言： 版主可以讲讲没部迭代的原理吗。</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

页: [1]

声振论坛's Archiver

单纯型法的c语言实现！

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