非齐次泊松过程课程设计.doc

课程名称：《随机过程》课程设计（论文）

题目: 非齐次泊松过程

在数控机床可靠

性建模中的应用

学院：理学院

专业：数学与应用数学

班级：数学12-1班

***名：***

学生学号： ********** ***师：***

2015 年 1月 3 日

随机过程课程设计

目录

任务书 (1)

摘要 (1)

前言 (2)

１非齐次泊松过程理论 (2)

1.1 非齐次泊松过程的基本理论简介 (2)

1.2 基于试验总时间法的趋势检验 (2)

2 数控机床的非齐次泊松过程可靠性建模 (3)

2.1强度函数的建立................................................. . (3)

2.2 Ｋ台数控机床强度函数的参数估计......................... (4)

2.3 非齐次泊松过程下的可靠性指标............................... ....... (5)

3实例分析 (5)

4结束语 (7)

5程序及结果 (8)

6参考文献 (9)

附录………………………………………………………………………………

评阅书……………………………………………………………………………

摘要

基于试验总时间法对多样本随机截尾的数控机床现场数据进行趋势检验，在故障过程为浴盆曲线的趋势条件下构建了数控机床的非齐次泊松过程的可靠性模型。本文使用极大似然估计法对非齐次泊松过程的强度函数进行参数估计得到了该模型的可靠性指标，以6台加工中心的现场数据为例建立了非齐次泊松过程的可靠性模型。再通过matlab曲线拟合，绘制出故障时间的曲线，通过曲线的拟合程度，可以确定非齐次泊松过程能够更恰当地表现故障的趋势。

关键词：数控机床可靠性非齐次泊松过程浴盆曲线

前言

数控机床是由数目众多的零部件组成的复杂机电液可修系统。在其可靠性研究中，需要考虑维修活动对其可靠度的影响【1】。以往的数控机床可靠性建模方法，是将故障间隔时间视为独立同分布来分析其寿命分布，即假设维修活动是“修复

如新”【2】

而在实际生产中数控机床的维修活动是以调整或者更换一部分零部件和

元器件为主的，对于复杂的系统来说这种维修活动只能使产品恢复到正常功能维修前后可靠度并没有很大改变，因此将数控机床的维修活动视为“修复如旧”更加合理。非齐次泊松过程经常被用于建立“修复如旧”的维修策略且维修时间可忽略的可修系统可靠性模型用于模拟出现故障间隔时间的趋势【3-4】。结合数控机床的维修特点，使用非齐次泊松过程建立的可靠性模型更能贴近于复杂系统的生产实际。

本文提出了非齐次泊松过程的数控机床可靠性建模方法，并结合数控机床的失效特点，建立故障率为浴盆曲线的非齐次泊松过程可靠性模型。同时，结合具有随机截尾特点的多样本数控机床现场试验故障数据，对数控机床的可靠性进行了深入分析。

１．非齐次泊松过程理论

1.1 非齐次泊松过程的基本理论简介：

非齐次泊松过程是随机点过程的一种典型类型，当可修系统的相邻故障间隔呈现某种趋势时可以使用这种方法来描述。 非齐次泊松过程的重要参数【4】：)(t ω为强度函数，是非负函数；其 累 积 故 障 强 度 函 数⎰=t

du u t W 0)()(ω，表示在[0,t]中的平均故障数,即E[N(t)]=W

（t ），N[t] 代表在[0,t]出现的故障次数ｔ表示机床从观测开始后的运行时间。

当强度函数为1)()(-=βα

αβωt

t 时，成为威布尔过程。其中，α、β>0，α为

尺度参数，β为形状参数。0<β<1, 表示不断改良的(好)系统; β>1, 表示不断恶化的(坏)系统；β=1，表示系统服从指数分布。

1.２基于试验总时间法的趋势检验：

本文采用基于试验总时间的方法，对具有多样本随机截尾现场数据的故障过程进行趋势检验【5】。

将在观测期间采集到的所有故障数据按照从大到小时间进行排序，得到t (i)

的时间序列。根据试验总时间的建模思想【6】，得到该序列的第ｉ个故障发生时的试验总时间：

⎰=)

(0)()()(i t i du u n t T （1）

=N ∑=k

i i n 1

（2）

式中：n(u)表示在ｕ时刻观察到的数控机床数量，N 表示在观测期间的故障数；i n 表示第ｉ台机床的故障数（共有Ｋ台机床），当时间序列)(i t 的最后一个时间是故障数据时，)(i n =i n -1；当时间序列)(i t 的最后一个时间不是故障数据而是截尾数据时)(i n =i n 。

在实际检验时同时使用Ｕ－检验、Ｊ－检验和Ｖ－检验等检验方法综合确定

有无趋势【7-8】。其中Ｖ－检验如下：

Ｈ0：齐次泊松过程； Ｈ1：具有非单调趋势；

()1

1()()

|()|24

~(0,1)()

48

N

i i T S T S T t N V N N T S =-

-⨯=

∑ (3)

22()1

22()()

|()|212()

180

N

i i T S T S T t N V N

T S =--⨯=

∑ ~(0,1)N (4)

()1

32

2log(()|2()()|)

(2)

N

i i T S T t T S V N χ=-=

∑~(0,1)N (5)

式中：Ｔ（Ｓ）表示总的观测时间。当｜Ｖ１｜＜ｚα／２，｜Ｖ２｜＜ｚα／２，Ｖ３＜χ２（２ｎ）时，接受Ｈ０。

一般情况下，Ｕ检验和Ｊ检验是检验具有单调趋势或齐次泊松过程和更新过程的故障数据的，而当故障数据具有非单调趋势时，则可以考虑Ｖ检验中的统计量，如表１所示。

表1 故障率和故障强度函数变化特性

２ 数控机床的非齐次泊松过程可靠性建模 2.1强度函数的建立

对于浴盆曲线趋势的故障过程，假设其故障强度函数由早期故障期和偶然故障期两部分组成，并且每一个阶段都是一个威布尔过程［９］，参数为尺度参数αｍ和形状参数βｍ（ｍ＝１，２）。结合以上假设和多重威布尔分布模型的性质，则