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

实验课程名称：测控系统仿真 实验项目名称 事件步长法仿真 实验成绩实验者雷阳专业班级测控1101班 组 别同组者 实验日期 2014年5月19 日第一部分 实验预习报告实验目的1、掌握事件步长法的原理及实现2、学会Matlab程序设计 实验内容对如下可靠性问题进行仿真: 一设备上三个相同的轴承，每个轴承正常工作寿命为随机变量，其概率分布如表所示在任何一个轴承损坏都可以使设备停止工作，从有轴承损坏，设备停止工作，到检修工到达开始更换部件为止，称为一个延迟时间延迟时间也是随机变量，其概率分布如下表所示寿命h1000 11001200 13001400 15001600 17001800 1900概率0.10 0.130.25 0.130.09 0.120.02 0.060.05 0.05延迟时间min5 10 15概率 0.6 0.3 0.1 设备停工时每分钟损失5元,检修工每小时工时费12元，轴承每个成本 16元更换一个轴承需要 20 min，同时更换两个轴承需要30min，同时更换三个轴承需要40min现在有两种方案：方案一是损坏一个更换一个；方案二是一旦有轴承损坏就全部更换试通过计算机仿真对这两种方案做出评价第二部分 实验过程记录（实验程序清单，截图，试验中遇到的问题等）一、程序如下方案一是损坏一个更换一个：clearcloseclct=0;%仿真时钟,成本一cost1=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.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(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);for j=1:3if tsj(j)=ti=j;break;endendcost1=cost1+(tyc(i)+20)*5+4+16;endt,cost1;ans = 100060 35185方案二是一旦有轴承损坏就全部更换：clearcloseclct=0;%仿真时钟,成本一 成本二cost2=0;tsj=0 0 0;while t=100000p1=rand(1,3);%产生三个轴承坏的时间概率p2=rand(1,3);%产生三个延迟时间概率 for i=1:3 if p1(i)0.1 ta(i)=1000; elseif p1(i)0.23 ta(i)=1100; elseif p1(i)0.48 ta(i)=1200; elseif p1(i)0.61 ta(i)=1300; elseif p1(i)0.7 ta(i)=1400; elseif p1(i)0.82 ta(i)=1500; elseif p1(i)0.84 ta(i)=1600; elseif p1(i)0.9 ta(i)=1700; elseif p1(i)0.95 ta(i)=1800; else ta(i)=1900; end if p2(i)0.6 tyc(i)=5; elseif p2(i)0.9 tyc(i)=10; else tyc(i)=15; end tsj(i)=tsj(i)+ta(i)+tyc(i); end t=min(tsj); for j=1:3 if tsj(j)=t i=j; break; end end cost2=cost2+(tyc(i)+40)*5+8+3*16;endt,cost2ans = 100185 22537第三部分 结果与讨论1、 结果分析： 由仿真结果可知