对如下可靠性题进行仿真

1、文档来源为 :从网络收集整理.word 版本可编辑.欢迎下载支持对如下可靠性问题进行仿真：一设备上三个相同的轴承，每个轴承正常工作寿命为随机变量，其概率分概率分布如下表所示.延时时间min51015概率0.60.30.1设备停工时每分钟损失5元，检修工每小时工时费12元，轴承每个成本16元.更 换一个轴承需要20 min ,同时更换两个轴承需要 30min,同时更换三个轴承需 要 40min.现在有两种方案：方案一是损坏一个更换一个；方案二是一旦有轴承损坏就全部 更换.试通过计算机仿真对这两种方案做出评价.方案1%事件步长法clearcloseclct=0;%仿真时钟，成本一成本二cost1=

2、0;cost2=0;%产生初始化事件表p1=rand(1,3);%产生三个轴承坏的时间概率p2=rand(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 p1(i)0.7ta(i)=1400;elseif p1(i)0.82ta(i)=1500;elseif p1(i)0.84ta(i)=1600;elseif p1(i)0.9ta(i)=1700;elseif p1(i)0

3、.95ta(i)=1800;else ta(i)=1900;endif p2(i)0.6tyc(i)=5;elseif p2(i)0.9tyc(i)=10;else tyc(i)=15;endtsj(i)=ta(i)+tyc(i);endt=min(tsj);for j=1:3if tsj(j)=ti=j;break;endendcost1=cost1+(tyc(i)+20)*5+4+16;%仿真while t=100000p1(i)=rand(1);p2(i)=rand(1);if p1(i)0.1ta(i)=1000;elseif p1(i)0.23ta(i)=1100;elseif p1

4、(i)0.48ta(i)=1200;elseif p1(i)0.61ta(i)=1300;elseif p1(i)0.7ta(i)=1400;elseif p1(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;endif p2(i)0.6tyc(i)=5;elseif p2(i)0.9tyc(i)=10;else tyc(i)=15;endtsj(i)=tsj(i)+ta(i)+tyc(i);t=min(tsj);fo

5、r j=1:3if tsj(j)=ti=j;break;endendcost1=cost1+(tyc(i)+20)*5+4+16;endt,cost1ans =10053034450方案 2%事件步长法clearcloseclct=0;%仿真时钟，成本一成本二cost2=0;tsj=0 0 0;while t=100000p1=rand(1,3);%产生三个轴承坏的时间概率 p2=rand(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;el

6、seif p1(i)0.61ta(i)=1300;elseif p1(i)0.7ta(i)=1400;elseif p1(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;endif p2(i)0.6tyc(i)=5;elseif p2(i)0.9tyc(i)=10;else tyc(i)=15;endtsj(i)=tsj(i)+ta(i)+tyc(i);endt=min(tsj);for j=1:3if tsj(j)=ti=j;break;endendcost2=cost2+(tyc(i)+40)*5+8+3*16;endt,cost2ans =10108022531第二种方案好实验四 蒙特卡洛法仿真实验实验目的了解蒙特卡洛法的原理学会运用matlab 编制m 文件学会运用matlab 进行蒙特卡法编程实验内容运用蒙特卡洛法求如下积分,并进行验证实验程序a=1;b=5;result=0;for j=0:1000if(j1000)gx=0;for i=0:1000i