单纯形法C程序源代码.doc

单纯形法C程序源代码2008-05-15 15:59#include#include#define m 3 /*定义约束条件方程组的个数*/#define n 5 /*定义未知量的个数*/float M=1000000.0;float Amn; /*用于记录方程组的数目和系数;*/float Cn; /*用于存储目标函数中各个变量的系数*/float bm; /*用于存储常约束条件中的常数*/float CBm; /*用于存储基变量的系数*/float setam; /*存放出基与入基的变化情况*/float deltan; /*存储检验数矩阵*/float xn;int numm; /*用于存放出基与进基变量的情况*/float ZB=0; /*记录目标函数值*/void input();void print();int danchunxing1();int danchunxing2(int a);void danchunxing3(int a,int b);int danchunxing1()int i,k=0;int flag=0;float min=0;for(i=0;i=0)flag=1;else flag=0;break;if(flag=1)return -1;for(i=0;ideltai) min=deltai;k=i;return k;int danchunxing2(int a)int i,k,j;int flag=0;float min;k=a;for(i=0;im;i+)if(Aik=0)flag=1;else flag=0;break;if(flag=1)printf(n该线性规划无最优解!n); return -1;for(i=0;i0) setai=bi/Aik; else setai=M;min=M;for(i=0;i=setai) min=setai;j=i;numj=k+1;CBj=Ck;return j;void danchunxing3(int p,int q)int i,j,c,l;float temp1,temp2,temp3;c=p;/*行号*/l=q;/*列号*/temp1=Acl;bc=bc/temp1;for(j=0;jn;j+)Acj=Acj/temp1; for(i=0;im;i+) if(i!=c) if(Ail!=0) temp2=Ail; bi=bi-bc*temp2; for(j=0;jn;j+) Aij=Aij-Acj*temp2; temp3=deltal;for(i=0;in;i+)deltai=deltai-Aci*temp3;void print()int i,j=0;printf(n-n);for(i=0;im;i+) printf(%8.2ftX(%d) %8.2f ,CBi,numi,bi); for(j=0;jn;j+) printf(%8.2f ,Aij); printf(n);printf(n-n);printf(ttt);for(i=0;in;i+)printf( %8.2f,deltai);printf(n-n);void input()int i,j; /*循环变量*/int k;printf(请输入方程组的系数矩阵A(%d行%d列):n,m,n);for(i=0;im;i+)for(j=0;jn;j+)scanf(%f,&Aij);printf(n请输入初始基变量的数字代码num矩阵:n);for(i=0;im;i+)scanf(%d,&numi);printf(n请输入方程组右边的值矩阵b:n);for(i=0;im;i+)scanf(%f,&bi);printf(n请输入目标函数各个变量的系数所构成的系数阵C:n);for(i=0;in;i+)scanf(%f,&Ci);for(i=0;in;i+)deltai=Ci;for(i=0;im;i+) k=numi-1; CBi=Ck;void main()int i,j=0;int p,q,temp;input();printf(n-n);printf( tCBtXBtbt);for(i=0;in;i+)printf( X(%d)t,i+1);for(i=0;in;i+)xi=0;printf(n);while(1) q=danchunxing1(); if(q=-1) print(); printf(n所得解已经是最优解!n); printf(n最优解为:n); for(j=0;jm;j+) temp=numj-1; xtemp=bj; for(i=0;in;i+) printf(x%d=%.2f ,i+1,xi); ZB=ZB-xi*Ci; printf(ZB=%.2f,ZB); break; print(); p=danchunxing2(q); printf(np=%d,q=%d,p,q); if(q=-1) break; danchunxing3(p,q);运行结果如下：请输入方程组的系数矩阵A(3行5列):1 2 1 0 04 0 0 1 00 4 0 0 1请输入初始基变量的数字代码num矩阵:3 4 5请输入方程组右边的值矩阵b:8 16 12请输入目标函数各个变量的系数所构成的系数阵C:-2 -3 0 0 0- CB XB b X(1) X(2) X(3) X(4) X(5)- 0.00 X(3) 8.00 1.00 2.00 1.00 0.00 0.00 0.00 X(4) 16.00 4.00 0.00 0.00 1.00 0.00 0.00 X(5) 12.00 0.00 4.00 0.00 0.00 1.00- -2.00 -3.00 0.00 0.00 0.00-p=2,q=1- 0.00 X(3) 2.00 1.00 0.00 1.00 0.00 -0.50 0.00 X(4) 16.00 4.00 0.00 0.00 1.00 0.00 -3.00 X(2) 3.00 0.00 1.00 0.00 0.00 0.25- -2.00 0.00 0.00 0.00 0.75-p=0,q=0- -2.00 X(1) 2.00 1.00 0.00 1.00 0.00 -0.50 0.00 X(4) 8.00 0.00 0.00 -4.00 1.00 2.00 -3.00 X(2) 3.00 0.00 1.00 0.00 0.00 0.25- 0.00 0.00 2.00 0.00 -0.25-p=1,q=4- -2.00 X(1) 4.00 1.00 0.00 0.00 0.25 0.00 0.00 X(5) 4.00 0.00 0.00 -2.00 0.50 1.00 -3.00 X(2