ReWiA 2.0 html manual
FIT dopasowywanie profili
Zobacz też : Dopasowywanie profili Gaussa i Lorentza
Zapis FIT
Komenda posiada własne menu i wyświetla własne parametry lokalne. Nie
akceptuje żadnych parametrów.
Opcje
Assoc
wybranie lub zmiana długości widma zawierającego tabelę profili. W
przypadku gdy widmo było puste zostanie utworzone widmo o 8 polach i
podanej długości w następującym formacie :
pole typ opis
0 strg*1 rodzaj linii, 'G' lub 'g' oznacza Gauss, 'L' lub 'l' oznacza
Lorentz, 'I' lub 'i' oznacza Gauss całkowy. Każda inna
wartość jest niedozwolona (EPTR=0).
1 real*8 położenie środka profilu x0 (XPTR=1). Przed dopasowaniem
należy wpisać przybliżoną wartość.
2 real*4 amplituda profilu A (YPTR=2). Przed dopasowaniem należy
wpisać przybliżoną wartość.
3 real*4 dyspersja profilu D. Przed dopasowaniem należy wpisać
przybliżoną wartość.
4 real*4 szerokość połówkowa (FWHM)
5 real*4 Pole profilu. Wartość ujemna oznacza linię absorpcyjną.
6 real*4 Szerokość równoważna. Wartość ujemna oznacza linię
absorpcyjną.
7 real*4 Logarytm wartości bezwzględnej szerokości równoważnej
Długość widma zawierającego tabelę profili może być większa od liczby
dopasowywanych profili. Z reguły zawiera ona informację o wszystkich
dopasowywanych profilach. Obecna wersja ReWii umożliwia dopasowywanie do
10 profili jednocześnie. Profile są określone następująco :
(x - x0)^2
Gauss(x) = A * exp(- ----------)
D^2
(x - x0)^2
Lorentz(x) = A / (1 + ----------)
D^2
CałkowyGauss(x) = erf(x + 0.5) - erf(x - 0.5)
erf(x) = A * całka (exp(-(((u-x0)/D)^2)), du) {-infty, x}
gdzie : x0 - położenie środka profilu (pole 1), A - głębokość profilu
(pole 2), D - dyspersja profilu (pole 3). Pola profili są obliczane
następująco :
area(Gauss) = area(CałkowyGauss) = A * D * sqrt(pi)
area(Lorentz) = A * D * pi
Gauss całkowy uwzględnia fakt, że punkt widma reprezentuje całkę z
wąskiego zakresu długości fali, i jest szczególnie przydatny przy
wąskich liniach otrzymywanych z CCD lub reticon'a. W porównaniu do
zwykłego Gaussa daje bardziej ostre profile. Niestety ma problemy ze
zbieżnością i warunkami początkowymi dopasowania.
Contin
ustalenie lokalnego kontinuum (linia prosta). Można je wprowadzić albo
podając numery 2 punktów przez które kontinuum ma przechodzić albo
podając bezpośrednio współczynniki prostej. Jeśli w opcji Options tak
ustawiono to przy wprowadzaniu kontinuum metodą 2 punktów te dwa punkty
automatycznie określają dopasowywany fragment widma.
Edit
edycja tabeli profili w trybie tekstowym. Umożliwia oglądanie wyników i
zmianę parametrów początkowych dopasowywanych linii.
Fit
dopasowanie profili do fragmentu widma (ustalonego opcją lImits).
Uwzględniane są tylko te punkty w wybranym fragmencie dla których E<>0.
Po zakończeniu wprowadzania parametrów linii opcją X należy przywrócić
właściwą linię bieżącą. Dopasowywanych będzie n linii (ustalone opcją
Lines) począwszy od linii bieżącej (parametry jej widoczne są w oknie
parametrów lokalnych).
lImits
ustalenie fragmentu widma do którego należy dopasować profile. Jeśli w
opcji Options tak ustawiono to przy wprowadzaniu kontinuum metodą 2
punktów te dwa punkty automatycznie określają dopasowywany fragment
widma.
Lines
wybór ilości dopasowywanych linii. Dopasowywane są linie od bieżącej
począwszy.
Mark
wybór bieżącej linii (marki). Można używać również klawiszy '-' i '+' po
prawej stronie klawiatury (szary kolor). UWAGA : po zakończeniu
wprowadzania parametrów linii opcją X należy przywrócić właściwą linię
bieżącą (pierwszą dopasowywaną linię).
Option
ustalenie sposobu działania opcji X (czy po wprowadzeniu położenia
środka profilu automatycnie obliczać i wstawiać głębokość, dyspersję
oraz kształt Gauss) oraz parametrów algorytmu dopasowującego : limitu
iteracji, limitu dokładności oraz współczynnika linearyzacji FLAM.
Algorytm zaczerpnięto z książki : Bevington "Data Reduction and Error
Analysis", procedura CURFIT. Inne komendy np. CURFIT, ENER, CALIB
również korzystają z tego algorytmu.
Quit
zakończenie komendy FIT
Subtr
odjęcie od widma wybranej linii w zakresie ustalonym ocją lImits.
View
wyświetlenie dopasowanych profili na tle widma. Można wyświetlić
wszystkie profile na raz, pojedynczy profil lub tylko sumę wszystkich
profili.
X
wprowadzenie przybliżonej pozycji środka profilu. Dodatkowo, o ile
ustawiono to w opcji Options, zostają automatycznie wprowadzone
początkowa głębokość i dyspersja linii, kształt Gauss a następnie numer
bieżącej linii jest zwiększany i opcja pyta się o pozycję następnej
linii. Tak więc inicjalizacji całej tabeli profili dla 2 linii,
dopasowania i wyświetlenia wyników dokonuje następująca sekwencja
klawiszy :
A 10
20
C 1
klawisze kursora Alt-spacja
klawisze kursora Alt-spacja
X
klawisze kursora Alt-X
klawisze kursora Alt-X
L 2
M 0
F V 1
Jeśli zablokowane jest automatyczne wprowadzanie głębokości, dyspersji i
kształtu profilu, to opcja X zapyta się o każdy w tych parametrów.
Naciśnięcie klawisza ESC przerywa wykonywanie się opcji X.