对如下可靠性题进行仿真

1、百度文库对如下可靠性问题进行仿真：一设备上三个相同的轴承，每个轴承正常工作寿命为随机变量，其概率分 布如表所示寿命h10001100120013001400M5001600170018001900概率0.100.130.250.130.090.120.020.060.050.05在任何一个轴承损坏都可以使设备停止工作， 从有轴承损坏，设备停止工作，到 检修工到达开始更换部件为止，称为一个延迟时间.延迟时间也是随机变量，其 概率分布如下表所示.延时时间min51015概率0.60.30.1设备停工时每分钟损失5元,检修工每小时工时费12元，轴承每个成本16元更 换一个轴承需要20 min，同时更

2、换两个轴承需要30min，同时更换三个轴承需要 40min.现在有两种方案：方案一是损坏一个更换一个；方案二是一旦有轴承损坏就全部 更换试通过计算机仿真对这两种方案做出评价.方案1%事件步长法clearcloseclct=0;%仿真时钟，成本一成本二cost1=0; cost2=0;%产生初始化事件表p仁ran d(1,3);%产生三个轴承坏的时间概率 p2=ra nd(1,3);%产生三个延迟时间概率for i=1:3if p1(i)<0.1ta(i)=1000;elseif p1(i)<0.23 ta(i)=1100;elseif p1(i)<0.48 ta(i)=120

3、0;elseif p1(i)<0.61ta(i)=1300;elseif p1(i)<0.7ta(i)=1400;elseif p1(i)<0.82 ta(i)=1500;elseif p1(i)<0.84ta(i)=1600;elseif p1(i)<0.9ta(i)=1700;elseif p1(i)<0.95ta(i)=1800;else ta(i)=1900; end if p2(i)<0.6/tyc(i)=5;/ elseif p2(i)<0.9tyc(i)=10;else tyc(i)=15; end tsj(i)=ta(i)+tyc

4、(i);endt=mi n(tsj);for j=1:3if tsj(j)=ti=j; break;endendcost 仁 cost1+(tyc(i)+20)*5+4+16;%仿真while t<=100000p1(i)=ra nd(1);p2(i)=ra nd(1);if p1(i)<0.1ta(i)=1000; elseif p1(i)<0.23 ta(i)=1100;elseif p1(i)<0.48ta(i)=1200;elseif p1(i)<0.61ta(i)=1300;elseif p1(i)<0.7ta(i)=1400;elseif p1(

5、i)<0.82ta(i)=1500;elseif p1(i)<0.84ta(i)=1600;elseif p1(i)<0.9ta(i)=1700;elseif p1(i)<0.95ta(i)=1800;else ta(i)=1900; end if p2(i)<0.6 tyc(i)=5;elseif p2(i)<0.9tyc(i)=10;else tyc(i)=15; end/(i)=tsj(i)+ta(i)+tyc(i);t=mi n(tsj);for j=1:3if tsj(j)=ti=j;break;endendcost 仁 cost1+(tyc(i)

6、+20)*5+4+16;endt,cost1ans =10053034450方案2%事件步长法clearcloseclct=0;%仿真时钟，成本一成本二cost2=0;tsj=0 0 0;while t<=100000p仁ran d(1,3);%产生三个轴承坏的时间概率 p2=ra nd(1,3);%产生三个延迟时间概率 for i=1:3if p1(i)<0.1ta(i)=1000;elseif p1(i)<0.23ta(i)=1100;elseif p1(i)<0.48ta(i)=1200;elseif p1(i)<0.61ta(i)=1300;elseif

7、p1(i)<0.7ta(i)=1400;elseif p1(i)<0.82ta(i)=1500;elseif p1(i)<0.84ta(i)=1600;elseif p1(i)<0.9/ta(i)=1700;/ elseif p1(i)<0.95ta(i)=1800;else ta(i)=1900; end if p2(i)<0.6 tyc(i)=5;elseif p2(i)<0.9tyc(i)=10;else tyc(i)=15; end(i)=tsj(i)+ta(i)+tyc(i);endt=mi n(tsj);for j=1:3if tsj(j)=ti=j;break;endend cost2=cost2+(tyc(i)+40)*5+8+3*16; endt,cost2ans =10108022531第二种方案好6实验四蒙特卡洛法仿真实验实验目的了解蒙特卡洛法的原理学会运用Matlab编制M文件学会运用Matlab进行蒙特卡法编程实验内容运用蒙特卡洛法求如下积分，并进行验证5 14八J 1 ( x x)dx实验程序a=1;b=5;result=0;for j=0:1000if(j<1000)gx=0;fo