《运筹学实验报告-》由会员分享,可在线阅读,更多相关《运筹学实验报告-(27页珍藏版)》请在金锄头文库上搜索。
1、运筹学实验报告运筹学实验报告专 业: 班 级: 姓 名: 学 号: 指导教师: 数学与应用数学专业2015-12-18实 验 目 录一、实验目的3二、实验要求3三、实验内容31、 线性规划32、 整数规划63、 非线性规划134、 动态规划145、 排队论19四、需用仪器设备26五、MATLAB优化工具箱使用方法简介26六、LINGO优化软件简介26七、实验总结27一、实验目的 1、 会利用适当的方法建立相关实际问题的数学模型;2、 会用数学规划思想及方法解决实际问题;3、 会用排队论思想及方法解决实际问题;4、 会用决策论思想及方法解决实际问题;5、 掌握MATLAB、LINGO等数学软件的
2、应用;2、 实验要求1、 七人一组每人至少完成一项实验内容;2、 每组上交一份实验报告;3、 每人进行12分钟实验演示;4、 实验成绩比例: 出 勤:40% 课堂提问:20% 实验报告:30% 实验演示:10%。3、 实验内容1、 线性规划例运筹学74页14题 Min z=-2x1-x2 s.t. 2x1+5x260 x1+x218 3x1+x244 X210 X1,x20用matlab运行后得到以下结果:the program is with the linear programmingPlease input the constraints number of the linear pro
3、gramming m=6m = 6Please input the variant number of the linear programming n=2n = 2Please input cost array of the objective function c(n)_T=-2,-1c = -2 -1Please input the coefficient matrix of the constraints A(m,n)=2,5;1,1;3,1;0,1;-1,0;0,-1A = 2 5 1 1 3 1 0 1 -1 0 0 -1Please input the resource arra
4、y of the program b(m)_T=60,18,44,10,0,0b = 60 18 44 10 0 0Optimization terminated.The optimization solution of the programming is:x = 13.0000 5.0000The optimization value of the programming is:opt_value = -31.0000LINDO程序在命令窗口键入以下内容:max -2x-ysubject to2x+5y=60x+y=183x+y=44y=10end按solve键在reports windo
5、w出现:Global optimal solution found. Objective value: 0.000000 Total solver iterations: 0 Variable Value Reduced Cost X 0.000000 2.000000 Y 0.000000 1.000000 Row Slack or Surplus Dual Price 1 0.000000 1.000000 2 60.00000 0.000000 3 18.00000 0.000000 4 44.00000 0.000000 5 10.00000 0.0000002、 整数规划课本第二章7
6、9页1题 Max z=100x1+180x2+70x3 s.t. 40x1+50x2+60x310000 3 x1+6x2+ 2x3600 x1130 X280 x3200 x1 x2 x30程序运行及结果:biprogramtheprogramiswiththebinarylinearprogrammingPleaseinputtheconstraintsnumberoftheprogrammingm=5m=5Pleaseinputthevariantnumberoftheprogrammingn=5n=5Pleaseinputcostarrayoftheobjectivefunctionc
7、(n)_T=100,180,70c=10018070PleaseinputthecoefficientmatrixoftheconstraintsA(m,n)=40,50,60;3,6,2;1,0,0;0,1,0;0,0,1A=405060362100010001Pleaseinputtheresourcearrayoftheprogramb(m)_T=10000;600;130;80;200b=1000060013080200Optimizationterminated.Theoptimizationsolutionoftheprogrammingis:x=000Theoptimizatio
8、nvalueoftheprogrammingis:opt_value= 0程序名:intprogram b程序说明:%theprogrammiswiththeintegerlinearprogrammingusebranchandboundmethod!%这个程序是用分支定界法解决整数规划问题%pleaseinputtheparametersinthemainfunctioninthecommandwinows%请在命令窗口输入这个主要定义函数的参数functionx,f=ILp(c,A,b,vlb,vub,x0,neqcstr,pre)%minf=c*x,s.t.A*x=b,vlb=x=vu
9、b%f的最小值等于c的转置乘以x,A乘以x小于等于b,x大于等于vlb小于等于vub%thevectorsofxisrequiredasintegersaswhole%x是整个的整数需要%x0istheinitialization,isalsook%x0是初始值,也可以是。%neqcstristhenumberofequationalconstraints,when0canbedelete%neqcstr是平均约束条件的数目,当0能删除时%preistheconciserate%pre是简明率%xistheintegeroptimizationandfistheoptimalvalue%x是整
10、数规划,f是最优值%ifnargin8,pre=0;%narginisthefactuallyinputvariantsnumber(这个参数是实际输入的变量个数)ifnargin7,neqcstr=0;ifnargin6,x0=;ifnargin5,vub=;ifnargin4,vlb=;endendendendend%settocolumnvectors%建立列向量x0=x0(:);c=c(:);b=b(:);vlb=vlb(:);vub=vub(:);mm=1;j=1;nvars=length(c);%numberofvariants(变量的个数)fvub=inf;xall=;fall=;x_f_b=;xtemp,ztemp,how=lp(c,A,b,vlb,vub,x0,neqcstr,-1);ftemp=c*xtemp;%ifstrcmp(how,ok)%comparebetweenhowandok(how和ok之间的比较)temp0=round(xtemp);%临时变量四舍五入temp1=floor(xtemp);%取比其小的整数temp