CN107220216A - 一种利用特征数的威布尔型备件需求量的近似计算方法 - Google Patents

一种利用特征数的威布尔型备件需求量的近似计算方法 Download PDF

Info

Publication number
CN107220216A
CN107220216A CN201710343600.1A CN201710343600A CN107220216A CN 107220216 A CN107220216 A CN 107220216A CN 201710343600 A CN201710343600 A CN 201710343600A CN 107220216 A CN107220216 A CN 107220216A
Authority
CN
China
Prior art keywords
mrow
mfrac
distribution
gamma
msup
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710343600.1A
Other languages
English (en)
Other versions
CN107220216B (zh
Inventor
邵松世
翟亚利
李华
王俊龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Naval University of Engineering PLA
Original Assignee
Naval University of Engineering PLA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Naval University of Engineering PLA filed Critical Naval University of Engineering PLA
Priority to CN201710343600.1A priority Critical patent/CN107220216B/zh
Publication of CN107220216A publication Critical patent/CN107220216A/zh
Application granted granted Critical
Publication of CN107220216B publication Critical patent/CN107220216B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Algebra (AREA)
  • Evolutionary Biology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开了一种利用特征数的威布尔型备件需求量的近似计算方法，该近似计算方法主要包括如下步骤：(1)利用所述威布尔分布参数α、b计算伽玛分布的参数αg、λ，计算正态分布的参数μ、σ；(2)计算偏度和峰度，依据所述伽玛分布的参数αg、λ及所述正态分布的参数μ、σ，按下式计算三种分布的所述特征数中的偏度和峰度；(3)比较所述特征数中的偏度和峰度，若所得伽玛分布与威布尔分布的偏度绝对差值较小，则按照一种方法计算备件保障概率；否则，按照另外的方法计算备件保障概率Ps。按照本发明实现的备件需求量的近似计算方法，能够简化计算过程以及提高近似计算的精度。

Description

一种利用特征数的威布尔型备件需求量的近似计算方法
技术领域
本发明属于备件需求量计算领域，特别是涉及一种利用特征数的威布尔型备件需求量的近似计算方法。
背景技术
威布尔分布常用来描述因逐渐老化导致故障的元器件寿命，具有这种威布尔分布的元器件为威布尔型单元。威布尔型单元主要适用于机电件，如：滚珠轴承、继电器、开关、断路器、某些电容器、电子管、磁控管、电位计、陀螺、电动机、航空发电机、蓄电池、液压泵、空气涡轮发动机、齿轮、活门、材料疲劳件等。
在上述单元使用于各类系统中时，需要对其备件的需求量进行预先评估计算，备件是在考虑备件寿命的情况下保障装备可持续工作的物质条件，在理论上，备件需求量计算涉及多重卷积。由于威布尔分布的多重卷积形式极为复杂，以致难以获得其多重卷积的数值积分结果。因此，在工程上，一般都采用近似方法来计算威布尔型备件需求量(例如指数近似、正态近似)，但目前在工程中使用的近似方法误差较大，其中指数近似是在威布尔形状参数接近于1的时候计算效果好，正态近似只在威布尔形状参数大于3，并且还要保证近似计算方法合理才能达到较好的计算效果，上述计算方式不仅计算过程复杂，并且不能有效覆盖形状参数的可能取值范围的所有情况，使得不能执行有效的备件需求量的计算。
发明内容
针对现有技术的以上缺陷或改进需求，本发明提供了一种利用特征数的威布尔型备件需求量的近似计算方法，利用威布尔型单元的寿命分布参数来计算伽玛分布和正态分布各自的参数，再分别计算这三种分布各自的偏度和峰度，从中选择和威布尔分布的偏度和峰度更接近的一种分布(伽玛或正态)，用于近似描述该威布尔型单元寿命，并以此计算备件需求量。
为实现上述目的，按照本发明，提供一种利用特征数的威布尔型备件需求量的近似计算方法，所述威布尔型备件的寿命服从威布尔分布W(α,b)，α、b为威布尔分布参数，α为尺度参数，b为形状参数；所述特征数为均值、方差、偏度和峰度，其特征在于，该计算方法包括如下步骤：
步骤一：利用所述威布尔分布参数α、b计算伽玛分布的参数αg、λ，
由威布尔分布的参数α、b，可得其均值为方差为其中Γ为伽玛函数；当伽玛分布的参数为αg、λ时，其均值为方差为按照所求伽玛分布的均值和方差，与威布尔分布的均值和方差相等的原则，计算出αg、λ：
计算正态分布的参数μ、σ，
当正态分布的参数为μ、σ时，其均值为μ，方差为σ2；按照所求正态分布的均值和方差，与威布尔分布的均值和方差相等的原则，计算出μ、σ，
步骤二：计算偏度和峰度，
依据所述伽玛分布的参数αg、λ及所述正态分布的参数μ、σ，按下式计算三种分布的所述特征数中的偏度和峰度：
伽玛分布，偏度为峰度为
正态分布，偏度为0，峰度为0；
威布尔分布，偏度为
峰度为
步骤三：按如下规则比较所述特征数中的偏度和峰度，
判断所述步骤二中的伽玛分布与正态分布中偏度和峰度与所述威布尔分布的偏度和峰度的绝对差值情况；
步骤三：计算备件需求量，
以所述偏度特征数的绝对差值为首要比较条件，相比正态分布，若所述步骤三所得伽玛分布与所述威布尔分布的偏度绝对差值较小，则按照下式计算备件保障概率：
否则，按照下式计算备件保障概率Ps：
其中，Tw为保障任务时间，所述保障任务时间为所述备件完成任务的预期累积工作时间；
设置所述备件保障概率阈值，令j从0开始逐一递增，使得所述保障概率Ps大于或等于所述概率阈值的j值即为计算出的备件需求量。
总体而言，通过本发明所构思的以上技术方案与现有技术相比，具有以下有益效果：
(1)现有技术中也有近似方法，但是由于存在近似程度无法预估的问题导致算法的精度无法达到有效的目标，本发明率先提出了采用两种分布来执行近似计算，并且计算出偏度和峰度，依据两个指标上挑选近似方案，并且尤其是偏度上来进行判断，这样可以比较出近似程度，从而选取相应的近似分布执行备件需求量的计算；
(2)在现有技术的结论上，在威布尔的形状参数取极端时，一个极端趋向于伽玛，另外一个极端趋向于正态，在这种情况下，利用偏度和峰度这两个指标来选择趋向的近似，能够显著地提高计算的精确程度；
(3)按照本发明的近似方法，在近似成伽玛和正态近似分布的时候，提出了参数转化公式，按照上述的参数转化公式直接执行计算，无需迭代过程，从而节约了复杂的迭代过程；
(4)本发明还提出了参数转化公式的近似原则，即不管是伽玛或者是正态近似，务必满足在均值和方差这两个指标上与威布尔分布相等，对任何一个分布，都可以使用均值、方差、偏度和峰度来描述，本发明的计算方法中采用均值和方差相等，利用偏度和峰度作为近似计算指标的约定，从而简化了计算过程并且提高了近似的精度。
附图说明
图1为威布尔分布的形状参数为1.3情况下的三种分布的概率密度曲线对比情况；
图2为威布尔分布的形状参数为1.9情况下的三种分布的概率密度曲线对比情况；
图3为威布尔分布的形状参数为2.5情况下的三种分布的概率密度曲线对比情况；
图4为威布尔分布的形状参数为3.5情况下的三种分布的概率密度曲线对比情况。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
记随机变量X服从威布尔分布W(α,b)，威布尔分布密度函数如式(1)，
其中α＞0为尺度参数，b＞0为形状参数。
记随机变量X服从伽玛分布Ga(αg,λ)，其中αg＞0为形状参数，λ＞0为尺度参数，伽玛分布密度函数如式(2)。
式(2)中Γ(αg)为伽玛函数，且
记随机变量X服从正态分布N(μ,σ2)，其中μ为位置参数，μ的物理含义为寿命均值；σ为尺度参数，σ2的物理含义为寿命的方差。正态分布密度函数如式(3)。
对于寿命服从威布尔分布W(α,b)的单元，本发明计算备件需求量的步骤如下：
1)计算伽玛分布的参数αg、λ
按照所求伽玛分布的均值和方差，与威布尔分布的均值和方差相等的原则，按式(4)计算αg、λ，
2)计算正态分布的参数μ、σ
按照所求正态分布的均值和方差，与威布尔分布的均值和方差相等的原则，按式(5)计算μ、σ
3)计算偏度和峰度
已知伽玛分布的参数αg、λ，正态分布的参数μ、σ和威布尔分布α、b参数，按表1计算三种分布的偏度和峰度。
表1偏度和峰度的计算式
4)计算备件需求量
4.1)在偏度和峰度这两个特征数上，相比正态分布，若步骤3)所得伽玛分布更接近威布尔分布的偏度和峰度(以偏度的误差绝对值指标为主)，则按照式(6)计算备件保障概率；否则，按照式(7)计算备件保障概率Ps
式(6)、(7)中，Tw为保障任务时间。
4.2)设置所述单元保障概率阈值，令j从0开始逐一递增，使得所述保障概率Ps大于或等于所述概率阈值的j值即为所计算出的备件需求量。
描述随机变量的常见工具除了分布函数、概率密度函数外，也可用均值、方差、偏度和峰度这4种特征数来进行描述。这4种特征数由该分布的1～4阶矩决定。在大部分情况下，只要知道1阶矩到4阶矩就已经足够描述分布的情况。因此，通过比较两种分布在均值、方差、偏度和峰度这4种特征数的差异程度，可以了解二者的相似程度。
由于本发明的方法采用“令伽玛/正态分布的均值和方差，与原威布尔分布的均值和方差都相等”的原则，因此首先保证了所求得的伽玛/正态分布与原威布尔分布具有一定的相似性，再通过进一步比较偏度和峰度，从中选择更为相似的近似分布结果。表3列出了威布尔形状参数1.1～4.1范围内，按照“令伽玛/正态分布的均值和方差，与原威布尔分布的均值和方差都相等”原则，计算伽玛和正态分布参数后，这三种分布的偏度和峰度情况。
表3偏度和峰度结果
从表2可以看出，当威布尔分布的形状参数在2.1以内时，伽玛分布比正态分布更接近威布尔分布。图1～图4展示了威布尔分布的形状参数4种典型取值时，三种分布的概率密度曲线对比情况，与上述结论相符。
为了解释上述算法的准确性，本实施例1执行上述方法来进行备件需求量的计算，并且利用以下备件保障仿真模型开展仿真验证。
对于某个不可修单元，配置n个备件，该类单元的寿命服从威布尔分布W(α,b)，保障任务时间记为Tw，则模拟一次备件保障的过程如下：
(1)产生1+n个随机数ti(1≤i≤1+n)，随机数ti服从威布尔分布W(α,b)；
(2)计算累积工作时间
(3)当simT≥Tw时，保障任务成功，输出结果flag＝1；否则保障任务失败，输出结果flag＝0。
多次重复运行上述备件保障仿真模型，对所有模拟结果flag进行统计，flag均值即备件保障概率。
实施例1：某单元寿命服从威布尔分布W(200,1.3)，保障任务时间为1000h，要求备件保障概率不小于0.8,计算备件需求量。
1)计算伽玛分布的参数αg、λ
按照所求伽玛分布的均值和方差，与威布尔分布的均值和方差相等的原则，计算αg、λ
2)计算正态分布的参数μ、σ
按照所求正态分布的均值和方差，与威布尔分布的均值和方差相等的原则，计算μ、σ
3)计算偏度和峰度
伽玛分布Ga(1.66,0.009)的偏度、峰度为1.55、3.61；
正态分布N(184.7,143.32)的偏度、峰度为0、0；
威布尔分布W(200,1.3)的偏度、峰度为1.35、2.43。
4)计算备件需求量
4.1)经比较三种分布的偏度和峰度，相比正态分布，伽玛分布Ga(1.66,0.009)与威布尔分布W(200,1.3)的相似程度更高，因此按下式计算备件保障概率：
4.2)令j从0开始逐一递增，使得所述保障概率Ps大于或等于所述概率阈值的j值即为所计算出的备件需求量。计算过程中的结果如表4.
表4计算过程的结果
从表4可知，算例1的备件需求量为7。
实施例2：某单元寿命服从威布尔分布W(200,2.9)，保障任务时间为1000h，要求备件保障概率不小于0.8,计算备件需求量。
1)计算伽玛分布的参数αg、λ
按照所求伽玛分布的均值和方差，与威布尔分布的均值和方差相等的原则，计算αg、λ
2)计算正态分布的参数μ、σ
按照所求正态分布的均值和方差，与威布尔分布的均值和方差相等的原则，计算μ、σ
3)计算偏度和峰度
伽玛分布Ga(7.12,0.04)的偏度、峰度为0.75、0.84；
正态分布N(178.3,66.82)的偏度、峰度为0、0；
威布尔分布W(200,2.9)的偏度、峰度为0.20、-0.26。
4)计算备件需求量
4.1)经比较三种分布的偏度和峰度，相比伽玛分布，正态分布N(178.3,66.82)与威布尔分布W(200,2.9)的相似程度更高，因此按下式计算备件保障概率：
4.2)令j从0开始逐一递增，使得所述保障概率Ps大于或等于所述概率阈值的j值即为所计算出的备件需求量。计算过程中的结果如表5.
表5计算过程的结果
从表5可知，算例2的备件需求量为6。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (1)

