W praktyce statystycznej, jednym z najważniejszych zagadnień jest określenie przedziałów ufności dla parametrów rozkładów. W kontekście rozkładu normalnego, istnieje kilka metod, które pozwalają oszacować przedział, w jakim znajduje się prawdziwa wartość średniej lub wariancji populacji. W tym rozdziale omówimy, jak określić przedział ufności dla średniej i wariancji rozkładu normalnego, zarówno wtedy, gdy wariancja jest znana, jak i gdy jest nieznana.

Zaczynając od przypadku, w którym wariancja populacji jest znana, obliczenie przedziału ufności dla średniej jest stosunkowo proste. Zakłada się, że mamy próbę o wielkości nn, a rozkład próby jest normalny z nieznanym parametrem średniej μ\mu i znaną wariancją σ2\sigma^2. Dla danego poziomu ufności γ\gamma, przedział ufności dla średniej jest obliczany według wzoru:

P(μkσnxˉμ+kσn)=γP\left( \mu - k \cdot \frac{\sigma}{\sqrt{n}} \leq \bar{x} \leq \mu + k \cdot \frac{\sigma}{\sqrt{n}} \right) = \gamma

gdzie kk jest wartością z rozkładu normalnego, zależną od poziomu ufności γ\gamma, a xˉ\bar{x} to średnia próby.

W praktyce najczęściej nie znamy wariancji populacji σ2\sigma^2, a zamiast tego mamy do dyspozycji próbkę z obliczoną estymatą wariancji próbki s2s^2. W takim przypadku zamiast wartości z rozkładu normalnego, używamy rozkładu t-Studenta, który uwzględnia zmienność wynikającą z faktu, że wariancja jest estymowana na podstawie próbki. Wzór na przedział ufności w tym przypadku ma postać:

P(μtn1,γsnxˉμ+tn1,γsn)=γP\left( \mu - t_{n-1, \gamma} \cdot \frac{s}{\sqrt{n}} \leq \bar{x} \leq \mu + t_{n-1, \gamma} \cdot \frac{s}{\sqrt{n}} \right) = \gamma

gdzie tn1,γt_{n-1, \gamma} to wartość krytyczna z rozkładu t-Studenta dla n1n-1 stopni swobody i poziomu ufności γ\gamma, a ss to odchylenie standardowe próbki.

Przykład 1: Załóżmy, że mamy próbkę pięciu niezależnych pomiarów punktu zapłonu oleju napędowego (D-2) o wartościach: 144, 147, 146, 142 i 144. Celem jest obliczenie 99% przedziału ufności dla średniej. Pierwszym krokiem jest wybranie poziomu ufności γ=0.99\gamma = 0.99. Z tabeli rozkładu t-Studenta, dla n1=4n-1 = 4 stopni swobody, otrzymujemy wartość t4,0.99=4.60t_{4, 0.99} = 4.60. Następnie obliczamy średnią próbki xˉ=144.6\bar{x} = 144.6 i wariancję s2=3.8s^2 = 3.8. Ostatecznie, obliczamy wartość k=t4,0.99sn=4.6013.85=4.01k = t_{4, 0.99} \cdot \frac{s}{\sqrt{n}} = 4.60 \cdot \frac{13.8}{\sqrt{5}} = 4.01. Przedział ufności to:

[144.64.01,144.6+4.01]=[140.5,148.7]\left[ 144.6 - 4.01, 144.6 + 4.01 \right] = \left[ 140.5, 148.7 \right]

Zauważmy, że przy założeniu, że wariancja populacji σ2\sigma^2 jest znana i wynosi 3.8, przedział ufności jest krótszy: [142.35,146.85]\left[ 142.35, 146.85 \right]. Dla małych prób różnica ta jest istotna.

