《对策论对策论.doc》由会员分享，可在线阅读，更多相关《对策论对策论.doc（7页珍藏版）》请在三一文库上搜索。
1、安徽师范大学数计 学院实验报告专业名称数学与应用数学实验室2号实验楼#201实验课程 Matlab实验名称 对策论姓名 周冈0学号 100701193同组人员无实验日期 2013.4.10 注：实验报告应包含（实验目的，实验原理，主要仪器设备和材料，实验过程 和步骤，实验原始数据记录和处理，实验结果和分析，成绩评定）等七项内容。 具体内容可根据专业特点和实验性质略作调整，页面不够可附页。一. 实验目的1、了解对策论建模的方法和模型的算法；2 、了解带线性规划的基本原理和解法；3、掌握Matlab优化工具箱求解线性规划的基本用法；二. 实验原理1.主要命令和注意事项线性规划模型：min z=cX
2、s.tAX £ bAeq ?X beqVLBC X< VUB在Matlab中可通过linprog函数来实现，其调用形式为：1 x=linprog（c，A，b，Aeq,beq, VLB，VUB2 x=linprog（c，A，b，Aeq,beq, VLB，VUB, X。）3 x,fval,exitflag,output,lambda=linprog（ c， A， b，Aeq,beq, VLB，VUB, X。options ）返回最优解X及X处的目标函数值 fval.注意：1若没有等式约束：Aeq?X beq,则令Aeq= , beq=.2 其中X0表示初始点三. 主要仪器设备和材料
3、计算机，matlab软件四. 实验原始数据记录和处理（1）求解线性规划min z = - 5X| + 4x2 + 2x3s.t 6 x1 - x2 + x3 810x-i + 2 x2 + 4x3-1 # Xi 3 0 # X2 2 x3 3 0.c=-5 4 2;A=6 -1 1;1 2 4;b=8;10;Aeq=; beq=;vlb=-1;0;0; vub=3;2;i nf;x,fval=li nprog(c,A,b,Aeq,beq,vlb,vub)Optimization terminated.x =1.33330.00000.0000fval =-6.6667(2)求解线性规划maxz
4、 = 6xi + 4x2 s.t2 x-i + 3x21004x-i + 2x2120Xj, x2 3 0解：c=-6 -4;A=2 3;4 2;b=100;120;Aeq=;beq=; vlb=O;O;vub=; x,g=li nprog(c,A,b,Aeq,beq,vlb,vub); x g=-gOptimization terminated.x =20.000020.0000g =200.0000O( 3)P94 ex6假设每吨混合药剂获得利益为单位1且让资源充分利用，即能生产400吨产品，总共利益为400.对成品的利润的分配，按各厂产品与合成产品的量来计 算。解：主程序：k=1;i=1
5、;phi=zeros(1,3);while k<4p=0;for l=0:1for m=0:1x=1 l m;if k=2z=x(1);x(1)=x(2);x(2)=z;elseif k=3z=x(1);x(1)=x (3);x(3)=z;ends=sum(x); y=x; y(k)=0;vs=tezhe ngha nshu(x)-tezhe ngha nshu(y);vs=vs*jieche ng(s-1)*jieche ng(3-s)/jieche ng(3);p=p+vs;endendPhi(i)=p；i=i+1;k=k+1;endphi调用特征函数程序：function v=te
6、zhenghanshu(x)s=sum(x);if s=0|s=1 v=-400; retur n; endif x=1,1,0 v=-400; return; endif x=1,0,1 v=-40; return; endif x=0,1,1 v=-120; return; endif s=3 v=0; end调用阶乘函数：fun cti on k=jieche ng(n)k=1;for i=1: nk=k*i;end求得结果为phi =112.5000 87.5000 200五. 实验结果和分析通过该实验，我们可以得到对策论可以帮助我们很好的解决一些问题。而且得到最佳问题的解决方案.8成绩评定:1、根据实验情况和实验报告质量作出写事性评价2、评分综合评分 折合成等级 指导教师签名：时间：年月日