Funkcja rozkładu zmiennej losowej. Ciągłe zmienne losowe. Znajdź funkcję rozkładu F(x)

Rozważana właściwość rozkładów dyskretnych stwarza znaczne trudności podczas zestawiania i stosowania takich rozkładów, ponieważ niemożliwe jest dokładne zachowanie typowych wartości liczbowych charakterystyk rozkładu. W szczególności dotyczy to wartości krytycznych i poziomów istotności nieparametrycznych testów statystycznych (patrz poniżej), ponieważ rozkłady statystyk tych testów są dyskretne.

Kolejność kwantylowa ma ogromne znaczenie w statystyce R= ½. Nazywa się to medianą (zmienna losowa X lub jego funkcje dystrybucyjne F(x)) i jest wyznaczony Ja(X). W geometrii istnieje pojęcie „środkowej” - linii prostej przechodzącej przez wierzchołek trójkąta i dzielącej jego przeciwny bok na pół. W statystyce matematycznej mediana dzieli na pół nie bok trójkąta, ale rozkład zmiennej losowej: równość F(x 0,5)= 0,5 oznacza prawdopodobieństwo dostania się w lewo x 0,5 i prawdopodobieństwo dotarcia w prawo x 0,5(lub bezpośrednio do x 0,5) są sobie równe i równe ½, tj.

P(X < X 0,5) = P(X > X 0,5) = ½.

Mediana wskazuje „centrum” rozkładu. Z punktu widzenia jednej ze współczesnych koncepcji – teorii stabilnych procedur statystycznych – mediana jest lepszą charakterystyką zmiennej losowej niż oczekiwanie matematyczne. Podczas przetwarzania wyników pomiarów na skali porządkowej (patrz rozdział o teorii pomiaru) można zastosować medianę, ale nie można zastosować oczekiwań matematycznych.

Cecha zmiennej losowej taka jak mod ma jasne znaczenie - wartość (lub wartości) zmiennej losowej odpowiadająca lokalnemu maksimum gęstości prawdopodobieństwa dla ciągłej zmiennej losowej lub lokalnemu maksimum prawdopodobieństwa dla dyskretnej zmiennej losowej .

Jeśli x 0– postać zmiennej losowej z gęstością f(x), wówczas, jak wiadomo z rachunku różniczkowego, .

Zmienna losowa może mieć wiele postaci. Zatem dla równomiernego rozkładu (1) każdy punkt X takie, że A< x < b , jest moda. Jest to jednak wyjątek. Większość zmiennych losowych stosowanych w probabilistycznych statystycznych metodach podejmowania decyzji i innych badaniach stosowanych ma jeden tryb. Zmienne losowe, gęstości i rozkłady posiadające jeden tryb nazywane są unimodalnymi.

Oczekiwanie matematyczne dla dyskretnych zmiennych losowych o skończonej liczbie wartości zostało omówione w rozdziale „Zdarzenia i prawdopodobieństwa”. Dla ciągłej zmiennej losowej X wartość oczekiwana M(X) spełnia równość

co jest analogią wzoru (5) ze stwierdzenia 2 rozdziału „Zdarzenia i prawdopodobieństwa”.

Przykład 5. Oczekiwanie na równomiernie rozłożoną zmienną losową X równa się

W przypadku zmiennych losowych rozpatrywanych w tym rozdziale wszystkie właściwości oczekiwań matematycznych i wariancji, które zostały wcześniej rozważone dla dyskretnych zmiennych losowych o skończonej liczbie wartości, są prawdziwe. Nie przedstawiamy jednak dowodu na te właściwości, ponieważ wymagają one zagłębienia się w subtelności matematyczne, co nie jest konieczne do zrozumienia i kwalifikowanego zastosowania probabilistyczno-statystycznych metod podejmowania decyzji.

Komentarz. W podręczniku celowo unika się subtelności matematycznych związanych w szczególności z pojęciami zbiorów mierzalnych i funkcji mierzalnych, algebry zdarzeń itp. Pragnący opanować te pojęcia powinni sięgnąć po literaturę specjalistyczną, w szczególności encyklopedię.

Każda z trzech cech – oczekiwanie matematyczne, mediana, moda – opisuje „środek” rozkładu prawdopodobieństwa. Pojęcie „centrum” można definiować na różne sposoby – stąd trzy różne cechy. Jednak w przypadku ważnej klasy rozkładów – symetrycznych unimodalnych – wszystkie trzy cechy są zbieżne.

Gęstość dystrybucji k(x)– gęstość rozkładu symetrycznego, jeśli istnieje liczba x 0 takie, że

. (3)

Równość (3) oznacza, że ​​wykres funkcji y = f(x) symetryczny względem linii pionowej przechodzącej przez środek symetrii X = X 0. Z (3) wynika, że ​​rozkład symetryczny spełnia zależność

(4)

W przypadku rozkładu symetrycznego z jedną modą oczekiwanie matematyczne, mediana i moda pokrywają się i są równe x 0.

Najważniejszym przypadkiem jest symetria wokół 0, tj. x 0= 0. Wtedy (3) i (4) stają się równościami

(6)

odpowiednio. Powyższe zależności pokazują, że nie ma potrzeby zestawiania rozkładów symetrycznych dla wszystkich X, wystarczy mieć stoliki przy X > x 0.

Zwróćmy uwagę na jeszcze jedną właściwość rozkładów symetrycznych, która jest stale wykorzystywana w probabilistyczno-statystycznych metodach podejmowania decyzji i innych badaniach stosowanych. Dla funkcji rozkładu ciągłego

P(|X| < a) = P(-a < X < a) = F(a) – F(-a),

Gdzie F– dystrybuanta zmiennej losowej X. Jeśli funkcja dystrybucji F jest symetryczny względem 0, tj. wówczas obowiązuje dla niego wzór (6).

