Co to jest rozkład Weibulla?
Rozkład Weibulla to rozkład ciągły używany przede wszystkim do modelowania czasu życia, czasu do awarii, trwałości elementów, niezawodności urządzeń i procesów, w których ryzyko zdarzenia zmienia się wraz z czasem. W przeciwieństwie do prostego rozkładu wykładniczego, Weibull pozwala opisać sytuacje, w których awaryjność maleje, pozostaje stała albo rośnie.
Rozkład zapisuje się najczęściej za pomocą parametrów k i λ. Parametr k jest parametrem kształtu, a λ parametrem skali. Dodatkowo kalkulator pozwala użyć przesunięcia γ, które określa minimalny punkt startu. Gdy γ=0, otrzymujemy najczęściej spotykaną wersję bez przesunięcia.
Weibull dobrze łączy się z kalkulatorami rozkładów czasu oczekiwania. Jeżeli zakładasz stałą intensywność zdarzeń, porównaj wynik z rozkładem wykładniczym. Jeżeli modelujesz sumę czasów albo bardziej ogólne czasy dodatnie, przydatny będzie też rozkład gamma.
Parametr k i trzy typy awaryjności
Najważniejszą informacją w rozkładzie Weibulla jest wartość parametru k. To ona mówi, jak zachowuje się hazard, czyli chwilowa intensywność wystąpienia zdarzenia pod warunkiem, że do tej pory zdarzenie jeszcze nie wystąpiło.
| Parametr k | Hazard | Praktyczna interpretacja |
|---|---|---|
| k < 1 | malejący | wczesne awarie, selekcja wadliwych sztuk, największe ryzyko na początku |
| k = 1 | stały | model taki jak rozkład wykładniczy, brak starzenia w czasie |
| k > 1 | rosnący | zużycie, starzenie, zmęczenie materiału, rosnące ryzyko awarii |
Dzięki temu rozkład Weibulla jest popularny w analizie niezawodności. Ten sam wzór może opisać bardzo różne mechanizmy: początkowe wady produkcyjne, przypadkowe awarie w trakcie pracy albo stopniowe zużywanie się elementu. Właśnie dlatego w praktyce parametr k bywa ważniejszy interpretacyjnie niż sama średnia.
Najważniejsze wzory dla rozkładu Weibulla
Dla rozkładu Weibulla z parametrem kształtu k, skalą λ i przesunięciem γ używamy zmiennej pomocniczej t=(x-γ)/λ. Dla x<γ dystrybuanta wynosi 0, bo czas lub wartość nie może pojawić się przed punktem startowym.
| Wielkość | Wzór | Znaczenie |
|---|---|---|
| dystrybuanta | F(x)=1-exp(-tk) | prawdopodobieństwo, że zdarzenie nastąpi do czasu x |
| niezawodność | R(x)=exp(-tk) | prawdopodobieństwo przetrwania powyżej x |
| gęstość | f(x)=(k/λ)tk-1exp(-tk) | kształt rozkładu czasu do zdarzenia |
| hazard | h(x)=(k/λ)tk-1 | chwilowa intensywność awarii |
| kwantyl | xp=γ+λ[-ln(1-p)]1/k | czas, do którego zdarzenie nastąpi z prawdopodobieństwem p |
Dystrybuanta i niezawodność - dwie strony tego samego wyniku
W rozkładzie Weibulla dystrybuanta F(x) i niezawodność R(x) są dopełnieniami. Jeśli F(x)=0,30, to oznacza, że do czasu x zdarzenie wystąpi z prawdopodobieństwem 30%. Jednocześnie R(x)=0,70 oznacza, że element przetrwa ponad x z prawdopodobieństwem 70%.
W analizie awarii często wygodniej mówić o niezawodności niż o dystrybuancie. Przykład: jeśli x oznacza 1000 godzin pracy, R(1000) informuje, jaka część elementów powinna działać dłużej niż 1000 godzin. Ten sposób interpretacji jest bardziej naturalny w utrzymaniu ruchu, testach trwałości, produkcji i analizie gwarancji.
Kalkulator pokazuje oba wyniki jednocześnie, bo pomylenie „awaria do czasu x” z „przetrwanie powyżej x” jest jednym z najczęstszych błędów przy Weibullu.
Kwantyle i czas B10, B50, B90
Kwantyl rozkładu Weibulla mówi, do jakiego czasu zdarzenie nastąpi z określonym prawdopodobieństwem. W niezawodności często spotyka się oznaczenia takie jak B10, B50 czy B90. B10 oznacza czas, do którego awarii ulegnie 10% populacji, a 90% nadal przetrwa. B50 jest medianą czasu do awarii.
W kalkulatorze tryb kwantyla pozwala wpisać dowolne p, np. 0,10, 0,50 albo 0,90. Wynik xp nie jest średnią. To granica, pod którą znajduje się p części rozkładu. Dla skośnych rozkładów Weibulla średnia, mediana i dominanta mogą znacząco się różnić.
Średnia, wariancja, mediana i dominanta
Parametry opisowe rozkładu Weibulla korzystają z funkcji gamma. Kalkulator liczy je automatycznie, ale warto znać ich interpretację. Średnia mówi o przeciętnym czasie lub wartości, mediana mówi o punkcie, przed którym pojawi się połowa zdarzeń, a wariancja i odchylenie standardowe opisują rozrzut.
E(X)=γ+λΓ(1+1/k) Var(X)=λ²[Γ(1+2/k)-Γ(1+1/k)²]
Dominanta, czyli najbardziej prawdopodobna okolica, istnieje wewnątrz rozkładu tylko dla k>1. Gdy k≤1, gęstość jest największa przy początku rozkładu. To dobrze pasuje do sytuacji, w której największe ryzyko występuje na starcie, np. przy wadach ujawniających się zaraz po rozpoczęciu pracy.
Przykład: czas do awarii elementu
Załóżmy, że czas pracy elementu opisuje rozkład Weibulla z k=1,8 i λ=1200 godzin. Wartość k>1 oznacza rosnący hazard, czyli ryzyko awarii wzrasta wraz z czasem. Chcemy policzyć, jaka jest szansa, że element ulegnie awarii do 1000 godzin.
Kalkulator wyznacza F(1000), czyli prawdopodobieństwo awarii do 1000 godzin, oraz R(1000), czyli prawdopodobieństwo przetrwania powyżej 1000 godzin. Dodatkowo pokazuje hazard w tym punkcie, co pozwala lepiej zrozumieć, czy ryzyko chwilowe jest wysokie względem wcześniejszych momentów.
Jeżeli interesuje Cię zakres, np. awaria między 800 a 1500 godzin, wybierz tryb przedziału. Wynik będzie równy F(1500)-F(800), czyli prawdopodobieństwu, że zdarzenie wystąpi właśnie w tym odcinku czasu.
Weibull w praktyce: produkcja, gwarancje i utrzymanie ruchu
Rozkład Weibulla jest użyteczny wtedy, gdy samo pytanie „jaka jest średnia trwałość?” nie wystarcza. W produkcji często ważniejsze jest, jaki odsetek elementów może zawieść przed określonym czasem granicznym. W serwisie i utrzymaniu ruchu istotne jest, kiedy ryzyko zaczyna rosnąć tak mocno, że opłaca się wymiana prewencyjna.
Przykład: dwa elementy mogą mieć podobną średnią trwałość, ale inny parametr k. Element z k bliskim 1 zachowuje się bardziej losowo, a element z dużym k ma bardziej wyraźny etap zużycia. To zmienia decyzje: przy rosnącym hazardzie sensowniejsze mogą być przeglądy po określonym czasie pracy, natomiast przy hazardzie stałym profilaktyczna wymiana może dawać mniejszą korzyść.
Kalkulator nie zastępuje pełnej analizy niezawodności na danych z testów, ale pomaga zrozumieć, jak parametry rozkładu wpływają na ryzyko, medianę, kwantyle i przedziały czasu życia.
Zadania przykładowe
Dla k=1, co przypomina rozkład Weibulla?
Co oznacza k>1 w analizie awarii?
Jak policzyć niezawodność R(x), gdy znasz F(x)?
Co oznacza kwantyl B10?
Kiedy dominanta leży wewnątrz rozkładu?
Jak policzyć P(a≤X≤b) w rozkładzie Weibulla?
Najczęstsze błędy przy rozkładzie Weibulla
- Mylenie dystrybuanty z niezawodnością - F(x) to awaria do czasu x, a R(x) to przetrwanie powyżej x.
- Ignorowanie parametru k - to on mówi, czy hazard maleje, jest stały, czy rośnie.
- Traktowanie λ jako średniej - skala λ nie jest zwykle równa średniej rozkładu.
- Porównywanie Weibulla z wykładniczym bez sprawdzenia k - wykładniczy odpowiada szczególnemu przypadkowi k=1.
- Używanie kwantyla jako średniej - B10, B50 i B90 są granicami procentowymi, nie przeciętnym czasem.
- Pomijanie jednostek - λ, x i γ muszą być w tych samych jednostkach, np. godzinach, dniach albo cyklach.
Sprawdź także
Do prostszego modelu czasu oczekiwania użyj rozkładu wykładniczego. Do dodatnich czasów i sum czasów sprawdź rozkład gamma. Jeśli analizujesz dane normalne albo standaryzację, przyda się rozkład normalny i standaryzacja z-score.
Wskazówka od KalkulatorXXL
Przy interpretacji Weibulla zacznij od parametru k. Dopiero potem patrz na średnią, kwantyle i niezawodność. Ten jeden parametr często od razu mówi, czy problem dotyczy wczesnych awarii, losowego czasu oczekiwania, czy zużycia w czasie.
Cały pakiet statystyki i rozkładów
Ten kalkulator jest częścią większego pakietu narzędzi do statystyki, prawdopodobieństwa, rozkładów i analizy danych. Jeśli chcesz przejść od pojedynczego wyniku do szerszego zestawu metod, sprawdź Statystyka, prawdopodobieństwo i rozkłady - akademia obliczeń. Znajdziesz tam ścieżki do rozkładów dyskretnych i ciągłych, testów statystycznych, z-score, wartości odstających, ANOVA oraz narzędzi do oceny modeli.