欢迎来到天天文库
浏览记录
ID:28713778
大小:32.00 KB
页数:6页
时间:2018-12-13
《测控系统仿真实验二报告材料书实用实用模板》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
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(
4、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)=ta(i)+tyc(i);endt=min(tsj)
5、;for精彩文档实用标准文案j=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;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
6、.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)=tsj(i)+ta(i)+tyc(i);t=min(tsj);forj=1:3if精彩文档实用标准文案tsj(j)==ti=j;break;endendcost1=c
7、ost1+(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)<0.23ta(i)=1100;elseifp1(i)<0.48ta(i)=1200;
8、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(
此文档下载收益归作者所有