1.一种利用特征数的威布尔型备件需求量的近似计算方法，所述威布尔型备件的寿命服从威布尔分布W(α,b)，α、b为威布尔分布参数，α为尺度参数，b为形状参数；所述特征数为均值、方差、偏度和峰度，其特征在于，该计算方法包括如下步骤：
步骤一：利用所述威布尔分布参数α、b计算伽玛分布的参数αg、λ，
由威布尔分布的参数α、b，可得其均值为方差为其中Γ为伽玛函数；当伽玛分布的参数为αg、λ时，其均值为方差为按照所求伽玛分布的均值和方差，与威布尔分布的均值和方差相等的原则，计算出αg、λ：
<mrow> <msub> <mi>&amp;alpha;</mi> <mi>g</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msup> <mi>&amp;Gamma;</mi> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mfrac> <mn>1</mn> <mi>b</mi> </mfrac> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mfrac> <mn>2</mn> <mi>b</mi> </mfrac> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>-</mo> <msup> <mi>&amp;Gamma;</mi> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mfrac> <mn>1</mn> <mi>b</mi> </mfrac> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>,</mo> <mi>&amp;lambda;</mi> <mo>=</mo> <mfrac> <mrow> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mfrac> <mn>1</mn> <mi>b</mi> </mfrac> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>a</mi> <mo>&amp;lsqb;</mo> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mfrac> <mn>2</mn> <mi>b</mi> </mfrac> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>-</mo> <msup> <mi>&amp;Gamma;</mi> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mfrac> <mn>1</mn> <mi>b</mi> </mfrac> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mfrac> <mo>;</mo> </mrow>
计算正态分布的参数μ、σ，
当正态分布的参数为μ、σ时，其均值为μ，方差为σ2；按照所求正态分布的均值和方差，与威布尔分布的均值和方差相等的原则，计算出μ、σ，
<mrow> <mi>&amp;mu;</mi> <mo>=</mo> <mi>&amp;alpha;</mi> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mfrac> <mn>1</mn> <mi>b</mi> </mfrac> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>,</mo> <mi>&amp;sigma;</mi> <mo>=</mo> <mi>a</mi> <msqrt> <mrow> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mfrac> <mn>2</mn> <mi>b</mi> </mfrac> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>-</mo> <msup> <mi>&amp;Gamma;</mi> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mfrac> <mn>1</mn> <mi>b</mi> </mfrac> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>;</mo> </mrow>
步骤二：计算偏度和峰度，
依据所述伽玛分布的参数αg、λ及所述正态分布的参数μ、σ，按下式计算三种分布的所述特征数中的偏度和峰度：
伽玛分布，偏度为峰度为
正态分布，偏度为0，峰度为0；
威布尔分布，偏度为
峰度为
步骤三：按如下规则比较所述特征数中的偏度和峰度，
判断所述步骤二中的伽玛分布与正态分布中偏度和峰度与所述威布尔分布的偏度和峰度的绝对差值情况；
步骤三：计算备件需求量，
以所述偏度特征数的绝对差值为首要比较条件，相比正态分布，若所述步骤三所得伽玛分布与所述威布尔分布的偏度绝对差值较小，则按照式(1)计算备件保障概率；否则，按照式(2)计算备件保障概率Ps：
<mrow> <msub> <mi>P</mi> <mi>s</mi> </msub> <mo>=</mo> <mi>P</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>&gt;</mo> <msub> <mi>T</mi> <mi>w</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mfrac> <msup> <mi>&amp;lambda;</mi> <msub> <mi>&amp;alpha;</mi> <mi>g</mi> </msub> </msup> <mrow> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <mi>j</mi> </mrow> <mo>)</mo> <msub> <mi>&amp;alpha;</mi> <mi>g</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <msub> <mi>T</mi> <mi>w</mi> </msub> </msubsup> <msup> <mi>x</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>j</mi> <mo>)</mo> <msub> <mi>&amp;alpha;</mi> <mi>g</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>&amp;lambda;</mi> <mi>x</mi> </mrow> </msup> <mi>d</mi> <mi>x</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>P</mi> <mi>s</mi> </msub> <mo>=</mo> <mi>P</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>&gt;</mo> <msub> <mi>T</mi> <mi>w</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mfrac> <mn>1</mn> <mrow> <mi>&amp;sigma;</mi> <msqrt> <mrow> <mn>2</mn> <mi>&amp;pi;</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>j</mi> <mo>)</mo> </mrow> </mrow> </msqrt> </mrow> </mfrac> <msubsup> <mo>&amp;Integral;</mo> <mrow> <mo>-</mo> <mi>&amp;infin;</mi> </mrow> <msub> <mi>T</mi> <mi>w</mi> </msub> </msubsup> <msup> <mi>e</mi> <mfrac> <mrow> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <mi>j</mi> </mrow> <mo>)</mo> <mi>&amp;mu;</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>j</mi> <mo>)</mo> </mrow> <msup> <mi>&amp;sigma;</mi> <mn>2</mn> </msup> </mrow> </mfrac> </msup> <mi>d</mi> <mi>x</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
式(1)、(2)中，Tw为保障任务时间，所述保障任务时间为所述备件完成任务的预期累积工作时间；
设置所述备件保障概率阈值，令j从0开始逐一递增，使得所述保障概率Ps大于或等于所述概率阈值的j值即为计算出的备件需求量。
CN201710343600.1A 2017-05-16 2017-05-16 一种利用特征数的威布尔型备件需求量的近似计算方法 Active CN107220216B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710343600.1A CN107220216B (zh) 2017-05-16 2017-05-16 一种利用特征数的威布尔型备件需求量的近似计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710343600.1A CN107220216B (zh) 2017-05-16 2017-05-16 一种利用特征数的威布尔型备件需求量的近似计算方法