P(|X| < a) = 2F(a) – 1.

Często używane jest inne sformułowanie omawianego stwierdzenia: jeśli

.

Jeżeli i są kwantylami rzędu i odpowiednio (patrz (2)) funkcji rozkładu symetrycznej wokół 0, to z (6) wynika, że

Z charakterystyki pozycji – oczekiwanie matematyczne, mediana, moda – przejdźmy do charakterystyki rozrzutu zmiennej losowej X: wariancja, średnia odchylenie kwadratowe i współczynnik zmienności w. Definicja i właściwości dyspersji dyskretnych zmiennych losowych zostały omówione w poprzednim rozdziale. Dla ciągłych zmiennych losowych

Odchylenie standardowe to nieujemna wartość pierwiastka kwadratowego wariancji:

Współczynnik zmienności to stosunek odchylenia standardowego do oczekiwań matematycznych:

Współczynnik zmienności stosuje się, gdy M(X)> 0. Mierzy rozrzut w jednostkach względnych, podczas gdy odchylenie standardowe jest wyrażone w jednostkach bezwzględnych.

Przykład 6. Dla równomiernie rozłożonej zmiennej losowej X Znajdźmy dyspersję, odchylenie standardowe i współczynnik zmienności. Wariancja wynosi:

Zmiana zmiennej umożliwia zapis:

Gdzie C = (B – A)/ 2. Zatem odchylenie standardowe jest równe, a współczynnik zmienności wynosi:

Dla każdej zmiennej losowej X określ jeszcze trzy wielkości - wyśrodkowane Y, znormalizowany V i dane U. Wyśrodkowana zmienna losowa Y jest różnicą między daną zmienną losową X i jego matematyczne oczekiwanie M(X), te. Y = X – M(X). Oczekiwanie wyśrodkowanej zmiennej losowej Y równa się 0, a wariancja jest wariancją danej zmiennej losowej: M(Y) = 0, D(Y) = D(X). Funkcja dystrybucyjna F Y(X) wyśrodkowana zmienna losowa Y związane z funkcją dystrybucji F(X) pierwotna zmienna losowa X stosunek:

F Y(X) = F(X + M(X)).

Gęstości tych zmiennych losowych zachowują następującą równość:

dla Y(X) = F(X + M(X)).

Znormalizowana zmienna losowa V jest stosunkiem danej zmiennej losowej X do odchylenia standardowego, tj. . Oczekiwanie i wariancja znormalizowanej zmiennej losowej V wyrażone poprzez cechy X Więc:

,

Gdzie w– współczynnik zmienności pierwotnej zmiennej losowej X. Dla funkcji dystrybucji F V(X) i gęstość f V(X) znormalizowana zmienna losowa V mamy:

Gdzie F(X) – rozkład pierwotnej zmiennej losowej X, A F(X) – jego gęstość prawdopodobieństwa.

Zredukowana zmienna losowa U jest wyśrodkowaną i znormalizowaną zmienną losową:

.

Dla danej zmiennej losowej

Znormalizowane, wyśrodkowane i zredukowane zmienne losowe są stale wykorzystywane zarówno w badaniach teoretycznych, jak i w algorytmach, produkty oprogramowania, dokumentacja normatywno-techniczna i instruktażowo-metodologiczna. W szczególności ze względu na równości umożliwiają uproszczenie uzasadniania metod, formułowania twierdzeń i wzorów obliczeniowych.

Stosowane są transformacje zmiennych losowych i nie tylko Ogólny plan. Więc jeśli Y = topór + B, Gdzie A I B– w takim razie kilka liczb

Przykład 7. Jeśli następnie Y jest zredukowaną zmienną losową, a wzory (8) zamieniają się we wzory (7).

Z każdą zmienną losową X możesz powiązać wiele zmiennych losowych Y, dane ze wzoru Y = topór + B w różnych A> 0 i B. Zestaw ten nazywa się rodzina z przesunięciem skali, generowane przez zmienną losową X. Funkcje dystrybucji F Y(X) stanowią rodzinę rozkładów z przesunięciem skali generowaną przez funkcję dystrybucji F(X). Zamiast Y = topór + B często korzystaj z nagrywania

Numer Z nazywa się parametrem przesunięcia i liczbą D- parametr skali. Pokazuje to wzór (9). X– wynik pomiaru określonej wielkości – trafia do U– wynik pomiaru tej samej wielkości w przypadku przesunięcia początku pomiaru do punktu Z, a następnie użyj nowej jednostki miary, in D razy większy od starego.

W przypadku rodziny przesunięć skali (9) rozkład X nazywany jest standardem. W probabilistycznych statystycznych metodach podejmowania decyzji i innych badaniach stosowanych stosuje się standardowy rozkład normalny, standardowy rozkład Weibulla-Gnedenko, standardowy rozkład gamma itp. (patrz poniżej).

Stosuje się także inne transformacje zmiennych losowych. Na przykład dla dodatniej zmiennej losowej X rozważają Y= log X, gdzie lg X– logarytm dziesiętny liczby X. Łańcuch równości

F Y (x) = P( lg X< x) = P(X < 10x) = F( 10X)

łączy funkcje dystrybucyjne X I Y.

Podczas przetwarzania danych wykorzystywane są następujące cechy zmiennej losowej X jako chwile porządku Q, tj. oczekiwania matematyczne wobec zmiennej losowej Xk, Q= 1, 2, ... Zatem samo oczekiwanie matematyczne jest momentem porządku 1. W przypadku dyskretnej zmiennej losowej momentem porządku Q można obliczyć jako

Dla ciągłej zmiennej losowej

Chwile porządku Q zwane także początkowymi momentami porządku Q, w przeciwieństwie do pokrewnych cech - centralne momenty porządku Q, podane przez wzór

