系统工程蒙特卡罗

《系统工程蒙特卡罗》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、系统工程用蒙特卡罗法对下面报童问题进行仿真。对新闻的重要度分布概率进行随机数编码。随机数采用两位数字（00~99），如下表1234重要程度概率累积概率随机编码重要0.350.3500~34一般0.450.8035~79平淡0.201.0089~99对报纸的销量分布概率进行随机数编码。随机数采用两位数字（00~99），如下表报纸需求量累积分布重大新闻随机编码一般新闻随机编码平淡新闻随机编码40.000.0300~020.1000~090.4400~4350.000.0803~070.2810~270.6644~6560.000.230

2、8~220.6828~670.8266~8170.000.4323~420.8868~870.9482~9380.000.7843~770.9688~951.0094~9990.000.9378~921.0096~991.00100.001.0093~991.001.00用excel随机产生100个数，模拟次数重要度报纸需求随机数程度随机数报纸需求量15重要8490241一般8470373一般3360453一般2250529重要14069重要7180767一般740869一般7370956一般45601093平淡86701161一般

3、2550123重要3501334重要47801455一般47601515重要17601651一般63601798平淡10401810重要14601952一般0402077一般91802125重要15602271一般70702388平淡20402457一般2402558一般94802666一般71702775一般18502828重要13602933重要61803019重要5580316重要68803273一般72703397平淡97803423重要68803561一般76703668一般81703715重要84903866一般6660

4、3958一般85704098平淡90704162一般74704224重要24704327重要48804462一般87704526重要2770464重要7504758一般42604813重要931004989平淡8405024重要53805120重要84905235一般70705348一般18505491平淡84705580平淡32405641一般13505740一般5405815重要0405988平淡85706048一般23506133重要18606288平淡9406324重要92906420重要42706570一般58606660

5、一般92806753一般1406875一般57606954一般94807064一般91807120重要3507294平淡6407393平淡63507441一般3660752重要9290764重要86907784平淡57507842一般35607912重要6780807重要54808133重要36708222重要51808396平淡14408454一般63608532重要8608616重要51808791平淡93708823重要45808917重要971009055一般6409152一般24509249一般12509318重要4680

6、9480平淡96809570一般24509670一般65609751一般77709820重要2770995重要377010027重要5180根据上述的结果得到报纸的需求量，如下表报纸需求量的个数个数概率值40150.1550140.1460180.1870230.2380220.229060.0610020.02总数100第二步计算报纸过多的概率和平均需求量。当需求点A=40，则不可能发生报纸多，过多的概率为0当需求点A=50，报纸需求量为40的情况下，就会发生多货，多货的概率从上表中可知为P（DDLT<50）=P（DDLT=40）

7、=0.15当需求点A=60，报纸需求量为40.50的情况下，就会发生多货，多货的概率从上表中可知为P（DDLT<60）=P（DDLT=40）+P（DDLT=50）=0.29当需求点A=70，报纸需求量为40.50.60的情况下，就会发生多货，、多货的概率从上表中可知为P（DDLT<70）=P（DDLT=40）+P（DDLT=50）+P（DDLT=60）=0.47当需求点A=80，报纸需求量为40.50.60.70的情况下，就会发生多货，多货的概率从上表中可知为P（DDLT<60）=P（DDLT=40）+P（DDLT=50）+P（DD

8、LT=60）+P（DDLT=70）=0.70当需求点A=90，报纸需求量为40.50.60.70.80的情况下，就会发生多货，多货的概率从上表中可知为P（DDLT<0）=P（DDLT=40）+P（DDLT=50）+P（DDLT=60）