测控系统仿真实验二报告书模板

《测控系统仿真实验二报告书模板》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、实验课程名称：测控系统仿真实验项目名称事件步长法仿真　　　　　　　实验成绩　　　　　　实　验　者　　　雷阳　　专业班级　　测控1101班　　组别　　　同　组　者　　　　　　　　　　　　　　　　　　　　实验日期2014年5月19日第一部分实验预习报告　实验目的1、掌握事件步长法的原理及实现2、学会Matlab程序设计实验内容对如下可靠性问题进行仿真:一设备上三个相同的轴承，每个轴承正常工作寿命为随机变量，其概率分布如表所示．在任何一个轴承损坏都可以使设备停止工作，从有轴承损坏，设备停止工作，到检修工到达开始更换部件为止，称为一个延迟时间．延迟时间也是随机变量，其概率分布如下表所示．寿命

2、h1000110012001300140015001600170018001900概率0.100.130.250.130.090.120.020.060.050.05延迟时间min51015概率0.60.30.1设备停工时每分钟损失5元,检修工每小时工时费12元，轴承每个成本16元．更换一个轴承需要20min，同时更换两个轴承需要30min，同时更换三个轴承需要40min．现在有两种方案：方案一是损坏一个更换一个；方案二是一旦有轴承损坏就全部更换．试通过计算机仿真对这两种方案做出评价．第二部分实验过程记录（实验程序清单，截图，试验中遇到的问题等）一、程序如下方案一是损坏一个更换一个：

3、clearcloseclct=0;%仿真时钟,成本一cost1=0;%产生初始化事件表p1=rand(1,3);%产生三个轴承坏的时间概率p2=rand(1,3);%产生三个延迟时间概率fori=1:3ifp1(i)<0.1ta(i)=1000;elseifp1(i)<0.23ta(i)=1100;elseifp1(i)<0.48ta(i)=1200;elseifp1(i)<0.61ta(i)=1300;elseifp1(i)<0.7ta(i)=1400;elseifp1(i)<0.82ta(i)=1500;elseifp1(i)<0.84ta(i)=1600;elseifp1(i)

4、<0.9ta(i)=1700;elseifp1(i)<0.95ta(i)=1800;elseta(i)=1900;endifp2(i)<0.6tyc(i)=5;elseifp2(i)<0.9tyc(i)=10;elsetyc(i)=15;endtsj(i)=ta(i)+tyc(i);endt=min(tsj);forj=1:3iftsj(j)==ti=j;break;endendcost1=cost1+(tyc(i)+20)*5+4+16;%仿真whilet<=100000p1(i)=rand(1);p2(i)=rand(1);ifp1(i)<0.1ta(i)=1000;elseif

5、p1(i)<0.23ta(i)=1100;elseifp1(i)<0.48ta(i)=1200;elseifp1(i)<0.61ta(i)=1300;elseifp1(i)<0.7ta(i)=1400;elseifp1(i)<0.82ta(i)=1500;elseifp1(i)<0.84ta(i)=1600;elseifp1(i)<0.9ta(i)=1700;elseifp1(i)<0.95ta(i)=1800;elseta(i)=1900;endifp2(i)<0.6tyc(i)=5;elseifp2(i)<0.9tyc(i)=10;elsetyc(i)=15;endtsj(i)=

6、tsj(i)+ta(i)+tyc(i);t=min(tsj);forj=1:3iftsj(j)==ti=j;break;endendcost1=cost1+(tyc(i)+20)*5+4+16;end[t,cost1];ans=10006035185方案二是一旦有轴承损坏就全部更换：clearcloseclct=0;%仿真时钟,成本一成本二cost2=0;tsj=[000];whilet<=100000p1=rand(1,3);%产生三个轴承坏的时间概率p2=rand(1,3);%产生三个延迟时间概率fori=1:3ifp1(i)<0.1ta(i)=1000;elseifp1(i)<

7、0.23ta(i)=1100;elseifp1(i)<0.48ta(i)=1200;elseifp1(i)<0.61ta(i)=1300;elseifp1(i)<0.7ta(i)=1400;elseifp1(i)<0.82ta(i)=1500;elseifp1(i)<0.84ta(i)=1600;elseifp1(i)<0.9ta(i)=1700;elseifp1(i)<0.95ta(i)=1800;elseta(i)=1900;endifp2(i)<0.