Zatem dyspersja jest centralnym momentem rzędu 2.

Rozkład normalny i centralne twierdzenie graniczne. W probabilistyczno-statystycznych metodach podejmowania decyzji często mówimy o rozkładzie normalnym. Czasami próbują go wykorzystać do modelowania rozkładu danych wyjściowych (próby te nie zawsze są uzasadnione - patrz poniżej). Co ważniejsze, wiele metod przetwarzania danych opiera się na tym, że obliczone wartości mają rozkłady zbliżone do normalnych.

Pozwalać X 1 , X 2 ,…, Xn M(X ja) = M i odchylenia D(X ja) = , I = 1, 2,…, N,... Jak wynika z wyników poprzedniego rozdziału,

Rozważ zredukowaną zmienną losową U n za kwotę , mianowicie,

Jak wynika ze wzorów (7), M(U n) = 0, D(U n) = 1.

(dla warunków o identycznym rozkładzie). Pozwalać X 1 , X 2 ,…, Xn, … – niezależne zmienne losowe o jednakowym rozkładzie z oczekiwaniami matematycznymi M(X ja) = M i odchylenia D(X ja) = , I = 1, 2,…, N,... Wtedy dla dowolnego x istnieje granica

Gdzie F(x)– funkcja standardowego rozkładu normalnego.

Więcej o tej funkcji F(x) – poniżej (czytaj „fi od x”, ponieważ F- Grecka wielka litera „phi”).

Centralne twierdzenie graniczne (CLT) ma swoją nazwę, ponieważ jest centralnym, najczęściej używanym matematycznym wynikiem teorii prawdopodobieństwa i statystyki matematycznej. Historia CLT trwa około 200 lat – od roku 1730, kiedy to angielski matematyk A. Moivre (1667-1754) opublikował pierwszy wynik związany z CLT (patrz niżej o twierdzeniu Moivre’a-Laplace’a), aż do lat dwudziestych i trzydziestych XX wieku XX wieku, kiedy Finn J.W. Lindeberg, Francuz Paul Levy (1886-1971), Jugosławia V. Feller (1906-1970), Rosjanin A.Ya. Khinchin (1894-1959) i inni naukowcy uzyskali warunki konieczne i wystarczające dla ważności klasycznego centralnego twierdzenia granicznego.

Na tym nie zakończył się rozwój rozważanego tematu - badano zmienne losowe, które nie mają dyspersji, tj. ci dla kogo

(akademik B.V. Gnedenko i inni), sytuacja, gdy sumowane są zmienne losowe (a dokładniej elementy losowe) o bardziej złożonym charakterze niż liczby (naukowcy Yu.V. Prochorow, A.A. Borovkov i ich współpracownicy) itp. .d.

Funkcja dystrybucyjna F(x) jest dana przez równość

,

gdzie jest gęstość standardowego rozkładu normalnego, który ma dość złożone wyrażenie:

.

Tutaj =3,1415925... jest liczbą znaną w geometrii, równą stosunkowi obwodu do średnicy, mi = 2,718281828... - podstawa logarytmów naturalnych (aby zapamiętać tę liczbę, należy pamiętać, że rok 1828 to rok urodzenia pisarza L.N. Tołstoja). Jak wiadomo z Analiza matematyczna,

Podczas przetwarzania wyników obserwacji funkcja rozkładu normalnego nie jest obliczana przy użyciu podanych wzorów, ale jest znajdowana przy użyciu specjalnych tabel lub programy komputerowe. Najlepsze „Tabele statystyki matematycznej” w języku rosyjskim zostały opracowane przez odpowiednich członków Akademii Nauk ZSRR L.N. Bolszewa i N.V. Smirnowa.

Postać gęstości standardowego rozkładu normalnego wynika z teorii matematycznej, której nie możemy tutaj rozważać, a także z dowodu CLT.

Dla ilustracji podajemy małe tabele funkcji rozkładu F(x)(Tabela 2) i jej kwantyle (Tabela 3). Funkcjonować F(x) symetryczny wokół 0, co przedstawiono w tabeli 2-3.

Tabela 2.

Standardowa funkcja rozkładu normalnego.

Jeśli zmienna losowa X ma funkcję dystrybucji F(x), To M(X) = 0, D(X) = 1. To stwierdzenie zostało udowodnione w teorii prawdopodobieństwa w oparciu o rodzaj gęstości prawdopodobieństwa. Jest to spójne z podobnym stwierdzeniem dotyczącym charakterystyki zredukowanej zmiennej losowej U n, co jest całkiem naturalne, gdyż CLT stwierdza, że ​​przy nieograniczonym wzroście liczby wyrazów funkcja dystrybucji U n zmierza do standardowej funkcji rozkładu normalnego F(x), i dla każdego X.

Tabela 3.

Kwantyle standardowego rozkładu normalnego.

Wprowadźmy pojęcie rodziny rozkładów normalnych. Z definicji rozkład normalny to rozkład zmiennej losowej X, dla którego rozkład zredukowanej zmiennej losowej wynosi F(x). Jak wynika z ogólnych właściwości rodzin rozkładów z przesunięciem skali (patrz wyżej), rozkład normalny jest rozkładem zmiennej losowej

Gdzie X– zmienna losowa z rozkładem F(X), I M = M(Y), = D(Y). Rozkład normalny z parametrami przesunięcia M i skala jest zwykle wskazana N(M, ) (czasami używana jest notacja N(M, ) ).

Jak wynika z (8), gęstość prawdopodobieństwa rozkładu normalnego N(M, ) Jest

Rozkłady normalne tworzą rodzinę z przesunięciem skali. W tym przypadku parametrem skali jest D= 1/ i parametr przesunięcia C = - M/ .

