实验一单纯形法.docVIP

实验一 单纯形法 I 实验目的： 1. 掌握单纯形法的基本思路和迭代步骤； 2 培养编程与上机调试程序的能力。 II算法思想： 从一个基本可行解出发，判断其是否为最优解，如果不是最优解，则转换到相邻的基本可行解，并使目标函数值不断增大，直到找出最优解或判断有无界解、无解为止。 III程序功能： 在输入化为标准型（含有初始基变量）后的约束条件系数矩阵、目标函数的系数及约束条件的右端项以后，该程序能找出初始基，能处理循环退化问题，也能判断出存在无界解和存在无穷多最优解的最大化问题。另外该程序只能计算最大化问题，对于最小化问题必须在转化为最大化问题后才能进行输入计算。若问题要加入人工变量，也必须在化为标准型（含有初始基变量）后才能进行输入计算。 IV算法整体计算过程： 输入各参数； 查找初始基，若找不到则返回； 若找到初始基，则计算初始检验数，并把初始基和初始检验数复制到用于比较的初始基和初始检验数组中； 检查所有的检验数是否都小于等于零，若不是则转到下一步，若是则判断非基变量检验数小于零，有则存在无穷多最优解，无则存在唯一最优解，然后调出循环； 若不是所有的检验数是否小于等于零，则判断大于零的检验数所对应的约束条件系数矩阵的列是否小于等于零，若有则存在无界解并调出循环，无则转到下一步； ]选入基变量和出基变量； 计算基变量并输出单纯形表； 比较检验数、入基变量和初始检验数、初始入基变量是否相等，若相等着设置循环退化标志并跳出循环，不相等，则转到（4）； 如存在循环退化，则转到进行以下几步： 检验所有的检验数是否都小于等于零，若不是则把初始基和初始检验数复制到用于比较的初始基和初始计算检验数数组中，然后转到下一步，若是则判断是否存在非基变量检验数小于零，有则存在无穷多最优解，无则存在唯一最优解，然后跳出循环； 若不是所有的检验数是否都小于等于零，则判断大于零的检验数所对应的约束条件系数矩阵的列是否都小于等于零，若有则存在无界解，并跳出循环，无则转到下一步； 用BLAND规则选入基变量和出基变量； 计算检验数并输出单纯形表，然后转到a。 注：程序流程图略 V运算实例与结果分析： Max z=2x1+x2 x2=15 6x1+2x2=24 s.t. x1+x2=5 xj=0 j=1,2 化为标准形式有： Max z=2x1+x2 x2+x3=15 6x1+2x2+x4=24 s.t. x1+x2+x5=5 xj=0 j=1,2,3,4,5 输入提示： equator Variable 变量个数； Restrictor 约束条件个数；（注：这里程序默认X(i)0；否则对每个X(i)我们均要添加一个约束条件X(i)0） b[] 约束条件右端项； Restrictor code 约束条件不等式符号编码;(0表示=；1表示=；2表示=) The XiShu 各约束条件中X(i)对应的系数 The Type Z的最优方向（0表示MIN;1表示MAX） The Z 目标函数的系数 VI计算结果如下： 最优解如下： X=(3.500,1.500,7.500,0.000,0.000)，其中： X(1)=3.500 X(2)=1.500 X(3)=7.500 最大化问题目标函数为：8.500 注：请检查所得的解中是否含有人工变量，若有则此问题无可行解。 VII 体会与总结: 通过这次实验使我真正知道，在清楚了一个算法思想后，但不能用语言描述出来，感到自己语言这方面的匮乏。更重要的是调试以及函数模块划分的重要性。总之，通过本次实验，基本掌握了单纯形法的算法思想以及语言描述。 VIII源程序： #includestdio.h #includemath.h #includeiostream.h float matrix[100][100],x[100]; /* 记录总方程的数组,解的数组 */ int a[100]; /* 记录基础,非基础的解的情况,0:非基础,1:基础 */ int m,n,s,type; /* 方程变量,约束数,求最大最小值的类型,0:最小 1:最大 */ int indexe,indexl,indexg; /* 剩余变量,松弛变量,人工变量 */ void Jckxj() { int i,j; for(i=0;in;i++) for(j=0;js;j++) if(matrix[i][j]==1a[j]==1) { x[j]=matrix[i][s]; j=s; } for(i=0;is;i++) if(a[i]==0) x[i]=0; } int