Metoda Iteracji Odwrotnej
function [u,lambda,it] =10 (A, xQ, e) %funkcja pobiera macierz A początkowy wektor przybliżeń x0 i ograniczenie dokładności rozwiązania(u i lambda) e
A=AA(- ); %tworzenie macierzy odwrotnej
c=sqrt(dot(x0,x0)); %cbliczenie długości wektora xQ
u0=x0/c; %norrr.alizacja wektora przybliżeń
it=G; %początkowa wartość licznika iteracji
el=l; %przypisanie początkowych wartości warunku zbieżności procesu dla lambda (el) i u (eu) gwarantujące wejście dc pętli eu=l;
lambda0= ; %przypisanie początkowej wartości lambdaO, które będzie słyżyło jako lambda z poprzeniej iteracji
Hwhile el>e && eu>e %pętla iteracyjna sprawdzająca warunki zbieżności dla lambda i u
x=A*uQ; %obliczanie wektora własnego macierzy A lambda=u0'*x; %obliczanie max wartości własnej macierzy A c=sqrt(dot(x,x)); %obliczenie długości wektora x u=x/c; %normalizacja wektora własnego
el=abs(lambda-lambdaO)/lambda; %obliczanie błędu względnego z lambdy (warunek zbieżności dla lambda) eu=sqrt(dot((u-uO),(u-uO))); %cbliczanie normy dla różnicy wektorów u i uO (warunek zbieżności dla lambda) it=it+L; %zwiek3zenie licznika iteracji
H if it>20C %sprawdzenie ograniczenia liczby iteracji disp('proces rozbieżny’); break; %przerwanie działania pętli
end
lambdaO=lambda; %przypisanie wartości lambda obliczonej w danej iteracji jako wartość z iteracji poprzedniej u0=u; %przypisanie wektora u obliczonego w danej iteracji jako wektora z iteracji poprzedniej
end
lambda=l/lambda; %lambda dla macierzy odwrotnej (ze wzoru) x=u; end
Wyszukiwarka
Podobne podstrony:
Metoda Potęgowa function [u,lambda,it]=MP(A,x0,e) %funJccja pobiera macierz A początkowy wektor przy
Metoda Gradientów Sprzężonych function x=gradienty(A,b,xQ,e) %funkcja pobiera macierz A, wektor b, x
Metoda Choleskiego Banachiewicza function x=ChB(A,b) %funkcja pobiera macierz A i wektor wyrazów wol
Metoda Gaussa Seidla (function x=GS(A,b,xO,e) %funkcja pobiera macierz A, wektor b, początkowy wekto
Metoda Jacobiego (function x=jacob(A,b,xO,e) %funkcja pobiera macierz A, wektor b, początkowy wektor
Podstawianie W Przód 1 2 3 4 5 6 7 8 9 10 11 function x=ForwardS (A,b) %funkcja pobiera rr.acierz A
CCF20121213 011 3. Metoda Obszarów Stetsona (WISC - R) • Wyróżnione obszary (10 sfer funkcjonowania
Eliminacja Gaussa function x=gauss(A,b) %funkcja pobiera macierz A i wektor wyrazów wolnych b, a zwr
Podstawianie Wstecz function x=3ackS(A,b) %funkcja pobiera macierz A i wektor wyrazów wolnych b, a z
iteracja Metoda iteracji. W przepadku gdy liczba niewiadomych układu równań jest duża, rozwiązanie t
gausa siedla Metoda Gaussa - Seidela jest metodą iteracyjną i pozwala nam obliczyć układ n równań z
2614. Betony przeciętnych i wysokich wytrzymałości projektowane metodą iteracji.
81746 vanna christmas traditions SnowflakesMATERIALS I >-15 yds whitc bcdspread weight cottou i
więcej podobnych podstron