Regularyzacja Tichonowa

Regularyzacja Tichonowa – metoda regularyzacji zagadnień nie postawionych poprawnie opracowana niezależnie przez Andrieja Tichonowa i Davida Phillipsa, jednak nazwana od nazwiska rosyjskiego matematyka.

Regularyzacja za pomocą parametru regularyzacyjnego

Rozważmy układ równań typu:

$A𝐱=𝐛,$

gdzie $A$ to macierz układu (macierz przekształcenia), $𝐱$ wektor niewiadomych, $𝐛$ wektor wyjściowy. Jeśli zagadnienie to jest źle postawione (np. rozwiązanie nie istnieje, lub jest niejednoznaczne) wtedy standardowym podejściem jest rozwiązanie w sensie najmniejszych kwadratów. Oznacza to poszukiwanie minimum następującego funkcjonału $\mathrm{\Pi }$:

$\mathrm{\Pi }=\{A𝐱-𝐛{\}}^T\{A𝐱-𝐛\}.$

Aby poprawić właściwości funkcjonału $\mathrm{\Pi }$ (np. żeby wyeliminować niefizyczne oscylacje w rozwiązaniu) do powyższej postaci dodaje się człon regularyzacyjny $\alpha {𝐱}^T𝐱:$

$\mathrm{\Pi }=\{A𝐱-𝐛{\}}^T\{A𝐱-𝐛\}+\alpha {𝐱}^T𝐱.$

Wyrażenie ${𝐱}^T𝐱$ to kwadrat długości wektora $𝐱,$ zaś $\alpha$ to parametr regularyzacyjny zwany również współczynnikiem tłumienia.

Szukanie minimum funkcjonału $\mathrm{\Pi }$ oznacza rozwiązanie równania:

$\frac{\partial \mathrm{\Pi }}{\partial x}=0,$

co daje:

$-2A^T\{𝐛-A𝐱\}+2\alpha 𝐱=0$

i po przekształceniu:

$A^T𝐛=[A^{T}A+\alpha 𝐈]𝐱,$

gdzie $𝐈$ to macierz jednostkowa.

Macierz stojąca po prawej stronie powyższego równania jest macierzą symetryczną, ponadto dzięki obecności członu $\alpha 𝐈$ jest również odwracalna. Poza tym parametr regularyzacyjny wpływa na określoność macierzy i poprawia jej uwarunkowanie. W związku z tym rozwiązanie wyjściowego równania może być wyznaczone za pomocą:

$𝐱={[A^{T}A+\alpha 𝐈]}^{-1}{A}^T𝐛.$

Dobór parametru regularyzacyjnego może być wykonany za pomocą tzw. metody L-curve.

Regularyzacja za pomocą macierzy regularyzacji

Jeśli dana jest informacja a priori na temat wektora niewiadomych $𝐱$ może ona zostać uwzględniona w funkcjonale $\mathrm{\Pi }$ w następujący sposób:

$\mathrm{\Pi }=\{A𝐱-𝐛{\}}^T\{A𝐱-𝐛\}+\alpha {\left[R𝐱\right]}^T\left[R𝐱\right],$

gdzie $R$ to macierz regularyzacji oraz $\alpha$ kontroluje stopień regularyzacji.

Pochodna funkcjonału $\mathrm{\Pi }$ względem wektora $𝐱$ wynosi:

$\frac{\partial \mathrm{\Pi }}{\partial x}=-2A^T\{𝐛-A𝐱\}+2\alpha R^{T}R𝐱.$

Co daje:

$A^T𝐛=[A^{T}A+\alpha R^{T}R]𝐱.$

Ponownie macierz stojąca po prawej stronie jest symetryczna i w większości przypadków odwaracalna. Postać macierzy regularyzacji zależy od rodzaju informacji, który jest wprowadzany do zagadnienia. Na przykład jeśli rozwiązaniem zagadnienia powinny być w przybliżeniu stałe wartości i-ty wiersz macierzy $R$ powinien wynosić:

$r_i=[0,\dots ,0,\underset{i-1}{\underbrace{-1}},\underset{i}{\underbrace{1}},0,\dots ,0].$

Inne przykłady sformułowania macierzy $R$ w zależności od informacji, która ma być w niej zawarta można odnaleźć w literaturze.

Przykład

Przykład zastosowania regularyzacji Tichonowa. Szczegóły w tekście.

Rozpatrzmy następujący przykład:

$A=\left[\begin{array}{cccc}1& 2& 3& 4\\ 1& 3& 1& 5\\ 1,99999999& 3,99999999& 6,00000001& 8,0000001\\ 7& 5& 3& 1\end{array}\right],$

$𝐱=\left[\begin{array}{c}1\\ 2\\ 2\\ 1\end{array}\right].$

Wtedy $𝐛$ wynosi:

$𝐛=A𝐱\approx \left[\begin{array}{c}15\\ 14\\ 30,00000009\\ 24\end{array}\right].$

Załóżmy teraz, że na podstawie $A$ i $𝐛$ chcemy wyznaczyć $𝐱.$ Ze względu na fakt, że wyznacznik macierzy jest bliski zera $(\det (A)=3,6\cdot 10^{-6})$ i wskaźnik uwarunkowania macierzy wynosi w przybliżeniu 909460299 wyznaczenie macierzy odwrotnej może nastręczać kłopotów, dlatego rozwiązanie poszukiwane jest za pomocą metody najmniejszych kwadratów. Na wykresie obok zielonym kolorem naniesiono przykładowy wynik rozwiązania tą metodą. Widać duże odstępstwa od założonych wcześniej wartości $𝐱$ oznaczonych za pomocą niebieskich kółek. Dodanie członu regularyzacyjnego zgodnie ze schematem przedstawionym we wcześniejszej części artykułu zdecydowanie poprawia uzyskane rozwiązanie (czerwona linia na wykresie obok). Przy czym parametr regularyzacyjny w tym przypadku wynosi $\alpha =0,00000007.$ Poniżej znajduje się kod źródłowy w pakiecie MATLAB z omawianym przykładem:

a = 0.00000007;% parametr regularyzacyjny
%macierz A
A=[1 2 3 4; 1 3 1 5; 1.99999999 3.99999999 6.00000001 8.0000001;7 5 3 1];
%wektor x
x = [1;2;2;1];
% wektor b
b = A*x;
% uchwyt do funkcji bez regularyzacji
fun = @(xx)(A*xx-b)*(A*xx-b)';
% uchwyt do funkcji z regularyzacją
fun2 = @(xx)(A*xx-b)*(A*xx-b)'+a*xx*xx';
%rozwiązanie metodą najmniejszych kwadratów
x1 = lsqnonlin(fun,0*x);
x2 = lsqnonlin(fun2,0*x);
hold all
plot(x,'o')
plot(x1)
plot(x2)
grid on
legend('szukane wartości','rozwiązanie bez regularyzacji',...
'rozwiązanie z regularyzacją')

Zobacz też

• problem odwrotny
• regularyzacja funkcją dzeta

Linki zewnętrzne

• Zastosowanie metod numerycznych