Metoda Potęgowa

c=sqrt(dot(x0,x0)); %cbliczenie długości wektora xQ u0=x0/c; %normalizacja wektora przybliżeń it= ; %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 do pętli eu=l ;

Iambda0= ; %przypisanie początkowej wartości lambdaO, które będzie słyżyło jako lambda z pcprzeniej iteracji

Rwhile el>e && eu>e %pętla iteracyjna sprawdzająca warunki zbieżności dla lambda i u

x=A*uO; %cbliczanie wektora własnego macierzy A Iambda=u0’*x; %obliczanie m.ax wartości własnej macierzy A c=sqrt(dot(x,x)); %obliczenie długości wektora x u=x/c; %ncrm.alizacja wektora własnego

el=abs(lambda-lambdaO)/lambda; %obiiczanie błędu względnego z lambdy (warunek zbieżności dla lambda) eu=sqrt(dot((u-uO),(u-uO))); %obliczanie normy dla różnicy wektorów u i uO (warunek zbieżności dla lambda) it=it+l; %zwiekszenie licznika iteracji

H if it>2C %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

end

Wyszukiwarka

Podobne podstrony:
Metoda Iteracji Odwrotnej function [u,lambda,it] =10 (A, xQ, e) %funkcja pobiera macierz A początkow
Metoda Gaussa Jordana function x=gaussjordan(A,b) %funJccja pobiera macierz A i wektor wyrazów wolny
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
6 „Budowanie społeczeństwa informacyjnego” metoda musi zakładać wykorzystanie IT do rozwiązania
Metoda Choleskiego Banachiewicza function x=ChB(A,b) %funkcja pobiera macierz A i wektor wyrazów wol
Metoda Gradientów Sprzężonych function x=gradienty(A,b,xQ,e) %funkcja pobiera macierz A, wektor b, x
DSC04631 •AZS •apm •AKZS p it^mp MApm) P*4**PZ P rjmnomp P
img083 83 6.5. Metoda funkcji nieliniowych dla każdego (dowolnie małego) e. Jak z tego wynika, począ
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
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
functionblock Usc thc pulldown menu :it the top of thc scrccn to sclcct a diffcrenł function bł

więcej podobnych podstron