PL EN DE RU UA
Medycyna nuklearna

Kalkulator

(wcześniej dostępny na www.nuk.bieganski.org)

Opcje ogólnomatematyczne i ogólnomedyczne

Średnia arytmetyczna i odchylenie standardowe
Regresja punktów do linii prostej, współczynnik korelacji
Znajdowanie naturalnych dzielników
Największy wspólny dzielnik
Wielomiany
Obliczanie różnicy czasu
Szacowanie niektórych parametrów biometrycznych człowieka
Przeliczanie stężeń
Farmakokinetyczne modele kompartmentowe

Obliczenia związane z fizyką nuklearną

Przykłady zastosowań poniższych programów znajdują się tutaj.
Rozpad promieniotwórczy prosty (1 nuklid promieniotwórczy: A → B)
Rozpad promieniotwórczy seryjny (2 nuklidy promieniotwórcze: A → B → C)
Rozpad promieniotwórczy seryjny (3 nuklidy promieniotwórcze: A → B → C → D)
Przeliczanie jednostek aktywności (tradycyjnych na pochodne układu SI i odwrotnie)
Przeliczanie masy na aktywność i odwrotnie (jednostki masy na jednostki aktywności i odwrotnie)

Obliczenia związane z medycyną nuklearną

Instrukcja do poniższych programów znajduje się tutaj.
(1.) Obliczanie objętości tarczycy na podstawie wymiarów płatów
(2.) Obliczanie jednorazowego wychwytu (radio)izotopu (na przykład radiojodu w tarczycy)
 (2a.) Obliczanie wychwytu radiojodu w tarczycy (Program można nagrać jako plik i uruchomić na innym komputerze, konieczna przeglądarka obsługująca HTML i JavaScript).
(3.) Wyliczanie dawki jodu promieniotwórczego (uproszczone)
(4.) Modelowanie kinetyczne I. (efektywny czas połowicznego zaniku, maksymalny wychwyt i in. - na podstawie serii pomiarów) - do leczenia izotopowego
(5.) Modelowanie kinetyczne II. (efektywny czas połowicznego zaniku, maksymalny wychwyt i in. - na podstawie trzech pomiarów) - do leczenia izotopowego
(6.) Obliczanie dawki jodu promieniotwórczego lub innego radionuklidu do leczenia (zmodyfikowany wzór Marinelli)
(7.) Dozymetria promieniowania α i β...
(8.) Natężenie i dawka promieniowania γ (w określonej odległości od źródła punktowego)

Wielomiany

Obliczanie miejsc zerowych równania liniowego, kwadratowego, sześciennego oraz równania czwartego (kwartycznego) i piątego stopnia (kwintycznego):
0= a5·x5+a4·x4+a3·x3+a2·x2+a1·x1+a0·x0
Podaj współczynniki wielomianu (a5, a4, a3, a2, a1, a0):
a5a4a3a2a1a0



Algorytmy

Algorytmy stosowane do obliczeń miejsc zerowych wielomianów stopni od 0 do 5 o wzorze ogólnym:
0= a5·x5+a4·x4+a3·x3+a2·x2+a1·x1+a0·x0

Rówanie liniowe 0= a·x+b

Miejsce zerowe: x0= -b/a.

Równanie kwadratowe 0= a·x2+b·x+c:

Wyznaczanie miejsc zerowych:
1. Obliczenie wyróżnika Δ= b2-4·a·c.
2. Trzy warianty:
- dla Δ>0: dwa pierwiastki rzeczywiste różne,
- dla Δ=0: pierwiastek podwójny (dwa równe),
- dla Δ<0: dwa pierwiastki zespolone.
3. Wyznaczanie wierzchołka funkcji (ekstremum) z pochodnej (2·a·x+b): x= -b/(2·a), y= -Δ/(4·a).
4. Wyznaczanie położenia ogniska (x= -b/(2·a), y= -Δ/(4·a)+1/(4·a)) oraz kierownicy (y= -Δ/(4·a)-1/(4·a)) paraboli.

Równanie sześcienne 0= a·x3+b·x2+c·x+d

1. Obliczanie wyróżnika Δ3= -27·a2·d2+18·a·b·c·d-4·a·c3-4·b3·d+b2·c2,
a także parametrów "postaci kanonicznej": 0= y3+p·x+q (po podstawieniu x=y-b/(3·a)):
Δ3b= Δ3/(-27·a4),
p= c/a-b2/(3·a2),
q= 2·b3/(27·a2)-b·c/(3·a2)+d/a,
przy czym: Δ3b= q2+4/27·p3.
2. Trzy warianty:
 I. Δ3<0 (istnieje jeden pierwiastek rzeczywisty r):
  r= {curt[-q/2+sqrt(Δ3b)/2]+curt[-q/2-sqrt(Δ3b)/2]}-b/(3·a)
  (przy czym curt() to pierwiastek sześcienny),
  rozkład na czynniki: 0= (x-r)·(aq·x2+bq·x+cq),
  gdzie: aq= a, bq= b+a·r, cq= c+b·r+a·r2.
 II. Δ3= 0 (istnieje pierwiastek wielokrotny):
  Dwa warianty:
   a) p= 0 i q= 0 (pierwiastek potrójny):
    x1,2,3= curt(d/a),
    rozkład na czynniki: 0= a·(x-x1,2,3)3;
   b) p<0 i q>0 (pierwiastek podwójny):
    x3= 2·curt(-q/2)-b/(3·a); x1,2= -curt(-q/2)-b/(3·a);
    rozkład na czynniki: 0= a·(x-x1,2)2·(x-x3)
 I. Δ3>0 (istnieją trzy różne pierwiastki rzeczywiste):
  Φ= atan2[sqrt(-Δ3b),-q];
  x1, x2, x3 = 2/3·sqrt(-3·p)·cos(Φ/3+2·k·π/3)-b/(3·a) dla k= 0, 1 i 2,
  rozkład na czynniki: 0= a·(x-x1)·(x-x2)·(x-x3)