Inny ważny aspekt dotyczy obliczania przedziałów ufności dla wariancji rozkładu normalnego, szczególnie gdy nie znamy średniej populacji. W takim przypadku wykorzystuje się rozkład chi-kwadrat, który jest odpowiedni do estymacji wariancji. Aby obliczyć przedział ufności dla wariancji, postępujemy zgodnie z następującymi krokami:

  1. Wybieramy poziom ufności γ\gamma (np. 95% lub 99%).

  2. Określamy wartości c1c_1 i c2c_2 z tabeli rozkładu chi-kwadrat dla n1n-1 stopni swobody.

  3. Obliczamy wariancję próbki s2s^2.

  4. Obliczamy granice przedziału ufności dla wariancji:

P((n1)s2c2σ2(n1)s2c1)=γP\left( \frac{(n-1)s^2}{c_2} \leq \sigma^2 \leq \frac{(n-1)s^2}{c_1} \right) = \gamma

Dla przykładu, rozważmy próbkę 13 wyników (np. wytrzymałość na rozciąganie blachy stalowej w kg/mm²): 89, 84, 87, 81, 89, 86, 91, 90, 78, 89, 87, 99, 83, 89. Chcemy obliczyć 95% przedział ufności dla wariancji. Z tabeli rozkładu chi-kwadrat dla n1=12n-1 = 12 stopni swobody, otrzymujemy wartości c1=5.01c_1 = 5.01 i c2=24.74c_2 = 24.74. Obliczamy wariancję próbki s2=9.05s^2 = 9.05. Ostatecznie przedział ufności dla wariancji wynosi:

P(129.0524.74σ2129.055.01)=[13.21,65.25]P\left( \frac{12 \cdot 9.05}{24.74} \leq \sigma^2 \leq \frac{12 \cdot 9.05}{5.01} \right) = \left[ 13.21, 65.25 \right]

Warto zauważyć, że przedział ufności dla wariancji jest dość szeroki, co sugeruje, że próbka nie jest wystarczająco duża, by uzyskać precyzyjne oszacowanie.

Wszystkie te obliczenia opierają się na założeniu, że dane pochodzą z rozkładu normalnego, a wyniki obliczeń mogą różnić się w zależności od wielkości próby i poziomu ufności. W praktyce, dla małych prób, wartość przedziału ufności może być znacznie szersza, a dla dużych prób staje się coraz bardziej precyzyjna.

Jak obliczyć rozwiązania równań różniczkowych przy użyciu przekształcenia Laplace'a?

Metoda przekształceń Laplace’a jest jednym z najpotężniejszych narzędzi do rozwiązywania równań różniczkowych. Pomaga ona sprowadzić trudne problemy do postaci algebraicznych, które są łatwiejsze do rozwiązania, a następnie przywrócić wyniki do przestrzeni czasowej. Główne zalety tej metody to możliwość łatwego uwzględnienia warunków początkowych oraz radzenie sobie z układami złożonymi, których prawe strony są funkcjami skokowymi lub impulsowymi. W tym rozdziale przedstawimy kilka przykładów, które ilustrują, jak przekształcenie Laplace’a może zostać zastosowane do rozwiązania równań różniczkowych z warunkami początkowymi.

Przykład 1: Rozwiązanie problemu początkowego z wykorzystaniem przekształcenia Laplace'a

Rozważmy zadanie:

yy=t,y(0)=1,y(0)=1y'' - y = t, \quad y(0) = 1, \quad y'(0) = 1
Kroki rozwiązania są następujące:

  1. Przekształcamy równanie różniczkowe przy użyciu przekształcenia Laplace'a. Z tabeli przekształceń (patrz tabela 6.1) wiemy, że:

    L{y(t)}=s2Y(s)sy(0)y(0)\mathcal{L}\{y''(t)\} = s^2 Y(s) - sy(0) - y'(0)

    Podstawiając wartości początkowe:

    (s21)Y(s)=s+1s2(s^2 - 1)Y(s) = \frac{s+1}{s^2}
  2. Po przekształceniu uzyskujemy funkcję przenoszenia:

    Q(s)=1s21Q(s) = \frac{1}{s^2 - 1}

    i rozwiązanie w przestrzeni ss-owej:

    Y(s)=s+1s2(s21)Y(s) = \frac{s + 1}{s^2 (s^2 - 1)}
  3. Przechodzimy teraz do przestrzeni czasowej, obliczając odwrotne przekształcenie Laplace'a. Na podstawie tabeli przekształceń uzyskujemy:

    y(t)=et+sinhtty(t) = e^t + \sinh t - t

    Wynik ten jest zgodny z rozwiązaniem, które uzyskujemy przy użyciu tradycyjnych metod rozwiązywania równań różniczkowych, ale w tym przypadku proces obliczeniowy jest znacznie uproszczony.