Publications (2)

Publication Number Publication Date
CN107220216A true CN107220216A (zh) 2017-09-29
CN107220216B CN107220216B (zh) 2020-07-28

Family

ID=59945251

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710343600.1A Active CN107220216B (zh) 2017-05-16 2017-05-16 一种利用特征数的威布尔型备件需求量的近似计算方法

Country Status (1)

Country Link
CN (1) CN107220216B (zh)

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108334720A (zh) * 2018-04-09 2018-07-27 中国人民解放军海军工程大学 一种贮存失效风险下的正态型单元备件需求量计算方法
CN108460495A (zh) * 2018-03-28 2018-08-28 联想(北京)有限公司 一种数据处理方法及装置
CN108470250A (zh) * 2018-03-25 2018-08-31 中国人民解放军海军工程大学 一种正态型串联部件的备件需求量计算方法
CN108509389A (zh) * 2018-03-25 2018-09-07 中国人民解放军海军工程大学 一种威布尔型串联部件的备件需求量计算方法
CN108509390A (zh) * 2018-03-25 2018-09-07 中国人民解放军海军工程大学 一种伽玛型串联部件的备件需求量计算方法
CN108520145A (zh) * 2018-04-09 2018-09-11 中国人民解放军海军工程大学 一种贮存失效风险下的威布尔型单元备件需求量计算方法
CN108564270A (zh) * 2018-04-09 2018-09-21 中国人民解放军海军工程大学 一种贮存失效风险下的伽玛型单元备件需求量计算方法
CN108932358A (zh) * 2018-04-16 2018-12-04 中国人民解放军海军工程大学 一种确定威布尔型单元后续备件需求量的方法
CN109145502A (zh) * 2018-09-17 2019-01-04 中国人民解放军海军工程大学 一种威布尔型单元寿命分布参数估计方法
CN109543276A (zh) * 2018-11-15 2019-03-29 青岛航讯网络技术服务有限公司 大型货轮长期保障任务多Gamma单元备件需求量确定方法
CN110704884A (zh) * 2019-09-30 2020-01-17 中国人民解放军海军工程大学 一种对数正态型通用件备件需求量计算方法及装置
CN110717265A (zh) * 2019-09-30 2020-01-21 中国人民解放军海军工程大学 一种对数正态型单元的备件需求量近似计算方法及装置
CN110737865A (zh) * 2019-09-09 2020-01-31 中国人民解放军海军工程大学 一种基于分段函数的威布尔型备件需求量确定方法及系统
EP3916495A1 (de) * 2020-05-28 2021-12-01 Carl Zeiss Industrielle Messtechnik GmbH Verfahren und anordnung zum auswerten oder zum vorbereiten einer auswertung einer stichprobe von messdaten aus einer vermessung einer vielzahl von werkstücken
CN114529018A (zh) * 2022-01-14 2022-05-24 中国人民解放军海军工程大学 一种基于Gamma分布的舰船备件需求近似计算方法
CN116579573A (zh) * 2023-05-23 2023-08-11 中国人民解放军海军工程大学 一种满足维修需求的机电件备件保障分析方法
CN116611654A (zh) * 2023-05-23 2023-08-18 中国人民解放军海军工程大学 机电通用件在维修工作下的备件需求量估计方法及系统
CN114529018B (zh) * 2022-01-14 2024-05-31 中国人民解放军海军工程大学 一种基于Gamma分布的舰船备件需求近似计算方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1363211A3 (en) * 2002-05-13 2004-06-23 Honda Giken Kogyo Kabushiki Kaisha System for predicting a demand for repair parts

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1363211A3 (en) * 2002-05-13 2004-06-23 Honda Giken Kogyo Kabushiki Kaisha System for predicting a demand for repair parts

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
刘天华 等: "威布尔型可修备件需求量的解析算法研究", 《系统工程与电子技术》 *
张建军 等: "备件保障度评估与备件需求量模型研究", 《电子产品可靠性与环境试验》 *
李能贵: "《电子元器件的可靠性》", 30 June 1990 *
白和丰 等: "可靠性试验数据分布类型的综合判别法", 《武汉工业大学学报》 *