3. Obliczanie extremów z pochodnej (3·a·x2+2·b·x+c) i punktu przegięcia z drugiej pochodnej (6·a·x+2·b).

Równanie czwartego stopnia 0= a4·x4+a3·x3+a2·x2+a1·x+a0

Wyróżnik: Δ4= 256·a43·a03 -192·a42·a3·a1·a02 -128·a42·a22·a02 +144·a42·a2·a12·a0 -27·a42·a14 +144·a4·a32·a2·a02 -6·a4·a32·a12·a0 -80·a4·a3·a22·a1·a0 +18·a4·a3·a2·a13 +16·a4·a24·a0 -4·a4·a23·a12 -27·a34·a02 +18·a33·a2·a1·a0 -4·a33·a13 -4·a32·a23·a0 +a32·a22·a12
1. Sprowadzanie do postaci unormowanej 0= a4·(x4+a·x3+b·x2+c·x+d),
gdzie a= a3/a4, b= a2/a4, c= a1/a4, d= a0/a4.
Potem następuje rozkład na czynniki (wg algorytmu S. Shmakova, IJPAM, 2011):
0= (x2+g1·x+h1)·(x2+g2·x+h2), gdzie:
  i:    a= g1+g2
  ii:    b= h1+h2+g1·g2
  iii:    c= g1·h2+g2·h1
  iiii:    d= h1·h2
2. Rozwiązanie pośredniego równania sześciennego: 0= y3+(-b)·y2+(a·c-4·d)·y+(4·b·d-c2-a2·d)
Następuje według algorytmu podanego powyżej. Potrzebne jest tylko jedno rozwiązanie, które należy wybrać.
3. Rozwiązywanie pośrednich (pomocniczych) równań kwadratowych G i H:
  G:  0= g2+(-a)·g+(b-y)
  H:  0= h2+(-y)·h+(d)
Następuje według algorytmu podanego powyżej. Otrzymujemy w ten sposób cztery liczby:
g1 i g2 (rozwiązanie równania G) oraz hk i hn (rozwiązanie równania H).
4. Sprawdzenie, w jaki sposób liczby hk oraz hn spełniają równanie iii. W ten sposób wyznaczamy, która z liczb hk oraz hn odpowiada h1 a która h2.
W ten sposób rozłożono równanie czwartego stopnia na czynniki kwadratowe.
5. Znalezienie pierwiastków powyższych równań kwadratowych (według algorytmu równania kwadratowego) jest równoznaczne ze znalezieniem pierwiastków początkowego równania czwartego stopnia.
6. Obliczanie extremów z pochodnej i punktu przegięcia z drugiej pochodnej.

Równanie piątego stopnia 0= a5·x5+a4·x4+a3·x3+a2·x2+a1·x+a0

1. Sprowadzanie do postaci unormowanej 0= a5·(x5+A·x4+B·x3+C·x2+D·x+E),
gdzie A= a4/a5, B= a3/a5, C= a2/a5, D= a1/a5, E= a0/a5.
2. Wyznaczanie przedziału, w którym znajdują się wszystkie pierwiastki, za pomocą twierdzenia Cauchy'ego:
  ρ2= 1+max(1,A,B,C,D,E).
  Wszystkie pierwiastki muszą być zawarte w przedziale od -ρ2 do +ρ2;
  ze względu na to, że mamy do czynienia z wielomianem unormowanym piątego stopnia, wartość funkcji dla -ρ2 musi być ujemna, a dla +ρ2 musi być dodatnia.
3. Iteracyjne przybliżanie wartości jednego z pierwiastków (x5) za pomocą metody siecznych (bisection method):
  3.1. Założenia wstępne: x51= ρ2, x52= -ρ2.
  3.2. Oliczenie x5: x5= (x51+x52)/2.
  3.3. Sprawdzanie wartości funkcji wielomianu unormowanego piątego stopnia dla znalezionego x5:
   3.3.1. Wartość funkcji = 0 →: opuszczenie pętli (znaleziono miejsce zerowe);
   3.3.2. Wartość funkcji > 0 →: x5→x51, idź do punktu 3.2.;
   3.3.3. Wartość funkcji < 0 →: x5→x52, idź do punktu 3.2.;
  w ten sposób jeden z pierwiastków jest stopniowo przybliżany aż do osiągnięcia satysfakcjonującej dokładności.
 Procedurę można uprościć, jeżeli E= 0. Wtedy x5= 0, zaś a= A, b= B, c= C oraz d= D.
4. Po wyznaczeniu x5 następuje rozłożenie wielomianu piątego stopnia na dwumian pierwszego stopnia i wielomian czwartego stopnia:
  0= a5·(x-x5)·(x4+a·x3+b·x2+c·x+d), gdzie:
  a= A+x5,
  b= B+A·x5+x52,
  c= C+B·x5+A·x52+x53,
  d= D+C·x5+B·x52+A·x53+x54;
5. Rozwiązanie wielomianu czwartego stopnia (celem rozłożenia na czynniki kwadratowe i ew. znalezienia kolejnych pierwiastków) następuje według schematu podanego powyżej.
6. Obliczanie extremów z pochodnej i punktu przegięcia z drugiej pochodnej.


©Autor: Cyprian Świętaszczyk, 2013; Ostatnia aktualizacja: 07.2020;