Dla momentów centralnych trzeciego i czwartego rzędu rozkładu normalnego obowiązują równości:

Równości te stanowią podstawę klasycznych metod sprawdzania, czy obserwacje mają rozkład normalny. Obecnie zwykle zaleca się badanie normalności za pomocą kryterium W Shapiro – Wilka. Problem testowania normalności omówiono poniżej.

Jeśli zmienne losowe X 1 I X2 mają funkcje dystrybucyjne N(M 1 , 1) I N(M 2 , 2) zatem odpowiednio X 1+ X2 ma dystrybucję Dlatego też, jeśli zmienne losowe X 1 , X 2 ,…, Xn N(M, ) , a następnie ich średnia arytmetyczna

ma dystrybucję N(M, ) . Te właściwości rozkładu normalnego są stale wykorzystywane w różnych probabilistycznych i statystycznych metodach podejmowania decyzji, w szczególności w statystycznej regulacji procesów technologicznych oraz w statystycznej kontroli akceptacji opartej na kryteriach ilościowych.

Korzystając z rozkładu normalnego, zdefiniowano trzy rozkłady, które są obecnie często stosowane w przetwarzaniu danych statystycznych.

Rozkład (chi – kwadrat) – rozkład zmiennej losowej

gdzie są zmienne losowe X 1 , X 2 ,…, Xn niezależne i mają ten sam rozkład N(0,1). W tym przypadku liczba terminów, tj. N, nazywana jest „liczbą stopni swobody” rozkładu chi-kwadrat.

Dystrybucja T T Studenta jest rozkładem zmiennej losowej

gdzie są zmienne losowe U I X niezależny, U ma standardowy rozkład normalny N(0,1) i X– rozkład chi – kwadrat c N stopnie swobody. W której N nazywana jest „liczbą stopni swobody” rozkładu Studenta. Podział ten wprowadził w 1908 roku angielski statystyk W. Gosset, który pracował w fabryce piwa. Do obliczeń ekonomicznych i statystycznych wykorzystano probabilistyczne metody statystyczne rozwiązania techniczne w tej fabryce, dlatego jej kierownictwo zabroniło V. Gossetowi publikować artykuły naukowe pod własnym nazwiskiem. W ten sposób był chroniony tajemnica handlowa, „know-how” w postaci metod probabilistycznych i statystycznych opracowanych przez V. Gosseta. Miał jednak okazję publikować pod pseudonimem „Student”. Historia Gosseta-Studenta pokazuje, że przez kolejne sto lat menadżerowie Wielkiej Brytanii byli świadomi tego, co wielkie wydajność ekonomiczna Probabilistyczne i statystyczne metody podejmowania decyzji.

Rozkład Fishera to rozkład zmiennej losowej

gdzie są zmienne losowe X 1 I X2 są niezależne i mają rozkład chi-kwadrat z liczbą stopni swobody k 1 I k 2 odpowiednio. W tym samym czasie para (k 1 , k 2 ) – parę „stopni swobody” rozkładu Fishera, a mianowicie: k 1 jest liczbą stopni swobody licznika, oraz k 2 – liczba stopni swobody mianownika. Rozkład zmiennej losowej F nosi imię wielkiego angielskiego statystyka R. Fishera (1890-1962), który aktywnie wykorzystywał go w swoich pracach.

Wyrażenia dla funkcji rozkładu chi-kwadrat, Studenta i Fishera, ich gęstości i charakterystyki, a także tabele można znaleźć w literaturze specjalistycznej (patrz na przykład).

Jak już wspomniano, rozkłady normalne są obecnie często stosowane w modelach probabilistycznych w różnych obszarach zastosowań. Jaki jest powód tak powszechnego rozpowszechnienia tej dwuparametrowej rodziny rozkładów? Wyjaśnia to następujące twierdzenie.

Centralne twierdzenie graniczne(dla terminów o różnym rozkładzie). Pozwalać X 1 , X 2 ,…, Xn,… - niezależne zmienne losowe z oczekiwaniami matematycznymi M(X 1 ), M(X 2 ),…, M(X n), ... i odchylenia D(X 1 ), D(X 2 ),…, D(X n), ... odpowiednio. Pozwalać

Następnie, jeśli spełnione są pewne warunki, które zapewniają niewielki udział któregokolwiek z terminów w U n,

dla kazdego X.

Nie będziemy tutaj formułować omawianych warunków. Można je znaleźć w literaturze specjalistycznej (patrz na przykład). „Wyjaśnienie warunków, w jakich działa CPT, jest zasługą wybitnych rosyjskich naukowców A.A. Markowa (1857–1922), a w szczególności A.M. Lapunowa (1857–1918).”

Centralne twierdzenie graniczne pokazuje, że w przypadku, gdy wynik pomiaru (obserwacja) kształtuje się pod wpływem wielu przyczyn, z których każda ma niewielki udział, i określa się wynik całkowity dodatkowo, tj.

dodając, to rozkład wyniku pomiaru (obserwacji) jest zbliżony do normalnego. X Czasami uważa się, że aby rozkład był normalny wystarczy, że wynik pomiaru (obserwacji) X powstaje pod wpływem wielu przyczyn, z których każdy ma niewielki wpływ. To jest źle. Ważne jest, jak te przyczyny działają. Jeśli dodatek, to ma w przybliżeniu rozkład normalny. Jeśli multiplikatywnie X(tj. działania poszczególnych przyczyn są mnożone, a nie dodawane), a następnie podział X zbliżone nie do normalnego, ale do tzw. logarytmicznie normalne, tj. Nie X, a log X ma w przybliżeniu rozkład normalny. Jeśli nie ma powodu sądzić, że działa jeden z tych dwóch mechanizmów kształtowania się wyniku końcowego (lub jakiś inny dobrze zdefiniowany mechanizm), to o dystrybucji

