C 单纯形法.docVIP

//在Visual C++控制台程序中编译执行 #includeiostream.h #includemath.h #define M 10000 //全局变量 float kernel[11][31];//核心矩阵表 int m=0,n=0,t=0;//m:结构向量的个数 //n:约束不等式个数 //t:目标函数类型：－1代表求求最小值，1代表求最大值 //输入接口函数 void input() { //读入所求问题的基本条件 cout－－－－－－－－－－参 数 输 入－－－－－－－－－－－endl; cout请按提示输入下列参数：endlendl; cout 结构向量数m: m= ; cinm; coutendl 约束不等式数n: n= ; cinn; int i,j; //初始化核心向量 for (i=0;i=n+1;i++) for (j=0;j=m+n+n;j++) kernel [i][j]=0; //读入约束条件 coutendl 约束方程矩阵的系数及不等式方向（1代表=,-1代表=）:endlendl ; for (i=1;i=m;i++) cout xi; cout 不等式方向 常数项endl; for (i=1;i=n;i++) { cout 不等式i ; for (j=1;j=m+2;j++) cinkernel [i][j]; } for (i=1;i=n;i++) { kernel [i][0]=kernel [i][m+2]; kernel [i][m+2]=0; } //读入目标条件 coutendlendl 目标函数的系数及类型（求最小值:1；求最大值:-1）:endlendl ; for(i=1;i=m;i++) coutxi ; cout类型endl ; cout 目标函数: ; for (i=1;i=m;i++) cinkernel [0][i]; cint; //矩阵调整 if(t==-1) for(i=1;i=m;i++) kernel [0][i]=(-1)*kernel [0][i]; for(i=1;i=n;i++) { kernel [i][m+i]=kernel [i][m+1]; if(i!=1) kernel [i][m+1]=0; } } //算法函数 void comput() { int i,j,flag,temp1,temp2,h,k=0,temp3[10]; float a,b[11],temp,temp4[11],temp5[11],f=0,aa,d,c; //初始化 for(i=1;i=n;i++) temp3[i]=0; for(i=0;i11;i++) { temp4[i]=0; temp5[i]=0; } for(i=1;i=n;i++) { if(kernel [i][m+i]==-1) { kernel [i][m+n+i]=1; kernel [0][m+n+i]=M; temp3[i]=m+n+i; } else temp3[i]=m+i; } for(i=1;i=n;i++) temp4[i]=kernel [0][temp3[i]]; //循环求解 do{ for(i=1;i=m+n+n;i++) { a=0; for(j=1;j=n;j++) a+=kernel [j][i]*temp4[j];