对偶单纯形法编程.docVIP

云南大学数学与统计学实验教学中心 实 验 报 告 课程名称：运筹学实验 学期：2012~2013学年第二学期 成绩： 指导教师：葛瑜 学生姓名：卢富毓 学生学号：20101910072 实验名称：对偶单纯形法编程 实验要求： 必做 实验学时：2学时 实验编号：02 实验日期： 2013/3/8 完成日期：2013/3/11 学院： 数学与统计学院 专业 ： 信息与计算科学 年级： 2010级 一、实验目的、环境 ，找到其所换入的列以保证的到的对偶问题仍有可行解。 然后再按照原单纯形算法经行求解。重复上述步骤。 B、编译运行程序，输入约束方程的系数矩阵A，常矩阵B,价值系数C,得到最优解 为了保证得到的算法既能对特殊情形（有bi0）的情况下能求解，也能对一般的问题求解。这里给出了一下两组数据。 1、上一个实验报告中的求解的数据。 B=[2 ; 14 ; 2]; X=[3 -4 2 0 -5 5 0 0] 2、参照课本P62例题6中的2-6表给出的数据 X=[-2 -3 -4 0 0] D、运行结果如下图： 1）下图表示的是第一个数据带入后得到的结果。表示对偶单纯形法对一般的数据成立 2）下图表示的对偶单纯形法对b中存在负值时纠正的结果： a、初始化后的结果： b、优化后的结果： c、最后得到的结果： 最后得到的结果与书上例题所有的结果完全吻合。所以所编写的算法为对偶单纯形法。 五、总结，b的变换有了更加清晰的理解。 以及在求解过程中，对线性约束的标准型也有了明确的掌握。 六、源代码 源代码如下： #includeiostream #includeiomanip using namespace std; #define M -1000 class DanChun{ private: double A[100][120],B[30]; //创建一个二维数组，用来存储xi值以及b的值 double C[110],CB[40];//定义一个C矩阵用来存储x 的价值系数 double CZ[110];//CZ用来存储cj-zj的值 double sita[30];//用来存储sita的值 int XB[30];//用来存XB的值 double X[100];//用来存储最后计算的X的值 double w[100]; //Cj-Zj/X[i]的替换。 public: int n,m; void init();//1.将不等式约束的各矩阵进行赋值 bool FindE(int i, int j);//2.寻找A矩阵中的基向量 void FindCB(); //3.寻找有用价值系数，Cj-Zj需要用到的 int Cj_Zj(); //4.1求出cj-zj的值 int C_sita(int k);//4.2计算sita的值 void Com_DC();//4.用来计算单纯形法的主函数 void Com_E(int l, int k);//4.3求其单位向量 bool P_CjZj();//4.4 int Bbool(double b[], int k); //判断b中是否有负值存在 int TH(int l); //b有负值是，进行替换操作 void Display(); void Display1(); }; void DanChun::init(){ int i=0,j=0; cout请输入 m,n ：; cinm; cinn; cout请输入线性约束矩阵：endl; for(i=0; im;i++){ for(j=0;jn+m; j++){ if(jn){ cinA[i][j]; } else { A[i][j]=0; if(j==n+i){ A[i][j]=1; } } } } coutendl; cout请输入B的值：endl; for(i=0;im;i++){ cinB[i]; CB[i]=0; } coutendl; cout请输入价值系数 C：endl; for(i=0;in+m;i++){ if(in){ cinC[i]; } else{ C[i] = M; } X[i]=0; } n = m+n; } void DanChun::Display(){ int i,j; cout价值系数 C：; for(i=0;in;i++){ coutC[i] ; } coutendl; coutB的值：; for(i=0;im