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
Przypisy
Linki zewnętrzne Zastosowanie metod numerycznych
Kategoria:Układy równań liniowych