运筹学实验报告单.doc

运筹学实验报告一实验目的与任务通过实验使学生进一步掌握运筹学相关问题的原理、方法和求解过程，加深对运筹学的有关理论、方法的理解，培养学生的分析能力、动手能力和创新意识以及解决问题的能力。为该专业学生学习其他相关专业课程提供有关系统决策和最优化问题的基础知识。同时也为学生今后从事科研和工程实践打下良好的基础。二实验内容根据给定的线性规划问题（Linear Programming），建立线性规划模型（Linear Programming Model）编写相应的单纯形法求解程序，完成题目的求解，并对结果进行分析和讨论。三问题的背景与生产问题（某木器制造厂年度产品生产计划的优化研究） 某木器制造厂是生产木器制品的专业厂家，有多年生产该类产品的丰富经验和悠久历史，无论从生产组织，产品研发，技术力量还是从经营销售等哪一方面来看，在当地同行中均占有相当的优势。但另一方面，也确实存在着管理方法陈旧，管理手段落后的实际问题，尤其是近年来，随着经济体制改革的逐步实施，整个市场情况发生了很大的变化。已计划经济为主，市场调节为辅的指导性生产计划代替了传统的统购包销的指令性计划。这种新形势的出现迫使企业必须提高自身对市场需求变化的适应能力，必须重视和加强对产品生产总体经济效益的分析与研究，努力提高企业的素质。四数学模型对应于产品的编号，记其年产量分别为并记一年的总收益为，则线性规划问题如下 ： 五程序代码（C语言版）#include#includeVoid main Float q,p1,A2131,A02131,b21,b021,x,c31,c031; Int n,i,j,k,m,p,J31,J031,s21; Printf(“Please input the size of the system of equations. the row is no more than 20 and the rank is no more than 30n”); Scanf(“%d,%d”,&m,&n); For(i=1;i=m;i+)Printf (“Please input NO.%d row of the coefficient array A:n”,i);For(j=1;j=n;j+) Scanaf(“%f”,&Aij); Printf(“Please input the resource coefficients b:n”); For(i=1;i=m;i+)Scanf(“%f”,&bi); Printf(“Please input the coefficient of these variables in the target function c:n”); For(j=1;j=n;js+)Scanf(“%f”,&cj); Printf(“Please input the initial basic variables. Distinguish them with 1 or 0:n”); For(j=1;j=n;j+)Scanf(“%d”,&Jj); For(j=1;j=n;j+)J0j=1-Jj; Printf(“Please input the subscripts of the initial basic variables:n”); For(i=1;i=m;i+)Scanf(“%d”,&si);D1:For(i=1;i=m;i+) For(j=1;j=n;j+) A0ij=(-1)*csi*Aij*J0j; For(j=1;j=n;j+)p1=0;For(i=1;i=m;i+) p1=p1+A0ij;c0j=p1+cj*J0j; D2:q=c01; p=1; For(j=2;j=n;j+)If(c0j=0)q=0;Printf(“the optimal solution of the linear programming is:n”);For(i=1;i=m;i+) Printf(“x%d=%fn”,si,bi); q=q+csi*bi; Printf(“the value of the target function is:n”);Return; ElseFor(i=1;i=m;i+) If(Ai,p=0) b0i=-1; Else if (b0i=0)&(Ai,p0) c0p=1; Goto D2; Else b0i=bi/Aip; q=b01; For(i=2;iq) q=b0i; If (q0) Printf(“there is no optimal solution about the linear programming”); Return; Else For(i=1;i0)&(b0i=q) q=b0i; k=i; Jp=1; Jsk=0; J0p=0; J0sk=1; sk=p; x=Akp; For(j=1;j=1;i-) x=Aip; For(j=1;j=n;j+) Aij=Aij-Akj*x; bi=bi-bk*x; For(i=k+1;i=m;i+) x=Aip; For(j=1;j=n;j+) Aij=Aij-Akj*x; bi=bi-bk*x Goto D1;六程序运行结果 七注意事项问题的代码为求min型的问题，故在输入时要对目标函数进行修正。八实验结果解读由于实验中变量下标大于26为松弛变量，故由此可知该工厂应该在接下来的一年中这样计划生产：在众多产品中，该工厂只需生产762步枪箱，B713卡板和