Questo metodo risulta molto efficace per il suo effetto di ``regolarizzazione''
della soluzione. Ad ogni passo viene trovata un'approssimazione che si avvicina
sempre più a quella ideale.
Il processo iterativo è composto dai seguenti passi:
(2.12)
dove:
-
è la trasposta della matrice
-
è l'operatore che garantisce il vincolo di positività, cioè
la proiezione convessa nello spazio chiuso e convesso di vettori non-negativi
definito da:
(2.13)
-
è il parametro di rilassamento che varia tra i valori:
(2.14)
dove è il più grande valore singolare della matrice
(Appendice A).
Per il nostro tipo di problema deve essere minore di 0.125
(nel codice si pone ).