管理实验三(答案).docVIP

鲁东大学信息与电气工程学院 管理决策模型与方法实验报告（三） 实验名称： 计算机模拟 专 业： 信息管理与信息系统 班 级： 2012级信管02班 姓 名： 邬阳阳 学 号： 20122213199 一、实验目的：学习计算机模拟的基本过程与方法。 二、实验要求：将相关程序（M文件）和运行结果粘贴到word文档中，说明操作过程中出现的问题及改正方法，以及后续学习应注意的问题。 三、实验内容： 1、产生相应的模拟随机数 （1）产生5×5阶［3,6］上均匀分布U（3,6）的随机数矩阵 unifrnd (3,6,5,5) （2）产生2×3阶［0,1］均匀分布的随机数矩阵 rand(2,3) （3）产生5×3阶均值为10，方差为25的正态分布的随机数矩阵 normrnd(10,25,5,3) （4）产生2×7阶期望值为0.25的指数分布的随机数矩阵 参数 exprnd(4,2,7) （5）产生4×3阶参数为5的泊松分布的随机数矩阵 poissrnd(5,4,3) 2、利用蒙特卡罗法求解 在MATLAB软件包中编程，共需3个Ｍ文件:randll.m, ff.m, st.m.主程序为randll.m. ff.m:建立目标函数 function y=ff(x) y=-x(1)-2*x(2)+0.5*x(1).^2+0.5*x(2).^2; end st.m:约束条件 function lpc=st(x) if 2*x(1)+3*x(2)-6=0 2*x(1)+3*x(2)-5=0 lpc=1; else lpc=0; end randll.m:主程序 function [sol,r1,r2]=randll(a,b,n) debug=1; a=0;%试验点下界 b=5;%试验点上界 n=10000;%试验点个数 r1=unifrnd(a,b,n,1); r2=unifrnd(a,b,n,1); sol=[r1(1) r2(1)]; z0=inf; for i=1:n x1=r1(i); x2=r2(i); lpc=st([x1 x2]); if lpc==1 z=ff([x1 x2]); if zz0 z0=z; sol=[x1 x2]; end end end z=z0 3、 某设备上安装有4只型号规格完全相同的电子管，已知电子管寿命服从1000～2000h之间的均匀分布．电子管损坏时有两种维修方案，一是每次更换损坏的那只；二是当其中1只损坏时4只同时更换．已知更换时间为换1只时需1h，4只同时换为2h．更换时机器因停止运转每小时的损失为20元，又每只电子管价格10元，试用模拟方法确定哪一个方案经济合理？ 方案一 money=[0,0,0,0]; s=0; for i=1:4 %四只电子管 r=unifrnd(1000,2000);% 电子管寿命服从1000～2000h之间的均匀分布 t=r; while t=10000 %模拟的电子管的总工作时间10000小时 money(i)=money(i)+30;%损失包括停止运转损失和电子管的价格 r=unifrnd(1000,2000);%维修后电子管的寿命 t=t+r;%维修前及维修后电子管可以工作的总时间 end end s=sum(money) %总的损失 方案二： money=0; r=unifrnd(1000,2000);% 电子管寿命服从1000～2000h之间的均匀分布 t=r; while t=10000 %模拟的电子管的总工作时间10000小时 money(i)=money(i)+200;%损失包括四只电子管停止运转损失和电子管的价格 r=unifrnd(1000,2000);%维修后电子管的寿命 t=t+r;%维修前及维修后电子管可以工作的总时间 end money 根据结果可以比较出两种方案。