Cited By (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108509390B (zh) * 2018-03-25 2021-11-02 中国人民解放军海军工程大学 一种伽玛型串联部件的备件需求量计算方法
CN108509390A (zh) * 2018-03-25 2018-09-07 中国人民解放军海军工程大学 一种伽玛型串联部件的备件需求量计算方法
CN108470250A (zh) * 2018-03-25 2018-08-31 中国人民解放军海军工程大学 一种正态型串联部件的备件需求量计算方法
CN108509389A (zh) * 2018-03-25 2018-09-07 中国人民解放军海军工程大学 一种威布尔型串联部件的备件需求量计算方法
CN108470250B (zh) * 2018-03-25 2021-11-02 中国人民解放军海军工程大学 一种正态型串联部件的备件需求量计算方法
CN108509389B (zh) * 2018-03-25 2022-01-14 中国人民解放军海军工程大学 一种威布尔型串联部件的备件需求量计算方法
CN108460495B (zh) * 2018-03-28 2020-12-18 联想(北京)有限公司 一种数据处理方法及装置
CN108460495A (zh) * 2018-03-28 2018-08-28 联想(北京)有限公司 一种数据处理方法及装置
CN108334720A (zh) * 2018-04-09 2018-07-27 中国人民解放军海军工程大学 一种贮存失效风险下的正态型单元备件需求量计算方法
CN108520145B (zh) * 2018-04-09 2022-04-05 中国人民解放军海军工程大学 一种贮存失效风险下的威布尔型单元备件需求量计算方法
CN108564270A (zh) * 2018-04-09 2018-09-21 中国人民解放军海军工程大学 一种贮存失效风险下的伽玛型单元备件需求量计算方法
CN108520145A (zh) * 2018-04-09 2018-09-11 中国人民解放军海军工程大学 一种贮存失效风险下的威布尔型单元备件需求量计算方法
CN108564270B (zh) * 2018-04-09 2021-11-02 中国人民解放军海军工程大学 一种贮存失效风险下的伽玛型单元备件需求量计算方法
CN108334720B (zh) * 2018-04-09 2021-08-31 中国人民解放军海军工程大学 一种贮存失效风险下的正态型单元备件需求量计算方法
CN108932358A (zh) * 2018-04-16 2018-12-04 中国人民解放军海军工程大学 一种确定威布尔型单元后续备件需求量的方法
CN108932358B (zh) * 2018-04-16 2021-08-13 中国人民解放军海军工程大学 一种确定威布尔型单元后续备件需求量的方法
CN109145502A (zh) * 2018-09-17 2019-01-04 中国人民解放军海军工程大学 一种威布尔型单元寿命分布参数估计方法
CN109543276A (zh) * 2018-11-15 2019-03-29 青岛航讯网络技术服务有限公司 大型货轮长期保障任务多Gamma单元备件需求量确定方法
CN109543276B (zh) * 2018-11-15 2022-11-01 青岛航讯网络技术服务有限公司 大型货轮长期保障任务多Gamma单元备件需求量确定方法
CN110737865B (zh) * 2019-09-09 2023-05-19 中国人民解放军海军工程大学 一种基于分段函数的威布尔型备件需求量确定方法及系统
CN110737865A (zh) * 2019-09-09 2020-01-31 中国人民解放军海军工程大学 一种基于分段函数的威布尔型备件需求量确定方法及系统
CN110704884B (zh) * 2019-09-30 2023-05-12 中国人民解放军海军工程大学 一种对数正态型通用件备件需求量计算方法及装置
CN110717265B (zh) * 2019-09-30 2022-06-14 中国人民解放军海军工程大学 一种对数正态型单元的备件需求量近似计算方法及装置
CN110704884A (zh) * 2019-09-30 2020-01-17 中国人民解放军海军工程大学 一种对数正态型通用件备件需求量计算方法及装置
CN110717265A (zh) * 2019-09-30 2020-01-21 中国人民解放军海军工程大学 一种对数正态型单元的备件需求量近似计算方法及装置
EP3916495A1 (de) * 2020-05-28 2021-12-01 Carl Zeiss Industrielle Messtechnik GmbH Verfahren und anordnung zum auswerten oder zum vorbereiten einer auswertung einer stichprobe von messdaten aus einer vermessung einer vielzahl von werkstücken
US11822321B2 (en) 2020-05-28 2023-11-21 Carl Zeiss Industrielle Messtechnik Gmbh Coordinate measuring machine measurement and analysis of multiple workpieces
CN114529018A (zh) * 2022-01-14 2022-05-24 中国人民解放军海军工程大学 一种基于Gamma分布的舰船备件需求近似计算方法
CN114529018B (zh) * 2022-01-14 2024-05-31 中国人民解放军海军工程大学 一种基于Gamma分布的舰船备件需求近似计算方法
CN116579573A (zh) * 2023-05-23 2023-08-11 中国人民解放军海军工程大学 一种满足维修需求的机电件备件保障分析方法
CN116611654A (zh) * 2023-05-23 2023-08-18 中国人民解放军海军工程大学 机电通用件在维修工作下的备件需求量估计方法及系统
CN116611654B (zh) * 2023-05-23 2024-02-09 中国人民解放军海军工程大学 机电通用件在维修工作下的备件需求量估计方法及系统
CN116579573B (zh) * 2023-05-23 2024-02-09 中国人民解放军海军工程大学 一种满足维修需求的机电件备件保障分析方法

Also Published As

Publication number Publication date
CN107220216B (zh) 2020-07-28

Similar Documents

Publication Publication Date Title
CN107220216A (zh) 一种利用特征数的威布尔型备件需求量的近似计算方法
CN108509389A (zh) 一种威布尔型串联部件的备件需求量计算方法
CN104236889B (zh) 橡胶密封圈的贮存寿命优化设计方法
CN112598248B (zh) 负荷预测方法、装置、计算机设备和存储介质
CN107230977A (zh) 基于误差修正和提升小波组合预测模型的风电功率预测方法
CN106780130A (zh) 一种含分布式光伏配电网的评价方法
Caldera et al. A comparison of methods of estimating missing daily rainfall data
CN106875047B (zh) 水库流域径流预报方法及系统
CN105046321A (zh) 一种基于相似性搜索的组合模型预测水位的方法
CN111784116A (zh) 一种多能耦合下配电网弹性能力评估方法及系统
Liu et al. Reliability mathematical models of repairable systems with uncertain lifetimes and repair times
CN107590570A (zh) 一种负载功率预测方法及系统
CN113987684A (zh) 基于sn曲线的随机应力载荷谱分级方法、系统和存储介质
CN108932358B (zh) 一种确定威布尔型单元后续备件需求量的方法
Latief et al. Preliminary cost estimation using regression analysis incorporated with adaptive neuro fuzzy inference system
Adnan et al. An integrated neural-fuzzy system of software reliability prediction
Zaitseva et al. Dynamic reliability indices for parallel, series and k-out-of-n multi-state systems
CN108051676A (zh) 一种雷电流幅值累积概率分布曲线拟合计算方法
CN112101797A (zh) 一种复杂工业系统动态故障诊断方法及系统
Wilkie et al. Fatigue reliability of offshore wind turbines using Gaussian processes
CN106875097B (zh) 基于混合总体的轨道交通产品现场可靠性评估方法及装置
Cluever et al. Estimation of Weibull Life Distributions From Expert Categorical Estimates of Failure Probabilities
Darby Tools for evaluating risk of terrorist acts using fuzzy sets and belief/plausibility
Li et al. Study on Perceptive Fuzzy Petri Net-based Autoloader Fault Analysis
Kordas et al. Matching using semiparametric propensity scores

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant