Elementy matematyki
ELEMENTY MATEMATYKI
I. OPERACJE PODSTAWOWE
Wartość bezwzględna (moduł) liczby
ae" 0
1. | a | = a, gdy ,
2. | a | = - a, gdy a < 0.
3. | a + b | d" | a | + | b |,
4. | a  b | e" | a | - | b |,
5. | a · b | = | a | · | b |.
Potęgowanie i pierwiastkowanie
1. a n = a · a ·.....· a (n czynników), n  liczba naturalna
2. a 0 = 1,
1
a-n = a`" 0
3. ,
an
1
4. m ,
a = am
n
5. m ,
am = an
6. an Å"am = an+m
,
7. an :am = an-m ,
(an)m = anÅ"m
8. ,
(aÅ" b)n = an Å" bn
9. ,
Logarytmowanie, funkcja wykładnicza
Funkcja wykładnicza jest zdefiniowana wzorem:
f ( x ) = a x
gdzie a jest dowolną liczbą rzeczywistą dodatnią. Jeżeli a x = N to funkcję logatytm o
podstawie a (a > 0 i a `" 1) z liczby N definiuje siÄ™ wzorem:
log a N = x
Funkcja logarytm jest zdefiniowana wyłącznie dla liczb dodatnich.
1. log a 1 = 0,
Zakład Biofizyki CMUJ 1
Elementy matematyki
2. log a a = 1,
3. log a (M · N) = log a M + log a N,
4. log a (N r) = r · log a N,
logb N
loga N =
5. .
logb a
6. Jeżeli a H" 2,7182818... (ta liczba niewymierna jest ściśle zdefiniowana i oznaczana
symbolem e) to wprowadza siÄ™ oznaczenie log e N = lnN i funkcjÄ™ takÄ… nazywamy
logarytmem naturalnym.
7. Jeżeli a=10 to wprowadza się oznaczenie log 10 N = logN i funkcję taką nazywamy
logarytmem dziesiętnym.
loga N -", gdy N 0,
8.
loga N ", gdy N "
9. ,
Funkcje trygonometryczne
Tabela 1 Miary kątów w stopniach i w
radianach
c
Miara kÄ…ta Ä… w Miara kÄ…ta Ä… w
a
stopniach radianach
Ä…
.
00 0
b
Ä„
300
6
Rys. 1 Schemat trójkąta prostokątnego
Ä„
450
przyjęty do zdefiniowana
4
funkcji trygonometrycznych.
Ä„
900
2
1800 Ä„
3600 2Ä„
Definicje funkcji trygonometrycznych na podstawie schematu trójkąta prostokątnego:
sin Ä… = a/c cos Ä… = b/c tg Ä… = a/b ctg Ä… = b/a
Zakład Biofizyki CMUJ 2
Elementy matematyki
sin( + 2Ä„) =sin(
x x)
1. ,
cos( + 2Ä„) = cos(
x x)
2. ,
tg( + 2Ä„) = tg(
x x)
3. ,
2
4. sin2 x + cos x = 1
,
sin(
x)
tg( =
x)
5. ,
cos(
x)
sin( Ä… y) = sin( y) Ä…sin( x)
x x)cos( y)cos(
6. ,
cos( Ä… y) = cos( y) mðsin( y)
x x)cos( x)sin(
7. ,
x + y x - y
sin( + sin( = 2sin cos
x) y)
8. ,
2 2
x + y x - y
cos( + cos( = 2cos cos
x) y)
9. ,
2 2
10. sin(-x)=-sin(x),
11. cos(-x)=cos(x),
Ä„
ëÅ‚x öÅ‚
sin Ä… = cos(
x)
ìÅ‚ ÷Å‚
12. ,
2Å‚Å‚
íÅ‚
Ä„
ëÅ‚x öÅ‚
cos Ä… = mðsin(
x)
ìÅ‚ ÷Å‚
13. .
2Å‚Å‚
íÅ‚
Tabela 2 Wartości funkcji trygonometrycznych dla wybranych kątów.
x sin(x) cos(x) tg(x)
00 0 1 0
1
1
3
300
2 3
2
1 1
450 1
2 2
1
3
600 3
2
2
900 1 0 Ä… "
Wykresy podstawowych funkcji
Zakład Biofizyki CMUJ 3
Elementy matematyki
4 10
8
3
6
2
y
y
y = x
4
y=x0.5
1
2
0 0
0 2 4 6 8 10
0 2 4 6 8 10
x
x
100 4
80
3
60
2
y y
y=x2 y=x-1
40
1
20
0 0
0 2 4 6 8 10 0 2 4 6 8 10
x x
Rys. 2 Wykres funkcji liniowej. Funkcja liniowa opisuje m.in. zwiÄ…zek drogi s
przebytej przez ciaÅ‚o i czasu t w ruchu jednostajnym: s = v · t, gdzie v to
prędkość ciała. Wykres przedstawiony na rysunku opisuje zatem np.
zależność drogi s od czasu t w ruchu jednostajnym z prędkością v=1m/s.
4 10
8
3
6
2
y
y=x
y=x0.5 y 4
1
2
0 0
0 2 4 6 8 10 0 2 4 6 8 10
x x
100 4
80
3
60
2
2
y
y = x
y
40 y=x-1
1
20
0 0
0 2 4 6 8 10
0 2 4 6 8 10
x
x
Rys. 3 Wykres funkcji kwadratowej. Funkcja kwadratowa opisuje m.in. zwiÄ…zek
drogi s przebytej przez ciało i czasu t w ruchu jednostajnie przyspieszonym:
aÅ" t2
s = v0 Å" t+ , gdzie a oznacza przyspieszenie ciaÅ‚a, a v to poczÄ…tkowa
0
2
prędkość ciała. Wykres przedstawiony na rysunku opisuje zatem np.
zależność drogi s od czasu t w ruchu jednostajnie przyspieszonym z
prędkością początkową v =0m/s, z przyspieszeniem a=1m/s2.
0
Zakład Biofizyki CMUJ 4
Elementy matematyki
4
10
8
3
6
2
y
y
4 y=x
1
2
0 0
0 2 4 6 8 10
0 2 4 6 8 10
x
x
100
4
1
2
y= x = x
80
3
60
2
y y=x2
y
y=x-1
40
1
20
0 0
0 2 4 6 8 10
0 2 4 6 8 10
x
x
Rys. 4 Wykres funkcji pierwiastek kwadratowy. Funkcja pierwiastek kwadratowy
opisuje m.in. związek prędkości v ciała spadającego swobodnie z wysokości
v = 2hg
h: , gdzie g oznacza przyspieszenie ziemskie.
4
10
8
3
6
2
y
y=x0.5
y
4 y=x
1
2
0
0
0 2 4 6 8 10
0 2 4 6 8 10
x
x
100 4
80
3
60
2
y y=x2 y
y = x -1
40
1
20
0 0
0 2 4 6 8 10 0 2 4 6 8 10
x x
Rys. 5 Wykres funkcji hiperbolicznej. Funkcja hiperboliczna opisuje m.in. zwiÄ…zek
ciśnienia p i objętości V gazu doskonałego, poddawanego przemianie
kRT
p=
izotermicznej: , gdzie k jest liczbą moli gazu, R jest stałą gazową,
V
a T oznacza temperaturÄ™ przemiany.
Zakład Biofizyki CMUJ 5
Elementy matematyki
25000 1.0
20000 0.8
15000 0.6
y = e x
y=e-x
y y
10000 0.4
0.2
5000
0.0
0
0 2 4 6 8 10 0 2 4 6 8 10
x
x
1.0
2
0.8
0
0.6 2
y=e-2(x-5)
y=ln(x)
y y
0.4 -2
0.2
-4
0.0
0 2 4 6 8 10 0 2 4 6 8 10
x x
Rys. 6 Wykres funkcji wykładniczej y = ex (oznacza się również e x = exp(x)).
Funkcja wykładnicza opisuje m.in. wzrost liczby zachorowań w począ-
tkowej fazie epidemii.
25000 1.0
20000 0.8
15000 0.6
y=ex
y = e - x
y
y
10000 0.4
0.2
5000
0.0
0
0 2 4 6 8 10 0 2 4 6 8 10
x
x
1.0
2
0.8
0
0.6 2
y=e-2(x-5)
y=ln(x)
y y
0.4 -2
0.2
-4
0.0
0 2 4 6 8 10 0 2 4 6 8 10
x
x
Rys. 7 Wykres funkcji wykładniczej (argument ujemny). Prezentowana funkcja
opisuje m.in. prawo rozpadu promieniotwórczego: N = N · e - t, gdzie N
0
oznacza liczbę jąder promieniotwórczych które nie uległy rozpadowi do
czasu t, N to początkowa liczba jąder, a stała , zwana stałą rozpadu, jest
0
odwrotnością czasu połowicznego rozpadu pierwiastka promienio-
twórczego.
Zakład Biofizyki CMUJ 6
Elementy matematyki
25000 1.0
20000 0.8
15000 0.6
y=ex
y=e-x
y y
10000 0.4
0.2
5000
0.0
0
0 2 4 6 8 10 0 2 4 6 8 10
x
x
2
1.0
0.8
0.6 2
0
y=e-2(x-5)
y = log(x)
y y
0.4
0.2
0.0
-2
0 2 4 6 8 10 0 2 4 6 8 10
x
x
Rys. 8 Funkcja logarytm dziesiętny. Funkcja ta opisuje m.in. związek natężenia
dz
wiÄ™ku I (wyrażonego w W·m-2) z wrażeniem sÅ‚uchowym ² (wyrażonym
² = logëÅ‚ I I0 öÅ‚
w belach): , gdzie I =10-12 W·m-2.
ìÅ‚ ÷Å‚ 0
íÅ‚ Å‚Å‚
25000 1.0
20000 0.8
15000 0.6
y=ex
y=e-x
y y
10000 0.4
0.2
5000
0.0
0
0 2 4 6 8 10 0 2 4 6 8 10
szerokość
x
x
połówkowa
2
1.0
2
y = e - 2 ( x - 5 )
0.8
wartość
0.6
0
oczekiwana
y=log(x)
y y
0.4
0.2
0.0
-2
0 2 4 6 8 10 0 2 4 6 8 10
x
x
Rys. 9 Funkcja Gaussa. Funkcja ta może opisywać m.in. rozkład częstości
występowania cech ilościowych w populacji. Jeżeli funkcję Gaussa
2
ëÅ‚
x- x0 öÅ‚
ëÅ‚ öÅ‚
ìÅ‚ ÷Å‚
f (x) = A Å"exp-ìÅ‚ ÷Å‚
zapiszemy w postaci , to wartość oczekiwana
ìÅ‚ ÷Å‚
s
íÅ‚ Å‚Å‚
íÅ‚ Å‚Å‚
jest równa x , a szerokość połówkowa jest równa 2Å" sÅ" ln2
.
0
Zakład Biofizyki CMUJ 7
Elementy matematyki
1.0 1.0
y=cos(x)
y=sin(x)
0.5 0.5
0.0 0.0
y y
-0.5 -0.5
-1.0 -1.0
0 2 4 6 8 10 12 0 2 4 6 8 10 12
x x
4
2
0
y
-2
y=tg(x)
-4
0 2 4 6 8 10 12
x
Rys. 10 Wykresy funkcji trygonometrycznych. Funkcje trygonometryczne sinus i
cosinus są używane m.in. w opisie zjawisk periodycznych.
II. ELEMENTY RACHUNKU RÓŻNICZKOWEGO I CAA
KOWEGO
Pochodna funkcji jednej zmiennej
dy
,
Pochodna funkcji jednej zmiennej y = f(x), oznaczana symbolami : y', f '(x),
dx
df(x)
jest to nowa funkcja zmiennej x, spełniająca następującą zależność:
dx
f(x + "x) - f(x)
f '(x) H" .
"x
Powyższe przybliżenie jest tym dokładniejsze, im mniejszy jest przyrost "x (Rysunek
11). Obliczanie pochodnej f '(x) nazywamy różniczkowaniem funkcji f(x).
Zakład Biofizyki CMUJ 8
Elementy matematyki
Interpretacja geometryczna pochodnej: Jeżeli wykresem funkcji y = f(x) w układzie
współrzędnych prostokątnych jest pewna krzywa, to wartość pochodnej f'(x) w punkcie o
współrzędnych (x,f(x)) równa się tgą, gdzie ą jest kątem między osią OX i styczną do
krzywej w punkcie (x,f(x)); kÄ…t ten liczy siÄ™ od dodatniego kierunku osi OX w kierunku
przeciwnym do ruchu wskazówek zegara.
dy "y
6.0
= tg( H"
Ä…)
15 "y<0
5.5
dx "x
5.0
"x>0
4.5
10
4.0
y
.
3.5
0.18 0.21 0.24 0.27
5
Ä…
0
0.0 0.2 0.4 0.6 0.8 1.0
x
Rys. 11 Geometryczna interpretacja pochodnej funkcji.
Szereg wielkości fizycznych jest zdefiniowana jako pochodne innych wielkości, np.
ds
v =
prędkość v jest pochodną po czasie t drogi s przebytej w trakcie ruchu tj. ,
dt
dv
a =
przyspieszenie a jest pochodną po czasie prędkości v, tj. , moc P jest pochodną po
dt
dW
P =
czasie pracy W wykonanej przez pewną siłę, tj. , prąd dyfuzyjny j w
dt
niejednorodnym roztworze pewnej substancji jest proporcjonalny do pochodnej stężenia Á tej
dÁ
j = D
substancji po współrzędnych przestrzennych x, tj. , gdzie D jest współczynnikiem
dx
dyfuzji.
Prawa fizyki formułuje się w postaci tzw. równań różniczkowych, tj. równań wiążących ze
sobą różne wielkości oraz ich pochodne, np. druga zasada dynamiki Newtona może
przyjmować postać:
ds
dëÅ‚ öÅ‚
ìÅ‚ ÷Å‚
dv dtłł
íÅ‚
F = ma= m = m
dt dt
Zakład Biofizyki CMUJ 9
Elementy matematyki
gdzie F to siła, a m masa ciała.
Zmiana stężenia c substancji A w przebiegu prostej reakcji typu A çÅ‚çÅ‚k B, gdzie k jest staÅ‚Ä…
A
szybkości reakcji opisana jest równaniem:
dcA
= -kcA
dt
Podstawowe reguły różniczkowania
Oznaczenia: u,v,w,...- funkcje zmiennej niezależnej x; u',v', w',...- pochodne tych
funkcji względem x.
1. Pochodna sumy algebraicznej dwóch lub kilku funkcji jest równa sumie algebraicznej
pochodnych każdej z tych funkcji:
(u + v  w + ...)' = u' + v'  w' + ... ;
2. Pochodna iloczynu dwóch funkcji jest równa sumie:
(uv)' = u'v + uv',
3. Stały czynnik można wynosić przed znak pochodnej
(cu)' = cu';
4. Pochodną ilorazu oblicza się według wzoru
'
u vu'-uv'
ëÅ‚ öÅ‚
= .
ìÅ‚ ÷Å‚
vłł v2
íÅ‚
6. Twierdzenie o pochodnej funkcji zÅ‚ożonej. Jeżeli y = f(u) i u = Õ(x), to
dy df(u) dÕ(x)
= Å"
;
dx du dx
Przykładowo wychylenie x ciała z położenia równowagi w trakcie ruchu drgającego
(np. ciało zaczepione na sprężynie, poruszające się bez tarcia w płaszczyz
nie poziomej) jest
opisane wzorem:
x = A sin( + Ć)
2Ä„½t
;
gdzie A jest amplitudÄ… ruchu, ½ - czÄ™stotliwoÅ›ciÄ… drgaÅ„, t  czasem, a Ć  kÄ…tem fazowym.
Prędkość v w ruchu drgającym dana jest wzorem (pochodne funkcji elementarnych są podane
w Tabeli 3):
dx
v = = 2Ä„½A cos( + Ć)
2Ä„½t
.
dt
Zakład Biofizyki CMUJ 10
Elementy matematyki
Tabela 3 Pochodne funkcji elementarnych.
Funkcja Pochodna
stała 0
x n n x n - 1
1 1
-
x
x2
1
x
2 x
ex ex
1
ln x
x
sin(x) cos(x)
cos(x) -sin(x)
1
tg(x)
2
cos x
Całka nieoznaczona i oznaczona funkcji
FunkcjÄ… pierwotnÄ… danej funkcji f(x) jednej zmiennej nazywamy takÄ… funkcjÄ™ F(x),
której pochodna jest równa f(x) tzn. F'(x)=f(x). Jeśli F(x) jest funkcją pierwotną funkcji f(x),
to każda funkcja F(x)+C, gdzie C jest dowolną stałą jest również funkcją pierwotną funkcji
f(x). Wyrażenie ogólne F(x)+C dla wszystkich funkcji pierwotnych danej funkcji f(x)
nazywamy całką nieoznaczoną funkcji f(x) i oznaczamy:
F(x) +C = f (
+"x)dx.
Jeśli F(x) jest funkcją pierwotną funkcji f(x) to całka oznaczona funkcji f(x) w
przedziale od a do b jest określona wzorem:
b
f
+"(x)dx= F(b) -F(a) .
a
Zakład Biofizyki CMUJ 11
Elementy matematyki
Y
b
+"f ( x ) d x
a
0
a b
X
Rys. 12 Geometryczna interpretacja całki oznaczonej.
Jeśli dana funkcja f(x) przedstawiona jest w postaci krzywej we współrzędnych XY, to
wartość całki oznaczonej funkcji f(x) w granicach od a do b jest równa polu ograniczonemu
krzywą y = f(x), osią 0X i dwiema rzędnym, stałą rzędną w punkcie x = a i stałą rzędną w
punkcie x = b (Rysunek 12).
Szereg wielkości fizycznych definiuje się jako całki z innych wielkości fizycznych,
W = F
+"Å"ds
np. praca W jest całką z siły F na drodze S, tj. , masa M substancji
S
M = Á
+"Å"dV
niejednorodnej zawartej w objÄ™toÅ›ci V jest caÅ‚kÄ… z gÄ™stoÅ›ci Á tej substancji, tj. ,
V
natężenie I promieniowania rentgenowskiego po przejściu przez ośrodek niejednorodny
ëÅ‚ öÅ‚
ìÅ‚- (
I =I Å"exp µ
wzdÅ‚uż drogi S jest równe 0 , gdzie µ(x) jest współczynnikiem
+"x)dx÷Å‚
ìÅ‚ ÷Å‚
íÅ‚ S Å‚Å‚
osłabienia w punkcie o współrzędnej x na drodze S, a I jest natężeniem promieniowania
0
przed wejściem do ośrodka pochłaniającego to promieniowanie.
III. TRANSFORMATA FOURIERA, SZEREGI FOURIERA.
Każdą funkcję okresową f(x), o okresie T (tzn. taką, że f(x)=f(x+T)), można przedstawić w
postaci tzw. szeregu Fouriera (pod warunkiem, że spełnia ona kilka dodatkowych warunków
 musi być ograniczona i przedziałami ciągła):
" "
1
f (x) = a0 + cos(Éx) + sin(
n nÉx)
,
"an "bn
2
n=1 n=1
"
gdzie É=2Ä„/T. Symbol , gdzie W jest dowolnym wyrażeniem zależnym od n,
"W n
n
n=1
oznacza nieskończoną sumę W +W +W +.... Współczynniki a i b wyznacza się ze wzorów:
1 2 3 n n
Zakład Biofizyki CMUJ 12
Elementy matematyki
T
2
an = f n n = 0, 1, 2, 3...
+"(x)cos(Éx)dx,
T
0
T
2
bn = f nÉx)dx. n = 1, 2, 3,...
+"(x)sin(
T
0
Szereg Fouriera można zapisać również w postaci:
"
1
f (x) = a0 + sin( + Õn)
nÉx
,
"A n
2
n=1
tg( = an / bn
Õn)
gdzie A = a2 + b2 , a . Znajdowanie współczynników szeregu Fouriera
n n n
danej funkcji f(x) nosi nazwÄ™ analizy harmonicznej.
Przykład 1
Na rysunku 13 przedstawiono wykresy prostych sygnałów sinusoidalnych, ich złożenie oraz
odpowiadające im widma fourierowskie. Widmo fourierowskie jest wykresem zależności
amplitudy od częstotliwości i stanowi graficzną metodę prezentacji szeregu Fouriera danej
funkcji. W przypadku funkcji okresowych widma fourierowskie mają postać dyskretną, tj.
składają się z przeliczalnej liczby punktów nie leżących na osi X. W przypadku funkcji
harmonicznych widmo składa się z dokładnie jednego punktu (Rysunek 13, pierwsze trzy
linie). W przypadku funkcji, których szereg fourierowski zawiera więcej niż jeden wyraz
widmo fourierowskie zbudowane jest odpowiednio z większej liczby punktów (Rysunek 13,
czwarta linia). Widma fourierowskie można tworzyć także dla funkcji nieokresowych. Szereg
Fouriera zastępowany jest wówczas całką Fouriera, a widmo fourierowskie przedstawia
funkcję ciągłą.
Przykład 2
W poprzednim przykładzie przedstawiono widma fourierowskie szeregów zbudowanych ze
skończonej liczby wyrazów. W przykładzie bieżącym zaprezentowany jest szereg nieskończony
i widmo funkcji okresowej o okresie T, określonej wzorem:
1 3
Å„Å‚
0< t< T, T< t< T
2 dla
ôÅ‚
4 4
,
f(t)=
òÅ‚0 dla
1 3
ôÅ‚
T< t< T
ół 4 4
Funkcję taką można rozwinąć w szereg Fouriera:
4 2Ä„t 1 2Ä„t 1 2Ä„t
ëÅ‚cos - cos + cos +...öÅ‚
f (t) =1+ 3 5
ìÅ‚ ÷Å‚
.
Ä„ T 3 T 5 T
íÅ‚ Å‚Å‚
Zakład Biofizyki CMUJ 13
Elementy matematyki
Na rysunku 14 przedstawiono wykres funkcji f(t) dla T=1 (górny wykres, lewa strona), widmo
fourierowskie funkcji f(t) (górny wykres, prawa strona) oraz wykresy sumy pierwszych pięciu
f (t) i dziesięciu f (t) wyrazów powyższego szeregu (dolne wykresy).
5 10
Przykład 3
Wykresy przedstawione na rysunku 15 wykonano w oparciu o cyfrową rejestrację dwóch
głosek "a" wypowiedzianych przez tą samą osobę głosem, którego wysokość różniła się o
oktawę. Bezpośrednio rejestrowaną wielkością jest wielkość wychylenia membrany
mikrofonu, przetworzona na sygnał elektryczny. Lewa strona, przedstawiająca wychylenie
membrany w zależności od czasu, ilustruje fakt, że dz
więk jest sygnałem periodycznym. Prawe
wykresy to widma fourierowskie zarejestrowanych dz
więków. Z porównania wykresów widać,
że słyszana oktawa różnicy wysokości dz
więków odpowiada podwojeniu częstotliwości.
4
10
5
3
0
2
-5
1
-10 y1=2sin(2Ä„x)
0
0 1 2 3 0 1 2 3 4
czas [s] częstotliwość [Hz]
10 4
5
3
0
2
-5
1
-10
y2=3sin(4Ä„x)
0
0 1 2 3 0 1 2 3 4
czas [s] częstotliwość [Hz]
Zakład Biofizyki CMUJ 14
wychylenie
amplituda
amplituda
wychylenie
Elementy matematyki
4
10
5
3
0
2
-5
1
-10 y3=4sin(6Ä„x)
0
0 1 2 3 0 1 2 3 4
czas [s]
częstotliwość [Hz]
10 4
5
3
0
2
-5
1
-10
y4=y1+y2+y3
0
0 1 2 3 0 1 2 3 4
czas [s] częstotliwość [Hz]
Rys. 13 Sygnały periodyczne (lewa kolumna) i ich
transformaty fourierowskie (prawa kolumna).
1.4
2.0
1.2
1.0
1.5
0.8
1.0
0.6
0.4
0.5
0.2
0.0
0.0
0 1 2 3 4 5 6 7 8
0.0 0.5 1.0 1.5 2.0 2.5 3.0
t częstotliwość [Hz]
2.5
2.5
2.0
2.0
1.5
1.5
1.0
1.0
0.5 0.5
0.0 0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0
t t
Rys. 14 Sumowanie składowych harmonicznych sygnału okresowego pozwala uzyskać tym
dokładniejsze przybliżenie sygnału oryginalnego, im większa ilość składowych
zostanie wysumowana. Na rysunku przedstawiono sygnał prostokątny i jego widmo
fourierowskie (górne wykresy) oraz wykresy sum pierwszych pięciu i pierwszych
dziesięciu składowych harmonicznych sygnału prostokątnego.
Zakład Biofizyki CMUJ 15
wychylenie
amplituda
amplituda
wychylenie
f(t)
amplituda
5
10
f (t)
f (t)
Elementy matematyki
20
200 "A" niskie
15
150
10
100
5
50
0
1.500 1.505 1.510 1.515 1.520 0 500 1000 1500 2000
ot Hz]
czas częst liwość [
200 20
15
150 "A" wysokie
10
100
5
50
0
1.500 1.505 1.510 1.515 1.520 0 500 1000 1500 2000
czas
częst
otliwość [Hz]
20
200 "A" niskie
15
150
10
100
5
50
0
1.500 1.505 1.510 1.515 1.520 0 500 1000 1500 2000
czas częstotliwość [Hz]
200 20
15
150 "A" wysokie
10
100
5
50
0
1.500 1.505 1.510 1.515 1.520 0 500 1000 1500 2000
czas
częstotliwość [Hz]
Rys. 15 Przykłady dz
więków mowy przedstawionych w dziedzinie czasu (lewa kolumna) i
w dziedzinie częstotliwości (prawa kolumna).
IV. ELEMENTY RACHUNKU BA

DÓW
y
ródła i rodzaje błędów
Celem pomiaru jest ustalenie wartości liczbowej miary danej wielkości fizycznej.
Każdy pomiar fizyczny obarczony jest niemożliwym do uniknięcia błędem pomiarowym.
Pomiary wielkości fizycznych są bezużyteczne, jeśli w ich wyniku nie zostanie oszacowana
wielkość błędu pomiarowego. Ze względu na sposób w jaki błędy wpływają na wynik
pomiaru dzielimy je na błędy systematyczne, przypadkowe i grube.
Błędy systematyczne zawsze w ten sam sposób wpływają na wyniki pomiarów wykonanych
za pomocą tej samej metody i aparatury pomiarowej. Przy zmianie warunków pomiaru
według określonej prawidłowości błąd systematyczny zachowuje stałą wartość lub zmienia
się w sposób prawidłowy. Minimalna wartość błędu systematycznego jest określona
dokładnością stosowanego przyrządu t.j. połową wartości najmniejszej działki skali przyrządu
"x. Wynik pojedynczego pomiaru wielkoÅ›ci fizycznej zapisujemy zatem w postaci x Ä… ½·"x.
Błędy przypadkowe ujawniają się najczęściej w przypadku wielokrotnego pomiaru danej
wielkości przyrządem, którego dokładność jest bardzo duża (a więc błąd systematyczny
pomiaru jest mały). W takim przypadku może się zdarzyć, że różnice między wynikami
kolejnych pomiarów będą znacznie przewyższać błąd systematyczny. Wiele przyczyn może
spowodować powstanie błędu przypadkowego. Może on wynikać z własności przedmiotu
mierzonego, np. przy pomiarze średnicy drutu może być efektem wahań średnicy. Innym jego
z
ródłem są własności samego przyrządu pomiarowego, którego wskazania zależą od
przypadkowych drgań budynku, ruchów powietrza, tarcia w łożyskach, docisku (np.
mikromierza). Błędy przypadkowe mogą mieć za przyczynę również podłoże fizjologiczne
Zakład Biofizyki CMUJ 16
amplituda
wychylenie
amplituda
wychylenie
amplituda
wychylenie
amplituda
wychylenie
Elementy matematyki
np. zdolność do określenia równości oświetlenia poszczególnych części pola widzenia lub
usłyszenia maksimum natężenia dz
więku. Błędów przypadkowych nie można wyeliminować,
lecz można określić ich wpływ na wynik ostateczny.
Błędy grube lub pomyłki wynikają z niestaranności eksperymentatora. Celem
wyeliminowania takich błędów należy powtórzyć pomiary. Efektem popełnienia błędu
grubego podczas ćwiczeń na Pracowni Biofizyki i umieszczenia wyników obarczonych
grubym błędem w sprawozdaniu może być brak zaliczenia ćwiczenia.
Wynik jakiegokolwiek pomiaru wielkości x podawany jest w następujący sposób:
(wartość zmierzona x) = x Ä…´x.
np
Zapis ten oznacza, że:
" w wyniku eksperymentu otrzymano liczbę x jako najlepsze przybliżenie wartości
np
mierzonej,
" z rozsądnym prawdopodobieństwem szukana wielkość należy do przedziału (x
np-´x,x +
np
´x).
Liczba ´x zwana jest niepewnoÅ›ciÄ… lub bÅ‚Ä™dem pomiaru x.
Przykład 1.
W wyniku pomiaru pewnej wielkości otrzymano liczbę 53,73 jako oszacowanie średniej i
liczbę 0.3 jako oszacowanie niepewności. Prawidłowy sposób przedstawienia wyniku to:
53,7 Ä… 0,3.
Jeśli niepewność równa jest 3, to ten sam rezultat należałoby zapisać jako
54 Ä… 3,
jeśli zaś niepewność wynosi 30, to odpowiedz
powinna brzmieć:
50 Ä… 30.
Jest jeden wyjątek od przedstawionej powyżej reguły. Jeśli pierwsza cyfra niepewności jest
mała (1 lub być może 2), to mogłoby być właściwe pozostawienie w odpowiedzi jeszcze jednej
cyfry znaczącej. Przykładowo, wynik taki jak
zmierzona długość = 27,6 ą 1 cm
jest zupełnie rozsądny. W tym przypadku należałoby się zgodzić, że jego zaokrąglenie do 28 ą
1 powodowałoby utratę informacji.
Liczby używane w obliczeniach powinny mieć jedną cyfrę znaczącą więcej niż te podawane
ostatecznie. Zmniejsza to niedokładności wprowadzone podczas zaokrąglania liczb. Końcowy
wynik powinien być zaokrąglony tak, aby usunąć tę dodatkową (i nieznaczącą) cyfrę.
Rozkłady empiryczne
Zakład Biofizyki CMUJ 17
Elementy matematyki
Rezultaty eksperymentu polegającego na wielokrotnym pomiarze pewnej wielkości
można przedstawić przy pomocy wykresu rozkładu wyników pomiaru, nazywanego
histogramem.
Przykład 2
W trakcie eksperymentu dokonano rejestracji sygnału EKG, a następnie wyznaczono czas R-
R między kolejnymi załamkami R sygnału. W wyniku analizy sygnału zbudowano następujący
rozkład wyników pomiaru czasu R-R:
Przedział wartości 0,95- 1,00-
0,75-0,8 0,8-0,85 0,85-0,9 0,9-0,95 1,05-1,1
czasu R-R (sekundy) 1,00 1,05
Liczba obserwacji w
3 38 55 174 92 29 4
przedziale
Wykres przedstawiający ilość zaobserwowanych wartości czasu R-R w kolejnych
przedziałach jest histogramem rozkładu wyników eksperymentu (Rysunek 16).
150
100
50
0
0.75 0.80 0.85 0.90 0.95 1.00 1.05 1.10
przedział wartości czasu R-R
Rys. 16 Typowy histogram czasów R-R.
Przykład 3
W trakcie eksperymentu dokonano rejestracji sygnału EKG osoby z mocno zaznaczoną
arytmią oddechową. Wyznaczono czas R-R między kolejnymi punktami R sygnału, wyniki
przedstawiono w postaci histogramu (Rysunek 17).
Zakład Biofizyki CMUJ 18
Liczba zliczeń
Elementy matematyki
20
10
0
0.6 0.8 1.0 1.2 1.4
przedział wartości czasu R-R
Rys. 17 Histogram R-R sygnału EKG, w którym zaznacza się
wyraz
nie arytmia oddechowa.
Przykład 4
W izbie przyjęć notowano ilość zachorowań tygodniowo na pewną chorobę. Obserwacje
prowadzone przez długi okres pozwoliły na stworzenie następującego wykresu rozkładu
częstości zachorowań. Histogram ten jest wyraz
nie asymetryczny (Rysunek 18).
20
15
10
5
0
1 2 3 4 5 6 7 8 9 10
liczba zachorowań
Rys. 18 Typowy histogram częstości zachorowań w jednostce czasu
(tu tygodniowo)  nie dotyczy przypadków epidemii.
Statystyka opisowa
Zakład Biofizyki CMUJ 19
liczba zliczeń
zanotowano
daną liczbę zachorowań
liczba tygodni w których
Elementy matematyki
Podstawowa ilościowa analiza wyników eksperymentu prowadzi do określenia liczb
dających opis zbiorowości, nazywanych parametrami statystycznymi. Parametry statystyczne
tak charakteryzują zbiorowość, że porównywanie różnych zbiorowości statystycznych można
sprowadzić do porównań parametrów statystycznych. Podstawowe zadania tych parametrów
to:
1. Określenie przeciętnego rozmiaru i rozmieszczenia wartości zmiennej. Dokonujemy tego
przez obliczenie miar położenia, przy czym najczęściej stosowaną miarą położenia są
wartości średnie. Średnią arytmetyczną definiujemy następującym wzorem:
n
1 x1 + x2 + x3 + ...+ xn
x = xi =
,
"
n n
i=1
gdzie x , x , ..., x to ciąg zmierzonych wartości zmiennej.
1 2 n
2. Określenie granic obszaru zmienności wartości zmiennej. Dokonujemy tego przez
obliczenie miar zmienności. W praktyce najczęściej stosowanymi miarami zmienności
są wariancja i odchylenie standardowe. Wariancją zmiennej X nazywamy średnią
arytmetyczną kwadratów odchyleń poszczególnych wartości zmiennej od średniej
arytmetycznej całej zbiorowości:
n
1
2
Ã2 = - x)
"(xi
n
i=1
Pierwiastek kwadratowy z wariancji nazywany jest odchyleniem standardowym i wyraża się
wzorem:
n
1
2
à = - x)
"(xi
n
i=1
Rozkłady modelowe
W miarę jak liczba wykonanych pomiarów zwiększa się, histogram może zbliżać się
do pewnej określonej ciągłej krzywej, nazywanej rozkładem granicznym. W przypadku
wielokrotnego pomiaru, obciążonego błędami przypadkowymi rozkładem granicznym jest
rozkład Gaussa (rozkład normalny) zdefiniowany wzorem:
ëÅ‚ - x)2 ÷Å‚
öÅ‚
1 (x
ìÅ‚
f (x) = exp-
ìÅ‚
à 2Ä„ 2Ã2 ÷Å‚
íÅ‚ Å‚Å‚
Zakład Biofizyki CMUJ 20
Elementy matematyki
gdzie à i x sÄ… parametrami rozkÅ‚adu. Współczynnik umieszczony po prawej stronie
powyższego równania przed funkcją exp() został wybrany w ten sposób, aby
prawdopodobieństwo P(x") otrzymania wyniku x pomiaru w przedziale
spełniało równość:
b
P(x") = f
+"(x)dx
a
Rozkładem Gaussa (rysunek 19) przybliża się szereg ilościowych cech
fizjologicznych, anatomicznych itp. na przykład rozkład wysokości lub masy ciała w
populacji, czy rozkład ilości punktów uzyskanych na egzaminie poprawkowym z biofizyki.
Cechą rozkładu Gaussa jest to, że wartość średnia wyników pomiarów, podlegających temu
rozkładowi zbliża się w miarę wzrostu liczby pomiarów do parametru rozkładu, a
x
odchylenie standardowe wyników zbliża siÄ™ do parametru à rozkÅ‚adu.
Rozkładem granicznym opisującym wyniki eksperymentu polegającego na zliczaniu
zdarzeń występujących losowo, ale z określoną przeciętną częstotliwością jest rozkład
Poissona. W odróżnieniu od rozkładu Gaussa, który jest funkcją określoną na zbiorze liczb
rzeczywistych, dziedziną rozkładu Poissona jest z definicji zbiór liczb naturalnych. Rozkład
Poissona jest określony wzorem:
N
P(N) = Å"exp(
-)
N!
gdzie  jest parametrem rozkładu. P(N) jest równa prawdopodobieństwu, że w zadanej
jednostce czasu zdarzenie zajdzie dokładnie N razy. Rozkład Poissona (Rysunek 19) opisuje
m.in. eksperyment polegający na zliczaniu liczby rozpadów promieniotwórczych w jednostce
czasu, liczby wyjazdów karetki pogotowia do pacjentów w ciągu jednego dnia lub ilość
podejść do egzaminu z biofizyki w ciągu semestru. Parametr  rozkładu Poissona jest
oczekiwaną średnią liczbą zliczeń w wybranym przedziale czasu. Rozkład Poissona ze

średnią liczbą zliczeń  ma odchylenie standardowe .
0.4
0.3
0.2
0.1
0.0
-6 -4 -2 0 2 4 6
x
Zakład Biofizyki CMUJ 21
f(x)
Elementy matematyki
Rys. 19 RozkÅ‚ad Gaussa dla parametrów x= 0 i à = 1.
0.20
0.15
0.10
0.05
0.00
0 2 4 6 8 10
N
Rys. 20 Rozkład Poissona z oczekiwaną liczbą zliczeń w
jednostce czasu równą 4. Odchylenie standardowe
à =  = 2
tego rozkładu jest równe .
Przenoszenie niepewności
Większości wielkości fizycznych nie da się określić na podstawie bezpośredniego
pomiaru. Zawsze jednak jest możliwe wyznaczenie wielkości interesującej ekperymentatora
poprzez bezpośredni pomiar innych wielkości, powiązanych związkiem funkcyjnym z
szukaną wielkością. Przykładowo aby znalez
ć powierzchnię prostokąta, trzeba zmierzyć
długości jego boków l oraz h, i na tej podstawie obliczyć powierzchnię A, korzystając ze
zwiÄ…zku A = l·h. Prawie wszystkie interesujÄ…ce doÅ›wiadczenia skÅ‚adajÄ… siÄ™ z dwóch etapów,
bezpośredniego pomiaru i następujących po nim obliczeń.
Kiedy pomiar dzieli się na dwa etapy, wówczas ocena niepewności również przebiega
dwuetapowo. Po pierwsze należy ocenić niepewności wielkości mierzonych bezpośrednio,
następnie zaś stwierdzić, w jaki sposób owe niepewności "przenoszą się" w trakcie obliczeń
na niepewność ostatecznego wyniku.
W przypadku gdy poszukiwana wielkość  oznaczana dalej symbolem q  jest funkcją
tylko jednej mierzonej wielkości  oznaczanej dalej x  tzn. q=q(x) (np. objętość sześcianu
jest równa trzeciej potędze długości boku sześcianu), to jeśli x jest najlepszym
np
oszacowaniem x, wówczas najlepszym oszacowaniem q jest q(x ). Jeśli ekstremalne (tj.
np
najmniejsza i najwiÄ™ksza) prawdopodobne wartoÅ›ci x równe sÄ… x +´x, to odpowiadajÄ…ce im
np.
ekstremalne wartoÅ›ci q równe sÄ… q(x Ä… ´x). KorzystajÄ…c z przybliżenia:
np.
dq
q(x + u) H" q(x) + u
,
dx
(dla dowolnego przyrostu u) można przepisać ekstremalne wartości jako
dq
q(xnp) Ä… ´x.
dx
Zakład Biofizyki CMUJ 22
P(N)
Elementy matematyki
W powyższym wzorze wartość bezwzględna pozwala na uwzględnienie sytuacji, gdy
dq/dx jest ujemne. Zaprezentowany wynik oznacza, że niepewność oszacowania q wynosi
´q H" ïÅ‚dq/dx ïÅ‚´x.
W przypadku, gdy q jest funkcją dwóch zmiennych tzn q = q(x, y) i jeśli x i y są
np np
najlepszymi przybliżeniami x i y, to spodziewamy się, że najlepsze przybliżenie q równe jest
q = q(x y ).
np np, np
Maksymalna niepewność wyznaczenia wartości q(x,y) równa jest:
" q
q "
´ H" ´ + ´
q x y.
" y
x "
gdzie ´x jest niepewnoÅ›ciÄ… wyznaczenia wielkoÅ›ci x, ´y jest niepewnoÅ›ciÄ… wyznaczenia
"
q
"
q
wielkości y, a symbole i (wprowadzane w przypadku funkcji wielu zmiennych)
"
y
"
x
oznaczajÄ… odpowiednio pochodnÄ… funkcji q(x,y) po zmiennej x i pochodnÄ… funkcji q(x,y) po
zmiennej y. (Pochodną funkcji q(x,y) po zmiennej x oblicza się zgodnie z regułami podanymi
w rozdziale Elementy rachunku różniczkowego i całkowego, traktując zmienną y tak jak stałą.
Podobnie postępuje się w celu obliczenia pochodnej funkcji q(x,y) po zmiennej y  wówczas
zmiennÄ… x traktuje siÄ™ tak jak staÅ‚Ä…: np. pochodna funkcji q(x,y) = x2·y po zmiennej x jest
"q
"q
= x2
= 2Å" xÅ" y
równa , a pochodna po y jest równa ).
"y
"x
JeÅ›li niepewnoÅ›ci ´x i ´y sÄ… niezależne i przypadkowe, to powyższa suma
zastępowana jest przez pierwiastek z sumy kwadratów. Jeśli funkcja q zależy od więcej niż
dwóch zmiennych, to po prostu dodajemy kolejne wyrażenie odpowiadające każdej następnej
zmiennej. Podsumowując, można sformułować następującą regułę:
JeÅ›li wielkoÅ›ci x, ...,z zmierzone z niepewnoÅ›ciami ´x, ...,´z sÅ‚użą do obliczenia
wartości funkcji q(x, ..., z) i jeśli niepewności wyznaczenia x, ...,z są niezależne i
przypadkowe, to niepewność wyznaczenia wartości funkcji q równa jest
2 2
"q "q
ëÅ‚ ëÅ‚
´q = ´xöÅ‚ + ...+ ´zöÅ‚ .
ìÅ‚ ÷Å‚ ìÅ‚ ÷Å‚
"x "z
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
W żadnym jednak wypadku niepewność ta nie jest większa niż zwykła suma:
" "
q q
´ d" ´ +... ´
q x z.
" "
x z
Przykład 5
W celu wyznaczenia gęstości odczynnika ustalono z pomocą wagi jego masę m=19,236g oraz
przy pomocy menzurki objÄ™tość V=23,5cm3. Wyznaczono gÄ™stość Á=m/V. Jako niepewnoÅ›ci
oszacowania przyjÄ™to minimalne dziaÅ‚ki skali przyrzÄ…dów: ´m=0,001g i ´V=0,1cm3.
Ponieważ błędy wynikające z dokładności skali są błędami systematycznymi, a nie
przypadkowymi, w celu oceny niepewności wyznaczenia gęstości odczynnika skorzystano ze
wzorów:
Zakład Biofizyki CMUJ 23
Elementy matematyki
"Á "Á
´Á= ´m+ ´V,
"m "V
"Á 1
= ,
"m V
"Á m
= .
2
"V V
W wyniku eksperymentu otrzymano wynik Á=0,819Ä…0,004 g/cm3.
Przykład 6
Celem eksperymentu jest wyznaczenie średniej powierzchni S pewnych komórek. Komórki te
mają w przybliżeniu kształt eliptyczny w związku z czym aby obliczyć ich pole należy najpierw
wyznaczyć ich największą i najmniejszą średnicę  oznaczane odpowiednio symbolami a i b 
oraz skorzystać ze związku S=Ąab. W wyniku wielokrotnego pomiaru wielkości a oszacowano
wartość Å›redniÄ…  0,0143mm i odchylenie standardowe ´a=0,0021mm. . W wyniku
wielokrotnego pomiaru wielkości b oszacowano wartość średnią  0,0318mm i odchylenie
standardowe ´b=0,0043mm. NiepewnoÅ›ci ´a i ´b wyznaczono w wyniku wielokrotnego
powtórzenia niezależnych pomiarów, są więc oszacowaniem błędów przypadkowych. Wobec
tego w celu określenia niepewności S należy skorzystać ze związków:
2 2
"S "S
´S = (´a)2 + (´b)2 ,
"a "b
"S
=Ä„
b,
"a
"S
=Ä„
a.
"b
W wyniku eksperymentu otrzymano zatem wynik S=(1,4Ä…0,3)10-3 mm2.
Metoda najmniejszych kwadratów
W praktyce codziennej mamy bardzo często do czynienia z sytuacjami, w których dla
poznania rzeczywistości konieczne jest badanie kilku wielkości równocześnie, ze zwróceniem
szczególnej uwagi na ich wzajemne powiązania ze sobą. Przykładowo można badać zależność
pomiędzy stanem pacjenta a ilością zaaplikowanego leku.
Metoda statystyczna zajmująca się opisem zależności pomiędzy dwoma lub więcej
zmiennymi nosi nazwę analizy korelacji i regresji. Korelacja oznacza fakt współzależności
zmiennych, czyli istnienie powiązania pomiędzy nimi. Siłę i kierunek powiązania określa się
za pomocą współczynnika korelacji. Termin regresja odnosi się natomiast do modelu
matematycznego (przedstawionego w postaci równania matematycznego lub wykresu)
opisującego wspomnianą współzależność.
Wykorzystywane zależności pomiędzy zmiennymi mającymi charakter losowy zwane
są niekiedy modelami matematycznymi lub zależnościami statystycznymi (stochastycznymi).
Przykładem takiej zależności jest np. związek masy ciała i wysokości ciała Zależności
funkcyjne (deterministyczne) pomiędzy mierzonymi wielkościami są natomiast konsekwencją
podstawowych praw fizyki. Przykładem takiej zależności jest związek wartości indukcji pola
magnetycznego i odległości od środka cewki, indukujacej pole magnetyczne. Element
losowości jest wprowadzany do takich zależności np. przez proces pomiaru.
Zakład Biofizyki CMUJ 24
Elementy matematyki
Wykresy, które reprezentują obrazowo związek pomiędzy zmiennymi, nazywane są
wykresami rozrzutu. Załóżmy, że próba jest badana ze względu na dwie zmienne X i Y. W
prostokątnym układzie współrzędnych na osi odciętych zaznaczamy zmienną niezależną X, a
na osi rzędnych wartości zmiennej zależnej Y. Punkty, odpowiadające poszczególnym
wartościom cech, tworzą korelacyjny wykres rozrzutu. Analizę zależności powinno się
rozpoczynać od jego sporządzenia. Wzrokowa ocena umożliwia często określenie siły i
rodzaju zależności. Rysunek 21 przedstawia wykresy rozrzutu w przypadku gdy między
zmiennymi występuje:
1. korelacja liniowa dodatnia;
2. korelacja liniowa ujemna;
3. brak korelacji;
4. korelacja krzywoliniowa.
Rzadko zdarza się, że zaznaczone punkty leżą dokładnie na linii prostej (pełna
korelacja), częściej spotykana konfiguracja składa się z wielu zaznaczonych punktów
leżących mniej więcej wzdłuż konkretnej krzywej. Gdy korelacja staje się coraz mniej
doskonała, wówczas punkty zaczynają się rozpraszać na płaszczyz
nie XY. Taka sytuacja
występuje w przypadku [3] na rysunku 21. Korelacja dodatnia występuje wtedy, gdy
wzrostowi wartości jednej cechy odpowiada wzrost średnich wartości drugiej cechy
(przypadek [1] na rysunku 21). Korelacja ujemna występuje wtedy, gdy wzrostowi jednej
cechy odpowiada spadek średnich wartości drugiej cechy (przypadek [2] na rysunku 21).
Natężenie współzależności dwóch zmiennych można wyrazić liczbowo za pomocą
współczynnika korelacji liniowej Pearsona. Współczynnik ten (oznaczony symbolem r i
xy
przyjmujący wartości z przedziału [-1, 1] jest miernikiem siły związku prostoliniowego
między dwiema cechami mierzalnymi. Wzór na współczynnik korelacji liniowej Pearsona ma
postać:
n
- x)(yi - y)
"(xi
i=1
rxy =
,
n n
- x) - y)
"(xi "(yi
i=1 i=1
x,y
gdzie - to średnie wartości cech x i y.
2
1
Zakład Biofizyki CMUJ 25
Elementy matematyki
3
4
Rys. 21 Wykresy rozrzutu dla przypadku korelacji liniowej dodatniej (1),
korelacji liniowej ujemnej (2), braku korelacji (3) i korelacji
krzywoliniowej (4).
Znak współczynnika korelacji informuje nas o kierunku korelacji, natomiast jego
bezwzględna wartość  o sile związku. Mamy oczywiście równość r = r . Gdy |r | jest
xy yx xy
bliskie jedności, to zależność korelacyjna przechodzi w zależność funkcyjną (funkcja
liniowa). Jeśli natomiast |r | jest bliskie zeru, oznacza to brak związku korelacyjnego między
xy
badanymi zmiennymi X i Y (przypadek [3] na rysunku 21).
Funkcja regresji to analityczny wyraz przyporządkowania średnich wartości zmiennej
zależnej Y konkretnym wartościom zmiennej niezależnej X. Najprostsze zależności między
zmiennymi to te, które mają postać liniową. Krzywe regresji będące liniami prostymi
nazywamy prostymi regresji. Mają one następującą postać:
y = ax+ b
.
Parametry równania regresji z próby szacuje się metodą najmniejszych kwadratów.
Metoda najmniejszych kwadratów polega na takim oszacowaniu parametrów wyżej
wymienionych funkcji, by dla danych z próby spełniony był warunek :
n n
2 2
wyrażenie - yi ) = - ax- b)
ma osiągnąć minimum,
"(y "(y
i i
i=1 i=1
yi
gdzie y oznaczają wartości empiryczne zmiennej Y, a wartości teoretyczne wyznaczone
i
y = ax+b
na podstawie równania .
Wartości współczynników a i b wyznacza się z wzorów:
n
"(x - x)(yi - y)
i
i=1
a =
n
2
- x)
"(xi
i=1
b = y -ax
.
Można udowodnić, że błędy  odchylenia standardowe wartości średnich wielkości a i b
wyrażają się wzorami :
Zakład Biofizyki CMUJ 26
Elementy matematyki
n
1 n
2
Ãa = µi
"
2
n n
n- 2
i=1
ëÅ‚ öÅ‚
2
n xi - ìÅ‚ xi ÷Å‚
" "
i=1 íÅ‚ i=1 Å‚Å‚
n
2
xi
"
n
1
2 i=1
à = µi ,
"
b
2
n n
n- 2
i=1
ëÅ‚ öÅ‚
2
n xi - ìÅ‚ xi ÷Å‚
" "
i=1 íÅ‚ i=1 Å‚Å‚
w których
n n n n
2 2
"µ = "y - a"x yi - b"y .
i i i i
i=1 i=1 i=1 i=1
Przykład 7
Wykonano pomiar gęstości kręgu lędz
wiowego L3 w grupie 16 osób przy pomocy dwóch
urządzeń - densytometru i tomografu rentgenowskiego. Wyniki przedstawiono na rysunku 22.
180
160
140
120
Y = a * X + b
100
80 a = 162 Ä… 34
b = -67 Ä…34
60
40
rxy2 = 0.783
20
0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5
BMD (g/cm2)
Rys. 22 Wykres zależności wyników pomiarów gęstości kości,
wykonanych przy pomocy różnych metod diagnostycznych.
Zakład Biofizyki CMUJ 27
3
BMD (mgHA/cm )
Elementy matematyki
Zakład Biofizyki CMUJ 28