[an error occurred while processing this directive]
Quantities in Probability & Statistics

Order Distribution from Distribution with Units

A system has three components, out of which one is primary and two are redundant. Each component has a lifetime following a WeibullDistribution, with shape parameter 2 and mean lifetime of 885 hours.

In[1]:=
Click for copyable input
{sol} = NSolve[ Mean[WeibullDistribution[2, b]] == Quantity[885, "Hours"], b]
Out[1]=
In[2]:=
Click for copyable input
\[ScriptD] = WeibullDistribution[2, b] /. sol
Out[2]=

The lifetime of this system can be described as the maximum lifetime of its components.

In[3]:=
Click for copyable input
\[ScriptCapitalD] = OrderDistribution[{\[ScriptD], 3}, 3]
Out[3]=

The mean lifetime of the system.

In[4]:=
Click for copyable input
Mean[\[ScriptCapitalD]]
Out[4]=

Compare with the computation using ReliabilityDistribution.

In[5]:=
Click for copyable input
Mean[ReliabilityDistribution[ comp1 \[Or] comp2 \[Or] comp3, {{comp1, \[ScriptD]}, {comp2, \[ScriptD]}, {comp3, \ \[ScriptD]}}]]
Out[5]=

Compare probability density function of the lifetime of a single component with that of the system.

show complete Wolfram Language input
In[6]:=
Click for copyable input
Plot[{PDF[\[ScriptD], Quantity[x, "Hours"]], PDF[\[ScriptCapitalD], Quantity[x, "Hours"]]}, {x, 0, 3000}, FrameLabel -> {"hr"}, PlotTheme -> "Detailed", PlotLegends -> {"single component", "system"}, PlotLabel -> "Lifetime Distribution Density"]
Out[6]=

Related Examples

[an error occurred while processing this directive]
de es fr ja ko pt-br ru zh