Wszechnica Wszechwiedzy - Baner

Rodzaje obciążeń w statyce – siły skupione, momenty i obciążenia rozłożone

Obciążeniem nazywamy zewnętrzne oddziaływanie na ciało lub konstrukcję. W zadaniach ze statyki może ono występować jako siła skupiona, moment skupiony albo obciążenie rozłożone o stałej lub zmiennej intensywności.

Rodzaj obciążenia decyduje o sposobie jego zapisu, jednostce oraz położeniu siły wypadkowej. Szczególnie ważne jest odróżnienie siły wyrażanej w niutonach od intensywności obciążenia podawanej w niutonach na metr oraz od momentu, którego jednostką jest niutonometr.

Podstawowe pojęcia związane z siłami, momentami i równowagą omówiliśmy wcześniej w artykule Czym zajmuje się statyka? Podstawowe pojęcia i zasady. W tym materiale skupimy się na klasyfikacji obciążeń oraz zastępowaniu obciążeń rozłożonych odpowiednimi siłami wypadkowymi.

Czym jest obciążenie?

Obciążenie jest zewnętrznym oddziaływaniem na rozpatrywane ciało albo układ mechaniczny. Może dążyć do jego przesunięcia, obrotu lub odkształcenia.

Do obciążeń zalicza się między innymi:

W zadaniach mechanicznych rzeczywiste oddziaływania zastępuje się odpowiednimi modelami. Ten sam rzeczywisty nacisk może zostać przedstawiony jako siła skupiona albo obciążenie rozłożone, zależnie od skali analizowanego zagadnienia i wymaganej dokładności.

Schematy siły skupionej, obciążenia równomiernego, obciążenia liniowo zmiennego i momentu skupionego
W płaskich układach statycznych stosuje się między innymi siły skupione, momenty skupione oraz obciążenia ciągłe o stałej lub zmiennej intensywności.

Siła skupiona

Siła skupiona jest idealizacją obciążenia przyłożonego w jednym punkcie. Na schemacie przedstawia się ją za pomocą pojedynczej strzałki.

W rzeczywistości siła zawsze działa na pewnej powierzchni. Można jednak potraktować ją jako skupioną, jeżeli obszar kontaktu jest niewielki w porównaniu z wymiarami analizowanego elementu.

Przykładami sił skupionych mogą być:

Siłę skupioną oznacza się często literami \(F\), \(P\) albo \(R\). Jej jednostką w układzie SI jest niuton:

\[ [F]=\mathrm{N}. \]

W zadaniach konstrukcyjnych najczęściej stosuje się kiloniutony:

\[ 1\ \mathrm{kN}=1000\ \mathrm{N}. \]

Moment skupiony

Moment skupiony opisuje oddziaływanie wywołujące obrót bez powstawania niezerowej siły wypadkowej. Na płaskim schemacie przedstawia się go zwykle za pomocą zakrzywionej strzałki.

Moment skupiony może powstać wskutek działania pary sił, czyli dwóch sił:

Suma tych sił jest równa zeru:

\[ \vec{F}+(-\vec{F})=\vec{0}, \]

ale pozostaje niezerowy moment:

\[ M=F a, \] 

gdzie \(a\) jest prostopadłą odległością między liniami działania sił.

W płaskim układzie obciążeń wektor momentu jest prostopadły do płaszczyzny rysunku. Symbol \(\odot\) oznacza wektor skierowany w stronę obserwatora, natomiast symbol \(\otimes\) — wektor skierowany w głąb rysunku.

Moment pary sił jest tak zwanym wektorem swobodnym. Można przenieść go do innego punktu ciała sztywnego bez zmiany zewnętrznego skutku mechanicznego.

\[ [M]=\mathrm{N\,m}. \]

Uwaga na jednostki: kiloniutonometr \(\mathrm{kN\,m}\) jest jednostką momentu. Nie należy go mylić z kiloniutonem na metr \(\mathrm{kN/m}\), który jest jednostką intensywności obciążenia liniowego.

Dwie równe i przeciwnie skierowane siły tworzące moment skupiony działający na belkę
Dwie równoległe siły o jednakowych wartościach i przeciwnych zwrotach tworzą parę sił o momencie \(M=F\cdot a\).

Obciążenia rozłożone

Obciążenie rozłożone działa na pewnej długości, powierzchni albo objętości. W przypadku belek najczęściej rozpatruje się obciążenie rozłożone wzdłuż odcinka.

Intensywność takiego obciążenia oznacza się zwykle symbolem \(q\) i podaje w jednostkach siły przypadającej na jednostkę długości:

\[ [q]=\frac{\mathrm{N}}{\mathrm{m}} \qquad\text{lub}\qquad [q]=\frac{\mathrm{kN}}{\mathrm{m}}. \]

Wartość \(q=5\ \mathrm{kN/m}\) oznacza, że na każdy metr długości elementu przypada obciążenie o wartości \(5\ \mathrm{kN}\).

Obciążenie liniowe może być:

Wypadkowa obciążenia rozłożonego jest równa polu pod wykresem jego intensywności. Jej linia działania przechodzi przez środek ciężkości tego pola.

Obciążenie ciągłe równomierne

Jeżeli intensywność obciążenia na całym odcinku jest jednakowa, otrzymujemy obciążenie ciągłe równomierne. Jego wykres ma kształt prostokąta.

Dla obciążenia o intensywności \(q\), działającego na długości \(l\), siła wypadkowa wynosi:

\[ R_q=q l. \] 

Wynik ma jednostkę siły:

\[ \frac{\mathrm{kN}}{\mathrm{m}}\cdot\mathrm{m} = \mathrm{kN}. \]

Wypadkowa działa w środku ciężkości prostokątnego wykresu obciążenia, czyli w połowie obciążonego odcinka:

\[ x_R=\frac{l}{2}. \]

Jeżeli obciążenie działa tylko na części belki, odległość \(l/2\) należy mierzyć od początku obciążonego odcinka, a nie od końca całej belki.

Obciążenie równomierne q zastąpione siłą wypadkową Rq działającą w połowie długości odcinka
Obciążenie równomierne o intensywności \(q\), działające na długości \(l\), zastępuje siła \(R_q=ql\) przyłożona w środku obciążonego odcinka.

Obciążenie trójkątne

Obciążenie trójkątne jest szczególnym przypadkiem obciążenia liniowo zmiennego. Jego intensywność zmienia się od zera na jednym końcu do wartości maksymalnej \(q_{\max}\) na drugim.

Siła wypadkowa jest równa polu trójkąta:

\[ R_q=\frac{q_{\max}l}{2}. \]

Wypadkowa przechodzi przez środek ciężkości trójkąta. Jej linia działania znajduje się:

Łatwy sposób zapamiętania: wypadkowa obciążenia trójkątnego leży bliżej tej strony, przy której strzałki obciążenia są dłuższe.

Obciążenie trójkątne zastąpione siłą wypadkową działającą jedną trzecią długości od strony większej intensywności
Wypadkowa obciążenia trójkątnego wynosi \(R_q=\frac{q_{\max}l}{2}\) i działa w odległości \(l/3\) od strony większej intensywności.

Obciążenie trapezowe

Jeżeli intensywność obciążenia zmienia się liniowo od wartości \(q_1\) do \(q_2\), a obie wartości są różne od zera, wykres obciążenia ma kształt trapezu.

Całkowita wypadkowa jest równa polu trapezu:

\[ R_q=\frac{q_1+q_2}{2}\,l. \]

Najwygodniej jest jednak rozłożyć obciążenie trapezowe na dwa prostsze obciążenia. Jeżeli \(q_2>q_1\), otrzymujemy:

Wypadkowa części prostokątnej wynosi:

\[ R_1=q_1l \]

i działa w odległości:

\[ x_1=\frac{l}{2} \]

od początku obciążonego odcinka.

Wypadkowa części trójkątnej wynosi:

\[ R_2=\frac{(q_2-q_1)l}{2} \]

i działa w odległości:

\[ x_2=\frac{2l}{3} \]

od strony intensywności \(q_1\), czyli w odległości \(l/3\) od strony większej intensywności \(q_2\).

Całkowita wypadkowa wynosi:

\[ R_q=R_1+R_2. \]

Jej położenie można wyznaczyć z warunku zgodności momentów:

\[ R_qx_R=R_1x_1+R_2x_2. \]

Stąd:

\[ x_R= \frac{R_1x_1+R_2x_2}{R_1+R_2}. \]

Po podstawieniu można otrzymać bezpośredni wzór na położenie wypadkowej, mierzone od strony \(q_1\):

\[ x_R= \frac{l(q_1+2q_2)} {3(q_1+q_2)}. \]

Jeżeli \(q_1>q_2\), rozkład można przeprowadzić od strony mniejszej wartości albo zastosować bezpośrednio wzór na położenie środka ciężkości trapezu.

Obciążenie trapezowe rozłożone na obciążenie prostokątne i trójkątne wraz z ich siłami wypadkowymi
Obciążenie trapezowe można rozłożyć na obciążenie równomierne oraz trójkątne, a następnie wyznaczyć wypadkowe obu części.

Dowolne obciążenie opisane funkcją

Jeżeli intensywność obciążenia zmienia się w dowolny sposób i jest opisana funkcją \(q(x)\), siłę wypadkową wyznacza się przez całkowanie:

\[ R_q=\int_a^b q(x)\,dx. \]

Całka ta oznacza pole pod wykresem intensywności obciążenia.

Położenie linii działania wypadkowej wyznacza się z warunku zgodności momentów:

\[ x_R= \frac{\displaystyle\int_a^b xq(x)\,dx} {\displaystyle\int_a^b q(x)\,dx}. \]

Jest to współrzędna środka ciężkości pola znajdującego się pod wykresem funkcji \(q(x)\).

Jeżeli obciążenia działają w przeciwnych kierunkach, należy przyjąć odpowiednią konwencję znaków i traktować pola algebraicznie.

Obciążenia powierzchniowe i objętościowe

Obciążenie może być rozłożone nie tylko wzdłuż linii, ale również na powierzchni lub w objętości ciała.

Obciążenie powierzchniowe

Intensywność obciążenia powierzchniowego oznacza się często symbolem \(p\). Jego jednostką jest:

\[ [p]=\frac{\mathrm{N}}{\mathrm{m}^2}. \]

Jednostka \(\mathrm{N/m^2}\) jest równa paskalowi:

\[ 1\ \mathrm{Pa}=1\ \frac{\mathrm{N}}{\mathrm{m}^2}. \]

Przykładami są parcie cieczy, napór wiatru oraz obciążenie śniegiem działające na powierzchnię dachu.

Jeżeli równomierne obciążenie powierzchniowe \(p\) działa na pas o szerokości \(b\), można przeliczyć je na obciążenie liniowe:

\[ q=pb. \]

Kontrola jednostek daje:

\[ \frac{\mathrm{kN}}{\mathrm{m}^2}\cdot\mathrm{m} = \frac{\mathrm{kN}}{\mathrm{m}}. \]

Obciążenie objętościowe

Obciążenie objętościowe działa w całej objętości materiału. Jego intensywność wyraża się w jednostkach:

\[ \frac{\mathrm{N}}{\mathrm{m}^3} \qquad\text{lub}\qquad \frac{\mathrm{kN}}{\mathrm{m}^3}. \]

Najczęstszym przykładem jest ciężar własny materiału.

Ciężar własny konstrukcji

Ciężar własny jest obciążeniem wynikającym z masy analizowanego elementu. Dla pręta o stałym polu przekroju \(A\), wykonanego z materiału o gęstości \(\rho\), intensywność liniowa ciężaru wynosi:

\[ q_g=\rho A g, \] 

gdzie \(g\) jest przyspieszeniem ziemskim.

Można również posługiwać się ciężarem objętościowym materiału:

\[ \gamma=\rho g. \] 

Wtedy:

\[ q_g=\gamma A. \] 

Jeżeli pole przekroju i gęstość są stałe, ciężar własny pręta można traktować jako obciążenie ciągłe równomierne.

Obciążenia ruchome

Obciążenie ruchome zmienia swoje położenie względem konstrukcji. Przykładem jest pojazd poruszający się po moście, suwnica albo wózek przemieszczający się po belce.

Obciążenie od pociągu można modelować jako układ sił skupionych odpowiadających naciskom poszczególnych osi. Jeżeli osi jest bardzo wiele i znajdują się blisko siebie, w uproszczonym modelu można zastąpić je obciążeniem rozłożonym.

Model zależy od celu obliczeń: lokomotywa i wagony nie stają się z natury różnymi rodzajami obciążenia. Zarówno naciski osi lokomotywy, jak i osi wagonów mogą być przedstawiane jako siły skupione. Obciążenie ciągłe jest jedynie możliwym uproszczeniem całego układu wielu nacisków.

Oddziaływania temperatury

Zmiana temperatury nie jest siłą ani momentem, ale może powodować zmianę wymiarów elementu. Swobodne wydłużenie pręta wynosi:

\[ \Delta L=\alpha\,\Delta T\,L, \]

gdzie:

Jeżeli element może wydłużyć się swobodnie, sama zmiana temperatury nie wywołuje reakcji podporowych. Jeżeli jednak jego przemieszczenie zostanie zablokowane, pojawią się siły i naprężenia termiczne.

Z tego powodu mosty i długie konstrukcje wyposaża się w łożyska przesuwne oraz szczeliny dylatacyjne. Umożliwiają one zmianę długości bez wywoływania nadmiernych sił wewnętrznych.

Jednostki i kontrola wymiarów

Kontrola jednostek jest jednym z najprostszych sposobów wykrywania błędów w obliczeniach ze statyki.

WielkośćTypowa jednostkaPrzykładowe oznaczenie
Siła skupiona\(\mathrm{N}\), \(\mathrm{kN}\)\(F\), \(P\), \(R\)
Moment skupiony\(\mathrm{N\,m}\), \(\mathrm{kN\,m}\)\(M\)
Obciążenie liniowe\(\mathrm{N/m}\), \(\mathrm{kN/m}\)\(q\)
Obciążenie powierzchniowe\(\mathrm{N/m^2}\), \(\mathrm{kN/m^2}\)\(p\)
Obciążenie objętościowe\(\mathrm{N/m^3}\), \(\mathrm{kN/m^3}\)\(\gamma\)

Obciążenie liniowe pomnożone przez długość daje siłę:

\[ q l \quad\Longrightarrow\quad \frac{\mathrm{kN}}{\mathrm{m}}\cdot\mathrm{m} = \mathrm{kN}. \]

Obciążenie liniowe pomnożone przez kwadrat długości ma jednostkę momentu:

\[ q l^2 \quad\Longrightarrow\quad \frac{\mathrm{kN}}{\mathrm{m}}\cdot\mathrm{m}^2 = \mathrm{kN\,m}. \]

Jeżeli w zadaniu intensywność obciążenia oznaczono symbolem \(q\), a długość symbolem \(l\), to:

W niektórych zadaniach intensywność obciążenia bywa oznaczana wielką literą \(Q\). Zawsze trzeba jednak sprawdzić definicję symbolu. Jeżeli \(Q\) jest już siłą, a nie siłą przypadającą na jednostkę długości, zależności wymiarowe będą inne.

Nie należy mylić:

  • \(\mathrm{kN/m}\) — kiloniuton na metr, jednostka obciążenia liniowego;
  • \(\mathrm{kN\,m}\) — kiloniutonometr, jednostka momentu.

W pierwszej jednostce długość znajduje się w mianowniku, a w drugiej jest mnożnikiem.

Kiedy można zastąpić obciążenie wypadkową?

Obciążenie rozłożone i odpowiadająca mu siła wypadkowa wywołują taki sam całkowity skutek mechaniczny względem analizowanego ciała sztywnego, jeżeli zachowane są:

Zastąpienie obciążenia wypadkową jest szczególnie wygodne podczas wyznaczania reakcji podporowych i sprawdzania równowagi całej konstrukcji.

Nie oznacza to jednak, że siła skupiona i obciążenie rozłożone są równoważne w każdej analizie. Powodują one inny przebieg sił wewnętrznych, ugięć i naprężeń wzdłuż elementu.

Przykład: przy wyznaczaniu reakcji belki można zastąpić obciążenie równomierne siłą \(R_q=ql\). Podczas wyznaczania siły tnącej i momentu zginającego w dowolnym przekroju trzeba jednak uwzględnić rzeczywiste obciążenie działające na odpowiednim fragmencie belki.

Dodawanie sił wypadkowych i momentów można wykonywać metodami rachunku wektorowego opisanymi w artykule Rachunek wektorowy w statyce – siły, składowe i momenty.

Najczęstsze błędy

Podsumowanie

W statyce obciążenia przedstawia się za pomocą uproszczonych modeli. Do najczęściej stosowanych należą siły skupione, momenty skupione oraz obciążenia rozłożone o stałej lub zmiennej intensywności.

Najważniejsze zależności dla obciążeń liniowych to:

\[ R_q=ql \qquad \text{dla obciążenia równomiernego}, \] \[ R_q=\frac{q_{\max}l}{2} \qquad \text{dla obciążenia trójkątnego}, \] \[ R_q=\frac{q_1+q_2}{2}\,l \qquad \text{dla obciążenia trapezowego}, \] \[ R_q=\int_a^b q(x)\,dx \qquad \text{dla dowolnej funkcji }q(x). \]

Wypadkowa obciążenia rozłożonego przechodzi przez środek ciężkości pola pod wykresem jego intensywności. Jej położenie musi zapewniać zachowanie tego samego momentu co pierwotne obciążenie.

Przy każdym obliczeniu należy kontrolować jednostki. Intensywność \(q\) ma jednostkę siły na długość, iloczyn \(ql\) jest siłą, natomiast \(ql^2\) ma jednostkę momentu.

Ruch prostoliniowy – położenie, prędkość, przyspieszenie i równania ruchu

Ruch prostoliniowy zachodzi wtedy, gdy ciało porusza się wzdłuż jednej prostej. Jego opis można sprowadzić do jednej współrzędnej \(x=x(t)\), a prędkość i przyspieszenie otrzymać przez kolejne różniczkowanie położenia względem czasu.

Choć ruch odbywa się tylko na jednej osi, poprawne rozwiązanie zadania wymaga uwzględnienia znaków położenia, prędkości i przyspieszenia. Szczególne znaczenie ma rozróżnienie między drogą a przemieszczeniem oraz między znakiem przyspieszenia a rzeczywistym wzrostem lub spadkiem szybkości.

Ogólne definicje prędkości, przyspieszenia i równania ruchu zostały omówione w artykule Kinematyka punktu materialnego – równania ruchu, prędkość, przyspieszenie i promień krzywizny. Tutaj zastosujemy je do ruchu wzdłuż jednej prostej.

Czym jest ruch prostoliniowy?

Ruch prostoliniowy jest ruchem, którego tor stanowi linia prosta. Po odpowiednim wyborze układu współrzędnych tor można utożsamić z jedną osią, najczęściej oznaczaną jako \(x\).

Położenie ciała opisuje wtedy jedna funkcja czasu:

\[ x=x(t). \]

Nie są potrzebne osobne równania dla współrzędnych \(y\) i \(z\), ponieważ podczas ruchu pozostają one stałe.

Ruchem prostoliniowym może być między innymi:

Opis ruchu na osi współrzędnych

Przed rozpoczęciem obliczeń należy określić:

Współrzędna \(x\) może być dodatnia, ujemna albo równa zeru. Jej znak informuje, po której stronie początku układu znajduje się ciało, ale nie określa kierunku jego ruchu.

Punkt poruszający się na osi x z zaznaczonym położeniem oraz wektorami prędkości i przyspieszenia
W ruchu prostoliniowym położenie określa współrzędna \(x=x(t)\), prędkość wynosi \(v=\frac{dx}{dt}\), a przyspieszenie \(a=\frac{dv}{dt}=\frac{d^2x}{dt^2}\).

Położenie i kierunek ruchu to dwie różne informacje. Ciało może znajdować się w obszarze \(x>0\), a jednocześnie poruszać się w lewo, czyli mieć prędkość \(v<0\).

Położenie, przemieszczenie i droga

Położenie ciała w danej chwili określa współrzędna \(x(t)\). Jeżeli w chwili \(t_1\) ciało znajduje się w położeniu \(x_1\), a w chwili \(t_2\) w położeniu \(x_2\), jego przemieszczenie wynosi:

\[ \Delta x = x_2-x_1. \] 

Przemieszczenie jest wielkością algebraiczną. Może być dodatnie, ujemne albo równe zeru.

Droga \(s\) jest całkowitą długością odcinka toru przebytego przez ciało. Jest zawsze nieujemna:

\[ s\geq 0. \] 

Jeżeli ciało nie zmienia zwrotu ruchu, droga jest równa wartości bezwzględnej przemieszczenia:

\[ s=|\Delta x|. \] 

Gdy ciało zawraca, droga jest większa od wartości przemieszczenia. Może się również zdarzyć, że ciało wróci do punktu początkowego. Wtedy:

\[ \Delta x=0, \qquad s>0. \]

Prędkość średnia i chwilowa

Prędkość średnia

Średnia prędkość w przedziale czasu od \(t_1\) do \(t_2\) jest ilorazem przemieszczenia i czasu:

\[ v_{\mathrm{śr}} = \frac{\Delta x}{\Delta t} = \frac{x_2-x_1}{t_2-t_1}. \]

Jej znak informuje o zwrocie wypadkowego przemieszczenia.

Średnia szybkość jest natomiast ilorazem drogi i czasu:

\[ v_{\mathrm{śr,szybkość}} = \frac{s}{\Delta t}. \] 

Średnia prędkość i średnia szybkość są sobie równe tylko wtedy, gdy ciało nie zmienia zwrotu i porusza się w dodatnim kierunku osi.

Prędkość chwilowa

Prędkość chwilowa jest granicą prędkości średniej dla coraz krótszego przedziału czasu:

\[ v(t) = \lim_{\Delta t\to 0} \frac{\Delta x}{\Delta t} = \frac{dx}{dt} = \dot x. \] 

Znak prędkości określa zwrot ruchu:

Wartość bezwzględna prędkości jest szybkością:

\[ \text{szybkość}=|v|. \]

Przyspieszenie średnie i chwilowe

Przyspieszenie średnie

Średnie przyspieszenie jest ilorazem zmiany prędkości i czasu:

\[ a_{\mathrm{śr}} = \frac{\Delta v}{\Delta t} = \frac{v_2-v_1}{t_2-t_1}. \]

Przyspieszenie chwilowe

Przyspieszenie chwilowe jest pochodną prędkości względem czasu:

\[ a(t) = \frac{dv}{dt} = \dot v. \] 

Ponieważ prędkość jest pochodną położenia, przyspieszenie jest jego drugą pochodną:

\[ a(t) = \frac{d^2x}{dt^2} = \ddot x. \] 

Dodatnie przyspieszenie oznacza, że prędkość algebraiczna rośnie. Ujemne przyspieszenie oznacza, że prędkość algebraiczna maleje. Nie jest to jeszcze jednoznaczne z przyspieszaniem albo zwalnianiem ciała.

Znaki prędkości i przyspieszenia

O tym, czy szybkość rośnie, czy maleje, decyduje wzajemny znak prędkości i przyspieszenia.

\[ va>0 \quad\Longrightarrow\quad |v|\ \text{rośnie}, \] \[ va<0 \quad\Longrightarrow\quad |v|\ \text{maleje}. \]

Jeżeli \(v\) i \(a\) mają ten sam znak, przyspieszenie jest skierowane zgodnie z ruchem i szybkość rośnie. Jeżeli mają przeciwne znaki, przyspieszenie jest skierowane przeciwnie do ruchu i szybkość maleje.

Cztery przypadki znaków prędkości i przyspieszenia pokazujące ruch przyspieszony oraz opóźniony
Jeżeli prędkość \(v\) i przyspieszenie \(a\) mają ten sam znak, wartość prędkości rośnie. Gdy mają przeciwne znaki, wartość prędkości maleje.
PrędkośćPrzyspieszenieZwrot ruchuZmiana szybkości
\(v>0\)\(a>0\)w praworośnie
\(v>0\)\(a<0\)w prawomaleje
\(v<0\)\(a<0\)w leworośnie
\(v<0\)\(a>0\)w lewomaleje

Ujemne przyspieszenie nie zawsze oznacza zwalnianie. Jeżeli ciało porusza się w kierunku ujemnym i ma \(v<0\) oraz \(a<0\), wartość jego prędkości rośnie.

Ruch prostoliniowy jednostajny

W ruchu prostoliniowym jednostajnym prędkość jest stała:

