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

实验课程名称：测控系统仿真

实验项目名称事件步长法仿真实验成绩

实验者雷阳专业班级测控1101班组别

同组者实验日期2014年5月19日

第一部分实验预习报告

实验目的

1、掌握事件步长法的原理及实现

2、学会Matlab程序设计

实验内容

对如下可靠性问题进行仿真:

一设备上三个相同的轴承，每个轴承正常工作寿命为随机变量，其概率分布如表所示．

在任何一个轴承损坏都可以使设备停止工作，从有轴承损坏，设备停止工作，到检修工到达开始更换部件为止，称为一个延迟时间．延迟时间也是随机变量，其概率分布如下表所示．

寿命h 10001100 12001300 14001500 16001700 18001900

概率 0.100.13 0.250.13 0.090.12 0.020.06 0.050.05

延迟时间min 5 10 15

概率 0.6 0.3 0.1

设备停工时每分钟损失5元,检修工每小时工时费12元，轴承每个成本16元．更换一个轴承需要20min，同时更换两个轴承需要30min，同时更换三个轴承需要40min．

现在有两种方案：方案一是损坏一个更换一个；方案二是一旦有轴承损坏就全部更换．试通过计算机仿真对这两种方案做出评价．

第二部分实验过程记录

（实验程序清单，截图，试验中遇到的问题等）

一、程序如下

方案一是损坏一个更换一个：

clear

close

clc

t=0;%仿真时钟,成本一

cost1=0;

%产生初始化事件表

p1=rand(1,3);%产生三个轴承坏的时间概率

p2=rand(1,3);%产生三个延迟时间概率

fori=1:3

ifp1(i)<0.1

ta(i)=1000;

elseifp1(i)<0.23

ta(i)=1100;

elseifp1(i)<0.48

ta(i)=1200;

elseifp1(i)<0.61

ta(i)=1300;

elseifp1(i)<0.7

ta(i)=1400;

elseifp1(i)<0.82

ta(i)=1500;

elseifp1(i)<0.84

ta(i)=1600;

elseifp1(i)<0.9

ta(i)=1700;

elseifp1(i)<0.95

ta(i)=1800;

elseta(i)=1900;

end

ifp2(i)<0.6

tyc(i)=5;

elseifp2(i)<0.9

tyc(i)=10;

elsetyc(i)=15;

end

tsj(i)=ta(i)+tyc(i);

end

t=min(tsj);

forj=1:3

iftsj(j)==t

i=j;

break;

end

end

cost1=cost1+(tyc(i)+20)*5+4+16;

%仿真

whilet<=100000

p1(i)=rand(1);

p2(i)=rand(1);

ifp1(i)<0.1

ta(i)=1000;

elseifp1(i)<0.23

ta(i)=1100;

elseifp1(i)<0.48

ta(i)=1200;

elseifp1(i)<0.61

ta(i)=1300;

elseifp1(i)<0.7

ta(i)=1400;

elseifp1(i)<0.82

ta(i)=1500;

elseifp1(i)<0.84

ta(i)=1600;

elseifp1(i)<0.9

ta(i)=1700;

elseifp1(i)<0.95

ta(i)=1800;

elseta(i)=1900;

end

ifp2(i)<0.6

tyc(i)=5;

elseifp2(i)<0.9

tyc(i)=10;

elsetyc(i)=15;

end

tsj(i)=tsj(i)+ta(i)+tyc(i);

t=min(tsj);

forj=1:3

iftsj(j)==t

i=j;

break;

end

end

cost1=cost1+(tyc(i)+20)*5+4+16;

end

[t,cost1];

ans=

10006035185

方案二是一旦有轴承损坏就全部更换：

clear

close

clc

t=0;%仿真时钟,成本一成本二

cost2=0;

tsj=[000];

whilet<=100000

p1=rand(1,3);%产生三个轴承坏的时间概率

p2=rand(1,3);%产生三个延迟时间概率

fori=1:3

if