单纯形法实验报告

《单纯形法实验报告》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、实验报告（2011—2012年度第一学期）名称：运筹学上机实验题目：用编程实现单纯形法院系：数理系班级：学号:学生姓名：指导教师：设计周数：1周成绩：口期：2011年11月20日一、实验的目的与要求1.培养信息类学生而向对象编程技术的基本能力。2.掌握单纯形法基本算法。3.了解而叫对象编程的基本概念，并用程序编写单纯形法，解决一定的应用问题。二、实验预期实验成果要求用程序编写简单的单纯形法计算程序，并能用其解决简单的单纯形法问题。分析实现过程1(1)先把方程化为标准型，列出增广矩阵，然后输入程序；(2)找到初始可行基，确定初始基可行解；(

2、3)检验各非基变單:Xj的检验数Efj]=SCiAij-Cj,若E[jl大于等于零，则以得到最优解，可停止计算；否则进入下一步；(4)根据min(E[jl<0)=E[j],确定Xj为换入变量，按()规则计算()=min(Bi/Aik)=B1/Alk可确定XI为换岀变量，转入下一步；(5)以Aik为主元素进行迭代，把Xk所对应的列向U变换为(O,O...l...,0)t,将Xb列屮的换为Xk,得到新的单纯形法。重复(2)〜(5),直到终止。系统实现(1)计算检验数publicvoidEB(double[,]JD,int[]XS,double

3、t]卜:，int[]B){for(intj=0;j〈c_l;j++){for(inti=0;i

4、;1++){if(E[l]<0){min=1;break;>}returnmin;(4)计算换出变量所对应的行下表publicintMin2(int[]B,doublet,]JD,doublet]E)intL=Mini(E);intmin=0;doubletemp=JD[0,c-1]/JD[0,L];for(int1=0;1

5、E,int[]B,doublet,]JD){intH,L;L=Minl(E);H=Min2(B,JD,E);B[H]=L+l;doublem=JD[H,L];for(intj=0;jfor(inti=II+1;i

6、le.WriteLine("谐输入矩阵:");for(inti=0;i

7、c.EB(JD,XS,E,B);while(dc.MinO(E)〈0)使用说明dc.HB(E,B,JD);E=newdouble[c-1];dc.EB(JD,XS,E,R);}for(inti=0;i〈r;i++){Console.WriteI」ne(〃x{0}={l}",B[i],JD[i,c-1]);Z=Z+XS[B[i]-l]*JD[i,c-1];}Console.WriteLine("其余X变量为零"):Console.WriteLineC目标函数:{0}"，Z):(1)输入增广矩阵fik:///D:/学习/C#/单注形iS/DC

8、XF/DCXF/bin/Debug/DCXF.EXE(2)请输入目标函数系数(3)请输入基变量下表(4)求得结果C:windowssystem32cmd.exe回22三、实验总结或结论1