\[ v=\operatorname{const}. \]

Przyspieszenie jest więc równe zeru:

\[ a=\frac{dv}{dt}=0. \]

Jeżeli w chwili początkowej \(t_0\) ciało znajduje się w położeniu \(x_0\), to jego położenie w chwili \(t\) wynosi:

\[ x(t) = x_0+v(t-t_0). \]

Dla często stosowanego wyboru \(t_0=0\):

\[ x=x_0+vt. \]

W równych odstępach czasu ciało pokonuje jednakowe przemieszczenia:

\[ \Delta x=v\Delta t. \] 

Jeżeli \(v>0\), położenie rośnie liniowo. Dla \(v<0\) położenie liniowo maleje.

Ruch prostoliniowy jednostajnie zmienny

Ruch prostoliniowy jednostajnie zmienny zachodzi wtedy, gdy przyspieszenie jest stałe:

\[ a=\operatorname{const}. \]

Z definicji przyspieszenia:

\[ a=\frac{dv}{dt}. \]

Po scałkowaniu od chwili \(t_0\) do chwili \(t\):

\[ v(t) = v_0+a(t-t_0). \]

Dla \(t_0=0\):

\[ v=v_0+at. \]

Ponieważ:

\[ v=\frac{dx}{dt}, \]

kolejne całkowanie prowadzi do równania położenia:

\[ x(t) = x_0 + v_0(t-t_0) + \frac{a(t-t_0)^2}{2}. \]

Dla \(t_0=0\):

\[ x = x_0+v_0t+\frac{at^2}{2}. \]
Porównanie równań i wykresów ruchu prostoliniowego jednostajnego oraz jednostajnie zmiennego
W ruchu jednostajnym \(v=\operatorname{const}\) i \(x=x_0+vt\), natomiast przy stałym przyspieszeniu \(v=v_0+at\) oraz \(x=x_0+v_0t+\frac{at^2}{2}\).

Ruch jednostajnie przyspieszony

Ruch jest przyspieszony, gdy szybkość rośnie. Zachodzi wtedy:

\[ va>0. \]

Ruch jednostajnie opóźniony

Ruch jest opóźniony, gdy szybkość maleje:

\[ va<0. \]

Określenia „przyspieszony” i „opóźniony” odnoszą się do wartości prędkości, a nie do samego znaku \(a\).

Równanie ruchu bez czasu

W zadaniach ze stałym przyspieszeniem czas można wyeliminować. Korzystamy z reguły łańcuchowej:

\[ a = \frac{dv}{dt} = \frac{dv}{dx}\frac{dx}{dt} = v\frac{dv}{dx}. \]

Dla stałego \(a\):

\[ v\,dv=a\,dx. \]

Całkując od położenia \(x_0\) i prędkości \(v_0\) do \(x\) i \(v\), otrzymujemy:

\[ \frac{v^2-v_0^2}{2} = a(x-x_0). \]

Po przekształceniu:

\[ v^2 = v_0^2+2a(x-x_0). \]

Wzór ten jest szczególnie wygodny, gdy znamy drogę lub przemieszczenie, ale nie znamy czasu ruchu.

Wykresy położenia, prędkości i przyspieszenia

Wykres położenia \(x(t)\)

Nachylenie wykresu położenia jest równe prędkości:

\[ v=\frac{dx}{dt}. \]

W ruchu jednostajnym wykres \(x(t)\) jest linią prostą. Przy stałym niezerowym przyspieszeniu ma kształt paraboli.

Wykres prędkości \(v(t)\)

Nachylenie wykresu prędkości jest równe przyspieszeniu:

\[ a=\frac{dv}{dt}. \]

Dla ruchu jednostajnego wykres \(v(t)\) jest linią poziomą. Dla stałego przyspieszenia jest linią prostą, której współczynnik nachylenia wynosi \(a\).

Wykres przyspieszenia \(a(t)\)

W ruchu jednostajnym:

\[ a(t)=0. \]

W ruchu jednostajnie zmiennym wykres przyspieszenia jest poziomą linią położoną powyżej albo poniżej osi czasu, zależnie od znaku \(a\).

Pole pod wykresem prędkości

Przemieszczenie w przedziale czasu od \(t_1\) do \(t_2\) jest równe całce z prędkości:

\[ \Delta x = x(t_2)-x(t_1) = \int_{t_1}^{t_2}v(t)\,dt. \]

Geometrycznie jest to pole algebraiczne między wykresem \(v(t)\) a osią czasu:

Droga jest całką z wartości bezwzględnej prędkości:

\[ s = \int_{t_1}^{t_2}|v(t)|\,dt. \]

Jeżeli dodatnia część pola ma wartość \(A_+\), a wartość bezwzględna części ujemnej wynosi \(A_-\), to:

\[ \Delta x=A_+-A_-, \] \[ s=A_++A_-. \]
Wykres prędkości w funkcji czasu z zaznaczonymi polami odpowiadającymi przemieszczeniu i drodze
Pole algebraiczne pod wykresem \(v(t)\) jest równe przemieszczeniu \(\Delta x=\int_{t_1}^{t_2}v(t)\,dt\), natomiast droga wynosi \(s=\int_{t_1}^{t_2}|v(t)|\,dt\).

Pole pod wykresem \(v(t)\) nie zawsze jest drogą. Bez uwzględnienia wartości bezwzględnej otrzymujemy przemieszczenie. Różnica pojawia się wtedy, gdy prędkość zmienia znak.

Pole pod wykresem przyspieszenia

Zmiana prędkości w przedziale czasu jest równa całce z przyspieszenia:

\[ \Delta v = v(t_2)-v(t_1) = \int_{t_1}^{t_2}a(t)\,dt. \]

Pole algebraiczne pod wykresem \(a(t)\) informuje więc o zmianie prędkości, a nie bezpośrednio o przemieszczeniu.

Dla stałego przyspieszenia:

\[ \Delta v=a\Delta t. \] 

Ruch ze zmiennym przyspieszeniem

Jeżeli przyspieszenie zależy od czasu, nie można bezpośrednio stosować wzorów właściwych dla \(a=\operatorname{const}\).

Prędkość wyznaczamy przez całkowanie przyspieszenia:

\[ v(t) = v(t_0) + \int_{t_0}^{t}a(\tau)\,d\tau. \]

Położenie otrzymujemy przez całkowanie prędkości:

\[ x(t) = x(t_0) + \int_{t_0}^{t}v(\tau)\,d\tau. \]

Symbol \(\tau\) jest pomocniczą zmienną całkowania i nie oznacza dodatkowej wielkości fizycznej.

Jeżeli natomiast znana jest zależność przyspieszenia od położenia \(a=a(x)\), można skorzystać z zależności:

\[ a = v\frac{dv}{dx}. \]

Zatrzymanie i zmiana zwrotu ruchu

Ciało zatrzymuje się w chwili, w której:

\[ v(t)=0. \]

Samo spełnienie tego warunku nie zawsze oznacza zmianę zwrotu. Należy sprawdzić znak prędkości przed i po rozpatrywanej chwili.

Przy stałym przyspieszeniu chwila zatrzymania spełnia:

\[ 0=v_0+at_{\mathrm{z}}, \]

czyli:

\[ t_{\mathrm{z}} = -\frac{v_0}{a}. \]

Rozwiązanie ma znaczenie fizyczne tylko wtedy, gdy otrzymany czas należy do analizowanego przedziału i jest zgodny z przyjętym początkiem pomiaru czasu.

Przykłady obliczeniowe

Przykład 1. Ruch jednostajny

Samochód znajduje się początkowo w położeniu \(x_0=50\ \mathrm{m}\) i porusza się ze stałą prędkością:

\[ v=20\ \mathrm{m/s}. \]

Wyznaczmy położenie po \(8\ \mathrm{s}\).

\[ x=x_0+vt, \] \[ x = 50+20\cdot 8 = 210\ \mathrm{m}. \]

Przemieszczenie wynosi:

\[ \Delta x = 210-50 = 160\ \mathrm{m}. \]

Przykład 2. Ruch ze stałym przyspieszeniem

Ciało rozpoczyna ruch w położeniu \(x_0=10\ \mathrm{m}\) z prędkością \(v_0=5\ \mathrm{m/s}\). Przyspieszenie wynosi:

\[ a=2\ \mathrm{m/s^2}. \]

Po \(4\ \mathrm{s}\) prędkość jest równa:

\[ v=v_0+at, \] \[ v = 5+2\cdot 4 = 13\ \mathrm{m/s}. \]

Położenie wynosi:

\[ x = x_0+v_0t+\frac{at^2}{2}, \] \[ x = 10+5\cdot 4+\frac{2\cdot 4^2}{2} = 46\ \mathrm{m}. \]

Przemieszczenie wyniosło:

\[ \Delta x=46-10=36\ \mathrm{m}. \]

Przykład 3. Hamowanie samochodu

Samochód porusza się z prędkością początkową:

\[ v_0=25\ \mathrm{m/s} \]

i hamuje ze stałym przyspieszeniem:

\[ a=-5\ \mathrm{m/s^2}. \]

Czas zatrzymania:

\[ 0=v_0+at, \] \[ t_{\mathrm{z}} = -\frac{v_0}{a} = -\frac{25}{-5} = 5\ \mathrm{s}. \]

Drogę hamowania wyznaczamy z równania bez czasu:

\[ v^2 = v_0^2+2a\Delta x. \] 

Dla \(v=0\):

\[ 0 = 25^2+2(-5)\Delta x, \] \[ \Delta x = 62{,}5\ \mathrm{m}. \]

Samochód zatrzyma się po \(5\ \mathrm{s}\), pokonując w tym czasie \(62{,}5\ \mathrm{m}\).

Przykład 4. Zmiana zwrotu ruchu

Prędkość ciała jest opisana równaniem:

\[ v(t)=-4+2t. \]

Chwila zatrzymania:

\[ -4+2t=0, \] \[ t=2\ \mathrm{s}. \]

Dla \(t<2\ \mathrm{s}\) prędkość jest ujemna, a dla \(t>2\ \mathrm{s}\) dodatnia. Ciało zatrzymuje się więc w chwili \(t=2\ \mathrm{s}\), a następnie zmienia zwrot ruchu.

Jednostki i kontrola wymiarów

WielkośćOznaczenieJednostka SI
Położenie\(x\)\(\mathrm{m}\)
Przemieszczenie\(\Delta x\)\(\mathrm{m}\)
Droga\(s\)\(\mathrm{m}\)
Czas\(t\)\(\mathrm{s}\)
Prędkość\(v\)\(\mathrm{m/s}\)
Przyspieszenie\(a\)\(\mathrm{m/s^2}\)

Kontrola wymiarów równania położenia:

\[ [x] = [x_0]+[v_0t]+\left[\frac{at^2}{2}\right], \] \[ \mathrm{m} = \mathrm{m} + \frac{\mathrm{m}}{\mathrm{s}}\mathrm{s} + \frac{\mathrm{m}}{\mathrm{s^2}}\mathrm{s^2}. \]

Każdy składnik ma jednostkę długości, więc równanie jest wymiarowo poprawne.

Najczęstsze błędy

Aby zamienić kilometry na godzinę na metry na sekundę, dzielimy wartość przez \(3{,}6\):

\[ v[\mathrm{m/s}] = \frac{v[\mathrm{km/h}]}{3{,}6}. \]

W przeciwną stronę mnożymy przez \(3{,}6\).

Podsumowanie

Ruch prostoliniowy można opisać za pomocą jednej współrzędnej:

\[ x=x(t). \]

Prędkość i przyspieszenie są kolejnymi pochodnymi położenia względem czasu:

\[ v = \frac{dx}{dt}, \qquad a = \frac{dv}{dt} = \frac{d^2x}{dt^2}. \]

W ruchu jednostajnym:

\[ v=\operatorname{const}, \qquad a=0, \qquad x=x_0+vt. \]

Dla stałego przyspieszenia:

\[ v=v_0+at, \] \[ x=x_0+v_0t+\frac{at^2}{2}, \] \[ v^2=v_0^2+2a(x-x_0). \]

Przemieszczenie jest polem algebraicznym pod wykresem prędkości:

\[ \Delta x = \int_{t_1}^{t_2}v(t)\,dt, \]

natomiast droga jest całką z wartości bezwzględnej prędkości:

\[ s = \int_{t_1}^{t_2}|v(t)|\,dt. \]

O tym, czy ciało przyspiesza, czy zwalnia, decydują łącznie znaki \(v\) i \(a\). Zgodne znaki oznaczają wzrost szybkości, a przeciwne — jej zmniejszanie.

Kinematyka punktu materialnego – równania ruchu, prędkość, przyspieszenie i promień krzywizny

Kinematyka punktu materialnego opisuje zmianę położenia punktu w czasie. Znając równanie ruchu, możemy wyznaczyć tor, prędkość i przyspieszenie, a w ruchu krzywoliniowym także składowe styczną i normalną przyspieszenia oraz promień krzywizny toru.

Punkt materialny jest uproszczonym modelem ciała, którego rozmiary i obrót nie mają znaczenia w rozpatrywanym zagadnieniu. W kinematyce nie analizujemy jeszcze sił powodujących ruch — interesuje nas wyłącznie jego geometryczny i czasowy opis.

Podstawowe pojęcia, takie jak układ odniesienia, tor, droga i przemieszczenie, zostały omówione w artykule Czym zajmuje się kinematyka? Podstawowe pojęcia opisu ruchu. W tym materiale przechodzimy do szczegółowego i rachunkowego opisu ruchu punktu.

Model punktu materialnego

Punkt materialny jest modelem ciała, którego wymiary geometryczne pomijamy. Całe ciało reprezentuje jeden punkt \(P\), a jego ruch opisujemy przez zmianę położenia tego punktu w czasie.

Model ten można stosować wtedy, gdy:

Na przykład satelitę poruszającego się po orbicie można w pierwszym przybliżeniu traktować jako punkt materialny, jeżeli interesuje nas jego trajektoria, a nie obrót paneli słonecznych czy orientacja przestrzenna całego obiektu.

Punkt poruszający się po zakrzywionym torze z wektorem położenia, prędkością, przyspieszeniem i promieniem krzywizny
W punkcie \(P\) prędkość \(\vec{v}\) jest styczna do toru, a przyspieszenie \(\vec{a}\) można rozłożyć na składową styczną \(\vec{a}_\tau\) i normalną \(\vec{a}_n\).

Wektorowe równanie ruchu

Położenie punktu \(P\) względem początku układu współrzędnych \(O\) opisuje wektor położenia, nazywany także wektorem wodzącym:

\[ \vec r=\overrightarrow{OP}. \]

Podczas ruchu wektor położenia zmienia się w czasie:

\[ \vec r=\vec r(t). \]

Zależność ta jest wektorowym równaniem ruchu punktu materialnego. Jeżeli znamy funkcję \(\vec r(t)\), możemy dla każdej chwili \(t\) wyznaczyć położenie punktu, a przez różniczkowanie także jego prędkość i przyspieszenie.

W prostokątnym układzie współrzędnych:

\[ \vec r(t) = x(t)\vec i + y(t)\vec j + z(t)\vec k, \] 

gdzie \(\vec i\), \(\vec j\) i \(\vec k\) są wersorami osi \(x\), \(y\) i \(z\).

Równania ruchu we współrzędnych

Wektorowe równanie ruchu można zastąpić układem równań skalarnych:

\[ x=x(t),\qquad y=y(t),\qquad z=z(t). \]

Są to parametryczne równania ruchu, w których czas pełni rolę parametru. W ruchu płaskim jedna ze współrzędnych jest stała, dlatego wystarczają dwa równania:

\[ x=x(t),\qquad y=y(t). \]

W ruchu prostoliniowym, jeżeli oś zostanie poprowadzona wzdłuż toru, położenie można opisać jednym równaniem:

\[ x=x(t). \]

Opis wektorowy i współrzędnościowy przedstawiają ten sam ruch. Różnią się jedynie sposobem zapisu, co szerzej omówiliśmy w sekcji Sposoby opisu ruchu.

Wyznaczanie toru ruchu

Parametryczne równania ruchu określają zarówno kształt toru, jak i położenie punktu w każdej chwili. Aby otrzymać samo równanie toru, należy wyeliminować czas.

Rozważmy przykład:

\[ x=2t,\qquad y=t^2. \]

Z pierwszego równania:

\[ t=\frac{x}{2}. \]

Po podstawieniu do równania współrzędnej \(y\):

\[ y= \left(\frac{x}{2}\right)^2 = \frac{x^2}{4}. \]

Tor punktu jest więc parabolą.

Wektorowe i współrzędnościowe równania ruchu oraz parabola otrzymana po wyeliminowaniu czasu
Dla równań parametrycznych \(x=2t\) i \(y=t^2\) wyeliminowanie czasu prowadzi do równania toru \(y=\frac{x^2}{4}\).

Równanie toru nie jest pełnym równaniem ruchu. Zależność \(y=y(x)\) opisuje kształt krzywej, ale nie informuje, w której chwili punkt znajduje się w określonym miejscu ani z jaką prędkością się porusza.

Oznaczenia pochodnych w kinematyce

W kinematyce spotyka się kilka równoważnych sposobów zapisywania pochodnych. Szczególnie często używane są zapis ilorazowy, notacja Newtona z kropką oraz notacja z primem.

Pochodne po czasie

Pierwszą pochodną wielkości \(x(t)\) względem czasu można zapisać jako:

\[ \dot x=\frac{dx}{dt}. \]

Druga pochodna jest oznaczana dwiema kropkami:

\[ \ddot x=\frac{d^2x}{dt^2}. \]

Dla wektora położenia:

\[ \dot{\vec r} = \frac{d\vec r}{dt} = \vec v, \] \[ \ddot{\vec r} = \frac{d^2\vec r}{dt^2} = \vec a. \] 

Notacja z primem

Prim nie określa samodzielnie, względem jakiej wielkości obliczana jest pochodna. Znaczenie zapisu wynika z wcześniej wskazanej zmiennej niezależnej.

Dla funkcji \(y=y(x)\):