nic konkretnego nie da się powiedzieć.

Z powyższego wynika, że ​​w konkretnym zastosowanym problemie normalności wyników pomiarów (obserwacji) z reguły nie można ustalić na podstawie rozważań ogólnych, należy to sprawdzić za pomocą kryteriów statystycznych; Lub użyj nieparametrycznych metod statystycznych, które nie opierają się na założeniach dotyczących przynależności funkcji rozkładu wyników pomiarów (obserwacji) do tej lub innej rodziny parametrycznej. Oprócz rodziny rozkładów normalnych z przesunięciem skali, szeroko stosuje się wiele innych rodzin rozkładów - rozkłady lognormalne, wykładnicze, Weibulla-Gnedenki, rozkłady gamma. Przyjrzyjmy się tym rodzinom.

Losowa wartość X ma rozkład lognormalny, jeśli zmienna losowa Y= log X ma rozkład normalny. Następnie Z= log X = 2,3026…Y ma również rozkład normalny N(A 1 ,σ 1), gdzie ln X - naturalny logarytm X. Gęstość rozkładu lognormalnego wynosi:

Z centralnego twierdzenia granicznego wynika, że ​​iloczyn X = X 1 X 2 … Xn niezależne dodatnie zmienne losowe X ja, I = 1, 2,…, N, ogólnie N można przybliżyć za pomocą rozkładu logarytmiczno-normalnego. W szczególności multiplikatywny model formacji wynagrodzenie lub dochodu prowadzi do zalecenia, aby logarytmicznie przybliżać rozkład wynagrodzeń i dochodów normalne prawa. Dla Rosji zalecenie to okazało się uzasadnione – potwierdzają to dane statystyczne.

Istnieją inne modele probabilistyczne, które prowadzą do prawa lognormalnego. Klasyczny przykład takiego modelu podał A.N. Kołmogorow, który na podstawie fizycznego systemu postulatów doszedł do wniosku, że wielkość cząstek podczas kruszenia kawałków rudy, węgla itp. w młynach kulowych mają rozkład logarytmiczno-normalny.

Przejdźmy do innej rodziny rozkładów, szeroko stosowanej w różnych probabilistyczno-statystycznych metodach podejmowania decyzji i innych badaniach stosowanych – rodziny rozkładów wykładniczych. Zacznijmy od modelu probabilistycznego prowadzącego do takich rozkładów. Aby to zrobić, należy wziąć pod uwagę „przebieg zdarzeń”, tj. sekwencja zdarzeń następujących po sobie w określonych momentach czasu. Przykłady obejmują: przepływ połączeń w centrali telefonicznej; przepływ awarii urządzeń w łańcuchu technologicznym; przepływ usterek produktów podczas testowania produktów; przepływ wniosków klientów do oddziału banku; przepływ nabywców ubiegających się o towary i usługi itp. W teorii przepływów zdarzeń obowiązuje twierdzenie podobne do centralnego twierdzenia granicznego, tyle że nie chodzi tu o sumowanie zmiennych losowych, ale o sumowanie przepływów zdarzeń. Rozważamy całkowity przepływ złożony z duża liczba przepływy niezależne, z których żaden nie ma dominującego wpływu na przepływ całkowity. Na przykład strumień połączeń przychodzących do centrali telefonicznej składa się z dużej liczby niezależnych strumieni połączeń pochodzących od poszczególnych abonentów. Udowodniono, że w przypadku, gdy charakterystyki przepływów nie zależą od czasu, przepływ całkowity opisuje się w całości jedną liczbą – natężeniem przepływu. Dla całkowitego przepływu rozważ zmienną losową X- długość odstępu czasu pomiędzy kolejnymi zdarzeniami. Jego funkcja rozkładu ma postać

(10)

Rozkład ten nazywa się rozkładem wykładniczym, ponieważ wzór (10) zawiera funkcję wykładniczą mi -λ X. Wartość 1/λ jest parametrem skali. Czasami wprowadzany jest również parametr przesunięcia Z, rozkład zmiennej losowej nazywany jest wykładniczym X + s, gdzie dystrybucja X wyraża się wzorem (10).

Rozkłady wykładnicze są szczególnym przypadkiem tzw. Rozkłady Weibulla-Gnedenki. Ich nazwy pochodzą od nazwisk inżyniera V. Weibulla, który wprowadził te rozkłady do praktyki analizy wyników badań zmęczeniowych, oraz matematyka B.V. Gnedenko (1912-1995), który otrzymał takie rozkłady jako granice przy badaniu maksimum wyniki testu. Pozwalać X- zmienna losowa charakteryzująca czas funkcjonowania produktu, złożonego systemu, elementu (tj. zasobu, czasu działania do stanu granicznego itp.), czasu funkcjonowania przedsiębiorstwa lub życia istoty żywej itp. Intensywność awarii odgrywa ważną rolę

(11)

Gdzie F(X) I F(X) - funkcja rozkładu i gęstość zmiennej losowej X.

Opiszmy typowe zachowanie wskaźnika awaryjności. Cały przedział czasowy można podzielić na trzy okresy. Na pierwszym z nich funkcja λ(x) To ma wysokie wartości i wyraźną tendencję spadkową (najczęściej maleje monotonicznie). Można to wytłumaczyć obecnością w danej partii wyrobów jednostek o wadach oczywistych i ukrytych, które prowadzą do stosunkowo szybkiej awarii tych jednostek wyrobów. Pierwszy okres nazywany jest „okresem docierania” (lub „okresem docierania”). Tyle zazwyczaj obejmuje okres gwarancji.

