運筹学实验单纯形法上机报告.docxVIP

单纯形法—大M法实验报告目录一、 实验目的3二、 单纯形法及大M法31.单纯形法（Simplex Method）3（1）单纯形法是解线性规划问题的一个重要方法。3（2）用程序进行运算前，要将目标函数及约束方程变成标准形式。3（3）对于标准形式的线性规划问题。用单纯形法计算步骤的框图3（4）在程序运算过程中，采用单纯形表显示运算过程。42.大M法4三、 数据结构及模块设计51.以下是程序中用到的数据结构：52.模块设计：5四、 详细设计61.文件格式定义62.void read()读取约束矩阵、目标函数63.void print()单纯形表显示函数74.void init_change()初始矩阵变换，加入松弛变量和人工变量115.void compute_value()计算检验数126.int best_Result()判断是否得到最优解，唯一最优1，无穷多最优2，无界3，无可行5，未得到返回4137.void in_base() 进基选子函数148.void out_base()出基选择子函数149.void row_change()行变换子函数15五、 程序测试及结果161.第1题16（1） 原题16（2） 文件存储16（3） 读取16（4） 初始变换17（5） 运算过程17（6） 运算结果172.第2题18（1） 原题18（2） 文件存储18（3） 读取18（4） 初始变换19（5） 运算过程19（6） 运算结果203.第3题21（1） 原题21（2） 文件存储21（3） 读取21（4） 初始变换22（5） 运算过程22（6） 运算结果224.第2题23（1） 原题23（2） 文件存储23（3） 读取23（4） 初始变换24（5） 运算过程24（6） 运算结果255.第2题25（1） 原题25（2） 文件存储26（3） 读取26（4） 初始变换26（5） 运算过程27（6） 运算结果27六、 工作分配及实验感想281.工作分配282.实验感想28实验目的使用目前熟悉的语言，实现所学的单纯形法之大M法，并正确运算测试结果。本组成员使用C语言实现。 单纯形法及大M法单纯形法（Simplex Method）单纯形法是解线性规划问题的一个重要方法。其原理的基本框架为：第一步：将LP线性规划变标准型，确定一个初始可行解（顶点）。第二步：对初始基可行解最优性判别，若最优，停止；否则转下一步。第三步：从初始基可行解向相邻的基可行解（顶点）转换，且使目标值有所改善—目标函数值增加，重复第二和第三步直到找到最优解。用程序进行运算前，要将目标函数及约束方程变成标准形式。对于非标准形式须作如下变换：目标函数为极小值min z=CX时，转换为max z=-CX形式；在约束方程中有 “≤”时，在加上一个松弛变量；在约束方程中有 “≥”时，采用减去一个松弛变量，再加上一个人工变量；在约束方程中有 “=”时，加上一个人工变量；所有的人工变量，松弛变量的目标函数系数置为0。对于标准形式的线性规划问题。用单纯形法计算步骤的框图在程序运算过程中，采用单纯形表显示运算过程。大M法方法：在约束条件中，加入人工变量后，要求目标函数不受影响，目标函数中人工变量的系数取 –M（M为系统所能表示范围内的一个非常大的值本程序取1000000），其运算过程同单纯形法。理由：目标函数实现最大化，就必须将人工变量从基变量中换出，否则目标函数就不可能取得最大化。数据结构及模块设计程序中用到的数据结构：#define M 20 //最大20个变量#define N 40 //40个约束方程#define Max 1000000//大Mdouble D[M][N];//系数矩阵double C[M];//目标函数系数double Cb[M];//基向量系数double B[M];//约束常数double Value[N];//检验数int Xb[M];//基向量double X0[M];//可行解int op[M];//约束方程符号0---、1---=、2---int m,n;//矩阵行数、列数int begin_n;//初始变量数int In_BaseX=-1;//进基变量int in_n=-1;//进基列标示int out_m=-1;//出基行标示int Out_BaseX=-1;//出基变量int best;//最优函数返回值char name[30];//文件名int ManX_num=0;//人工变量数目int ManX_list[M];//人工变量存放数组模块设计：void read();//读取方程子函数void print();//显示单纯表子函数void init_change();//初始变换子函数void compute_value();//计算检验数子函数