\[ y'=\frac{dy}{dx}, \qquad y''=\frac{d^2y}{dx^2}. \]

Jeżeli natomiast funkcja zależy od współrzędnej naturalnej \(s\), można przyjąć:

\[ f'=\frac{df}{ds}, \qquad f''=\frac{d^2f}{ds^2}. \]
Porównanie notacji z kropką, zapisu d przez dt oraz notacji z primem stosowanych w kinematyce
Kropka oznacza zwykle pochodną po czasie, np. \(\dot{x}=\frac{dx}{dt}\), natomiast prim oznacza pochodną po przyjętej zmiennej niezależnej, np. \(y'=\frac{dy}{dx}\) lub \(f'=\frac{df}{ds}\).

Reguła łańcuchowa

Jeżeli:

\[ f=f(s), \qquad s=s(t), \]

to pochodna funkcji \(f\) względem czasu wynosi:

\[ \dot f = \frac{df}{dt} = \frac{df}{ds}\frac{ds}{dt} = f'\dot s. \] 

Dla drugiej pochodnej:

\[ \ddot f = f''\dot s^{\,2} + f'\ddot s. \] 

Najbardziej jednoznaczny jest zapis ilorazowy. Symbole \(\frac{df}{dt}\), \(\frac{df}{dx}\) i \(\frac{df}{ds}\) od razu wskazują zmienną, względem której wykonujemy różniczkowanie.

Prędkość punktu materialnego

Prędkość chwilowa jest pochodną wektora położenia względem czasu:

\[ \vec v = \frac{d\vec r}{dt} = \dot{\vec r}. \] 

W kartezjańskim układzie współrzędnych:

\[ \vec v = \dot x\vec i + \dot y\vec j + \dot z\vec k. \] 

Składowe prędkości wynoszą:

\[ v_x=\dot x=\frac{dx}{dt}, \qquad v_y=\dot y=\frac{dy}{dt}, \qquad v_z=\dot z=\frac{dz}{dt}. \]

Wartość wektora prędkości, czyli szybkość, jest równa:

\[ v=|\vec v| = \sqrt{\dot x^2+\dot y^2+\dot z^2}. \]

Wektor prędkości chwilowej jest styczny do toru i skierowany zgodnie z chwilowym zwrotem ruchu.

Przyspieszenie punktu materialnego

Przyspieszenie chwilowe jest pochodną prędkości względem czasu albo drugą pochodną położenia:

\[ \vec a = \frac{d\vec v}{dt} = \frac{d^2\vec r}{dt^2} = \ddot{\vec r}. \] 

W kartezjańskim układzie współrzędnych:

\[ \vec a = \ddot x\vec i + \ddot y\vec j + \ddot z\vec k. \] 

Składowe przyspieszenia wynoszą:

\[ a_x=\ddot x, \qquad a_y=\ddot y, \qquad a_z=\ddot z. \] 

Wartość całkowitego przyspieszenia:

\[ a=|\vec a| = \sqrt{\ddot x^2+\ddot y^2+\ddot z^2}. \]

Przyspieszenie opisuje każdą zmianę wektora prędkości. Może być wynikiem zmiany jego wartości, kierunku albo obu tych cech jednocześnie.

Opis naturalny ruchu

W opisie naturalnym wykorzystuje się lokalny układ związany z torem ruchu. Podstawowymi kierunkami są:

Wektor prędkości można zapisać jako:

\[ \vec v=v\vec\tau. \]

Jeżeli \(s\) oznacza długość drogi mierzoną wzdłuż toru, to:

\[ v=\frac{ds}{dt}=\dot s. \] 

W ruchu płaskim zmiana wersora stycznego wraz z położeniem na torze spełnia zależność:

\[ \frac{d\vec\tau}{ds} = \frac{1}{\rho}\vec n, \] 

gdzie \(\rho\) jest promieniem krzywizny toru.

Przyspieszenie styczne i normalne

Różniczkując wektor prędkości:

\[ \vec a = \frac{d}{dt}\left(v\vec\tau\right), \]

otrzymujemy:

\[ \vec a = \frac{dv}{dt}\vec\tau + v\frac{d\vec\tau}{dt}. \]

Z reguły łańcuchowej:

\[ \frac{d\vec\tau}{dt} = \frac{d\vec\tau}{ds}\frac{ds}{dt} = \frac{v}{\rho}\vec n. \] 

Ostatecznie:

\[ \vec a = a_\tau\vec\tau + a_n\vec n, \] 

gdzie:

\[ a_\tau = \frac{dv}{dt}, \qquad a_n = \frac{v^2}{\rho}. \]

Składowa styczna

Składowa styczna odpowiada za zmianę wartości prędkości:

Składowa normalna

Składowa normalna odpowiada za zmianę kierunku wektora prędkości. Jest prostopadła do prędkości i skierowana ku środkowi krzywizny toru.

Dla ruchu prostoliniowego:

\[ \rho\to\infty, \qquad a_n=0. \]

W ruchu jednostajnym po okręgu szybkość jest stała, dlatego \(a_\tau=0\), ale występuje niezerowe przyspieszenie normalne.

Prędkość styczna do toru oraz przyspieszenie rozłożone na składową styczną i normalną skierowaną do środka krzywizny
Przyspieszenie punktu ma składowe \(a_\tau=\frac{dv}{dt}\) i \(a_n=\frac{v^2}{\rho}\), dlatego jego wartość wynosi \(a=\sqrt{a_\tau^2+a_n^2}\).

Ponieważ kierunki styczny i normalny są prostopadłe, wartość przyspieszenia całkowitego obliczamy z twierdzenia Pitagorasa:

\[ a = \sqrt{a_\tau^2+a_n^2}. \]

Składową styczną można również wyrazić za pomocą pochodnej szybkości względem drogi:

\[ a_\tau = \frac{dv}{dt} = \frac{dv}{ds}\frac{ds}{dt} = v\frac{dv}{ds}. \]

Przy umowie \(v'=dv/ds\) otrzymujemy:

\[ \dot v=vv'. \]

Promień krzywizny toru

Promień krzywizny \(\rho\) określa stopień zakrzywienia toru w danym punkcie. Jest promieniem okręgu oskulacyjnego, czyli okręgu najlepiej dopasowanego do toru w bezpośrednim otoczeniu rozpatrywanego punktu.

Z zależności na przyspieszenie normalne wynika:

\[ \rho = \frac{v^2}{a_n}, \]

o ile \(a_n>0\).

Promień krzywizny jest wielkością lokalną. Na jednym torze może przyjmować różne wartości w różnych punktach. Nie należy utożsamiać go z odległością punktu od początku układu współrzędnych.

Wyznaczanie promienia krzywizny

Na podstawie prędkości i przyspieszenia

Dla dowolnego ruchu przestrzennego krzywiznę toru można zapisać jako:

\[ \kappa = \frac{|\vec v\times\vec a|}{v^3}. \]

Promień krzywizny jest odwrotnością krzywizny:

\[ \rho = \frac{1}{\kappa} = \frac{v^3}{|\vec v\times\vec a|}. \] 

Iloczyn wektorowy oraz jego własności zostały szerzej omówione w artykule Rachunek wektorowy w statyce – siły, składowe i momenty.

Dla równań parametrycznych \(x(t)\) i \(y(t)\)

W ruchu płaskim:

\[ v = \sqrt{\dot x^2+\dot y^2}. \]

Promień krzywizny wynosi:

\[ \rho = \frac{\left(\dot x^2+\dot y^2\right)^{3/2}} {\left|\dot x\ddot y-\dot y\ddot x\right|}. \]

W mianowniku występuje wartość bezwzględna, ponieważ promień krzywizny jest dodatnią długością.

Dla toru \(y=y(x)\)

Jeżeli tor jest podany bezpośrednio jako funkcja współrzędnej \(x\), stosujemy:

\[ \rho = \frac{\left[1+\left(y'\right)^2\right]^{3/2}} {|y''|}, \]

gdzie:

\[ y'=\frac{dy}{dx}, \qquad y''=\frac{d^2y}{dx^2}. \]

Jest to przykład sytuacji, w której primy oznaczają pochodne po współrzędnej \(x\), a nie po czasie.

Dwa rozkłady wektora przyspieszenia

Ten sam wektor przyspieszenia można rozłożyć w różnych bazach. Najczęściej stosuje się bazę kartezjańską oraz bazę styczno-normalną.

Ten sam wektor przyspieszenia rozłożony na składowe x i y oraz na składowe styczną i normalną
Ten sam wektor przyspieszenia można zapisać jako \(\vec{a}=a_x\vec{i}+a_y\vec{j}\) albo \(\vec{a}=a_\tau\vec{\tau}+a_n\vec{n}\), zależnie od wybranej bazy.

Rozkład kartezjański

\[ \vec a = a_x\vec i+a_y\vec j, \] \[ a = \sqrt{a_x^2+a_y^2}. \]

Kierunki \(\vec i\) i \(\vec j\) są związane ze stałymi, nieruchomymi osiami układu współrzędnych.

Rozkład styczno-normalny

\[ \vec a = a_\tau\vec\tau+a_n\vec n, \] \[ a = \sqrt{a_\tau^2+a_n^2}. \]

Kierunki \(\vec\tau\) i \(\vec n\) są związane z lokalną geometrią toru i zmieniają orientację wraz z ruchem punktu.

Nie należy utożsamiać składowych:

\[ a_x\neq a_\tau, \qquad a_y\neq a_n \] 

w ogólnym przypadku. Obie pary są składowymi tego samego wektora, ale wyznaczonymi w dwóch różnych bazach ortogonalnych. Mogą być sobie równe jedynie w szczególnym położeniu, gdy odpowiednie kierunki obu baz chwilowo się pokrywają.

Przykład obliczeniowy

Ruch punktu w płaszczyźnie opisują równania:

\[ x=2t, \qquad y=t^2. \]

Wektor położenia

\[ \vec r(t) = 2t\,\vec i+t^2\vec j. \] 

Prędkość

\[ \vec v(t) = \dot{\vec r} = 2\vec i+2t\vec j. \] 

Szybkość wynosi:

\[ v = \sqrt{2^2+(2t)^2} = 2\sqrt{1+t^2}. \]

Przyspieszenie

\[ \vec a(t) = \dot{\vec v} = 2\vec j. \] 

Jego wartość jest stała:

\[ a=2. \]

Składowa styczna

\[ a_\tau = \frac{dv}{dt} = \frac{2t}{\sqrt{1+t^2}}. \]

Składowa normalna

Ponieważ:

\[ a^2=a_\tau^2+a_n^2, \]

otrzymujemy:

\[ a_n = \sqrt{a^2-a_\tau^2} = \frac{2}{\sqrt{1+t^2}}. \]

Promień krzywizny

\[ \rho = \frac{v^2}{a_n} = 2\left(1+t^2\right)^{3/2}. \]

Dla chwili \(t=1\):

\[ \vec r(1)=2\vec i+\vec j, \] \[ \vec v(1)=2\vec i+2\vec j, \qquad v(1)=2\sqrt 2, \] \[ \vec a(1)=2\vec j, \] \[ a_\tau(1)=\sqrt 2, \qquad a_n(1)=\sqrt 2, \] \[ \rho(1)=4\sqrt 2. \] 

Wektor przyspieszenia ma w tym przykładzie stały kierunek kartezjański, ale jego składowe styczna i normalna zmieniają się wraz z położeniem punktu na torze. Dobrze pokazuje to różnicę pomiędzy oboma sposobami rozkładu.

Jednostki i kontrola wymiarów

WielkośćOznaczenieJednostka SI
Położenie\(\vec r\), \(x\), \(y\), \(z\)\(\mathrm{m}\)
Czas\(t\)\(\mathrm{s}\)
Prędkość\(\vec v\)\(\mathrm{m/s}\)
Przyspieszenie\(\vec a\)\(\mathrm{m/s^2}\)
Krzywizna\(\kappa\)\(\mathrm{1/m}\)
Promień krzywizny\(\rho\)\(\mathrm{m}\)

Dla prędkości:

\[ [\vec v] = \frac{[\vec r]}{[t]} = \frac{\mathrm{m}}{\mathrm{s}}. \]

Dla przyspieszenia normalnego:

\[ \left[\frac{v^2}{\rho}\right] = \frac{\mathrm{m^2/s^2}}{\mathrm{m}} = \frac{\mathrm{m}}{\mathrm{s^2}}. \]

Dla promienia krzywizny:

\[ \left[\frac{v^2}{a_n}\right] = \frac{\mathrm{m^2/s^2}}{\mathrm{m/s^2}} = \mathrm{m}. \]

Najczęstsze błędy

Podsumowanie

Ruch punktu materialnego opisuje zależny od czasu wektor położenia:

\[ \vec r=\vec r(t). \]

Prędkość i przyspieszenie otrzymujemy przez kolejne różniczkowanie po czasie:

\[ \vec v = \dot{\vec r} = \frac{d\vec r}{dt}, \] \[ \vec a = \ddot{\vec r} = \frac{d^2\vec r}{dt^2}. \]

W opisie naturalnym przyspieszenie rozkłada się na wzajemnie prostopadłe składowe:

\[ \vec a = a_\tau\vec\tau+a_n\vec n, \] \[ a_\tau=\frac{dv}{dt}, \qquad a_n=\frac{v^2}{\rho}. \]

Składowa styczna odpowiada za zmianę szybkości, a normalna za zmianę kierunku ruchu. Promień krzywizny można wyznaczyć z zależności:

\[ \rho=\frac{v^2}{a_n}, \]

albo bezpośrednio z parametrów ruchu:

\[ \rho = \frac{\left(\dot x^2+\dot y^2\right)^{3/2}} {\left|\dot x\ddot y-\dot y\ddot x\right|}. \]

Rozkład na składowe kartezjańskie i rozkład styczno-normalny opisują ten sam wektor przyspieszenia w dwóch różnych bazach ortogonalnych. Wybór sposobu opisu zależy od danych zadania i geometrii toru.

Czym zajmuje się kinematyka? Podstawowe pojęcia opisu ruchu

Kinematyka jest działem mechaniki zajmującym się opisem ruchu ciał bez analizowania sił i innych przyczyn, które ten ruch wywołują. Odpowiada przede wszystkim na pytania: gdzie znajduje się ciało, po jakim torze się porusza, jak szybko zmienia swoje położenie i jakie ma przyspieszenie.

W kinematyce badamy więc położenie, drogę, przemieszczenie, prędkość i przyspieszenie jako funkcje czasu. Dopiero dynamika wyjaśnia, dlaczego ruch przebiega w określony sposób i jakie oddziaływania są za niego odpowiedzialne.

Dla porównania, statyka zajmuje się przede wszystkim warunkami równowagi ciał. Kinematyka opisuje natomiast zarówno ciała pozostające w spoczynku, jak i poruszające się w dowolny sposób.

Czym zajmuje się kinematyka?

Kinematyka opisuje ruch ciał w przestrzeni i czasie. Nie bada przy tym masy ciała ani sił powodujących zmianę jego ruchu.

Podstawowymi zagadnieniami kinematyki są:

Znając równanie ruchu, można określić położenie ciała w dowolnej chwili, a następnie wyznaczyć jego prędkość i przyspieszenie.

Kinematyka a dynamika

Kinematyka odpowiada na pytanie:

Jak porusza się ciało?

Dynamika odpowiada natomiast na pytanie:

Dlaczego ciało porusza się w taki sposób?

Przykładowo, obserwując rowerzystę poruszającego się po zakręcie, w kinematyce możemy opisać jego tor, prędkość oraz przyspieszenie. Nie analizujemy jednak jeszcze siły tarcia, siły nacisku ani innych oddziaływań umożliwiających wykonanie zakrętu.

Kinematyka nie oznacza braku przyczyn ruchu. Przyczyny te istnieją, lecz na danym etapie analizy nie są rozpatrywane.

Układ odniesienia

Ruch zawsze opisuje się względem określonego układu odniesienia. Aby określić położenie ciała, potrzebujemy:

Najczęściej stosuje się prostokątny kartezjański układ współrzędnych \(Oxyz\). W zagadnieniach płaskich wystarczają dwie osie \(Ox\) i \(Oy\), natomiast w ruchu prostoliniowym często używa się tylko jednej osi.

To samo ciało może pozostawać w spoczynku względem jednego układu odniesienia i jednocześnie poruszać się względem innego. Pasażer siedzący w jadącym pociągu nie porusza się względem wagonu, lecz zmienia swoje położenie względem powierzchni Ziemi.

Punkt materialny i bryła sztywna

W najprostszym modelu ciało zastępuje się punktem materialnym. Jest to punkt posiadający masę, lecz pozbawiony wymiarów geometrycznych.

Takie uproszczenie można zastosować, gdy rozmiary i obrót ciała nie mają znaczenia dla badanego zagadnienia. Samochód jadący z miasta do miasta może być traktowany jako punkt materialny, jeżeli interesuje nas jedynie jego położenie na trasie.

Jeżeli trzeba uwzględnić rozmiary, kształt i obrót ciała, stosuje się model bryły sztywnej. W bryle sztywnej odległości pomiędzy dowolnymi punktami pozostają stałe.

Kinematykę dzieli się więc między innymi na:

Położenie punktu materialnego

Położenie punktu \(P\) względem początku układu współrzędnych \(O\) opisuje wektor położenia, nazywany także wektorem wodzącym:

\[ \vec r=\overrightarrow{OP}. \]

Jeżeli punkt porusza się, jego położenie zmienia się w czasie:

\[ \vec r=\vec r(t). \]

W kartezjańskim układzie współrzędnych:

\[ \vec r(t) = x(t)\vec i + y(t)\vec j + z(t)\vec k. \] 

Funkcje \(x(t)\), \(y(t)\) i \(z(t)\) nazywamy równaniami ruchu punktu w postaci parametrycznej. Znając czas \(t\), możemy obliczyć aktualne współrzędne punktu.

Sposoby opisu ruchu

Ruch punktu materialnego można opisać na kilka równoważnych sposobów. Wybór metody zależy od kształtu toru i rodzaju rozwiązywanego zagadnienia.

Porównanie wektorowego, współrzędnościowego i naturalnego sposobu opisu ruchu punktu materialnego
Ruch punktu materialnego można opisać za pomocą wektora położenia \(\vec r=\vec r(t)\), współrzędnych \(x=x(t)\), \(y=y(t)\), \(z=z(t)\) albo współrzędnej naturalnej \(s=s(t)\) mierzonej wzdłuż toru.

Opis wektorowy

W opisie wektorowym położenie punktu określa wektor:

\[ \vec r=\vec r(t). \]

Jest to zapis zwięzły i niezależny od konkretnego sposobu rozłożenia ruchu na osie współrzędnych.

Opis za pomocą współrzędnych

W opisie współrzędnościowym podaje się każdą współrzędną punktu jako funkcję czasu:

\[ x=x(t),\qquad y=y(t),\qquad z=z(t). \]

W zagadnieniu płaskim wystarczają zwykle równania \(x=x(t)\) oraz \(y=y(t)\). Eliminując parametr \(t\), można niekiedy otrzymać równanie toru w postaci:

\[ y=y(x). \]

Opis naturalny

Jeżeli tor ruchu jest znany, położenie punktu można określić za pomocą współrzędnej naturalnej \(s\), mierzonej wzdłuż toru od przyjętego punktu początkowego:

\[ s=s(t). \]

Należy przy tym ustalić początek pomiaru oraz dodatni zwrot współrzędnej \(s\). Opis naturalny jest szczególnie wygodny w ruchu po znanym, zakrzywionym torze.

Ważne: współrzędna naturalna \(s(t)\) nie zawsze jest tym samym co całkowita droga przebyta od rozpoczęcia obserwacji. Jeżeli punkt zmieni zwrot ruchu, współrzędna \(s\) może maleć, natomiast całkowita przebyta droga nadal rośnie.

Tor ruchu

Torem ruchu nazywamy linię utworzoną przez kolejne położenia poruszającego się punktu.

W zależności od kształtu toru ruch może być:

Tor jest pojęciem geometrycznym. Nie informuje samodzielnie, jak szybko punkt porusza się po danej linii ani w której chwili znajduje się w określonym miejscu.

Dwa punkty mogą poruszać się po takim samym torze, lecz mieć zupełnie inne prędkości i inne zależności położenia od czasu.

Droga i przemieszczenie

Drogaprzemieszczenie opisują zmianę położenia ciała, ale są różnymi wielkościami.

Punkt poruszający się po zakrzywionym torze od A do B z zaznaczoną drogą i wektorem przemieszczenia
Droga \(s\) jest długością części toru przebytej przez ciało, natomiast przemieszczenie \(\Delta\vec r=\vec r_B-\vec r_A\) jest wektorem łączącym położenie początkowe z końcowym. Zachodzi \(s\geq|\Delta\vec r|\).

Droga

Droga \(s\) jest długością fragmentu toru przebytego przez ciało. Jest wielkością skalarną i nie przyjmuje wartości ujemnych.

Jeżeli ciało porusza się, całkowita droga przebyta od chwili początkowej nie maleje, nawet gdy ciało zawraca.

Przemieszczenie

Przemieszczenie jest zmianą wektora położenia:

\[ \Delta\vec r = \vec r_2-\vec r_1. \] 

Jest to wektor skierowany od położenia początkowego do położenia końcowego. Zależy jedynie od tych dwóch położeń, a nie od kształtu przebytego toru.

Długość przemieszczenia nie może być większa od drogi:

\[ |\Delta\vec r|\leq s. \] 

Równość zachodzi w ruchu po odcinku prostym bez zmiany zwrotu.

Jeżeli ciało wróci do miejsca początkowego, jego przemieszczenie wynosi zero:

\[ \Delta\vec r=\vec 0, \] 

ale przebyta droga może być znacznie większa od zera.

Prędkość i szybkość

Prędkość średnia

Wektor prędkości średniej jest równy ilorazowi przemieszczenia i czasu, w którym ono nastąpiło:

\[ \vec v_{\mathrm{śr}} = \frac{\Delta\vec r}{\Delta t}. \] 

Kierunek prędkości średniej jest zgodny z kierunkiem wektora przemieszczenia, a niekoniecznie z przebiegiem toru.

Szybkość średnia

Szybkość średnia jest ilorazem przebytej drogi i czasu:

\[ v_{\mathrm{śr,droga}} = \frac{\Delta s}{\Delta t}. \] 

Jest to wielkość skalarna. Prędkość średnia i szybkość średnia nie są więc tym samym.

Prędkość chwilowa

Prędkość chwilową otrzymujemy jako pochodną wektora położenia względem czasu:

\[ \vec v = \frac{d\vec r}{dt}. \]

W kartezjańskim układzie współrzędnych:

\[ \vec v = \frac{dx}{dt}\vec i + \frac{dy}{dt}\vec j + \frac{dz}{dt}\vec k. \] 

Wektor prędkości chwilowej jest styczny do toru ruchu i skierowany zgodnie z chwilowym zwrotem ruchu.

Szybkość jest wartością wektora prędkości:

\[ v=|\vec v|. \] 

W opisie naturalnym szybkość można zapisać jako:

\[ v=\frac{ds}{dt}, \]

jeżeli \(s\) oznacza rosnącą długość drogi mierzonej wzdłuż toru.

Prędkość jest wektorem, a szybkość skalarem. Sam zapis „samochód jedzie z prędkością 50 km/h” podaje w rzeczywistości jedynie wartość prędkości, czyli szybkość. Pełny wektor wymaga także kierunku i zwrotu.

Przyspieszenie

Przyspieszenie opisuje zmianę wektora prędkości w czasie. Przyspieszenie średnie wynosi:

\[ \vec a_{\mathrm{śr}} = \frac{\Delta\vec v}{\Delta t}. \] 

Przyspieszenie chwilowe jest pochodną prędkości względem czasu:

\[ \vec a = \frac{d\vec v}{dt} = \frac{d^2\vec r}{dt^2}. \]

W kartezjańskim układzie współrzędnych:

\[ \vec a = \frac{d^2x}{dt^2}\vec i + \frac{d^2y}{dt^2}\vec j + \frac{d^2z}{dt^2}\vec k. \] 

Przyspieszenie może wynikać ze:

Oznacza to, że ciało może poruszać się ze stałą szybkością i mimo to mieć niezerowe przyspieszenie. Dzieje się tak na przykład w ruchu jednostajnym po okręgu, ponieważ kierunek wektora prędkości stale się zmienia.

Składowa styczna i normalna przyspieszenia

W opisie naturalnym przyspieszenie można rozłożyć na składową styczną i normalną:

\[ \vec a = a_{\tau}\vec\tau + a_n\vec n. \] 

Składowa styczna odpowiada za zmianę szybkości:

\[ a_{\tau} = \frac{dv}{dt}. \]

Składowa normalna odpowiada za zmianę kierunku ruchu:

\[ a_n = \frac{v^2}{\rho}, \]

gdzie \(\rho\) jest promieniem krzywizny toru. Wektor składowej normalnej jest skierowany ku środkowi krzywizny toru.

Szczegółowe wyznaczanie prędkości i przyspieszenia omówimy w osobnym artykule poświęconym kinematyce punktu materialnego.

Podstawowa klasyfikacja ruchów

Ruch można klasyfikować według różnych kryteriów.

Ze względu na kształt toru

Ze względu na zmianę szybkości

Ze względu na ruch bryły

Ruch jednostajny nie musi być ruchem prostoliniowym. W ruchu jednostajnym po okręgu szybkość jest stała, ale wektor prędkości zmienia kierunek.

Ruch względny

Opis ruchu zależy od wybranego układu odniesienia. Jeżeli dwa punkty \(A\) i \(B\) poruszają się względem tego samego układu, położenie punktu \(B\) względem punktu \(A\) wynosi:

\[ \vec r_{B/A} = \vec r_B-\vec r_A. \] 

Prędkość względna jest równa:

\[ \vec v_{B/A} = \vec v_B-\vec v_A. \] 

Przykładem może być pasażer idący w jadącym pociągu. Jego prędkość względem wagonu jest inna niż prędkość względem ziemi.

Szczegółowy opis ruchu względnego i złożonego wymaga uwzględnienia sposobu poruszania się układów odniesienia i zostanie rozwinięty w osobnym materiale.

Wykresy wielkości kinematycznych

Ruch można przedstawiać za pomocą wykresów pokazujących zależność wielkości kinematycznych od czasu, na przykład:

Nachylenie wykresu położenia względem czasu odpowiada prędkości:

\[ v_x = \frac{dx}{dt}. \]

Nachylenie wykresu prędkości odpowiada przyspieszeniu:

\[ a_x = \frac{dv_x}{dt}. \]

Pole algebraiczne pod wykresem prędkości w przedziale czasu odpowiada zmianie położenia:

\[ \Delta x = \int_{t_1}^{t_2}v_x(t)\,dt. \]

Pole pod wykresem wartości szybkości odpowiada natomiast przebytej drodze.

Wykres \(x(t)\) nie jest zwykle torem ruchu. Pokazuje on zmianę jednej współrzędnej w czasie. Tor przedstawia natomiast położenie ciała w przestrzeni, na przykład zależność \(y(x)\).

Jednostki i kontrola wymiarów

WielkośćOznaczenieJednostka SI
Czas\(t\)\(\mathrm{s}\)
Położenie\(\vec r\), \(x\), \(y\), \(z\)\(\mathrm{m}\)
Droga\(s\)\(\mathrm{m}\)
Przemieszczenie\(\Delta\vec r\)\(\mathrm{m}\)
Prędkość\(\vec v\)\(\mathrm{m/s}\)
Przyspieszenie\(\vec a\)\(\mathrm{m/s^2}\)
Promień krzywizny\(\rho\)\(\mathrm{m}\)

Zależności kinematyczne można sprawdzać za pomocą analizy jednostek. Dla prędkości:

\[ [\vec v] = \frac{[\vec r]}{[t]} = \frac{\mathrm{m}}{\mathrm{s}}. \]

Dla przyspieszenia:

\[ [\vec a] = \frac{[\vec v]}{[t]} = \frac{\mathrm{m}}{\mathrm{s^2}}. \]

Dla przyspieszenia normalnego:

\[ \left[\frac{v^2}{\rho}\right] = \frac{\mathrm{m^2/s^2}}{\mathrm{m}} = \frac{\mathrm{m}}{\mathrm{s^2}}. \]

Przeliczając kilometry na godzinę na metry na sekundę, korzystamy z zależności:

\[ 1\ \frac{\mathrm{km}}{\mathrm{h}} = \frac{1}{3{,}6}\ \frac{\mathrm{m}}{\mathrm{s}}. \]

Przykładowo:

\[ 72\ \frac{\mathrm{km}}{\mathrm{h}} = 20\ \frac{\mathrm{m}}{\mathrm{s}}. \]

Najczęstsze błędy

Podsumowanie

Kinematyka opisuje ruch bez analizowania jego przyczyn. Do pełnego opisu należy określić układ odniesienia oraz zależność położenia od czasu.

Położenie punktu można przedstawić wektorowo:

\[ \vec r=\vec r(t), \]

za pomocą współrzędnych:

\[ x=x(t),\qquad y=y(t),\qquad z=z(t), \]

albo naturalnie, za pomocą współrzędnej mierzonej wzdłuż znanego toru:

\[ s=s(t). \]

Prędkość jest pochodną położenia względem czasu:

\[ \vec v=\frac{d\vec r}{dt}, \]

a przyspieszenie jest pochodną prędkości:

\[ \vec a = \frac{d\vec v}{dt} = \frac{d^2\vec r}{dt^2}. \]

Droga jest długością przebytego toru, natomiast przemieszczenie jest wektorem łączącym położenie początkowe z końcowym. Prędkość określa zarówno tempo, jak i kierunek zmiany położenia, a przyspieszenie opisuje każdą zmianę wektora prędkości — także samą zmianę jego kierunku.

Prognozowanie – podstawowe pojęcia, metody i błędy prognoz

Prognozowanie polega na formułowaniu przewidywań dotyczących przyszłych wartości lub stanów badanego zjawiska. Prognozy mogą dotyczyć między innymi sprzedaży, cen, wielkości produkcji, zużycia energii, liczby ludności, poziomu inflacji, kursów walutowych albo zapotrzebowania na określone produkty.

Przyszłości nie można poznać ze stuprocentową pewnością. Można jednak wykorzystać dostępne dane, wiedzę o badanym procesie oraz odpowiednio dobrane metody, aby ograniczyć niepewność i przygotować możliwie wiarygodne przewidywania.

W artykule przedstawiamy podstawowe pojęcia związane z prognozowaniem, najważniejsze rodzaje metod prognostycznych, prognozy punktowe i przedziałowe, błędy ex ante i ex post oraz popularne miary służące do oceny dokładności prognoz.


Spis treści

Czym jest prognozowanie?

Prognozowanie jest procesem formułowania sądów dotyczących przyszłego przebiegu badanego zjawiska. Rezultatem tego procesu jest prognoza, czyli przewidywana wartość albo przewidywany stan zjawiska w określonym momencie lub okresie przyszłości.

Prognoza nie jest zwykłym przypuszczeniem. Powinna wynikać z określonych przesłanek, takich jak:

Przedsiębiorstwo może na przykład prognozować sprzedaż na kolejny miesiąc, aby odpowiednio zaplanować produkcję i wielkość zapasów. Operator systemu elektroenergetycznego przewiduje przyszłe zapotrzebowanie na energię, a instytucja publiczna może sporządzać prognozy liczby ludności lub wpływów podatkowych.

Prognozowanie nie usuwa niepewności dotyczącej przyszłości, ale pozwala ją opisać, ograniczyć i uwzględnić podczas podejmowania decyzji.

Podstawowe pojęcia prognostyczne

Zmienna prognozowana

Zmienna prognozowana jest wielkością, której przyszłe wartości chcemy przewidzieć. Może nią być na przykład miesięczna sprzedaż przedsiębiorstwa, cena produktu, liczba pasażerów, wielkość zużycia energii albo stopa inflacji.

Rzeczywistą wartość zmiennej w okresie \(t\) oznacza się najczęściej symbolem \(y_t\). Prognozowaną wartość tej zmiennej można natomiast zapisać jako \(\widehat{y}_t\).

Jeżeli prognoza dotyczy przyszłego okresu \(T+h\) i została sporządzona przy wykorzystaniu informacji dostępnych do okresu \(T\), można zastosować dokładniejsze oznaczenie:

\[ \widehat{y}_{T+h\mid T} \]

Symbol \(h\) oznacza liczbę okresów, o które prognoza wybiega w przyszłość.

Moment sporządzenia prognozy

Moment sporządzenia prognozy jest chwilą, w której na podstawie dostępnych informacji przewidujemy przyszłe wartości zmiennej. W tym momencie znamy już obserwacje pochodzące z przeszłości, ale nie znamy jeszcze wartości, które wystąpią w kolejnych okresach.

Jeżeli znamy miesięczną sprzedaż przedsiębiorstwa do końca czerwca i na tej podstawie przewidujemy sprzedaż lipcową, koniec czerwca jest momentem sporządzenia prognozy.

Horyzont prognozy

Horyzont prognozy określa, jak daleko w przyszłość sięga przewidywanie. Jeżeli dysponujemy danymi do czerwca, prognoza na lipiec ma horyzont jednego okresu, a prognoza na grudzień — sześciu okresów.

Prognoza na jeden okres naprzód może zostać zapisana jako:

\[ \widehat{y}_{T+1\mid T} \]

Prognoza na \(h\) okresów naprzód ma natomiast postać:

\[ \widehat{y}_{T+h\mid T} \]

Prognozy o dłuższym horyzoncie są zazwyczaj bardziej niepewne. Im dalej wybiegamy w przyszłość, tym większe znaczenie mogą mieć zdarzenia i zmiany, których nie dało się uwzględnić w momencie budowy modelu.

Okres prognozy

Okres prognozy obejmuje wszystkie przyszłe okresy, dla których wyznaczamy prognozowane wartości. Jeżeli na podstawie danych do 2025 roku przewidujemy wartości dla lat 2026 – 2030, okresem prognozy są lata od 2026 do 2030.

Rodzaje prognoz

Prognozy można klasyfikować według różnych kryteriów. Najczęściej bierze się pod uwagę długość horyzontu, sposób przedstawienia wyniku, charakter wykorzystywanych informacji oraz cel prognozy.

Prognozy krótko-, średnio- i długookresowe

Rodzaj prognozyCharakterystykaPrzykład
krótkookresowaDotyczy najbliższych okresów i zwykle opiera się przede wszystkim na aktualnych danych.Prognoza zużycia energii na następny dzień.
średniookresowaObejmuje okres pozwalający planować działania operacyjne lub inwestycyjne.Prognoza sprzedaży na następny rok.
długookresowaDotyczy odległej przyszłości i jest związana z dużą niepewnością.Prognoza liczby ludności za 20 lat.

Granice pomiędzy poszczególnymi rodzajami prognoz zależą od charakteru badanego zjawiska. Prognoza roczna może być bardzo długa w przypadku cen energii na rynku dnia następnego, ale stosunkowo krótka w analizach demograficznych.

Prognozy ilościowe i jakościowe

Prognozy ilościowe mają postać liczbową. Możemy na przykład przewidywać, że sprzedaż wyniesie 1200 sztuk albo że stopa inflacji osiągnie 3,5%.

Prognozy jakościowe opisują przyszły stan zjawiska bez konieczności wskazywania jednej konkretnej wartości liczbowej. Można na przykład przewidywać wzrost popytu, pogorszenie koniunktury albo pojawienie się określonego rozwiązania technologicznego.

Prognozy punktowe i przedziałowe

Prognoza punktowa wskazuje jedną przewidywaną wartość zmiennej. Prognoza przedziałowa określa natomiast zakres wartości, w którym z ustalonym prawdopodobieństwem albo przy określonym poziomie ufności powinna znaleźć się przyszła obserwacja.

Prognozy warunkowe i bezwarunkowe

Prognoza warunkowa jest formułowana przy założeniu spełnienia określonych warunków. Możemy na przykład przewidywać sprzedaż przy założeniu, że cena produktu nie ulegnie zmianie, a wydatki na reklamę pozostaną na dotychczasowym poziomie.

Prognoza bezwarunkowa nie jest bezpośrednio przedstawiana jako wynik zależny od jednego wyraźnie wskazanego scenariusza. W praktyce również ona opiera się jednak na pewnych założeniach dotyczących trwałości zależności występujących w danych.

Metody prognozowania

Metody prognozowania można podzielić na dwie główne grupy: metody jakościowe oraz metody ilościowe.

Metody jakościowe

Metody jakościowe wykorzystują wiedzę, doświadczenie i oceny ludzi. Stosuje się je szczególnie wtedy, gdy:

Do metod jakościowych należą między innymi:

Metody ilościowe

Metody ilościowe wykorzystują dane liczbowe oraz modele matematyczne i statystyczne. Można je podzielić na metody oparte na szeregach czasowych oraz metody przyczynowo-skutkowe.

Metody szeregów czasowych

W metodach szeregów czasowych prognoza powstaje przede wszystkim na podstawie wcześniejszych wartości tej samej zmiennej. Poszukuje się prawidłowości występujących w jej zmianach w czasie.

Do tej grupy należą między innymi:

Metody przyczynowo-skutkowe

Metody przyczynowo-skutkowe zakładają, że prognozowana zmienna zależy od innych czynników. Sprzedaż określonego produktu może na przykład zależeć od jego ceny, dochodów konsumentów, temperatury, wydatków na reklamę i aktywności konkurencji.

Zależności te mogą być przedstawiane za pomocą modeli regresji lub bardziej rozbudowanych modeli ekonometrycznych. Przykładowy model liniowy można zapisać jako:

\[ y_t=\beta_0+\beta_1x_{1t}+\beta_2x_{2t}+\ldots+\beta_kx_{kt}+\varepsilon_t \]

W równaniu:

Aby wykorzystać taki model do prognozowania, należy znać albo również prognozować przyszłe wartości zmiennych objaśniających.

Szereg czasowy i jego składniki

Szereg czasowy jest zbiorem obserwacji uporządkowanych według czasu. Może przedstawiać na przykład miesięczną sprzedaż, dzienne ceny energii, kwartalny produkt krajowy brutto albo roczną liczbę mieszkańców.

MiesiącSprzedaż
Styczeń100
Luty104
Marzec109
Kwiecień107
Maj115

W szeregu czasowym mogą występować następujące składniki:

Poziom

Poziom określa typową lub przeciętną wartość, wokół której zmienia się badane zjawisko.

Trend

Trend jest długookresowym kierunkiem zmian. Może mieć charakter rosnący, malejący albo w przybliżeniu stały.

Przykładowy liniowy model trendu ma postać:

\[ y_t=\alpha_0+\alpha_1t+\varepsilon_t \]

Parametr \(\alpha_1\) określa przeciętną zmianę wartości zmiennej przypadającą na jeden okres.

Sezonowość

Sezonowość oznacza regularne wahania powtarzające się w określonych odstępach. Przykładem może być wzrost sprzedaży lodów w miesiącach letnich albo większe zużycie energii elektrycznej w określonych porach doby.

Wahania cykliczne

Wahania cykliczne są zmianami trwającymi zwykle dłużej niż jeden sezon. Mogą być związane na przykład z kolejnymi fazami cyklu koniunkturalnego. W przeciwieństwie do sezonowości nie muszą powtarzać się w idealnie regularnych odstępach.

Składnik losowy

Składnik losowy obejmuje nieregularne zmiany, których model nie potrafi wyjaśnić. Mogą one wynikać z przypadkowych zdarzeń, błędów pomiaru albo działania czynników nieuwzględnionych w analizie.

Dobór metody prognozowania powinien uwzględniać budowę szeregu. Innego modelu potrzebuje szereg o stałym poziomie, innego szereg z trendem, a jeszcze innego dane zawierające jednocześnie trend i sezonowość.

Prognoza punktowa i przedziałowa

Prognoza punktowa

Prognoza punktowa wskazuje jedną konkretną wartość, która według modelu najlepiej opisuje przyszły poziom badanego zjawiska.

Jeżeli przedsiębiorstwo przewiduje, że w następnym miesiącu sprzeda 1200 sztuk produktu, wartość 1200 jest prognozą punktową:

\[ \widehat{y}_{T+1\mid T}=1200 \]

Prognoza punktowa jest łatwa do wykorzystania w planowaniu, ale sama nie pokazuje niepewności przewidywania.

Prognoza przedziałowa

Prognoza przedziałowa wskazuje zakres, w którym przy określonym poziomie ufności lub prawdopodobieństwa powinna znaleźć się przyszła wartość zmiennej.

Możemy na przykład przewidywać, że sprzedaż wyniesie 1200 sztuk, a przedział prognozy obejmuje wartości od 1080 do 1320 sztuk:

\[ 1080\leq y_{T+1}\leq1320 \]

Przedział prognozy jest zazwyczaj tym szerszy, im:

Warto odróżnić przedział prognozy dla przyszłej obserwacji od przedziału ufności wyznaczanego dla przeciętnego poziomu zjawiska. Przedział prognozy pojedynczej przyszłej wartości jest zazwyczaj szerszy, ponieważ oprócz niepewności oszacowania modelu uwzględnia również losową zmienność przyszłej obserwacji.

Prognoza punktowa, średni błąd ex ante, przedział prognozy oraz błąd ex post na wykresie szeregu czasowego
Prognoza punktowa jest wyznaczana przed poznaniem przyszłej wartości. Średni błąd ex ante opisuje jej przewidywaną niepewność, natomiast błąd ex post można obliczyć dopiero po poznaniu wartości zrealizowanej.

Średni błąd ex ante i błąd ex post

W prognozowaniu należy odróżnić ocenę niepewności dokonywaną przed poznaniem przyszłej wartości od błędu, który można obliczyć dopiero po jej zrealizowaniu.

Średni błąd prognozy ex ante

Średni błąd prognozy ex ante jest oszacowaniem niepewności prognozy dokonywanym w momencie jej sporządzania. Przyszła wartość zmiennej nie jest jeszcze znana, dlatego nie można wtedy obliczyć rzeczywistej różnicy pomiędzy prognozą a wartością zrealizowaną.

W zależności od zastosowanego modelu oblicza się średni błąd prognozy ex ante, standardowy błąd prognozy albo inną miarę opisującą przewidywaną skalę możliwego odchylenia przyszłej wartości od prognozy punktowej.

Jeżeli średni błąd ex ante dla prognozy na okres \(T+h\) oznaczymy symbolem \(s_{T+h}^{*}\), zakres odpowiadający jednemu takiemu błędowi można przedstawić jako:

\[ \widehat{y}_{T+h\mid T}\pm s_{T+h}^{*} \]

Sam średni błąd ex ante nie musi jeszcze wyznaczać właściwego przedziału prognozy. Przy budowie przedziału wykorzystuje się odpowiedni kwantyl rozkładu, zależny między innymi od przyjętego poziomu ufności.

Ogólną postać przedziału prognozy można zapisać jako:

\[ \widehat{y}_{T+h\mid T} \pm q_{1-\alpha/2}\,s_{T+h}^{*} \]

Symbol \(q_{1-\alpha/2}\) oznacza odpowiedni kwantyl przyjętego rozkładu. W modelach, w których wykorzystuje się rozkład Studenta, przedział może mieć postać:

\[ \widehat{y}_{T+h\mid T} \pm t_{1-\alpha/2,\nu}\,s_{T+h}^{*} \]

Symbol \(\nu\) oznacza liczbę stopni swobody. Dla powszechnie stosowanych poziomów ufności, takich jak 90%, 95% lub 99%, wartość odpowiedniego kwantyla jest większa od jedności. Z tego powodu przedział prognozy jest szerszy od zakresu wyznaczonego jedynie przez dodanie i odjęcie pojedynczego średniego błędu ex ante.

Średni błąd ex ante nie mówi, jaki dokładnie błąd zostanie popełniony. Określa jedynie przewidywaną skalę niepewności prognozy przed poznaniem wartości rzeczywistej.

Błąd ex post

Błąd ex post można obliczyć dopiero wtedy, gdy znana jest rzeczywista wartość prognozowanej zmiennej. Jest on różnicą pomiędzy wartością zrealizowaną a wartością prognozowaną.

\[ e_{T+h}=y_{T+h}-\widehat{y}_{T+h\mid T} \]

Jeżeli prognozowana sprzedaż wynosiła 1200 sztuk, a rzeczywista sprzedaż osiągnęła 1150 sztuk, błąd prognozy wynosi:

\[ e=1150-1200=-50 \]

Przy zastosowanej tutaj konwencji ujemny błąd oznacza, że prognoza była zbyt wysoka. Dodatni błąd oznacza natomiast, że wartość rzeczywista okazała się wyższa od prognozowanej.

W niektórych opracowaniach błąd definiuje się odwrotnie:

\[ e_{T+h}=\widehat{y}_{T+h\mid T}-y_{T+h} \]

Dlatego przed interpretacją znaku błędu zawsze należy sprawdzić, którą konwencję przyjęto.

Najważniejsze miary błędów prognoz

Pojedynczy błąd ex post informuje o dokładności jednej prognozy. Aby ocenić model na podstawie większej liczby okresów, stosuje się miary podsumowujące zbiór błędów.

Załóżmy, że oceniane są prognozy dla \(n\) okresów, a błędy mają postać:

\[ e_t=y_t-\widehat{y}_t \]

Średni błąd — ME

Średni błąd, oznaczany jako ME, obliczamy według wzoru:

\[ ME=\frac{1}{n}\sum_{t=1}^{n}e_t \]

Miara ta pomaga ocenić, czy prognozy są systematycznie zawyżane lub zaniżane. Wartość dodatnia oznacza przy przyjętej konwencji, że wartości rzeczywiste były przeciętnie wyższe od prognozowanych.

Wartość \(ME\) bliska zeru nie musi jednak oznaczać dużej dokładności. Błędy dodatnie i ujemne mogą się wzajemnie znosić.

Średni błąd bezwzględny — MAE

Średni błąd bezwzględny ma postać:

\[ MAE=\frac{1}{n}\sum_{t=1}^{n}|e_t| \]

Wskaźnik \(MAE\) informuje, o ile jednostek prognoza różni się przeciętnie od wartości rzeczywistej. Jeżeli \(MAE=5\), prognozy odbiegały od wartości rzeczywistych przeciętnie o 5 jednostek.

Średni błąd kwadratowy — MSE

Średni błąd kwadratowy obliczamy jako:

\[ MSE=\frac{1}{n}\sum_{t=1}^{n}e_t^2 \]

Podnoszenie błędów do kwadratu sprawia, że duże odchylenia mają szczególnie silny wpływ na wartość wskaźnika. \(MSE\) jest więc użyteczny wtedy, gdy duże błędy prognoz są szczególnie niepożądane.

Pierwiastek średniego błędu kwadratowego — RMSE

Pierwiastek średniego błędu kwadratowego ma postać:

\[ RMSE=\sqrt{\frac{1}{n}\sum_{t=1}^{n}e_t^2} \]

W przeciwieństwie do \(MSE\), wskaźnik \(RMSE\) jest wyrażony w tych samych jednostkach co prognozowana zmienna. Dzięki temu jego wartość jest łatwiejsza do interpretacji.

Średni bezwzględny błąd procentowy — MAPE

Średni bezwzględny błąd procentowy obliczamy jako:

\[ MAPE= \frac{100\%}{n} \sum_{t=1}^{n} \left| \frac{e_t}{y_t} \right| \]

Jeżeli \(MAPE=8\%\), prognozy różniły się od wartości rzeczywistych przeciętnie o około 8%.

Miara ta nie powinna być stosowana bezkrytycznie. Gdy wartość rzeczywista \(y_t\) wynosi zero, błąd procentowy jest nieokreślony. Gdy \(y_t\) jest bardzo bliskie zeru, wskaźnik może przyjmować bardzo wysokie i mało użyteczne wartości.

Przykład obliczenia błędów

OkresWartość rzeczywista \(y_t\)Prognoza \(\widehat{y}_t\)Błąd \(e_t\)\(|e_t|\)
11009822
211010644
3105108-33
411511233

Dla powyższych danych średni błąd wynosi:

\[ ME=\frac{2+4-3+3}{4}=1{,}5 \]

Średni błąd bezwzględny jest równy:

\[ MAE=\frac{2+4+3+3}{4}=3 \]

Prognozy różniły się zatem od wartości rzeczywistych przeciętnie o 3 jednostki, a dodatnia wartość \(ME\) wskazuje, że były przeciętnie nieco zaniżone.

Etapy procesu prognozowania

Poprawne prognozowanie nie sprowadza się do wybrania wzoru i wykonania obliczeń. Jest procesem składającym się z kilku etapów.

1. Określenie celu prognozy

Najpierw należy ustalić, co ma być prognozowane, na jaki okres oraz do jakiej decyzji prognoza będzie wykorzystywana.

Innego modelu może wymagać prognoza sprzedaży na następny dzień, a innego planowanie rozwoju przedsiębiorstwa na kilka lat.

2. Zebranie danych

Dane powinny być kompletne, wiarygodne i porównywalne. Należy sprawdzić między innymi:

3. Wstępna analiza danych

Przed wyborem modelu warto sporządzić wykres szeregu czasowego i sprawdzić, czy występują w nim:

4. Wybór metody prognozowania

Metoda powinna odpowiadać charakterowi danych i celowi prognozy. Szereg bez trendu i sezonowości może być prognozowany za pomocą prostej metody naiwnej lub średniej. Szereg zawierający trend wymaga modelu uwzględniającego kierunek zmian, a dane sezonowe — modelu uwzględniającego powtarzalne wahania.

5. Oszacowanie parametrów modelu

Na tym etapie na podstawie danych historycznych wyznacza się parametry modelu, na przykład współczynniki trendu, parametry regresji albo współczynniki wygładzania.

6. Ocena modelu

Model należy sprawdzić nie tylko pod względem dopasowania do danych historycznych, ale przede wszystkim pod względem zdolności do przewidywania nowych obserwacji.

7. Sporządzenie prognozy

Po zaakceptowaniu modelu wyznacza się prognozę punktową oraz — jeżeli jest to możliwe — przedział prognozy i średni błąd ex ante.

8. Kontrola i aktualizacja prognoz

Po pojawieniu się nowych danych prognozę należy porównać z wartością rzeczywistą. Pozwala to obliczyć błędy ex post i sprawdzić, czy model nadal prawidłowo opisuje badane zjawisko.

Modele powinny być okresowo aktualizowane, szczególnie gdy zmieniają się warunki gospodarcze, technologiczne, społeczne lub prawne.

Jak oceniać i porównywać modele?

Model najlepiej dopasowany do danych historycznych nie musi tworzyć najdokładniejszych prognoz. Bardzo rozbudowany model może niemal idealnie odtwarzać wcześniejsze obserwacje, ale słabo przewidywać przyszłe wartości.

Zjawisko nadmiernego dostosowania modelu do danych historycznych nazywa się przeuczeniem albo nadmiernym dopasowaniem.

Zbiór uczący i zbiór testowy

Dane można podzielić na:

W przypadku szeregów czasowych należy zachować kolejność obserwacji. Nie można losowo mieszać danych z różnych okresów, ponieważ mogłoby to prowadzić do wykorzystania informacji z przyszłości podczas prognozowania przeszłości.

Walidacja krocząca

Jedną z metod oceny prognoz jest walidacja krocząca. Polega ona na wielokrotnym powtarzaniu następujących czynności:

  1. oszacowanie modelu na podstawie dostępnych obserwacji,
  2. sporządzenie prognozy na kolejny okres,
  3. porównanie prognozy z wartością rzeczywistą,
  4. dołączenie nowej obserwacji do zbioru danych,
  5. ponowne oszacowanie lub aktualizacja modelu.

Taki sposób postępowania przypomina rzeczywiste wykorzystywanie modelu, ponieważ każda prognoza powstaje wyłącznie na podstawie informacji dostępnych w chwili jej sporządzania.

Prognozy aktualizowane i prognozy wielookresowe

Podczas historycznej weryfikacji modelu należy rozróżnić prognozy jednookresowe stale aktualizowane od prognozy wielookresowej wyznaczonej ze stałego momentu. Obie procedury mogą prowadzić do obliczenia błędów ex post, ale odpowiadają innym sposobom praktycznego wykorzystywania prognoz.

Prognozy jednookresowe stale aktualizowane

W pierwszym podejściu w każdym kolejnym okresie sporządzamy prognozę tylko na jeden krok naprzód:

\[ \widehat{y}_{t+1\mid t} \]

Po pojawieniu się rzeczywistej wartości \(y_{t+1}\) dołączamy ją do zbioru informacji i na jej podstawie wyznaczamy kolejną prognozę:

\[ \widehat{y}_{t+2\mid t+1} \]

Początek prognozy przesuwa się więc wraz z każdą nową obserwacją. Model jest na bieżąco korygowany, dzięki czemu korzysta z najświeższych dostępnych danych.

Takie podejście jest właściwe wtedy, gdy w praktyce prognoza jest regularnie aktualizowana i potrzebujemy przede wszystkim przewidywania na jeden najbliższy okres. Może ono jednak sprawiać wrażenie, że prosty model jest bardzo dokładny, ponieważ każda następna prognoza korzysta już z ostatniej znanej wartości rzeczywistej.

Prognoza wielookresowa ze stałego momentu

W drugim podejściu wybieramy jeden moment \(T\), w którym kończy się zbiór informacji wykorzystywany do sporządzenia prognozy. Następnie wyznaczamy całą ścieżkę przyszłych wartości:

\[ \widehat{y}_{T+1\mid T}, \quad \widehat{y}_{T+2\mid T}, \quad \ldots, \quad \widehat{y}_{T+H\mid T} \]

Podczas obliczania tych prognoz nie wykorzystujemy rzeczywistych wartości \(y_{T+1},y_{T+2},\ldots,y_{T+H}\). Jeżeli wykonujemy analizę na danych historycznych, wartości te są wprawdzie znane badaczowi, ale traktujemy je tak, jakby nie były dostępne w momencie sporządzania prognozy.

Dopiero po wyznaczeniu całej ścieżki porównujemy prognozy z wartościami rzeczywistymi i obliczamy błędy ex post dla kolejnych horyzontów:

\[ e_{T+h\mid T} = y_{T+h} - \widehat{y}_{T+h\mid T} \]

Takie postępowanie pozwala sprawdzić, jak model zachowałby się w rzeczywistej sytuacji, w której w momencie \(T\) trzeba przygotować prognozę na kilka kolejnych okresów bez możliwości bieżącego korygowania jej za pomocą nowych obserwacji.

Wraz z wydłużaniem horyzontu prognoza może coraz bardziej oddalać się od rzeczywistego przebiegu zjawiska. Jest to szczególnie widoczne w prostych metodach, które utrzymują stały poziom, ostatni przyrost albo ostatnie tempo zmian.

Porównanie prognoz jednookresowych stale aktualizowanych z prognozą wielookresową wyznaczoną ze stałego momentu
Prognozy jednookresowe są aktualizowane po każdej nowej obserwacji. Prognoza wielookresowa powstaje w jednym momencie i rozwija się bez wykorzystywania późniejszych danych rzeczywistych.

Na lewym wykresie każda prognoza dotyczy tylko jednego okresu naprzód i zostaje zaktualizowana po poznaniu kolejnej wartości rzeczywistej. Z tego powodu linia prognoz może pozostawać stosunkowo blisko obserwowanego szeregu.

Na prawym wykresie prognoza zostaje sporządzona tylko raz, w momencie \(T\), a następnie rozwija się przez kilka okresów bez korzystania z nowych danych. Porównanie jej z rzeczywistymi wartościami pokazuje, jak wraz z horyzontem może narastać rozbieżność między prognozą a faktycznym przebiegiem zjawiska.

Niewielkie błędy prognoz jednookresowych nie oznaczają automatycznie, że model będzie równie dokładny przy prognozowaniu kilku okresów naprzód bez bieżącej aktualizacji.

Prognozy wygasłe i wartości dopasowane

Prognozę wyznaczoną dla okresu historycznego, którego rzeczywista wartość jest już znana, można określić jako prognozę wygasłą. Pozwala ona obliczyć błąd ex post i sprawdzić, jak model poradziłby sobie z przewidywaniem obserwacji, której — zgodnie z przyjętą procedurą — nie wykorzystywano podczas sporządzania prognozy.

Prognoz wygasłych nie należy utożsamiać ze zwykłymi wartościami dopasowanymi modelu. Jeżeli model oszacowano na podstawie całego szeregu, a następnie porównano jego wartości teoretyczne z obserwacjami wykorzystanymi wcześniej do estymacji, oceniamy przede wszystkim dopasowanie modelu w próbie.

Wartość dopasowana dla danego okresu może pośrednio wykorzystywać informacje pochodzące z późniejszych obserwacji, jeżeli obserwacje te posłużyły do oszacowania parametrów modelu. Nie jest więc zawsze równoważna prognozie, którą można byłoby rzeczywiście sporządzić w analizowanym momencie.

Ocena błędów według horyzontu prognozy

Najpełniejszą ocenę modelu można uzyskać, powtarzając prognozowanie wielookresowe dla wielu historycznych momentów rozpoczęcia prognozy. Dla każdego z nich wyznacza się prognozy na jeden, dwa, trzy lub więcej okresów naprzód.

Następnie błędy grupuje się według horyzontu i oblicza osobne miary dokładności:

\[
MAE_1,\quad MAE_2,\quad MAE_3,\quad \ldots,\quad MAE_H
\]

Dzięki temu można sprawdzić, czy model dobrze prognozuje jedynie najbliższy okres, czy zachowuje użyteczną dokładność również przy dłuższym horyzoncie.

Porównanie z prostą metodą odniesienia

Bardziej zaawansowany model powinien być porównany z prostą metodą odniesienia, na przykład metodą naiwną prognozowania. Jeżeli rozbudowany model nie tworzy prognoz dokładniejszych niż założenie, że przyszła wartość będzie równa ostatniej obserwacji, jego stosowanie może nie być uzasadnione.

Bardziej zaawansowany model powinien być porównany z prostą metodą odniesienia, na przykład metodą naiwną. Jeżeli rozbudowany model nie tworzy prognoz dokładniejszych niż założenie, że przyszła wartość będzie równa ostatniej obserwacji, jego stosowanie może nie być uzasadnione.

Model powinien być na tyle rozbudowany, na ile wymaga tego badane zjawisko, ale nie bardziej skomplikowany, niż jest to potrzebne do uzyskania wiarygodnych prognoz.

Cechy dobrej prognozy

Dobra prognoza powinna być:

Nie istnieje metoda, która byłaby najlepsza dla wszystkich danych i wszystkich zastosowań. Skuteczność modelu zależy od charakteru badanego zjawiska, jakości informacji, długości horyzontu oraz celu prognozy.

Podsumowanie

Prognozowanie polega na przewidywaniu przyszłych wartości lub stanów zjawisk na podstawie dostępnych danych, wiedzy i rozpoznanych prawidłowości. Rezultatem procesu jest prognoza, która może mieć charakter punktowy albo przedziałowy.

Metody prognozowania dzielimy między innymi na jakościowe i ilościowe. Metody ilościowe mogą opierać się na analizie szeregów czasowych albo na zależnościach przyczynowo-skutkowych pomiędzy zmiennymi.

Średni błąd ex ante opisuje przewidywaną niepewność przed poznaniem przyszłej wartości. Nie jest on rzeczywistym błędem konkretnej prognozy. Błąd ex post można obliczyć dopiero po pojawieniu się wartości zrealizowanej:

\[
e_{T+h}=y_{T+h}-\widehat{y}_{T+h\mid T}
\]

Do oceny większej liczby prognoz wykorzystuje się między innymi wskaźniki \(ME\), \(MAE\), \(MSE\), \(RMSE\) i \(MAPE\). Modele powinny być porównywane na tych samych danych testowych oraz z odpowiednio dobraną prostą metodą odniesienia.

Najbardziej skomplikowany model nie zawsze daje najlepsze wyniki. Dobra metoda prognozowania powinna być dopasowana do charakteru danych, dostatecznie dokładna, możliwa do interpretacji oraz użyteczna z punktu widzenia podejmowanej decyzji.



Rzut poziomy – tor ruchu, czas lotu, zasięg i prędkość

Rzut poziomy jest ruchem ciała, któremu nadano początkową prędkość skierowaną poziomo. Przy pominięciu oporu powietrza ruch ten można rozłożyć na dwa niezależne ruchy: jednostajny ruch w poziomie oraz swobodne spadanie w pionie. Połączenie obu ruchów sprawia, że ciało porusza się po torze parabolicznym.

Opis rzutu poziomego wykorzystuje wiadomości przedstawione w artykułach Ruch prostoliniowy – położenie, prędkość, przyspieszenie i równania ruchu oraz Swobodne spadanie i rzut pionowy.

Czym jest rzut poziomy?

Rzut poziomy zachodzi wtedy, gdy ciało otrzymuje początkową prędkość skierowaną poziomo, a następnie porusza się pod wpływem grawitacji.

W chwili wyrzutu prędkość ma wyłącznie składową poziomą:

\[ v_{x0}=v_0, \qquad v_{y0}=0. \]

Natychmiast po rozpoczęciu ruchu grawitacja powoduje narastanie pionowej składowej prędkości. Ciało przemieszcza się więc jednocześnie w prawo i w dół.

Prędkość jest dokładnie pozioma tylko w chwili wyrzutu. W każdej późniejszej chwili pojawia się składowa pionowa, dlatego wektor prędkości jest nachylony w dół.

Założenia modelu

W podstawowym modelu rzutu poziomego przyjmujemy, że:

Wartość przyspieszenia ziemskiego przyjmujemy jako:

\[ g\approx 9{,}81\ \mathrm{m/s^2}. \]

W prostszych zadaniach może być stosowane przybliżenie:

\[ g\approx 10\ \mathrm{m/s^2}. \]

Układ współrzędnych

W tym artykule przyjmujemy początek układu współrzędnych w punkcie wyrzutu. Oś \(x\) jest skierowana poziomo w prawo, a oś \(y\) pionowo w dół.

Przy takim wyborze osi:

\[ a_x=0, \qquad a_y=g. \]

Zarówno pionowa współrzędna \(y\), jak i pionowa składowa prędkości \(v_y\), mają podczas spadania wartości dodatnie.

Wektor przyspieszenia można zapisać jako:

\[ \vec a = g\vec e_y, \] 

gdzie \(\vec e_y\) jest wersorem osi skierowanej pionowo w dół.

Można również przyjąć oś pionową skierowaną w górę. Wtedy \(a_y=-g\), a pionowe położenie ciała podczas spadania maleje. Oba układy są poprawne, ale nie wolno mieszać wzorów zapisanych dla różnych zwrotów osi.

Niezależność ruchu poziomego i pionowego

Rzut poziomy można rozpatrywać jako złożenie dwóch niezależnych ruchów odbywających się w tym samym czasie.

Ruch w kierunku poziomym

W kierunku poziomym nie działa żadne przyspieszenie:

\[ a_x=0. \]

Składowa pozioma prędkości jest więc stała:

\[ v_x=v_0. \]

Położenie poziome rośnie liniowo z czasem:

\[ x=v_0t. \]

Ruch w kierunku pionowym

W chwili wyrzutu pionowa składowa prędkości jest równa zeru:

\[ v_{y0}=0. \]

Ruch pionowy jest więc swobodnym spadaniem:

\[ a_y=g, \] \[ v_y=gt, \] \[ y=\frac{gt^2}{2}. \]

Grawitacja nie zmienia składowej poziomej prędkości, a nadana poziomo prędkość początkowa nie wpływa na ruch pionowy.

Dwa ciała puszczone jednocześnie z tej samej wysokości spadną w tym samym czasie, nawet jeżeli jedno zostanie puszczone bez prędkości początkowej, a drugie wyrzucone poziomo. Jest to prawdziwe przy pominięciu oporu powietrza.

Równania rzutu poziomego

Dla początku układu współrzędnych umieszczonego w punkcie wyrzutu równania położenia mają postać:

\[ x(t)=v_0t, \] \[ y(t)=\frac{gt^2}{2}. \]

Składowe prędkości:

\[ v_x(t)=v_0, \] \[ v_y(t)=gt. \]

Składowe przyspieszenia:

\[ a_x=0, \] \[ a_y=g. \]

W zapisie wektorowym:

\[ \vec r(t) = v_0t\,\vec e_x + \frac{gt^2}{2}\,\vec e_y, \] \[ \vec v(t) = v_0\,\vec e_x + gt\,\vec e_y, \] \[ \vec a(t) = g\,\vec e_y. \] 
Schemat rzutu poziomego z torem parabolicznym oraz składowymi prędkości poziomej i pionowej
W rzucie poziomym \(v_x=v_0\) pozostaje stałe, natomiast \(v_y=gt\) rośnie wraz z czasem. Wektor prędkości wypadkowej jest styczny do toru ruchu.

Równanie toru

Równanie toru otrzymujemy przez wyeliminowanie czasu z równań położenia.

Z równania ruchu poziomego:

\[ x=v_0t \]

wyznaczamy czas:

\[ t=\frac{x}{v_0}. \]

Podstawiamy tę zależność do równania ruchu pionowego:

\[ y = \frac{g}{2} \left(\frac{x}{v_0}\right)^2. \]

Otrzymujemy:

\[ \boxed{ y = \frac{g}{2v_0^2}x^2 } \]

Jest to równanie paraboli. W przyjętym układzie współrzędna \(y\) rośnie w dół, dlatego tor na rysunku odchyla się coraz bardziej ku powierzchni ziemi.

Większa prędkość początkowa \(v_0\) powoduje, że dla tej samej wartości \(x\) pionowe odchylenie \(y\) jest mniejsze. Tor staje się bardziej płaski.

Składowe prędkości

Wektor prędkości jest sumą składowej poziomej i pionowej:

\[ \vec v = \vec v_x+\vec v_y. \] 

Składowa pozioma pozostaje stała:

\[ v_x=v_0. \]

Składowa pionowa rośnie liniowo wraz z czasem:

\[ v_y=gt. \]

Położenie ciała na torze zmienia się więc w taki sposób, że wektor prędkości jest coraz bardziej nachylony w dół.

Rozkład prędkości ciała w rzucie poziomym na składową poziomą, pionową i prędkość wypadkową
Prędkość wypadkowa \(\vec v\) jest sumą składowej poziomej \(\vec v_x\) i pionowej \(\vec v_y\). Składowa \(v_x\) jest stała, a \(v_y\) rośnie wraz z czasem.

Wartość i kierunek prędkości

Wartość prędkości, czyli szybkość ciała, obliczamy z twierdzenia Pitagorasa:

\[ v = |\vec v| = \sqrt{v_x^2+v_y^2}. \]

Po podstawieniu składowych:

\[ \boxed{ v = \sqrt{v_0^2+g^2t^2} } \]

Niech \(\alpha\) oznacza kąt między wektorem prędkości a kierunkiem poziomym. Wtedy:

\[ \tan\alpha = \frac{v_y}{v_x}. \]

Zatem:

\[ \boxed{ \tan\alpha = \frac{gt}{v_0} } \]

Korzystając z zależności \(t=x/v_0\), można również zapisać:

\[ \tan\alpha = \frac{gx}{v_0^2}. \]

Kąt \(\alpha\) rośnie wraz z czasem. Oznacza to, że kierunek ruchu staje się coraz bardziej pionowy.

Czas lotu

Załóżmy, że ciało zostało wyrzucone z wysokości \(h\) nad poziomym podłożem. W chwili uderzenia pionowe przemieszczenie wynosi:

\[ y=h. \]

Z równania ruchu pionowego:

\[ h = \frac{gt_{\mathrm{lotu}}^2}{2}. \]

Po przekształceniu otrzymujemy:

\[ \boxed{ t_{\mathrm{lotu}} = \sqrt{\frac{2h}{g}} } \]

Czas lotu zależy wyłącznie od wysokości \(h\) i przyspieszenia ziemskiego \(g\). Nie występuje w nim pozioma prędkość początkowa \(v_0\).

Zwiększenie prędkości poziomej nie wydłuża czasu spadania. Ciało poruszy się dalej w poziomie, ale uderzy w podłoże po takim samym czasie jak ciało wyrzucone wolniej z tej samej wysokości.

Zasięg rzutu poziomego

Zasięg \(L\) jest poziomą odległością między punktem wyrzutu a punktem uderzenia.

W kierunku poziomym ciało porusza się ze stałą prędkością \(v_0\), dlatego:

\[ L = v_0t_{\mathrm{lotu}}. \]

Po podstawieniu czasu lotu:

\[ \boxed{ L = v_0\sqrt{\frac{2h}{g}} } \]

Zasięg jest wprost proporcjonalny do prędkości początkowej:

\[ L\sim v_0. \] 

Jeżeli prędkość początkową zwiększymy dwukrotnie, zasięg także zwiększy się dwukrotnie.

Zasięg jest proporcjonalny do pierwiastka z wysokości:

\[ L\sim\sqrt{h}. \]

Aby dwukrotnie zwiększyć zasięg przez samą zmianę wysokości, trzeba zwiększyć wysokość czterokrotnie.

Schemat rzutu poziomego z zaznaczoną wysokością, czasem lotu i zasięgiem
Dla rzutu poziomego z wysokości \(h\) czas lotu wynosi \(t_{\mathrm{lotu}}=\sqrt{\frac{2h}{g}}\), a zasięg \(L=v_0\sqrt{\frac{2h}{g}}\).

Prędkość w chwili uderzenia

W chwili uderzenia składowa pozioma nadal wynosi:

\[ v_x=v_0. \]

Składowa pionowa jest równa:

\[ v_y = gt_{\mathrm{lotu}}. \]

Po podstawieniu czasu lotu:

\[ v_y = g\sqrt{\frac{2h}{g}} = \sqrt{2gh}. \]

Wartość prędkości w chwili uderzenia wynosi:

\[ v_{\mathrm{uderzenia}} = \sqrt{v_0^2+v_y^2}. \]

Zatem:

\[ \boxed{ v_{\mathrm{uderzenia}} = \sqrt{v_0^2+2gh} } \]

Kąt uderzenia spełnia zależność:

\[ \boxed{ \tan\alpha_{\mathrm{uderzenia}} = \frac{\sqrt{2gh}}{v_0} } \]

Wykresy ruchu

Położenie poziome \(x(t)\)

Położenie poziome opisuje funkcja liniowa:

\[ x(t)=v_0t. \]

Nachylenie wykresu jest stałe i równe prędkości poziomej \(v_0\).

Położenie pionowe \(y(t)\)

Położenie pionowe opisuje funkcja kwadratowa:

\[ y(t)=\frac{gt^2}{2}. \]

Przy osi skierowanej dodatnio w dół wykres rozpoczyna się w punkcie \(y=0\) i odchyla w dół coraz szybciej.

Prędkość pozioma \(v_x(t)\)

Składowa pozioma jest stała:

\[ v_x(t)=v_0. \]

Jej wykres jest poziomą linią.

Prędkość pionowa \(v_y(t)\)

Składowa pionowa rośnie liniowo:

\[ v_y(t)=gt. \]

Nachylenie wykresu jest równe przyspieszeniu ziemskiemu \(g\).

Wykresy położenia poziomego i pionowego oraz składowych prędkości w rzucie poziomym
Przy osi \(y\) skierowanej w dół położenie \(y(t)=\frac{gt^2}{2}\) i prędkość \(v_y(t)=gt\) mają wartości dodatnie i rosną wraz z czasem.

Rzut na inny poziom końcowy

Wzór:

\[ t_{\mathrm{lotu}} = \sqrt{\frac{2h}{g}} \]

obowiązuje wtedy, gdy \(h\) jest pionową odległością między punktem wyrzutu a poziomem uderzenia, a pionowa prędkość początkowa wynosi zero.

W bardziej ogólnym przypadku, gdy początek osi nie znajduje się w punkcie wyrzutu, stosujemy:

\[ x(t)=x_0+v_0t, \] \[ y(t)=y_0+\frac{gt^2}{2}. \]

Jeżeli znamy współrzędną poziomu końcowego \(y_k\), czas uderzenia wyznaczamy z warunku:

\[ y_k = y_0+\frac{gt^2}{2}. \]

Po znalezieniu czasu zasięg obliczamy z równania ruchu poziomego.

Przykłady obliczeniowe

Przykład 1. Czas lotu i zasięg

Ciało wyrzucono poziomo z wysokości:

\[ h=20\ \mathrm{m} \]

z prędkością:

\[ v_0=8\ \mathrm{m/s}. \]

Czas lotu:

\[ t_{\mathrm{lotu}} = \sqrt{\frac{2h}{g}} = \sqrt{\frac{2\cdot20}{9{,}81}} \approx 2{,}02\ \mathrm{s}. \]

Zasięg:

\[ L = v_0t_{\mathrm{lotu}} = 8\cdot2{,}02 \approx 16{,}2\ \mathrm{m}. \]

Ciało spadnie po około \(2{,}02\ \mathrm{s}\) w odległości około \(16{,}2\ \mathrm{m}\) od pionu przechodzącego przez punkt wyrzutu.

Przykład 2. Prędkość w chwili uderzenia

Dla danych z poprzedniego przykładu pionowa składowa prędkości w chwili uderzenia wynosi:

\[ v_y = gt_{\mathrm{lotu}} \approx 9{,}81\cdot2{,}02 \approx 19{,}8\ \mathrm{m/s}. \]

Wartość prędkości wypadkowej:

\[ v = \sqrt{v_0^2+v_y^2} = \sqrt{8^2+19{,}8^2} \approx 21{,}4\ \mathrm{m/s}. \]

Kąt nachylenia prędkości:

\[ \tan\alpha = \frac{19{,}8}{8} \approx 2{,}48. \]

Stąd:

\[ \alpha\approx68^\circ. \]

Przykład 3. Wyznaczenie wysokości

Ciało wyrzucone poziomo z prędkością:

\[ v_0=10\ \mathrm{m/s} \]

spadło w odległości:

\[ L=30\ \mathrm{m}. \]

Czas lotu wynika z ruchu poziomego:

\[ t_{\mathrm{lotu}} = \frac{L}{v_0} = \frac{30}{10} = 3\ \mathrm{s}. \]

Wysokość:

\[ h = \frac{gt_{\mathrm{lotu}}^2}{2} = \frac{9{,}81\cdot3^2}{2} \approx 44{,}1\ \mathrm{m}. \]

Przykład 4. Porównanie dwóch prędkości początkowych

Z tej samej wysokości wyrzucono poziomo dwa ciała. Pierwsze otrzymało prędkość \(v_0\), a drugie \(2v_0\).

Czasy lotu są jednakowe:

\[ t_1=t_2=\sqrt{\frac{2h}{g}}. \]

Zasięg pierwszego ciała:

\[ L_1=v_0t_{\mathrm{lotu}}. \]

Zasięg drugiego ciała:

\[ L_2 = 2v_0t_{\mathrm{lotu}} = 2L_1. \] 

Dwukrotnie większa prędkość pozioma daje dwukrotnie większy zasięg, ale nie zmienia czasu lotu.

Jednostki i kontrola wymiarów

WielkośćOznaczenieJednostka SI
Położenie poziome\(x\)\(\mathrm{m}\)
Położenie pionowe\(y\)\(\mathrm{m}\)
Wysokość\(h\)\(\mathrm{m}\)
Zasięg\(L\)\(\mathrm{m}\)
Czas\(t\)\(\mathrm{s}\)
Prędkość\(v\), \(v_x\), \(v_y\)\(\mathrm{m/s}\)
Przyspieszenie\(g\)\(\mathrm{m/s^2}\)

Kontrola wymiarów wzoru na czas lotu:

\[ \left[ \sqrt{\frac{2h}{g}} \right] = \sqrt{ \frac{\mathrm{m}} {\mathrm{m/s^2}} } = \sqrt{\mathrm{s^2}} = \mathrm{s}. \]

Kontrola wzoru na zasięg:

\[ [L] = \left[ v_0\sqrt{\frac{2h}{g}} \right] = \frac{\mathrm{m}}{\mathrm{s}}\cdot\mathrm{s} = \mathrm{m}. \]

Kontrola równania toru:

\[ \left[ \frac{g}{2v_0^2}x^2 \right] = \frac{\mathrm{m/s^2}}{\mathrm{m^2/s^2}} \cdot \mathrm{m^2} = \mathrm{m}. \]

Najczęstsze błędy

Składowych prędkości nie dodajemy algebraicznie. Ponieważ są prostopadłe, wartość prędkości wypadkowej wynosi:

\[ v=\sqrt{v_x^2+v_y^2}, \]

a nie \(v=v_x+v_y\).

Podsumowanie

Rzut poziomy jest złożeniem jednostajnego ruchu poziomego i swobodnego spadania.

Dla osi \(x\) skierowanej w prawo i osi \(y\) skierowanej w dół:

\[ x=v_0t, \qquad y=\frac{gt^2}{2}, \] \[ v_x=v_0, \qquad v_y=gt, \] \[ a_x=0, \qquad a_y=g. \]

Tor ruchu jest parabolą:

\[ y=\frac{g}{2v_0^2}x^2. \]

Wartość prędkości wynosi:

\[ v=\sqrt{v_0^2+g^2t^2}. \]

Dla rzutu z wysokości \(h\) czas lotu i zasięg są równe:

\[ t_{\mathrm{lotu}} = \sqrt{\frac{2h}{g}}, \] \[ L = v_0\sqrt{\frac{2h}{g}}. \]

Czas lotu nie zależy od poziomej prędkości początkowej. Zwiększenie \(v_0\) powoduje wzrost zasięgu, ale nie zmienia czasu spadania z tej samej wysokości.

Testowanie hipotez dotyczących wskaźnika struktury w dwóch populacjach

Testowanie różnicy dwóch wskaźników struktury pozwala sprawdzić, czy udział określonej cechy różni się między dwiema niezależnymi populacjami. W artykule omawiamy test dla dwóch proporcji: hipotezy jedno- i dwustronne, proporcję łączną, warunki stosowania przybliżenia normalnego, wartość p, przedział ufności dla różnicy proporcji oraz najczęstsze błędy.

W wielu badaniach nie interesuje nas już tylko odpowiedź na pytanie, czy pojedynczy wskaźnik struktury jest równy wartości referencyjnej. Chcemy porównać dwa odsetki: skuteczność dwóch metod leczenia, udział produktów wadliwych na dwóch liniach produkcyjnych, poparcie dla danego rozwiązania w dwóch grupach respondentów albo odsetek zdających egzamin w dwóch rocznikach.

Jeżeli obie grupy są niezależne, podstawowym narzędziem jest test dla różnicy dwóch proporcji. Pozwala on ocenić, czy zaobserwowana różnica między wskaźnikami z prób może wynikać z losowego doboru obserwacji, czy też stanowi podstawę do uznania, że proporcje w populacjach są różne.

Artykuł jest rozwinięciem zagadnienia testowania hipotez dotyczących wskaźnika struktury w jednej populacji. Ogólne zasady interpretacji hipotezy zerowej, poziomu istotności i wartości p omówiono natomiast w artykule Testowanie hipotez statystycznych.

Co testujemy?

Niech:

\[ p_1 \] 

oznacza wskaźnik struktury w pierwszej populacji, a:

\[ p_2 \] 

— wskaźnik struktury w drugiej populacji.

Interesuje nas różnica:

\[ \delta=p_1-p_2 \]

Najczęściej sprawdzamy, czy różnica ta wynosi zero:

\[ \delta=0 \]

co jest równoważne hipotezie o równości wskaźników struktury:

\[ p_1=p_2 \]

W praktyce możemy badać na przykład, czy skuteczność nowej metody jest większa niż metody dotychczasowej, czy odsetek wadliwych produktów jest różny na dwóch liniach produkcyjnych albo czy deklarowane poparcie dla danego rozwiązania różni się między dwiema grupami społecznymi.

Dwie proporcje w dwóch populacjach

Załóżmy, że z pierwszej populacji pobrano próbę o liczebności \(n_1\), w której zaobserwowano \(m_1\) sukcesów. Z drugiej populacji pobrano niezależną próbę o liczebności \(n_2\), w której zaobserwowano \(m_2\) sukcesów.

Wskaźniki struktury z prób wynoszą:

\[ \hat p_1= \frac{m_1}{n_1} \]
\[ \hat p_2= \frac{m_2}{n_2} \]

Naturalnym estymatorem różnicy proporcji jest:

\[ \hat p_1-\hat p_2 \] 

Jeżeli wynik jest dodatni, wskaźnik z pierwszej próby jest większy. Jeżeli jest ujemny, większy wskaźnik uzyskano w drugiej próbie.

Sama różnica z prób nie przesądza jednak jeszcze o różnicy w populacjach. Jej interpretacja wymaga uwzględnienia liczebności prób oraz losowej zmienności wskaźników struktury.

Próby niezależne a próby zależne

Klasyczny test dla dwóch wskaźników struktury dotyczy prób niezależnych. Oznacza to, że wynik uzyskany w pierwszej grupie nie jest bezpośrednio powiązany z wynikiem uzyskanym w drugiej grupie.

Inaczej wygląda sytuacja, gdy te same osoby odpowiadają przed i po interwencji albo gdy obserwacje są zestawione w pary. Wtedy próby są zależne i nie należy stosować zwykłego testu dwóch niezależnych proporcji. W takich przypadkach odpowiednią procedurą jest zazwyczaj test McNemara.

Nie stosuj testu dla prób niezależnych do tych samych osób badanych dwa razy

Jeżeli porównujemy odpowiedzi tych samych osób „przed” i „po”, wyniki nie są niezależne. Test dla dwóch niezależnych wskaźników struktury nie uwzględnia tej zależności i może prowadzić do błędnych wniosków.

Hipotezy jedno- i dwustronne

Najczęściej formułujemy hipotezę zerową:

\[ H_0:p_1=p_2 \]

czyli równoważnie:

\[ H_0:p_1-p_2=0 \]

Test dwustronny

\[ H_0:p_1=p_2 \] \[ H_1:p_1\neq p_2 \] 

Stosujemy go, gdy interesuje nas każda różnica między proporcjami, niezależnie od jej kierunku.

Test prawostronny

\[ H_0:p_1=p_2 \] \[ H_1:p_1>p_2 \]

Stosujemy go na przykład wtedy, gdy pierwsza metoda ma być skuteczniejsza od drugiej albo gdy oczekujemy większego odsetka sukcesów w pierwszej populacji.

Test lewostronny

\[ H_0:p_1=p_2 \] \[ H_1:p_1

Ten wariant wybieramy wtedy, gdy interesuje nas wyłącznie niższy wskaźnik w pierwszej populacji.

Kierunek testu wybieramy przed analizą danych

Nie należy wybierać testu prawostronnego tylko dlatego, że po zebraniu danych otrzymaliśmy \(\hat p_1>\hat p_2\). Hipoteza alternatywna powinna wynikać z pytania badawczego, projektu badania albo wcześniej przyjętego założenia.

Warunki stosowania testu normalnego

Klasyczny test z dla dwóch proporcji korzysta z przybliżenia normalnego. Jest ono odpowiednie wtedy, gdy w obu próbach oczekiwana liczba sukcesów i porażek przy założeniu prawdziwości \(H_0\) jest dostatecznie duża.

W teście równości proporcji pod hipotezą zerową zakładamy:

\[ p_1=p_2=p \]

Nieznaną wspólną proporcję oszacujemy za pomocą proporcji łącznej \(\hat p\), którą omówimy dokładniej za chwilę. Warunki dużej próby zapisujemy wtedy jako:

\[ n_1\hat p \] 

oraz:

\[ n_1(1-\hat p) \] 

dla pierwszej próby, a także:

\[ n_2\hat p \] 

oraz:

\[ n_2(1-\hat p) \] 

dla drugiej próby.

W praktyce często wymaga się, aby każda z tych wartości była nie mniejsza niż 5 lub — ostrożniej — niż 10. Są to reguły praktyczne, które mają zapewnić dostatecznie dobre przybliżenie normalne.

Statystyka testowa dla dwóch proporcji

W klasycznym teście równości dwóch proporcji statystyka ma postać:

\[ Z= \frac{ \hat p_1-\hat p_2 }{ \sqrt{ \hat p(1-\hat p) \left( \frac{1}{n_1}+\frac{1}{n_2} \right) } } \]

gdzie \(\hat p\) jest proporcją łączną obliczoną z obu prób:

\[ \hat p= \frac{ m_1+m_2 }{ n_1+n_2 } \]

Przy prawdziwości hipotezy zerowej statystyka \(Z\) ma w przybliżeniu rozkład normalny standaryzowany:

\[ Z\approx N(0,1) \]

Jeżeli testujemy ogólnie różnicę równą ustalonej wartości \(\delta_0\), statystyka ma postać:

\[ Z= \frac{ (\hat p_1-\hat p_2)-\delta_0 }{ SE_0 } \] 

W typowym teście równości proporcji przyjmujemy jednak:

\[ \delta_0=0 \]

Dlaczego w teście stosujemy proporcję łączną?

To jeden z najważniejszych elementów testu dla dwóch proporcji. Jeżeli hipoteza zerowa mówi, że:

\[ H_0:p_1=p_2 \]

to pod hipotezą zerową obie próby mają wspólną proporcję sukcesów. Nie znamy jej, więc estymujemy ją, łącząc informacje z obu prób:

\[ \hat p= \frac{ m_1+m_2 }{ n_1+n_2 } \]

Proporcja łączna jest zatem oszacowaniem wspólnego parametru obowiązującego przy prawdziwości \(H_0\).

Dlatego w mianowniku statystyki testowej występuje:

\[ \sqrt{ \hat p(1-\hat p) \left( \frac{1}{n_1}+\frac{1}{n_2} \right) } \]

Nie używaj proporcji łącznej automatycznie w każdym wzorze

Proporcję łączną \(\hat p\) stosujemy w klasycznym teście równości dwóch proporcji, ponieważ wariancję różnicy wyznaczamy przy założeniu prawdziwości \(H_0:p_1=p_2\).

W przedziale ufności dla różnicy proporcji używamy natomiast osobno \(\hat p_1\) oraz \(\hat p_2\), ponieważ przedział ma opisywać nieznaną rzeczywistą różnicę na podstawie danych, a nie zakładać jej równość zero.

Infografika przedstawiająca dwie niezależne próby, proporcje z prób, proporcję łączną oraz wzór statystyki Z dla testu równości dwóch proporcji
W klasycznym teście równości proporcji proporcję łączną stosujemy pod hipotezą zerową, że oba wskaźniki struktury są równe.

Obszar krytyczny i wartość p

Po obliczeniu wartości \(z_{\text{obs}}\) podejmujemy decyzję przez porównanie jej z odpowiednią wartością krytyczną rozkładu normalnego albo na podstawie wartości p.

Test dwustronny

\[ H_1:p_1\neq p_2 \] 

Odrzucamy \(H_0\), gdy:

\[ |z_{\text{obs}}|> u_{\alpha/2} \]

Wartość p wynosi:

\[ p= 2P(Z\geq |z_{\text{obs}}|) \]

Test prawostronny

\[ H_1:p_1>p_2 \]

Odrzucamy \(H_0\), gdy:

\[ z_{\text{obs}}> u_{\alpha} \]

Wartość p wynosi:

\[ p= P(Z\geq z_{\text{obs}}) \]

Test lewostronny

\[ H_1:p_1

Odrzucamy \(H_0\), gdy:

\[ z_{\text{obs}} -u_{\alpha} \]

Wartość p wynosi:

\[ p= P(Z\leq z_{\text{obs}}) \]

Wartości krytyczne można odczytać z tablicy rozkładu normalnego i kalkulatora.

Przedział ufności dla różnicy proporcji

Testowanie hipotezy o równości proporcji warto uzupełnić przedziałem ufności dla różnicy:

\[ p_1-p_2 \] 

Klasyczny przybliżony dwustronny przedział ufności \(1-\alpha\) ma postać:

\[ \left( \hat p_1-\hat p_2 - u_{\alpha/2} \sqrt{ \frac{ \hat p_1(1-\hat p_1) }{ n_1 } + \frac{ \hat p_2(1-\hat p_2) }{ n_2 } }, \; \hat p_1-\hat p_2 + u_{\alpha/2} \sqrt{ \frac{ \hat p_1(1-\hat p_1) }{ n_1 } + \frac{ \hat p_2(1-\hat p_2) }{ n_2 } } \right) \]

W tym wzorze używamy osobnych proporcji z prób \(\hat p_1\) i \(\hat p_2\), a nie proporcji łącznej. Przedział ufności nie zakłada bowiem prawdziwości hipotezy o równości proporcji.

W teście dwustronnym na poziomie istotności \(\alpha\) klasyczny przedział ufności \(1-\alpha\) i test prowadzą do zgodnych wniosków:

Infografika przedstawiająca test dwustronny, prawostronny i lewostronny dla dwóch proporcji oraz przedział ufności dla różnicy wskaźników struktury
W teście dwustronnym liczba zero należąca do przedziału ufności dla różnicy proporcji oznacza brak podstaw do odrzucenia hipotezy o równości wskaźników.

Przykład: porównanie dwóch wskaźników struktury

Porównujemy skuteczność dwóch niezależnych metod przygotowania do egzaminu. W pierwszej grupie 84 spośród 120 osób zdały egzamin. W drugiej grupie zdało 63 spośród 120 osób.

Chcemy sprawdzić, czy wskaźniki zdawalności w populacjach są różne.

Krok 1. Wskaźniki z prób

\[ n_1=120 \]
\[ m_1=84 \]
\[ \hat p_1= \frac{84}{120} = 0{,}70 \]
\[ n_2=120 \]
\[ m_2=63 \]
\[ \hat p_2= \frac{63}{120} = 0{,}525 \]

Różnica wskaźników z prób wynosi:

\[ \hat p_1-\hat p_2= 0{,}70-0{,}525= 0{,}175 \]

Krok 2. Hipotezy

\[ H_0:p_1=p_2 \]
\[ H_1:p_1\neq p_2 \] 

Krok 3. Proporcja łączna

\[ \hat p= \frac{ 84+63 }{ 120+120 } = \frac{147}{240} = 0{,}6125 \]

Krok 4. Warunki dużej próby

\[ 120\cdot0{,}6125=73{,}5 \]
\[ 120\cdot(1-0{,}6125)=46{,}5 \]

W obu próbach oczekiwane liczby sukcesów i porażek są duże, dlatego przybliżenie normalne jest uzasadnione.

Krok 5. Statystyka testowa

\[ z_{\text{obs}}= \frac{ 0{,}70-0{,}525 }{ \sqrt{ 0{,}6125(1-0{,}6125) \left( \frac{1}{120}+\frac{1}{120} \right) } } \]
\[ z_{\text{obs}} \approx 2{,}801 \]

Krok 6. Decyzja

Przyjmujemy poziom istotności:

\[ \alpha=0{,}05 \]

Dla testu dwustronnego wartość krytyczna wynosi:

\[ u_{0{,}025}=1{,}96 \]

Ponieważ:

\[ |2{,}801|>1{,}96 \]

odrzucamy hipotezę zerową o równości wskaźników struktury.

Krok 7. Wartość p

\[ p\approx0{,}005 \]

Wartość p jest mniejsza od 0,05, więc potwierdza decyzję o odrzuceniu \(H_0\).

Krok 8. Przedział ufności dla różnicy

\[ 0{,}175 \pm 1{,}96 \sqrt{ \frac{ 0{,}70(1-0{,}70) }{ 120 } + \frac{ 0{,}525(1-0{,}525) }{ 120 } } \] 
\[ (0{,}053;\;0{,}297) \]

Liczba \(0\) nie należy do przedziału ufności, co jest zgodne z wynikiem testu.

Wniosek: dane dostarczają podstaw do stwierdzenia, że wskaźnik zdawalności w pierwszej populacji jest różny od wskaźnika zdawalności w drugiej populacji. Ponieważ różnica \(\hat p_1-\hat p_2\) jest dodatnia, wyniki wskazują na wyższy wskaźnik zdawalności w pierwszej grupie.

Małe próby i metody dokładne

Gdy liczebności prób są małe albo liczba sukcesów lub porażek jest niewielka, klasyczne przybliżenie normalne może być niedokładne. Dotyczy to zwłaszcza sytuacji, w których proporcje są bliskie 0 lub 1.

W takich przypadkach można stosować dokładne procedury oparte na rozkładzie hipergeometrycznym, z których najbardziej znanym przykładem jest dokładny test Fishera dla tablicy kontyngencji \(2\times2\).

Test Fishera jest szczególnie przydatny przy małych liczebnościach w kontroli jakości, badaniach medycznych, analizie rzadkich zdarzeń oraz wszędzie tam, gdzie klasyczny warunek dużej próby nie jest spełniony.

Małe próby: nie mieszaj procedur

Test z dla dwóch proporcji jest metodą przybliżoną opartą na rozkładzie normalnym. Dokładny test Fishera jest inną procedurą, opartą na rozkładzie hipergeometrycznym. Nie należy zastępować jednego testu drugim wyłącznie przez zmianę wzoru na błąd standardowy.

Kalkulatory i programy statystyczne

W klasycznych zadaniach wartości krytyczne można odczytywać z tablic rozkładu normalnego. W praktyce wygodniejsze jest korzystanie z kalkulatorów i programów statystycznych, które wyznaczają statystykę testową, wartość p oraz przedział ufności dla różnicy proporcji.

Do kontroli wartości krytycznych przyda się tablica rozkładu normalnego i kalkulator. Test dla dwóch proporcji i dokładny test Fishera są dostępne między innymi w Excelu, SPSS, Statistica, Gretl, R oraz innych pakietach statystycznych.

Najczęstsze błędy

Podsumowanie

Testowanie różnicy dwóch wskaźników struktury pozwala sprawdzić, czy dwie niezależne populacje różnią się udziałem określonej cechy.

Testowanie hipotez dotyczących wskaźnika struktury w jednej populacji

Testowanie hipotez dotyczących wskaźnika struktury pozwala sprawdzić, czy udział określonej cechy w populacji jest zgodny z wartością deklarowaną, wymaganą przez normę albo przyjętą w badaniu jako punkt odniesienia. W artykule omawiamy test dla jednej proporcji: hipotezy jedno- i dwustronne, warunki stosowania przybliżenia normalnego, statystykę testową, wartość p oraz dokładny test dwumianowy dla małych prób.

Wskaźnik struktury opisuje udział jednostek posiadających określoną cechę w całej zbiorowości. Może oznaczać odsetek klientów zadowolonych z usługi, udział produktów wadliwych, procent osób popierających danego kandydata, udział zdających egzamin albo częstość występowania określonego zdarzenia.

W praktyce sama wartość wskaźnika obliczona w próbie nie wystarcza. Jeżeli w próbie 59% respondentów deklaruje określoną opinię, pojawia się pytanie: czy rzeczywisty udział w całej populacji różni się od 50%, czy zaobserwowana różnica może być wyłącznie skutkiem losowego doboru próby?

Do odpowiedzi na takie pytania służy test dotyczący wskaźnika struktury w jednej populacji. Najczęściej opiera się on na przybliżeniu normalnym rozkładu proporcji z próby. W przypadku małych prób lub bardzo małych albo bardzo dużych proporcji stosuje się dokładne metody dwumianowe.

Ogólne pojęcia, takie jak hipoteza zerowa, poziom istotności, wartość p, błędy pierwszego i drugiego rodzaju oraz moc testu, zostały omówione w artykule Testowanie hipotez statystycznych.

Co testujemy?

Niech \(X\) oznacza liczbę jednostek posiadających badaną cechę w próbie o liczebności \(n\). Jeżeli każda obserwacja może przyjąć jedną z dwóch wartości, na przykład „tak” albo „nie”, „wadliwy” albo „niewadliwy”, „zdany” albo „niezdany”, w populacji interesuje nas prawdopodobieństwo sukcesu:

\[ p=P(\text{sukces}) \]

Parametr \(p\) jest jednocześnie populacyjnym wskaźnikiem struktury, czyli udziałem jednostek posiadających badaną cechę.

Jeżeli w próbie zaobserwowano \(m\) sukcesów, wskaźnik struktury z próby wynosi:

\[ \hat p= \frac{m}{n} \]

Test polega na porównaniu nieznanego wskaźnika populacyjnego \(p\) z wartością hipotetyczną:

\[ p_0 \] 

Wartość \(p_0\) może wynikać z deklaracji producenta, normy jakościowej, wcześniejszego badania, wymaganego poziomu skuteczności albo założenia przyjętego przez badacza.

Wskaźnik struktury i proporcja populacyjna

W języku statystyki opisowej mówi się zwykle o wskaźniku struktury, czyli udziale jednostek określonego typu w całej zbiorowości. W języku rachunku prawdopodobieństwa i statystyki matematycznej ten sam parametr oznacza się najczęściej literą \(p\) i nazywa proporcją populacyjną.

Przykłady:

W dalszej części symbolem \(p\) oznaczamy wskaźnik struktury w populacji, a symbolem \(\hat p\) — wskaźnik obliczony w próbie.

Hipotezy jedno- i dwustronne

Najczęściej hipoteza zerowa ma postać:

\[ H_0:p=p_0 \]

Postać hipotezy alternatywnej zależy od pytania badawczego.

Test dwustronny

\[ H_0:p=p_0 \] \[ H_1:p\neq p_0 \] 

Stosujemy go wtedy, gdy interesuje nas każda różnica względem wartości referencyjnej. Przykład: producent deklaruje, że odsetek wadliwych produktów wynosi 3%, a kontrola jakości sprawdza, czy rzeczywisty odsetek różni się od tej wartości w dowolnym kierunku.

Test prawostronny

\[ H_0:p=p_0 \] \[ H_1:p>p_0 \]

Test prawostronny ma zastosowanie na przykład wtedy, gdy chcemy sprawdzić, czy odsetek produktów wadliwych przekracza dopuszczalny poziom albo czy poparcie dla określonego rozwiązania wzrosło ponad wartość referencyjną.

Test lewostronny

\[ H_0:p=p_0 \] \[ H_1:p

Test lewostronny stosujemy wtedy, gdy interesuje nas spadek proporcji, na przykład zmniejszenie udziału produktów wadliwych poniżej wcześniejszego poziomu.

Kierunek testu wybieramy przed analizą danych

W testach jednostronnych kierunek hipotezy alternatywnej powinien wynikać z celu badania. Nie należy wybierać testu prawostronnego albo lewostronnego dopiero po zobaczeniu, czy wskaźnik z próby \(\hat p\) okazał się większy czy mniejszy od \(p_0\).

Kiedy można stosować test normalny?

Klasyczny test dla jednej proporcji opiera się na przybliżeniu normalnym rozkładu proporcji z próby. Jest ono wiarygodne wtedy, gdy oczekiwana liczba sukcesów i oczekiwana liczba porażek przy założeniu prawdziwości hipotezy zerowej są dostatecznie duże:

\[ np_0 \] 

oraz:

\[ n(1-p_0) \]

W podręcznikach spotyka się różne praktyczne progi, na przykład co najmniej 5 albo co najmniej 10 dla obu tych wartości. Nie jest to bezwzględne prawo matematyczne, lecz reguła oceny jakości przybliżenia normalnego.

Ważne: warunek sprawdzamy względem wartości hipotetycznej \(p_0\), ponieważ to właśnie przy założeniu \(H_0\) wyznaczamy rozkład statystyki testowej.

Przykładowo, gdy testujemy hipotezę \(H_0:p=0{,}03\) na próbie \(n=50\), otrzymujemy:

\[ np_0=50\cdot0{,}03=1{,}5 \]

Przybliżenie normalne jest wtedy słabe, nawet gdy wskaźnik z próby przypadkowo okaże się bliski wartości 0,50. W takim przypadku właściwszy będzie dokładny test dwumianowy.

Statystyka testowa dla jednej proporcji

Jeżeli warunki przybliżenia normalnego są spełnione, statystyka testowa ma postać:

\[ Z= \frac{ \hat p-p_0 }{ \sqrt{ \frac{ p_0(1-p_0) }{ n } } } \] 

Przy prawdziwości hipotezy zerowej statystyka \(Z\) ma w przybliżeniu rozkład normalny standaryzowany:

\[ Z\approx N(0,1) \]

Po podstawieniu danych z próby otrzymujemy wartość obserwowaną:

\[ z_{\text{obs}}= \frac{ \hat p-p_0 }{ \sqrt{ \frac{ p_0(1-p_0) }{ n } } } \] 

Równoważnie można zapisać statystykę przez liczbę sukcesów \(m\):

\[ Z= \frac{ m-np_0 }{ \sqrt{ np_0(1-p_0) } } \]

Obie postaci prowadzą do dokładnie tej samej wartości statystyki testowej.

Dlaczego we wzorze występuje \(p_0\), a nie \(\hat p\)?

To jeden z najczęstszych błędów przy testowaniu wskaźnika struktury. We wzorze na statystykę testową występuje:

\[ \sqrt{ \frac{ p_0(1-p_0) }{ n } } \] 

a nie:

\[ \sqrt{ \frac{ \hat p(1-\hat p) }{ n } } \] 

Powód jest prosty: w teście sprawdzamy, jak bardzo wynik próby odbiega od tego, czego należałoby oczekiwać, gdyby hipoteza zerowa była prawdziwa. Jeżeli \(H_0:p=p_0\) jest prawdziwa, wariancja wskaźnika z próby wynosi:

\[ Var(\hat p)= \frac{ p_0(1-p_0) }{ n } \] 

Dlatego błąd standardowy w teście obliczamy właśnie z użyciem \(p_0\).

Test i przedział ufności mogą wykorzystywać różne błędy standardowe

W klasycznym teście dla proporcji stosujemy błąd standardowy liczony pod hipotezą zerową, czyli z \(p_0\). Natomiast w typowym przedziale ufności dla proporcji wykorzystuje się zwykle wartość z próby \(\hat p\), ponieważ przedział ma opisywać nieznany parametr na podstawie zaobserwowanych danych.

Nie wolno automatycznie przenosić jednego wzoru do drugiej procedury.

Schemat porównujący proporcję referencyjną p zero i proporcję z próby p z daszkiem we wzorze na błąd standardowy testu dla jednej proporcji
Wartość p₀ opisuje rozkład statystyki przy prawdziwości hipotezy zerowej, dlatego występuje w mianowniku testu.

Obszar krytyczny i wartość krytyczna

Po obliczeniu wartości statystyki porównujemy ją z odpowiednią wartością krytyczną rozkładu normalnego.

Test dwustronny

Odrzucamy hipotezę zerową, gdy:

\[ |z_{\text{obs}}|> u_{\alpha/2} \]

Test prawostronny

Odrzucamy hipotezę zerową, gdy:

\[ z_{\text{obs}}> u_{\alpha} \]

Test lewostronny

Odrzucamy hipotezę zerową, gdy:

\[ z_{\text{obs}} -u_{\alpha} \]

Dla poziomu istotności \(\alpha=0{,}05\) wartości krytyczne wynoszą:

\[ u_{0{,}025}\approx1{,}96 \]

dla testu dwustronnego oraz:

\[ u_{0{,}05}\approx1{,}645 \]

dla testu jednostronnego. Wartości można odczytać z tablicy rozkładu normalnego i kalkulatora.

Infografika przedstawiająca test dla jednej proporcji, wzór statystyki Z z proporcją referencyjną p zero, warunek dużej próby oraz obszary krytyczne testu dwustronnego i jednostronnego
W teście dla proporcji błąd standardowy obliczamy pod hipotezą zerową, dlatego we wzorze występuje wartość referencyjna p₀, a nie wskaźnik z próby p̂.

Wartość p

Zamiast korzystać z wartości krytycznej można obliczyć wartość p.

Test prawostronny

\[ p= P(Z\geq z_{\text{obs}}) \]

Test lewostronny

\[ p= P(Z\leq z_{\text{obs}}) \]

Test dwustronny

\[ p= 2P(Z\geq|z_{\text{obs}}|) \]

Reguła decyzji jest standardowa:

\[ p\leq\alpha \quad\Longrightarrow\quad \text{odrzucamy }H_0 \]
\[ p>\alpha \quad\Longrightarrow\quad \text{brak podstaw do odrzucenia }H_0 \]

Wartość p nie jest prawdopodobieństwem, że hipoteza zerowa jest prawdziwa. Oznacza ona, jak skrajny byłby wynik z próby przy założeniu prawdziwości \(H_0\).

Test a przedział ufności dla proporcji

W klasycznym podejściu przybliżony dwustronny przedział ufności dla proporcji ma postać:

\[ \left( \hat p - u_{\alpha/2} \sqrt{ \frac{ \hat p(1-\hat p) }{ n } }, \; \hat p + u_{\alpha/2} \sqrt{ \frac{ \hat p(1-\hat p) }{ n } } \right) \]

W odróżnieniu od statystyki testowej w tym wzorze występuje \(\hat p\), czyli proporcja wyznaczona z próby. Dzieje się tak dlatego, że budujemy przedział dla nieznanego parametru na podstawie danych, a nie opisujemy rozkład statystyki przy założeniu prawdziwości konkretnej hipotezy.

Dla dużych prób test dwustronny na poziomie istotności \(\alpha\) i klasyczny dwustronny przedział ufności \(1-\alpha\) zwykle prowadzą do zgodnych wniosków:

Przy małych próbach oraz proporcjach bliskich 0 lub 1 należy zachować ostrożność. Klasyczny przedział normalny i klasyczny test z mogą wtedy mieć słabe własności przybliżone. W praktyce warto stosować metody dokładne albo bardziej odporne przedziały dla proporcji.

Małe próby i dokładny test dwumianowy

Jeżeli warunki przybliżenia normalnego nie są spełnione, nie należy mechanicznie stosować statystyki \(Z\). Dotyczy to przede wszystkim małych prób oraz sytuacji, w których \(p_0\) jest bliskie 0 albo 1.

Wtedy można zastosować dokładny test dwumianowy. Zakładamy, że liczba sukcesów w próbie ma przy prawdziwości \(H_0\) rozkład:

\[ X\sim Bin(n,p_0) \]

W teście dokładnym nie zastępujemy rozkładu dwumianowego rozkładem normalnym. Wartość p obliczamy bezpośrednio z prawdopodobieństw odpowiednich wyników dwumianowych.

W teście prawostronnym sumujemy prawdopodobieństwa wyników co najmniej tak dużych jak zaobserwowany. W teście lewostronnym sumujemy prawdopodobieństwa wyników co najmniej tak małych. W przypadku testu dwustronnego szczegółowa definicja dokładnej wartości p może nieznacznie różnić się między programami statystycznymi, dlatego warto podać zastosowaną metodę obliczeń.

Małe próby w kontroli jakości

Dokładny test dwumianowy ma szczególne znaczenie między innymi w kontroli jakości. Jeżeli z małej partii wybieramy niewielką próbę i liczymy liczbę wadliwych produktów, normalne przybliżenie może być niedokładne, zwłaszcza gdy dopuszczalny odsetek wad jest niski.

W typowych zadaniach akademickich dotyczących dużych prób najczęściej stosuje się jednak klasyczny test z dla jednej proporcji.

Przykład: test dwustronny dla wskaźnika struktury

W badaniu opinii publicznej przeprowadzono ankietę wśród 200 osób. Za określonym rozwiązaniem opowiedziało się 118 respondentów.

Chcemy sprawdzić, czy w populacji wskaźnik poparcia różni się od 50%.

Krok 1. Wskaźnik z próby

\[ n=200 \]
\[ m=118 \]
\[ \hat p= \frac{118}{200} = 0{,}59 \]

Krok 2. Hipotezy

\[ H_0:p=0{,}50 \]
\[ H_1:p\neq0{,}50 \]

Krok 3. Warunki przybliżenia normalnego

\[ np_0=200\cdot0{,}50=100 \]
\[ n(1-p_0)=200\cdot0{,}50=100 \]

Obie wartości są duże, więc przybliżenie normalne jest uzasadnione.

Krok 4. Statystyka testowa

\[ z_{\text{obs}}= \frac{ 0{,}59-0{,}50 }{ \sqrt{ \frac{ 0{,}50(1-0{,}50) }{ 200 } } } \] 
\[ z_{\text{obs}} = \frac{0{,}09}{0{,}03536} \approx 2{,}546 \]

Krok 5. Obszar krytyczny

Przyjmujemy poziom istotności:

\[ \alpha=0{,}05 \]

Dla testu dwustronnego wartość krytyczna wynosi:

\[ u_{0{,}025}=1{,}96 \]

Ponieważ:

\[ |2{,}546|>1{,}96 \]

odrzucamy hipotezę zerową na poziomie istotności 0,05.

Krok 6. Wartość p

\[ p\approx0{,}011 \]

Ponieważ wartość p jest mniejsza od 0,05, metoda wartości p prowadzi do tej samej decyzji.

Krok 7. Kontrola przez przedział ufności

\[ 0{,}59 \pm 1{,}96 \sqrt{ \frac{ 0{,}59(1-0{,}59) }{ 200 } } \] 
\[ (0{,}522;\;0{,}658) \]

Wartość hipotetyczna \(p_0=0{,}50\) nie należy do przedziału ufności. Potwierdza to decyzję o odrzuceniu hipotezy zerowej.

Wniosek: dane z próby dostarczają podstaw do stwierdzenia, że wskaźnik poparcia w populacji różni się od 50%.

Kalkulatory i programy statystyczne

W klasycznych zadaniach wartość krytyczną można odczytać z tablic rozkładu normalnego. W praktyce wygodniejsze jest korzystanie z kalkulatorów i programów statystycznych, które podają statystykę testową, wartość p oraz przedział ufności dla proporcji.

Do kontroli wartości krytycznych przyda się tablica rozkładu normalnego i kalkulator. Test dla jednej proporcji oraz dokładny test dwumianowy oferują również między innymi Excel, SPSS, Statistica, Gretl, R i inne pakiety statystyczne.

Najczęstsze błędy

Podsumowanie

Test dla jednego wskaźnika struktury pozwala sprawdzić, czy udział określonej cechy w populacji jest zgodny z wartością referencyjną.

Ruch po okręgu – prędkość kątowa, przyspieszenie dośrodkowe, okres i częstotliwość

Ruch po okręgu jest szczególnym przypadkiem ruchu krzywoliniowego, w którym punkt pozostaje w stałej odległości od środka toru. Do jego opisu wykorzystuje się zarówno wielkości liniowe, takie jak droga, prędkość i przyspieszenie, jak i wielkości kątowe: położenie kątowe, prędkość kątową, przyspieszenie kątowe, okres oraz częstotliwość.

Nawet w jednostajnym ruchu po okręgu, gdy szybkość punktu pozostaje stała, jego wektor prędkości nieustannie zmienia kierunek. Z tego powodu występuje przyspieszenie normalne, nazywane w tym przypadku przyspieszeniem dośrodkowym.

Ogólny rozkład przyspieszenia na składową styczną i normalną został przedstawiony w artykule Kinematyka punktu materialnego – przyspieszenie styczne i normalne. W ruchu po okręgu promień krzywizny toru jest stały i równy promieniowi okręgu \(R\).

Czym jest ruch po okręgu?

Ruch po okręgu zachodzi wtedy, gdy odległość poruszającego się punktu \(P\) od ustalonego środka \(O\) jest stała:

\[ OP=R=\operatorname{const}. \]

Tor punktu jest wówczas okręgiem o promieniu \(R\). Przykładami ruchu po okręgu lub ruchu zbliżonego do kołowego są:

W kinematycznym opisie ruchu po okręgu nie analizujemy jeszcze sił powodujących ten ruch. Pojęcie siły dośrodkowej należy do dynamiki, natomiast w kinematyce wyznaczamy przede wszystkim prędkość i przyspieszenie punktu.

Geometria ruchu po okręgu

Położenie punktu na okręgu można określić za pomocą kąta \(\varphi\), jaki promień wodzący \(\overrightarrow{OP}\) tworzy z przyjętym kierunkiem odniesienia.

Jeżeli punkt przemieszcza się po łuku okręgu o długości \(s\), zachodzi zależność:

\[ s=R\varphi, \]

pod warunkiem, że kąt \(\varphi\) jest wyrażony w radianach.

Okrąg z zaznaczonym promieniem, kątem środkowym i długością łuku
Długość łuku \(s\) zakreślonego przez punkt poruszający się po okręgu wynosi \(s=R\varphi\), jeżeli kąt \(\varphi\) jest wyrażony w radianach.

Dla małych przyrostów kąta i długości łuku:

\[ ds=R\,d\varphi. \]

Zależność ta stanowi podstawę powiązania prędkości liniowej z prędkością kątową.

Miara łukowa kąta – radiany

Jeden radian jest kątem środkowym, któremu odpowiada łuk o długości równej promieniowi okręgu:

\[ s=R \quad\Longrightarrow\quad \varphi=1\ \mathrm{rad}. \]

Pełnemu obrotowi odpowiada kąt:

\[ 2\pi\ \mathrm{rad}=360^\circ. \]

Stąd:

\[ \pi\ \mathrm{rad}=180^\circ, \qquad \frac{\pi}{2}\ \mathrm{rad}=90^\circ. \]

Zamiana stopni na radiany:

\[ \varphi[\mathrm{rad}] = \varphi[^\circ]\frac{\pi}{180^\circ}. \]

Zamiana radianów na stopnie:

\[ \varphi[^\circ] = \varphi[\mathrm{rad}] \frac{180^\circ}{\pi}. \]

We wzorach różniczkowych i zależności \(s=R\varphi\) kąt musi być podany w radianach. Bezpośrednie podstawienie kąta wyrażonego w stopniach prowadzi do błędnego wyniku.

Położenie i przemieszczenie kątowe

Położenie punktu na okręgu opisuje funkcja:

\[ \varphi=\varphi(t). \]

Jeżeli w chwili \(t_1\) punkt ma położenie kątowe \(\varphi_1\), a w chwili \(t_2\) położenie \(\varphi_2\), jego przemieszczenie kątowe wynosi:

\[ \Delta\varphi = \varphi_2-\varphi_1. \]

Najczęściej za dodatni przyjmuje się zwrot przeciwny do ruchu wskazówek zegara. Wtedy:

Liczba wykonanych obrotów jest związana z przemieszczeniem kątowym zależnością:

\[ n = \frac{\Delta\varphi}{2\pi}. \]

Przemieszczenie kątowe może być większe od \(2\pi\). Na przykład wartość \(6\pi\) oznacza trzy pełne obroty w dodatnim zwrocie.

Prędkość kątowa

Średnia prędkość kątowa

Średnia prędkość kątowa w przedziale czasu wynosi:

\[ \omega_{\mathrm{śr}} = \frac{\Delta\varphi}{\Delta t}. \] 

Chwilowa prędkość kątowa

Chwilowa prędkość kątowa jest pochodną położenia kątowego względem czasu:

\[ \omega = \frac{d\varphi}{dt} = \dot\varphi. \]

Jednostką prędkości kątowej jest:

\[ [\omega] = \mathrm{\frac{rad}{s}}. \]

Znak \(\omega\) informuje o zwrocie obrotu. Wartość bezwzględna \(|\omega|\) określa szybkość zmian położenia kątowego.

Przyspieszenie kątowe

Średnie przyspieszenie kątowe jest ilorazem zmiany prędkości kątowej i czasu:

\[ \varepsilon_{\mathrm{śr}} = \frac{\Delta\omega}{\Delta t}. \] 

Chwilowe przyspieszenie kątowe definiujemy jako:

\[ \varepsilon = \frac{d\omega}{dt} = \frac{d^2\varphi}{dt^2} = \ddot\varphi. \]

Jego jednostką jest:

\[ [\varepsilon] = \mathrm{\frac{rad}{s^2}}. \]

Podobnie jak w ruchu prostoliniowym, sam znak przyspieszenia nie wystarcza do stwierdzenia, czy ruch przyspiesza.

\[ \omega\varepsilon>0 \quad\Longrightarrow\quad |\omega|\ \text{rośnie}, \] \[ \omega\varepsilon<0 \quad\Longrightarrow\quad |\omega|\ \text{maleje}. \]

Związki wielkości liniowych i kątowych

Z zależności:

\[ s=R\varphi \]

po zróżniczkowaniu względem czasu otrzymujemy:

\[ \frac{ds}{dt} = R\frac{d\varphi}{dt}. \]

Ponieważ:

\[ v=\frac{ds}{dt}, \qquad \omega=\frac{d\varphi}{dt}, \]

otrzymujemy podstawowy związek:

\[ v=\omega R. \] 

Po ponownym zróżniczkowaniu wartości prędkości:

\[ a_t = \frac{dv}{dt} = R\frac{d\omega}{dt}, \]

czyli:

\[ a_t=\varepsilon R. \] 
Wielkość kątowaWielkość liniowaZależność
kąt \(\varphi\)długość łuku \(s\)\(s=R\varphi\)
prędkość kątowa \(\omega\)prędkość liniowa \(v\)\(v=\omega R\)
przyspieszenie kątowe \(\varepsilon\)przyspieszenie styczne \(a_t\)\(a_t=\varepsilon R\)

Wektor prędkości

Wektor prędkości punktu poruszającego się po okręgu jest w każdej chwili styczny do toru. Jest więc prostopadły do promienia poprowadzonego ze środka okręgu do aktualnego położenia punktu:

\[ \vec v\perp\overrightarrow{OP}. \]

Jego wartość wynosi:

\[ v=|\vec v|=|\omega|R. \]

W ruchu jednostajnym po okręgu wartość \(v\) jest stała, ale kierunek wektora \(\vec v\) stale się zmienia. Oznacza to, że punkt ma niezerowe przyspieszenie mimo stałej szybkości.

Przyspieszenie dośrodkowe

Składowa normalna przyspieszenia jest skierowana od punktu \(P\) do środka okręgu \(O\). W ruchu po okręgu nazywamy ją przyspieszeniem dośrodkowym.

\[ a_n = \frac{v^2}{R}. \]

Po wykorzystaniu zależności \(v=\omega R\):

\[ a_n = \frac{(\omega R)^2}{R} = \omega^2R. \]

Można również zapisać:

\[ a_n=v\omega. \]

Kierunek przyspieszenia normalnego nie zależy od tego, czy punkt porusza się zgodnie, czy przeciwnie do ruchu wskazówek zegara. W obu przypadkach wektor \(\vec a_n\) jest skierowany do środka okręgu.

Przedłużenie wektora przyspieszenia normalnego zawsze przechodzi przez środek okręgu. Wektor ten nie może być styczny do toru ani skierowany na zewnątrz.

Skąd bierze się wzór na przyspieszenie dośrodkowe?

Dla małego przemieszczenia po okręgu kierunek wektora prędkości zmienia się o niewielki kąt \(\Delta\varphi\). Dla małych kątów zmiana wartości wektora prędkości spełnia przybliżenie:

\[ \Delta v\approx v\Delta\varphi. \]

Po podzieleniu przez \(\Delta t\) i przejściu do granicy:

\[ a_n = \lim_{\Delta t\to 0} \frac{\Delta v}{\Delta t} = v\lim_{\Delta t\to 0} \frac{\Delta\varphi}{\Delta t} = v\omega. \]

Ponieważ \(v=\omega R\), otrzymujemy:

\[ a_n = v\omega = \frac{v^2}{R} = \omega^2R. \]

Przyspieszenie styczne i całkowite

Jeżeli wartość prędkości zmienia się w czasie, oprócz przyspieszenia normalnego występuje także składowa styczna:

\[ a_t = \frac{dv}{dt} = \varepsilon R. \] 

Przyspieszenie styczne jest styczne do okręgu:

Przyspieszenie całkowite jest sumą wektorową obu prostopadłych składowych:

\[ \vec a = \vec a_t+\vec a_n. \] 

Jego wartość wynosi:

\[ a = \sqrt{a_t^2+a_n^2}. \]

Po podstawieniu zależności kątowych:

\[ a = \sqrt{(\varepsilon R)^2+(\omega^2R)^2}, \] \[ a = R\sqrt{\varepsilon^2+\omega^4}. \]
Punkt na okręgu z prędkością styczną, przyspieszeniem dośrodkowym, przyspieszeniem stycznym i przyspieszeniem wypadkowym
Prędkość \(\vec v\) jest styczna do okręgu, przyspieszenie normalne \(\vec a_n\) jest skierowane do jego środka, a przyspieszenie styczne \(\vec a_t\) odpowiada za zmianę wartości prędkości.

W ruchu niejednostajnym całkowite przyspieszenie zazwyczaj nie jest skierowane dokładnie do środka okręgu. Do środka skierowana jest wyłącznie jego składowa normalna.

Ruch jednostajny po okręgu

Ruch jednostajny po okręgu zachodzi wtedy, gdy prędkość kątowa jest stała:

\[ \omega=\operatorname{const}. \]

Wówczas:

\[ \varepsilon=0, \qquad a_t=0. \]

Położenie kątowe zmienia się liniowo z czasem:

\[ \varphi(t) = \varphi_0+\omega(t-t_0). \]

Dla \(t_0=0\):

\[ \varphi = \varphi_0+\omega t. \] 

Wartość prędkości liniowej jest stała:

\[ v=|\omega|R=\operatorname{const}. \]

Przyspieszenie jest jednak niezerowe:

\[ a=a_n=\frac{v^2}{R}=\omega^2R. \]

Ruch jednostajny po okręgu nie jest ruchem bez przyspieszenia. Stała jest tylko wartość prędkości. Jej kierunek zmienia się w każdej chwili, dlatego występuje przyspieszenie dośrodkowe.

Okres i częstotliwość

Okres \(T\) jest czasem wykonania jednego pełnego obrotu.

Częstotliwość \(f\) określa liczbę pełnych obrotów wykonywanych w jednostce czasu.

\[ f=\frac{1}{T}, \qquad T=\frac{1}{f}. \]

Podczas jednego pełnego obrotu:

\[ \Delta\varphi=2\pi. \]

W jednostajnym ruchu po okręgu:

\[ \omega = \frac{2\pi}{T} = 2\pi f. \] 

Prędkość liniowa wynosi zatem:

\[ v = \omega R = \frac{2\pi R}{T} = 2\pi Rf. \] 

Przyspieszenie dośrodkowe można zapisać jako:

\[ a_n = \frac{4\pi^2R}{T^2}, \]

albo:

\[ a_n = 4\pi^2Rf^2. \]
Schemat pełnego obrotu punktu po okręgu z wzorami na okres, częstotliwość, prędkość kątową i liczbę obrotów
Okres \(T\) jest czasem jednego pełnego obrotu, częstotliwość wynosi \(f=\frac{1}{T}\), a prędkość kątowa w ruchu jednostajnym \(\omega=\frac{2\pi}{T}=2\pi f\).

Liczba obrotów wykonanych w czasie \(t\) wynosi:

\[ n = \frac{t}{T} = ft. \] 

Droga przebyta podczas \(n\) pełnych obrotów jest równa:

\[ s=2\pi Rn. \] 

Obroty na minutę

Prędkość obrotowa urządzeń technicznych jest często podawana w obrotach na minutę, oznaczanych jako \(\mathrm{obr/min}\) albo rpm.

Jeżeli urządzenie wykonuje \(N\) obrotów na minutę, częstotliwość w hercach wynosi:

\[ f=\frac{N}{60}. \]

Prędkość kątowa:

\[ \omega = 2\pi f = \frac{2\pi N}{60}. \]

W drugą stronę:

\[ N=60f. \]

Ruch niejednostajny po okręgu

W ruchu niejednostajnym prędkość kątowa zależy od czasu:

\[ \omega=\omega(t). \]

Występuje wtedy niezerowe przyspieszenie kątowe:

\[ \varepsilon = \frac{d\omega}{dt}. \]

Przyspieszenie punktu ma dwie składowe:

\[ a_t=\varepsilon R, \qquad a_n=\omega^2R. \]

Składowa styczna zmienia wartość prędkości, natomiast składowa normalna zmienia jej kierunek. Ponieważ \(a_n\) zależy od kwadratu prędkości kątowej, może zmieniać się nawet przy stałym przyspieszeniu kątowym.

Ruch jednostajnie zmienny kątowo

Ruch jest jednostajnie zmienny kątowo, gdy przyspieszenie kątowe jest stałe:

\[ \varepsilon=\operatorname{const}. \]

Równanie prędkości kątowej:

\[ \omega(t) = \omega_0+\varepsilon(t-t_0). \]

Dla \(t_0=0\):

\[ \omega = \omega_0+\varepsilon t. \] 

Równanie położenia kątowego:

\[ \varphi(t) = \varphi_0 + \omega_0(t-t_0) + \frac{\varepsilon(t-t_0)^2}{2}. \]

Dla \(t_0=0\):

\[ \varphi = \varphi_0+\omega_0t+\frac{\varepsilon t^2}{2}. \]

Po wyeliminowaniu czasu otrzymujemy odpowiednik równania znanego z ruchu prostoliniowego:

\[ \omega^2 = \omega_0^2 + 2\varepsilon(\varphi-\varphi_0). \]
Porównanie ruchu jednostajnego po okręgu z ruchem o stałym przyspieszeniu kątowym
W ruchu jednostajnym po okręgu \(\omega=\operatorname{const}\) i \(a_t=0\), natomiast przy stałym przyspieszeniu kątowym \(\varepsilon=\operatorname{const}\) prędkość kątowa zmienia się liniowo z czasem.

Warto podkreślić, że nawet przy stałym \(\varepsilon\) całkowite przyspieszenie punktu nie jest zazwyczaj stałe. Składowa styczna ma stałą wartość:

\[ a_t=\varepsilon R=\operatorname{const}, \]

ale składowa normalna zmienia się zgodnie z:

\[ a_n=\omega^2R. \]

Opis wektorowy w układzie kartezjańskim

Jeżeli środek okręgu znajduje się w początku układu współrzędnych, położenie punktu można zapisać jako:

\[ \vec r = R\cos\varphi\,\vec i + R\sin\varphi\,\vec j. \] 

Po zróżniczkowaniu względem czasu otrzymujemy wektor prędkości:

\[ \vec v = -R\omega\sin\varphi\,\vec i + R\omega\cos\varphi\,\vec j. \] 

Jego wartość:

\[ v = \sqrt{R^2\omega^2\sin^2\varphi+ R^2\omega^2\cos^2\varphi} = |\omega|R. \]

Wektor przyspieszenia wynosi:

\[ \vec a = -R\omega^2 \left( \cos\varphi\,\vec i+\sin\varphi\,\vec j \right) + R\varepsilon \left( -\sin\varphi\,\vec i+\cos\varphi\,\vec j \right). \]

Pierwszy składnik jest skierowany do środka okręgu i odpowiada przyspieszeniu normalnemu, natomiast drugi jest styczny do toru.

Przy użyciu wersorów radialnego \(\vec e_r\) i stycznego \(\vec e_t\) zapis jest krótszy:

\[ \vec v = R\omega\vec e_t, \] \[ \vec a = -R\omega^2\vec e_r + R\varepsilon\vec e_t. \] 

Minus przy składowej radialnej oznacza, że przyspieszenie normalne jest skierowane przeciwnie do wersora \(\vec e_r\), czyli do środka okręgu.

Wykresy wielkości kątowych

Wykres położenia kątowego \(\varphi(t)\)

Nachylenie wykresu położenia kątowego jest równe prędkości kątowej:

\[ \omega=\frac{d\varphi}{dt}. \]

Wykres prędkości kątowej \(\omega(t)\)

Nachylenie wykresu prędkości kątowej jest równe przyspieszeniu kątowemu:

\[ \varepsilon=\frac{d\omega}{dt}. \]

Pole algebraiczne pod wykresem \(\omega(t)\) jest równe przemieszczeniu kątowemu:

\[ \Delta\varphi = \int_{t_1}^{t_2}\omega(t)\,dt. \]

Wykres przyspieszenia kątowego \(\varepsilon(t)\)

Pole algebraiczne pod wykresem przyspieszenia kątowego jest równe zmianie prędkości kątowej:

\[ \Delta\omega = \int_{t_1}^{t_2}\varepsilon(t)\,dt. \]

Przykłady obliczeniowe

Przykład 1. Samochód na okrągłym torze

Samochód porusza się ze stałą szybkością:

\[ v=10\ \mathrm{m/s} \]

po okręgu o promieniu:

\[ R=50\ \mathrm{m}. \]

Prędkość kątowa:

\[ \omega = \frac{v}{R} = \frac{10}{50} = 0{,}2\ \mathrm{rad/s}. \]

Okres jednego okrążenia:

\[ T = \frac{2\pi}{\omega} = \frac{2\pi}{0{,}2} = 10\pi\ \mathrm{s} \approx 31{,}4\ \mathrm{s}. \]

Przyspieszenie dośrodkowe:

\[ a_n = \frac{v^2}{R} = \frac{10^2}{50} = 2\ \mathrm{m/s^2}. \]

Przykład 2. Wirnik wykonujący 1200 obrotów na minutę

Wirnik obraca się z prędkością:

\[ N=1200\ \mathrm{obr/min}. \]

Częstotliwość:

\[ f = \frac{1200}{60} = 20\ \mathrm{Hz}. \]

Okres:

\[ T = \frac{1}{f} = \frac{1}{20} = 0{,}05\ \mathrm{s}. \]

Prędkość kątowa:

\[ \omega = 2\pi f = 40\pi\ \mathrm{rad/s} \approx 125{,}7\ \mathrm{rad/s}. \]

Jeżeli rozpatrywany punkt znajduje się w odległości \(R=0{,}08\ \mathrm{m}\) od osi, jego prędkość liniowa wynosi:

\[ v = \omega R = 40\pi\cdot 0{,}08 = 3{,}2\pi\ \mathrm{m/s} \approx 10{,}1\ \mathrm{m/s}. \]

Przykład 3. Ruch ze stałym przyspieszeniem kątowym

Koło o promieniu:

\[ R=0{,}5\ \mathrm{m} \]

ma początkową prędkość kątową:

\[ \omega_0=2\ \mathrm{rad/s} \]

i przyspieszenie kątowe:

\[ \varepsilon=0{,}5\ \mathrm{rad/s^2}. \]

Po \(6\ \mathrm{s}\) prędkość kątowa wynosi:

\[ \omega = \omega_0+\varepsilon t = 2+0{,}5\cdot 6 = 5\ \mathrm{rad/s}. \]

Przemieszczenie kątowe:

\[ \Delta\varphi = \omega_0t+\frac{\varepsilon t^2}{2}, \] \[ \Delta\varphi = 2\cdot 6+ \frac{0{,}5\cdot 6^2}{2} = 21\ \mathrm{rad}. \]

Przyspieszenie styczne:

\[ a_t = \varepsilon R = 0{,}5\cdot 0{,}5 = 0{,}25\ \mathrm{m/s^2}. \]

Przyspieszenie normalne:

\[ a_n = \omega^2R = 5^2\cdot 0{,}5 = 12{,}5\ \mathrm{m/s^2}. \]

Przyspieszenie całkowite:

\[ a = \sqrt{0{,}25^2+12{,}5^2} \approx 12{,}50\ \mathrm{m/s^2}. \]

Przykład 4. Hamowanie ruchu obrotowego

Tarcza ma początkową prędkość kątową:

\[ \omega_0=6\ \mathrm{rad/s} \]

i stałe przyspieszenie kątowe:

\[ \varepsilon=-2\ \mathrm{rad/s^2}. \]

Chwila zatrzymania spełnia:

\[ 0=\omega_0+\varepsilon t_z, \] \[ t_z = -\frac{\omega_0}{\varepsilon} = -\frac{6}{-2} = 3\ \mathrm{s}. \]

Do chwili zatrzymania tarcza obróci się o kąt:

\[ \Delta\varphi = \omega_0t_z+\frac{\varepsilon t_z^2}{2}, \] \[ \Delta\varphi = 6\cdot 3+ \frac{-2\cdot 3^2}{2} = 9\ \mathrm{rad}. \]

Liczba obrotów:

\[ n = \frac{9}{2\pi} \approx 1{,}43. \]

Jeżeli ujemne przyspieszenie kątowe działałoby nadal, po zatrzymaniu tarcza zaczęłaby obracać się w przeciwnym kierunku.

Jednostki i kontrola wymiarów

WielkośćOznaczenieJednostka SI
Promień\(R\)\(\mathrm{m}\)
Kąt\(\varphi\)\(\mathrm{rad}\)
Długość łuku\(s\)\(\mathrm{m}\)
Prędkość kątowa\(\omega\)\(\mathrm{rad/s}\)
Przyspieszenie kątowe\(\varepsilon\)\(\mathrm{rad/s^2}\)
Prędkość liniowa\(v\)\(\mathrm{m/s}\)
Przyspieszenie\(a\)\(\mathrm{m/s^2}\)
Okres\(T\)\(\mathrm{s}\)
Częstotliwość\(f\)\(\mathrm{Hz}=\mathrm{s^{-1}}\)

Kontrola wymiarów prędkości liniowej:

\[ [v] = [\omega][R] = \frac{1}{\mathrm{s}}\cdot\mathrm{m} = \frac{\mathrm{m}}{\mathrm{s}}. \]

Radian jest formalnie jednostką bezwymiarową, dlatego nie wpływa na końcowy wymiar fizyczny.

Kontrola wymiarów przyspieszenia normalnego:

\[ [a_n] = \frac{[v]^2}{[R]} = \frac{\mathrm{m^2/s^2}}{\mathrm{m}} = \frac{\mathrm{m}}{\mathrm{s^2}}. \]

Kontrola wymiarów przyspieszenia stycznego:

\[ [a_t] = [\varepsilon][R] = \frac{1}{\mathrm{s^2}}\cdot\mathrm{m} = \frac{\mathrm{m}}{\mathrm{s^2}}. \]

Najczęstsze błędy

W chwili, w której \(\omega=0\), przyspieszenie normalne jest równe zeru:

\[ a_n=\omega^2R=0. \]

Przyspieszenie styczne może jednak nadal być niezerowe, jeżeli \(\varepsilon\neq 0\). Punkt może więc chwilowo się zatrzymać, a następnie rozpocząć ruch w przeciwnym kierunku.

Podsumowanie

Położenie punktu na okręgu opisuje kąt:

\[ \varphi=\varphi(t). \]

Długość łuku i kąt są związane zależnością:

\[ s=R\varphi. \]

Prędkość i przyspieszenie kątowe wynoszą:

\[ \omega = \frac{d\varphi}{dt}, \qquad \varepsilon = \frac{d\omega}{dt}. \]

Wielkości liniowe i kątowe łączą zależności:

\[ v=\omega R, \qquad a_t=\varepsilon R. \] 

Przyspieszenie normalne jest skierowane do środka okręgu:

\[ a_n = \frac{v^2}{R} = \omega^2R. \]

W ruchu niejednostajnym całkowite przyspieszenie wynosi:

\[ a = \sqrt{a_t^2+a_n^2}. \]

W jednostajnym ruchu po okręgu:

\[ \omega=\operatorname{const}, \qquad \varepsilon=0, \qquad a_t=0, \] \[ a=a_n=\omega^2R. \]

Okres, częstotliwość i prędkość kątową łączą wzory:

\[ f=\frac{1}{T}, \qquad \omega=\frac{2\pi}{T}=2\pi f. \] 

Ruch po okręgu pokazuje, że przyspieszenie nie musi oznaczać wyłącznie zmiany wartości prędkości. Może wynikać również ze zmiany jej kierunku, nawet gdy szybkość pozostaje stała.

Testowanie równości dwóch wariancji — test F-Snedecora i homogeniczność wariancji

Testowanie równości dwóch wariancji pozwala sprawdzić, czy dwie populacje różnią się zmiennością badanej cechy. W artykule omawiamy klasyczny test F-Snedecora dla dwóch wariancji: próby jedno- i dwustronne, wybór licznika i mianownika, obszary krytyczne, wartość p, przedział ufności dla ilorazu wariancji oraz znaczenie tego testu przy porównywaniu średnich.

Porównując dwie populacje, często skupiamy się przede wszystkim na ich średnich. Równie ważne może być jednak pytanie o zmienność. Dwa procesy technologiczne mogą dawać produkty o takiej samej średniej masie, ale jeden z nich może być wyraźnie mniej stabilny. Dwie grupy uczące się różnymi metodami mogą osiągać podobne przeciętne wyniki, lecz w jednej grupie rezultaty mogą być znacznie bardziej zróżnicowane.

Do badania tego problemu służy test równości dwóch wariancji, oparty na rozkładzie F-Snedecora. Pozwala on ocenić, czy różnica między wariancjami z prób jest na tyle duża, że stanowi podstawę do odrzucenia hipotezy o równej zmienności w populacjach.

Wcześniej omówiliśmy już testowanie hipotez dotyczących jednej wariancji oraz testowanie różnicy średnich w dwóch populacjach. Test F stanowi naturalne rozwinięcie obu tych zagadnień.

Co testujemy?

Niech \(\sigma_1^2\) oraz \(\sigma_2^2\) oznaczają wariancje w dwóch populacjach. W najczęstszym przypadku chcemy sprawdzić, czy są one równe:

\[ \sigma_1^2=\sigma_2^2 \]

Równoważnie możemy badać iloraz wariancji:

\[ \theta= \frac{\sigma_1^2}{\sigma_2^2} \]

Hipoteza o równości wariancji jest wtedy hipotezą, że:

\[ \theta=1 \]

Test może służyć zarówno do bezpośredniego badania zmienności dwóch procesów, jak i do oceny założenia pomocniczego potrzebnego w niektórych klasycznych modelach statystycznych.

Wariancja a odchylenie standardowe

Test F konstruuje się formalnie dla wariancji. Ponieważ jednak odchylenie standardowe jest dodatnim pierwiastkiem z wariancji, hipoteza o równości wariancji jest równoważna hipotezie o równości odchyleń standardowych:

\[ \sigma_1^2=\sigma_2^2 \quad\Longleftrightarrow\quad \sigma_1=\sigma_2 \]

W obliczeniach stosujemy wariancje z prób, ponieważ ich iloraz ma rozkład F-Snedecora. Wnioski można jednak opisywać także prostszym językiem: „zmienność w pierwszej populacji jest większa”, „odchylenia standardowe są zgodne” albo „nie ma podstaw do stwierdzenia różnicy rozproszenia”.

Założenia klasycznego testu F

Klasyczny test równości dwóch wariancji wymaga spełnienia kilku istotnych założeń:

Normalność obu populacji jest szczególnie ważna. Rozkład F ilorazu wariancji z prób wynika dokładnie z założenia normalności. Nawet umiarkowane odstępstwa od normalności mogą istotnie zmienić rzeczywisty poziom istotności testu.

Test F jest czuły na brak normalności

Klasyczny test F warto traktować jako właściwy przede wszystkim wtedy, gdy normalność populacji ma uzasadnienie merytoryczne albo została wiarygodnie potwierdzona diagnostycznie. Przy danych silnie asymetrycznych lub zawierających obserwacje odstające bezpieczniejsze mogą być metody bardziej odporne, takie jak test Levene’a lub test Browna – Forsythe’a.

Statystyka testowa i rozkład F-Snedecora

Niech \(S_1^2\) oraz \(S_2^2\) oznaczają skorygowane wariancje z dwóch niezależnych prób:

\[ S_1^2= \frac{1}{n_1-1} \sum_{i=1}^{n_1} (x_{1i}-\overline{x}_1)^2 \]
\[ S_2^2= \frac{1}{n_2-1} \sum_{i=1}^{n_2} (x_{2i}-\overline{x}_2)^2 \]

Statystyka testowa ma postać ilorazu:

\[ F= \frac{S_1^2}{S_2^2} \]

Jeżeli obie populacje są normalne i hipoteza zerowa o równości wariancji jest prawdziwa, statystyka \(F\) ma rozkład F-Snedecora z liczbą stopni swobody:

\[ \nu_1=n_1-1 \]
\[ \nu_2=n_2-1 \]
\[ F\sim F_{\nu_1,\nu_2} \]

Pierwsza liczba stopni swobody należy do wariancji umieszczonej w liczniku, a druga — do wariancji umieszczonej w mianowniku. Kolejność ma więc znaczenie.

Infografika przedstawiająca dwie niezależne próby, wariancje z prób, iloraz F oraz stopnie swobody w teście F-Snedecora
Statystyka F jest ilorazem wariancji z dwóch niezależnych prób.

Hipotezy jedno- i dwustronne

Podobnie jak w innych testach, rodzaj hipotezy alternatywnej powinien wynikać z celu analizy i zostać wybrany przed obejrzeniem wyników.

Test dwustronny

\[ H_0:\sigma_1^2=\sigma_2^2 \]
\[ H_1:\sigma_1^2\neq\sigma_2^2 \]

Ten wariant stosujemy wtedy, gdy interesuje nas każda różnica zmienności: zarówno sytuacja, w której pierwsza populacja jest bardziej zróżnicowana, jak i sytuacja odwrotna.

Test prawostronny

\[ H_0:\sigma_1^2=\sigma_2^2 \]
\[ H_1:\sigma_1^2>\sigma_2^2 \]

Test prawostronny ma zastosowanie na przykład wtedy, gdy chcemy sprawdzić, czy nowy proces produkcyjny zwiększył zmienność wyników względem procesu dotychczasowego.

Test lewostronny

\[ H_0:\sigma_1^2=\sigma_2^2 \]
\[ H_1:\sigma_1^2\sigma_2^2 \]

Test lewostronny wybieramy wtedy, gdy interesuje nas wyłącznie zmniejszenie zmienności w pierwszej populacji względem drugiej.

Którą wariancję umieścić w liczniku?

W testach jednostronnych kolejność wariancji w ilorazie musi odpowiadać kierunkowi hipotezy alternatywnej.

Jeżeli testujemy:

\[ H_1:\sigma_1^2>\sigma_2^2 \]

naturalnie stosujemy statystykę:

\[ F= \frac{S_1^2}{S_2^2} \]

i szukamy dużych wartości w prawym ogonie rozkładu F.

Jeżeli natomiast testujemy:

\[ H_1:\sigma_1^2\sigma_2^2 \]

możemy pozostawić iloraz \(S_1^2/S_2^2\) i porównywać go z dolnym kwantylem rozkładu F albo odwrócić kolejność wariancji:

\[ F^\ast= \frac{S_2^2}{S_1^2} \]

oraz przeprowadzić równoważny test prawostronny. Po odwróceniu ilorazu trzeba jednocześnie zamienić kolejność stopni swobody.

Nie zamieniaj licznika z mianownikiem po to, aby „wyszło większe F”

W teście jednostronnym samowolna zamiana wariancji w ilorazie zmienia kierunek hipotezy alternatywnej. Można ją wykonać wyłącznie wtedy, gdy świadomie przepisujemy test na równoważny wariant dotyczący odwrotnego ilorazu i równocześnie zamieniamy stopnie swobody.

W teście dwustronnym często stosuje się wygodną konwencję:

\[ F^\ast= \frac{ \max(S_1^2,S_2^2) }{ \min(S_1^2,S_2^2) } \]

Wtedy zawsze:

\[ F^\ast\geq1 \]

i wystarczy porównywać statystykę z prawostronną wartością krytyczną. Liczba stopni swobody w liczniku należy wtedy do grupy o większej wariancji z próby, a w mianowniku — do grupy o mniejszej wariancji z próby.

Obszary krytyczne i wartości krytyczne

W dalszej części przyjmujemy oznaczenie lewostronnych kwantyli rozkładu F:

\[ P\left( F_{\nu_1,\nu_2} \leq F_{p;\nu_1,\nu_2} \right) = p \] 

Wartość \(F_{p;\nu_1,\nu_2}\) oznacza więc kwantyl rzędu \(p\). Trzeba pamiętać, że różne tablice i programy mogą stosować inną konwencję oznaczeń, zwłaszcza dla wartości „z prawego ogona”.

Test prawostronny

Dla hipotezy:

\[ H_1:\sigma_1^2>\sigma_2^2 \]

odrzucamy \(H_0\), gdy:

\[ f_{\text{obs}} > F_{1-\alpha;\nu_1,\nu_2} \]

Test lewostronny

Dla hipotezy:

\[ H_1:\sigma_1^2\sigma_2^2 \]

odrzucamy \(H_0\), gdy:

\[ f_{\text{obs}} F_{\alpha;\nu_1,\nu_2} \]

Test dwustronny — zapis bez zmiany kolejności wariancji

Dla hipotezy:

\[ H_1:\sigma_1^2\neq\sigma_2^2 \]

odrzucamy \(H_0\), gdy iloraz wariancji znajduje się w jednym z dwóch ogonów rozkładu:

\[ f_{\text{obs}} F_{\alpha/2;\nu_1,\nu_2} \]

albo gdy:

\[ f_{\text{obs}} > F_{1-\alpha/2;\nu_1,\nu_2} \]

Test dwustronny — większa wariancja w liczniku

W wielu zadaniach najpierw umieszcza się większą wariancję z próby w liczniku. Wtedy wystarczy badać prawy ogon rozkładu F:

\[ F^\ast= \frac{ S_{\max}^2 }{ S_{\min}^2 } \]

Odrzucamy hipotezę o równości wariancji na poziomie istotności \(\alpha\), gdy:

\[ F^\ast > F_{1-\alpha/2;\nu_{\max},\nu_{\min}} \]

Ten zapis jest równoważny badaniu obu ogonów. Zawiera on jednak ważny warunek: stopnie swobody muszą odpowiadać dokładnie tej grupie, której wariancja znalazła się w liczniku i mianowniku.

Trzy wykresy rozkładu F-Snedecora przedstawiające obszary krytyczne dla testu prawostronnego, lewostronnego i dwustronnego
Położenie obszaru krytycznego zależy od kierunku hipotezy alternatywnej oraz od kolejności wariancji w ilorazie F

Wartość p

Alternatywą dla porównywania z wartością krytyczną jest obliczenie wartości p. Wartość p określa, jak skrajny jest obserwowany iloraz wariancji przy założeniu prawdziwości hipotezy zerowej.

Test prawostronny

\[ p= P\left( F_{\nu_1,\nu_2} \geq f_{\text{obs}} \right) \]

Test lewostronny

\[ p= P\left( F_{\nu_1,\nu_2} \leq f_{\text{obs}} \right) \]

Test dwustronny

W klasycznym dwustronnym teście F wartość p można wyznaczyć jako podwojone mniejsze z prawdopodobieństw ogonowych:

\[ p= 2 \min \left\{ P(F_{\nu_1,\nu_2}\leq f_{\text{obs}}), \; P(F_{\nu_1,\nu_2}\geq f_{\text{obs}}) \right\} \]

Jeżeli wcześniej umieścimy większą wariancję w liczniku, można wyznaczyć prawostronne prawdopodobieństwo dla otrzymanego \(F^\ast\), a następnie podwoić je dla testu dwustronnego.

Reguła decyzji jest standardowa:

\[ p\leq\alpha \quad\Longrightarrow\quad \text{odrzucamy }H_0 \]
\[ p>\alpha \quad\Longrightarrow\quad \text{brak podstaw do odrzucenia }H_0 \]

Przedział ufności dla ilorazu wariancji

Testowanie równości wariancji można połączyć z estymacją przedziałową ilorazu:

\[ \frac{\sigma_1^2}{\sigma_2^2} \]

Dwustronny przedział ufności \(1-\alpha\) dla ilorazu wariancji ma postać:

\[ \left( \frac{ f_{\text{obs}} }{ F_{1-\alpha/2;\nu_1,\nu_2} }, \; \frac{ f_{\text{obs}} }{ F_{\alpha/2;\nu_1,\nu_2} } \right) \]

W teście dwustronnym na poziomie istotności \(\alpha\) obowiązuje zasada:

Tak jak w innych testach, przedział ufności daje dodatkową informację o kierunku i skali możliwej różnicy. Test odpowiada przede wszystkim na pytanie decyzyjne, a przedział pokazuje, jakie wartości ilorazu wariancji są zgodne z danymi.

Związek z testowaniem różnicy średnich

Test równości wariancji jest często omawiany jako test pomocniczy przed zastosowaniem klasycznego testu t-Studenta dla dwóch niezależnych średnich.

W modelu ze wspólną wariancją zakładamy:

\[ \sigma_1^2=\sigma_2^2 \]

i wykorzystujemy połączoną wariancję z próby. Szczegóły tej procedury znajdują się w części nieznane, ale równe wariancje.

Jeżeli założenie równości wariancji nie jest wiarygodne, stosujemy test Welcha, który nie wymaga homogeniczności wariancji. Zobacz część nieznane i nierówne wariancje — test Welcha.

Test F nie powinien być jedyną „bramką” do wyboru testu średnich

Brak podstaw do odrzucenia równości wariancji nie dowodzi, że wariancje są równe. Z kolei odrzucenie hipotezy może wynikać z dużej liczebności próby albo wrażliwości testu F na brak normalności.

W praktycznych analizach często stosuje się test Welcha jako domyślną, bardziej odporną procedurę porównywania średnich. Klasyczny test t ze wspólną wariancją jest uzasadniony wtedy, gdy równość wariancji wynika z projektu badania lub mocnych przesłanek merytorycznych.

Schemat pokazujący wpływ wyniku testu F-Snedecora na wybór klasycznego testu t-Studenta ze wspólną wariancją lub testu Welcha dla dwóch średnich
Równość wariancji ma znaczenie przy klasycznym teście t ze wspólną wariancją, natomiast test Welcha nie wymaga tego założenia.

Przykład: test dwustronny równości wariancji

Porównujemy zmienność wyników uzyskanych w dwóch niezależnych grupach. Załóżmy, że obie populacje można uznać za normalne. Otrzymano:

\[ n_1=16 \]
\[ S_1^2=36 \]
\[ n_2=12 \]
\[ S_2^2=16 \]

Chcemy sprawdzić, czy wariancje w populacjach są równe. Formułujemy test dwustronny:

\[ H_0:\sigma_1^2=\sigma_2^2 \]
\[ H_1:\sigma_1^2\neq\sigma_2^2 \]

Przyjmujemy poziom istotności:

\[ \alpha=0{,}05 \]

Większą wariancję umieszczamy w liczniku:

\[ f_{\text{obs}} = \frac{36}{16} = 2{,}25 \]

Liczba stopni swobody wynosi:

\[ \nu_1=16-1=15 \]
\[ \nu_2=12-1=11 \]

Dla testu dwustronnego stosujemy prawostronną wartość krytyczną odpowiadającą \(\alpha/2=0{,}025\):

\[ F_{0{,}975;15,11} \approx 3{,}12 \]

Ponieważ:

\[ 2{,}253{,}12 \]

nie ma podstaw do odrzucenia hipotezy o równości wariancji na poziomie istotności 0,05.

Wniosek należy sformułować ostrożnie: dane z prób nie dostarczają wystarczających podstaw do stwierdzenia, że wariancje w obu populacjach są różne. Nie oznacza to, że udowodniliśmy ich dokładną równość.

Przy wykorzystaniu kalkulatora można dodatkowo wyznaczyć wartość p oraz 95-procentowy przedział ufności dla ilorazu \(\sigma_1^2/\sigma_2^2\). Jeżeli przedział obejmuje wartość 1, otrzymujemy ten sam wniosek co w teście dwustronnym.

Ograniczenia testu F i alternatywy

Klasyczny test F-Snedecora jest bardzo użyteczny w zadaniach opartych na założeniu normalności populacji. Nie należy jednak traktować go jako uniwersalnego testu zmienności dla dowolnych danych.

Przy wyraźnym braku normalności, asymetrii rozkładu lub obecności obserwacji odstających warto rozważyć bardziej odporne procedury. Do najczęściej stosowanych należą test Levene’a oraz jego modyfikacja Browna – Forsythe’a.

Gdy porównujemy więcej niż dwie grupy, problem homogeniczności wariancji można badać między innymi testem Levene’a albo testem Bartletta. Test Bartletta, podobnie jak test F, jest silnie związany z założeniem normalności.

W dalszych materiałach dotyczących testów nieparametrycznych i diagnostyki warto będzie szerzej omówić test Levene’a, test Browna – Forsythe’a oraz ich zastosowanie przed analizą wariancji.

Kalkulatory i programy statystyczne

Wartości krytyczne rozkładu F-Snedecora można odczytywać z tablic, ale wygodniejsze jest korzystanie z kalkulatorów, ponieważ rozkład F zależy od dwóch liczb stopni swobody i łatwo pomylić ich kolejność.

Do obliczeń przyda się tablica rozkładu F-Snedecora i kalkulator. Wartości p, przedziały ufności i testy równości wariancji oferują również między innymi Excel, SPSS, Statistica, Gretl oraz R.

Najczęstsze błędy

Podsumowanie

Test F-Snedecora pozwala sprawdzić, czy dwie normalne populacje mają taką samą wariancję. Jest podstawowym klasycznym narzędziem do badania homogeniczności wariancji w dwóch populacjach.