Następnie następuje okres normalnej pracy, charakteryzujący się w przybliżeniu stałą i stosunkowo niską awaryjnością. Charakter awarii w tym okresie jest nagły (wypadki, błędy obsługi itp.) i nie jest zależny od czasu eksploatacji urządzenia.

Wreszcie ostatnim okresem eksploatacji jest okres starzenia i zużycia. Awarie występujące w tym okresie polegają na nieodwracalnych zmianach fizycznych, mechanicznych i chemicznych w materiałach, prowadzących do postępującego pogorszenia jakości jednostki produktu i jego ostatecznej awarii.

Każdy okres ma swój własny typ funkcji λ(x). Rozważmy klasę zależności mocy

λ(x) = λ 0bx b -1 , (12)

Gdzie λ 0 > 0 i B> 0 - niektóre parametry numeryczne. Wartości B < 1, B= 0 i B> 1 odpowiadają rodzajowi awaryjności odpowiednio w okresach docierania, normalnej pracy i starzenia.

Zależność (11) przy zadanym współczynniku awaryjności λ(x)- równanie różniczkowe funkcji F(X). Z teorii równania różniczkowe wynika z tego

(13)

Podstawiając (12) do (13) otrzymujemy to

(14)

Rozkład dany wzorem (14) nazywany jest rozkładem Weibulla – Gnedenko. Ponieważ

to ze wzoru (14) wynika, że ​​ilość A, dany wzorem (15), jest parametrem skali. Czasami wprowadzany jest także parametr shift, tj. Nazywa się dystrybuanty Weibulla-Gnedenki F(X - C), Gdzie F(X) jest określony wzorem (14) dla niektórych λ 0 i B.

Gęstość rozkładu Weibulla-Gnedenko ma postać

(16)

Gdzie A> 0 - parametr skali, B> 0 - parametr formularza, Z- parametr zmiany. W tym przypadku parametr A ze wzoru (16) jest powiązany z parametrem λ 0 ze wzoru (14) przez zależność określoną we wzorze (15).

Rozkład wykładniczy jest bardzo szczególnym przypadkiem rozkładu Weibulla-Gnedenki, odpowiadającym wartości parametru kształtu B = 1.

Rozkład Weibulla-Gnedenki wykorzystuje się także przy konstruowaniu probabilistycznych modeli sytuacji, w których o zachowaniu obiektu decyduje „najsłabsze ogniwo”. Istnieje analogia z łańcuchem, którego bezpieczeństwo zależy od ogniwa o najmniejszej wytrzymałości. Inaczej mówiąc, niech X 1 , X 2 ,…, Xn- niezależne zmienne losowe o jednakowym rozkładzie,

X(1)=min( X 1, X 2,…, X rz), X(n)=maks.( X 1, X 2,…, X rz).

W wielu stosowanych problemach odgrywają one ważną rolę X(1) I X(N) w szczególności przy badaniu maksymalnych możliwych wartości („zapisów”) niektórych wartości, na przykład płatności ubezpieczeniowych lub strat wynikających z ryzyka handlowego, podczas badania granic elastyczności i wytrzymałości stali, szeregu cech niezawodności itp. . Pokazano, że dla dużych n rozkładów X(1) I X(N) , są z reguły dobrze opisane rozkładami Weibulla-Gnedenki. Zasadniczy wkład w badanie rozkładów X(1) I X(N) wniesiony przez radzieckiego matematyka B.V. Gnedenko. Prace V. Weibulla, E. Gumbela, V.B. poświęcone są wykorzystaniu wyników uzyskanych w ekonomii, zarządzaniu, technologii i innych dziedzinach. Nevzorova, E.M. Kudlaev i wielu innych specjalistów.

Przejdźmy do rodziny rozkładów gamma. Znajdują szerokie zastosowanie w ekonomii i zarządzaniu, teorii i praktyce niezawodności oraz testowaniu, w różnych dziedzinach techniki, meteorologii itp. W szczególności w wielu sytuacjach rozkład gamma zależy od takich wielkości, jak całkowity okres użytkowania produktu, długość łańcucha przewodzących cząstek pyłu, czas, w którym produkt osiąga stan graniczny podczas korozji, czas pracy do k-ta odmowa, k= 1, 2, … itd. Oczekiwana długość życia pacjentów choroby przewlekłe, czas do osiągnięcia określonego efektu podczas leczenia w niektórych przypadkach ma rozkład gamma. Rozkład ten jest najwłaściwszy do opisu popytu w ekonomicznych i matematycznych modelach zarządzania zapasami (logistyka).

Gęstość rozkładu gamma ma postać

(17)

Gęstość prawdopodobieństwa we wzorze (17) wyznaczają trzy parametry A, B, C, Gdzie A>0, B>0. W której A jest parametrem formularza, B- parametr skali i Z- parametr zmiany. Czynnik 1/Γ(а) normalizuje się, do tego wprowadzono

Tutaj Γ(a)- jeden z tych stosowanych w matematyce funkcje specjalne, tzw. „funkcja gamma”, od której nazwany jest rozkład określony wzorem (17),

Na stałe A wzór (17) określa rodzinę rozkładów z przesunięciem skali generowaną przez rozkład z gęstością

(18)

Rozkład postaci (18) nazywany jest standardowym rozkładem gamma. Otrzymuje się go ze wzoru (17) w B= 1 i Z= 0.

