jlmao运筹学-上机报告-v2.0.docVIP

上 机 实 验 报 告 班级： 121 专业/方向：测控 姓名： 李志荣 实 验 成 绩 (10分制) 学号： 2108上机实验名称：编程实现Simplex算法求解LP问题 上机时间：2014-5-28 上机地点：信自楼234自动化系机房 一、实验目的和要求 1、目的： 掌握单纯形算法的计算步骤，并能熟练使用该方法求解线性规划问题。 了解算法(程序实现的过程和方法。 2、要求： 使用熟悉的编程语言编制单纯形算法的程序。 独立编程，完成实验，撰写实验报告并总结。 二、实验内容和结果 1、单纯形算法的步骤及程序流程图。 单纯形算法的步骤 对于一般的标准形式线性规划问题(求极小问题),首先给定一个初始基本可行解。设初始基为B,然后执行如下步骤： (1).解,求得， (2).计算单纯形乘子w, ,得到,对于非基变量，计算判别数,可直接计算令 ,R为非基变量集合 若判别数 ,则得到一个最优基本可行解，运算结束；否则，转到下一步 (3).解,得到；若,即的每个分量均非正数， 则停止计算，问题不存在有限最优解，否则，进行步骤(4).确定下标r,使单纯形算法的程序流程图2、单纯形算法程序的规格说明（程序以附录形式给出）。 输入：（运行程序所需输入的参数） 2 2 1 0 0 4 0 0 1 0 0 5 0 0 1 2 3 0 0 0 12 16 15 0 0 0 输出：（程序执行之后的结果） 功能描述：求解单纯形法 3、使用所编程序求解如下LP问题并给出结果。 例：编程求解LP问题： 输入： 2 2 1 0 0 0 1 2 0 1 0 0 1 0 0 0 1 0 0 1 0 0 0 1 2 3 0 0 0 0 12 8 16 12 0 0 0 0结果： 三、实验总结 1、上机过程中出现的问题及解决方法。 输入会出现输错的情况。 解决：重新输入。 2、本次实验的体会。 深深地感受到计算机的强大，手工要很长时间甚至无法完成的计算，计算机瞬间就可以计算出结果。 【附录】单纯形算法程序（C、C++、matlab语言或其它，请注明！） /* HELLO.C -- Hello, world */ #includestdio.h #includemath.h #define m 3 #define n 5 /*定义未知量的个数*/ float M=1000000.0; float A[m][n]; /*用于记录方程组的数目和系数;*/ float C[n];/*用于存储目标函数中各个变量的系数*/ float b[m];/*用于存储常约束条件中的常数*/ float CB[m];/*用于存储基变量的系数*/ float seta[m]; /*存放出基与入基的变化情况*/ float delta[n]; /*存储检验数矩阵*/ float x[n]; int num[m];/*用于存放出基与进基变量的情况*/ float ZB=0; /*记录目标函数值*/ void input(); void print(); int danchunxing1(); int danchunxing2(int a); void danchunxing3(int a,int b); /* 检查校验数delta 如果所有校验数delta=0，当前解为最优解 否则找出最大的校验数delta及所在的列号，即确定入基变量 */ int danchunxing1(){int i=0;/*定义循环变量i*/int k=0;/*定义最大校验数下表变量k*/int flag=1;/*定义检查校验数标志flag，初始为1*/float max=0;/*从第1列开始检查校验数，如果校验数delta存在大于0的，flag=0，并跳出循环*/for(i=0;in;i++)if(delta[i]0){flag=0;break;}/*flag=1当前解为最优解，返回主函数-1*/if(flag==1)return -1;/*从第1列的校验数开始检查，寻找最大的校验数delta并返回校验数最大的列号*/for(i=0;in;i++){if(delta[i]max){max=delta[i];k=i;}}return k; } /* 确定入基变量 input：int a入基变量对应的列号 output：如果所有校验数delta=0，当前解为最优解；否则找出最小的seta及所在的行号 */ int danchunxing2(int a){int i,k,j;int flag=0;float min;k=a;/*针对入基变量对应的列从第1行开始检查A[i][k]，如果所有的A[i][k]=0则该线性规划