Przykład 2: Porównanie z tradycyjną metodą

Rozważmy inny problem początkowy:

y+y+9y=0,y(0)=0.16,y(0)=0y'' + y' + 9y = 0, \quad y(0) = 0.16, \quad y'(0) = 0
  1. Po przekształceniu przy użyciu Laplace'a otrzymujemy:

    (s2+s+9)Y(s)=0.16s(s^2 + s + 9)Y(s) = 0.16s
  2. Rozwiązanie algebraiczne daje:

    Y(s)=0.16(s1)(s2+s+9)Y(s) = \frac{0.16(s - 1)}{(s^2 + s + 9)}
  3. Następnie obliczamy odwrotne przekształcenie Laplace'a:

    y(t)=e0.5t(0.16cos(2.96t)+0.027sin(2.96t))y(t) = e^{ -0.5t}(0.16\cos(2.96t) + 0.027\sin(2.96t))

    To rozwiązanie jest zgodne z tym, które uzyskujemy za pomocą klasycznych metod, ale w tym przypadku, obliczenia zostały znacznie uproszczone.

Zaletami metody przekształcenia Laplace'a są:

  1. Rozwiązywanie równań niejednorodnych: Przekształcenie Laplace'a pozwala na rozwiązanie równań różniczkowych bez konieczności wcześniejszego rozwiązywania odpowiednich równań jednorodnych, co jest typowe dla tradycyjnych metod.

  2. Bezproblemowe uwzględnienie warunków początkowych: Warunki początkowe są automatycznie uwzględniane podczas przekształcania równania, co upraszcza proces obliczeniowy.

  3. Radzenie sobie z bardziej skomplikowanymi funkcjami wymuszeń: Metoda przekształcenia Laplace'a pozwala na łatwe rozwiązywanie równań z trudnymi funkcjami wymuszenia, jak np. impulsy, funkcje skokowe, czy funkcje okresowe bardziej złożone niż tylko sinus i cosinus.

Przykład 3: Problemy z danymi przesuniętymi w czasie

Kiedy warunki początkowe są określone w czasie różnym od zera, czyli dla t00t_0 \neq 0, zamiast klasycznego przekształcenia Laplace'a stosuje się tzw. przesunięcie w czasie. Rozważmy przykład:

y+y=2t,y(1)=1,y(1)=2y'' + y = 2t, \quad y(1) = 1, \quad y'(1) = 2

Aby rozwiązać ten problem, przesuwamy układ do czasu t0t_0. Dla t0=1t_0 = 1 ustawiamy t=t+1t = t' + 1, gdzie tt' jest nową zmienną czasową. Tak przekształcone równanie daje nowy problem początkowy, który możemy rozwiązać przy użyciu przekształcenia Laplace'a, uzyskując rozwiązanie w formie:

y(t)=2tsint+costy(t) = 2t - \sin t + \cos t

Ostatecznie po powrocie do oryginalnej zmiennej czasowej otrzymujemy wynik:

y(t)=2t+sin(t1)+cos(t1)y(t) = 2t + \sin(t - 1) + \cos(t - 1)

Metoda przekształceń Laplace’a znajduje szerokie zastosowanie w rozwiązywaniu równań różniczkowych z warunkami początkowymi i przesunięciami w czasie, szczególnie w inżynierii i naukach stosowanych, gdzie wiele problemów związanych z siłami wymuszającymi jest opisywanych przez funkcje skokowe lub impulsowe.