Szczególny przypadek rozkładów gamma dla A= 1 to rozkłady wykładnicze (z λ = 1/B). Z naturalnym A I Z Rozkłady gamma =0 nazywane są rozkładami Erlanga. Z prac duńskiego naukowca K.A. Erlanga (1878–1929), pracownika kopenhaskiej firmy telefonicznej, który studiował w latach 1908–1922. funkcjonowania sieci telefonicznych rozpoczął się rozwój teorii kolejkowania. Teoria ta zajmuje się probabilistycznym i statystycznym modelowaniem systemów, w których obsługiwany jest przepływ żądań w celu podejmowania optymalnych decyzji. Rozkłady Erlanga są używane w tych samych obszarach zastosowań, w których używane są rozkłady wykładnicze. Opiera się to na następującym fakcie matematycznym: suma k niezależnych zmiennych losowych o rozkładzie wykładniczym o tych samych parametrach λ i Z, ma rozkład gamma z parametrem kształtu a =k, parametr skali B= 1/λ i parametr przesunięcia kc. Na Z= 0 otrzymujemy rozkład Erlanga.

Jeśli zmienna losowa X ma rozkład gamma z parametrem kształtu A takie, że D = 2 A- liczba całkowita, B= 1 i Z= 0, następnie 2 X ma rozkład chi-kwadrat z D stopnie swobody.

Losowa wartość X z dystrybucją gvmma ma następujące cechy:

Wartość oczekiwana M(X) =ok + C,

Zmienność D(X) = σ 2 = ok 2 ,

Współczynnik zmienności

Asymetria

Nadmiar

Rozkład normalny jest skrajnym przypadkiem rozkładu gamma. Dokładniej, niech Z będzie zmienną losową o standardowym rozkładzie gamma określonym wzorem (18). Następnie

dla dowolnej liczby rzeczywistej X, Gdzie F(x)- standardowa funkcja rozkładu normalnego N(0,1).

W badaniach stosowanych wykorzystuje się także inne parametryczne rodziny rozkładów, z których najbardziej znane to układ krzywych Pearsona, szeregi Edgewortha i Charliera. Nie są one tutaj brane pod uwagę.

Oddzielny rozkłady stosowane w probabilistycznych i statystycznych metodach podejmowania decyzji. Najczęściej stosowane są trzy rodziny rozkładów dyskretnych – dwumianowy, hipergeometryczny i Poissona, a także kilka innych rodzin – geometryczna, ujemna dwumianowa, wielomianowa, ujemna hipergeometryczna itp.

Jak już wspomniano, rozkład dwumianowy występuje w niezależnych próbach, z których w każdej z prawdopodobieństwem R pojawia się wydarzenie A. Jeśli Łączna testy N podane, następnie liczba testów Y, w którym pojawiło się wydarzenie A, ma rozkład dwumianowy. W przypadku rozkładu dwumianowego prawdopodobieństwo uznania go za zmienną losową wynosi Y wartości y określa się na podstawie wzoru

Liczba kombinacji N elementy wg y, znane z kombinatoryki. Dla wszystkich y, z wyjątkiem 0, 1, 2, …, N, mamy P(Y= y)= 0. Rozkład dwumianowy ze stałą wielkością próby N jest określony przez parametr P, tj. rozkłady dwumianowe tworzą rodzinę jednoparametrową. Wykorzystuje się je w analizie danych z badań reprezentacyjnych, w szczególności w badaniu preferencji konsumentów, selektywnej kontroli jakości produktu według jednostopniowych planów kontroli, przy badaniu populacji jednostek w demografii, socjologii, medycynie, biologii itp. .

Jeśli Y 1 I Y 2 - niezależne dwumianowe zmienne losowe o tym samym parametrze P 0 , określone na podstawie próbek o objętościach N 1 I N 2 zatem odpowiednio Y 1 + Y 2 - dwumianowa zmienna losowa o rozkładzie (19) z R = P 0 I N = N 1 + N 2 . Uwaga ta rozszerza stosowalność rozkładu dwumianowego, umożliwiając łączenie wyników kilku grup testów, gdy istnieją podstawy, aby sądzić, że ten sam parametr odpowiada wszystkim tym grupom.

Charakterystyki rozkładu dwumianowego obliczono wcześniej:

M(Y) = n.p., D(Y) = n.p.( 1- P).

W części „Zdarzenia i prawdopodobieństwa” udowodniono prawo wielkich liczb dla dwumianowej zmiennej losowej:

dla kazdego . Korzystając z centralnego twierdzenia granicznego, prawo wielkich liczb można udoskonalić, wskazując, ile Y/ N różni się od R.

Twierdzenie de Moivre’a-Laplace’a. Dla dowolnych liczb a i B, A< B, mamy

Gdzie F(X) jest funkcją standardowego rozkładu normalnego z oczekiwaniem matematycznym 0 i wariancją 1.

Aby to udowodnić, wystarczy skorzystać z reprezentacji Y w postaci sumy niezależnych zmiennych losowych odpowiadających wynikom poszczególnych testów, wzory na M(Y) I D(Y) oraz centralne twierdzenie graniczne.

To twierdzenie dotyczy przypadku R= ½ udowodnił angielski matematyk A. Moivre (1667-1754) w 1730 r. W powyższym sformułowaniu udowodnił to w 1810 r. francuski matematyk Pierre Simon Laplace (1749 - 1827).

Rozkład hipergeometryczny zachodzi podczas selektywnego sterowania skończonym zbiorem obiektów o objętości N według alternatywnego kryterium. Każdy kontrolowany obiekt jest klasyfikowany jako posiadający dany atrybut A lub jako nieposiadający tej cechy. Rozkład hipergeometryczny ma zmienną losową Y, równa liczbie obiekty posiadające daną cechę A w losowej próbce objętości N, Gdzie N< N. Na przykład liczba Y wadliwe jednostki produktu w losowej próbce objętościowej N od objętości partii N ma rozkład hipergeometryczny jeśli N< N. Innym przykładem jest loteria. Niech znak A bilet jest oznaką „bycia zwycięzcą”. Niech całkowita liczba biletów N, i pewna osoba nabyta N z nich. Wówczas liczba zwycięskich kuponów dla tej osoby ma rozkład hipergeometryczny.

