运筹学上机报告.docVIP

上 机 实 验 报 告 班级：测控121 专业/方向：测控技术与仪器 姓名： 实 验 成 绩 (10分制) 学号： 上机实验名称：编程实现Simplex算法求解LP问题 上机时间：2014年5月21日 上机地点：信自(程序实现的过程和方法。 2、要求： 使用熟悉的编程语言编制单纯形算法的程序。 独立编程，完成实验，撰写实验报告并总结。 二、实验内容和结果 1、单纯形算法的步骤及程序流程图。 单纯形法的步骤： （1）将线性规划化为标准形。 （2）用最快的方法确定一个初始基本可行解X（0）。当s·t均为“≤”形式时，以松驰变量做初始基本变量最快。 （3）求X（0）中非基本变量xj的检验数σj。若，则停止运算，X（0）=X*（表示最优解），否则转下一步。 （4） ①由确定xk进基； ②由确定xl出基，其中alk称为主元素； ③利用初等变换将alk化为1，并利用alk将同列中其它元素化为0，得新解X（1）。 （5）返回（3），直至求得最优解为止。 程序流程图： 2、单纯形算法程序的规格说明（程序以附录形式给出）。 输入：（运行程序所需输入的参数） 系数矩阵A 初始基变量的数字代码num矩阵 方程组右边的值矩阵b 目标函数各个变量的系数所构成的系数阵C 输出：（程序执行之后的结果） 功能描述： 用单纯型法优P26例5 输入： 结果： 三、实验总结 1、上机过程中出现的问题及解决方法。 程序提示错误,但是无法找出错误通过找同学帮忙解决警告---函数未定义返回值，添加返回值返回值后结果不个松弛变量代码编写要严谨心细若所研究的线性规划模型的约束条件全是小于类型，引入松弛变量。 #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); int danchunxing1() { int i,k=0; int flag=0; float min=0; for(i=0;in;i++) if(delta[i]=0) flag=1; else {flag=0;break;} if(flag==1) return -1; for(i=0;in;i++) { if(mindelta[i]) {min=delta[i];k=i;} } return k; } int danchunxing2(int a) { int i,k,j; int flag=0; float min; k=a; for(i=0;im;i++) if(A[i][k]=0) flag=1; else {flag=0;break;} if(flag==1) {printf(\n该线性规划无最优解!\n); return -1;} for(i=0;im;i++) { if(A[i][k]0) seta[i]=b[i]/A[i][k]; else seta[i]=M; } min=M; for(i=0;im;i++) { if(min=seta[i]) {min=seta[i];j=i;} } num[j]=k+1; CB[j]=C[k]; return j; } void danchunxing3(int p,int q) { int i,j,c,l; float temp1,temp2,temp3; c=p;/*行号*/ l=q;/*列号*/ temp1=A[c][l]; b[c]=b[c]/temp1; for(j=0;jn;j++) A[c][j]=A[c][j]/temp1; for(i=0;im;i++) { if(i!=c) if(A[i][l]!=0) {