单纯形法实验报告.doc

实验报告 ( 2011 -- 2012 年度第 一 学期) 名 称： 运筹学上机实验 题 目： 用编程实现单纯形法 院 系： 数理系 班 级： 学 号： 学生姓名： 指导教师： 设计周数： 1 周 成 绩： 日期：2011年11月20日 一、实验的目的与要求 培养信息类学生面向对象编程技术的基本能力。 掌握单纯形法基本算法。 了解面向对象编程的基本概念，并用程序编写单纯形法，解决一定的应用问题。 二、实验 预期实验成果要求 用程序编写简单的单纯形法计算程序，并能用其解决简单的单纯形法问题。 分析实现过程 1 先把方程化为标准型，列出增广矩阵，然后输入程序； 找到初始可行基，确定初始基可行解； 检验各非基变量Xj的检验数 E[j]= ∑CiAij-Cj，若E[j]大于等于零， 则以得到最优解，可停止计算；否则进入下一步； 根据min（E[j]0）=E[j]，确定Xj为换入变量，按θ规则计算 θ=min（Bi/Aik）=Bl/Alk 可确定Xl为换出变量，转入下一步； 以Alk为主元素进行迭代，把Xk所对应的列向量 变换为（0,0…1…，0）t，将Xb列中的换为Xk，得到新的单纯形法。 重复（2）~（5），直到终止。 系统实现 计算检验数 public void EB(double[,] JD, int[] XS, double[] E,int []B) { for (int j = 0; j c-1; j++) { for (int i = 0; i r ; i++) { E[j] += JD[i, j] * XS[B[i]-1]; } E[j] = E[j] - XS[j]; } } 计算检验数中最小值 public double Min0(double[] E) { double temp = E[0]; for (int l = 0; l c-1; l++) if (E[l] temp) temp = E[l]; return temp ; } 计算检验数最小值所对应的列下表 public int Min1(double[] E) { int min = 0; for (int l = 0; l c-1; l++) { if (E[l] 0) { min = l; break; } } return min; } 计算换出变量所对应的行下表 public int Min2(int[] B, double[,] JD, double[] E) { int L = Min1(E); int min=0; double temp = JD[0, c-1] / JD[0,L]; for (int l = 0; l r; l++) { if (JD[l , c - 1] / JD[l , L]