Dla rozkładu hipergeometrycznego prawdopodobieństwo, że zmienna losowa Y przyjmie wartość y, ma postać

(20)

Gdzie D– liczba obiektów posiadających dany atrybut A, w rozważanym zestawie objętości N. W której y przyjmuje wartości od max(0, N - (N - D)) do min( N, D), inne rzeczy y prawdopodobieństwo we wzorze (20) jest równe 0. Zatem rozkład hipergeometryczny wyznaczają trzy parametry - objętość populacji N, liczba obiektów D w nim, posiadająca daną cechę A i wielkość próbki N.

Proste losowe próbkowanie objętości N od całkowitej objętości N to próbka otrzymana w wyniku losowego doboru, w którym dowolny ze zbiorów N obiekty mają takie samo prawdopodobieństwo wybrania. Metody losowego doboru próby respondentów (rozmówców) lub jednostek towarów na sztuki omówione są w dokumentach instruktażowych, metodologicznych i regulacyjnych. Jedna z metod selekcji jest następująca: obiekty wybiera się jeden po drugim i na każdym etapie każdy z pozostałych obiektów w zestawie ma taką samą szansę na wybranie. W literaturze na określenie rodzaju rozpatrywanych próbek stosowane są także określenia „próbka losowa” i „próbka losowa bez zwrotu”.

Ponieważ wielkość populacji (partia) N i próbki N są zwykle znane, wówczas parametrem rozkładu hipergeometrycznego, który należy oszacować, jest: D. W statystycznych metodach zarządzania jakością produktów D– zwykle liczba wadliwych jednostek w partii. Interesująca jest także charakterystyka dystrybucji D/ N– poziom usterek.

Dla rozkładu hipergeometrycznego

Ostatni czynnik w wyrażeniu wariancji jest bliski 1, jeśli N>10 N. Jeśli dokonasz wymiany P = D/ N, wówczas wyrażenia matematycznego oczekiwania i wariancji rozkładu hipergeometrycznego zamienią się w wyrażenia matematycznego oczekiwania i wariancji rozkładu dwumianowego. To nie przypadek. Można to wykazać

Na N>10 N, Gdzie P = D/ N. Stosunek ograniczający obowiązuje

i tę ograniczającą relację można zastosować kiedy N>10 N.

Trzecią powszechnie stosowaną dystrybucją dyskretną jest rozkład Poissona. Zmienna losowa Y ma rozkład Poissona jeśli

,

gdzie λ jest parametrem rozkładu Poissona, oraz P(Y= y)= 0 dla wszystkich pozostałych y(dla y=0 oznacza się 0! =1). Dla rozkładu Poissona

M(Y) = λ, D(Y) = λ.

Rozkład ten został nazwany na cześć francuskiego matematyka S. D. Poissona (1781-1840), który uzyskał go po raz pierwszy w 1837 r. Rozkład Poissona jest granicznym przypadkiem rozkładu dwumianowego, gdy prawdopodobieństwo R realizacja zdarzenia jest niewielka, ale ilość testów Nświetnie, i n.p.= λ. Dokładniej, obowiązuje relacja graniczna

Dlatego rozkład Poissona (w starej terminologii „prawo dystrybucji”) często nazywany jest także „prawem rzadkich zdarzeń”.

Rozkład Poissona pojawia się w teorii strumieni zdarzeń (patrz wyżej). Udowodniono, że dla najprostszego przepływu o stałym natężeniu Λ liczba zdarzeń (wywołań), które wystąpiły w czasie T, ma rozkład Poissona z parametrem λ = Λ T. Dlatego prawdopodobieństwo, że w tym czasie Tżadne zdarzenie nie nastąpi, równe mi - Λ T, tj. funkcja rozkładu długości odstępu między zdarzeniami jest wykładnicza.

Rozkład Poissona wykorzystuje się do analizy wyników przykładowych badań marketingowych konsumentów, obliczania charakterystyki operacyjnej statystycznych planów kontroli akceptacji w przypadku małych wartości poziomu akceptacji wad, do opisu liczby awarii kontrolowanego statystycznie procesu technologicznego w jednostce czasu, liczby „zapotrzebowań na usługi” otrzymanych w jednostce czasu w systemie kolejkowym, statystycznych wzorców wypadków i rzadkie choroby itp.

Opis innych rodzin parametrycznych rozkładów dyskretnych i ich możliwości praktyczne użycie są rozpatrywane w literaturze.

W niektórych przypadkach, na przykład podczas badania cen, wielkości produkcji lub całkowitego czasu między awariami w problemach z niezawodnością, funkcje rozkładu są stałe w pewnych przedziałach, w których wartości badanych zmiennych losowych nie mogą spaść.

Dystrybuantą zmiennej losowej X jest funkcja F(x), która dla każdego x wyraża prawdopodobieństwo, że zmienna losowa X przyjmie wartość, mniejszy x

Przykład 2.5. Biorąc pod uwagę szereg rozkładów zmiennej losowej

Znajdź i graficznie przedstaw jego funkcję rozkładu. Rozwiązanie. Zgodnie z definicją

F(jc) = 0 w X X

F(x) = 0,4 + 0,1 = 0,5 przy 4 F(x) = 0,5 + 0,5 = 1 przy X > 5.

Zatem (patrz ryc. 2.1):

Własności funkcji rozkładu:

1. Dystrybucja zmiennej losowej jest funkcją nieujemną mieszczącą się w przedziale od zera do jedynki:

2. Funkcja rozkładu zmiennej losowej jest funkcją niemalejącą na całej osi liczbowej, tj. Na X 2 >x

3. Przy minus nieskończoności funkcja rozkładu jest równa zeru, przy plus nieskończoności jest równa jeden, tj.