Budownictwo 2006 Wydział Inżynierii Lądowej Politechniki Krakowskiej

rzabaa

jaki byla tresc zadania?? wszyscy mieli to samo???

Szczur

ma ktoś może już to napisane?

ewelka87

może by ktoś napisał już coś co było?

Szczur

trzeba było obliczyć wektor i wartości własne macierzy 3 stopnia, za pomocą metody potęgowej albo metody iteracji odwrotnej.

metoda potęgowa:

Woodchack

Na mailu naszym gr.4 jest ten program poprawiony działający w pełni...chwila to wkleje...

A=[-2 5 9;4 -13 11;1 0 -18]

function [iter,lambda,x]=pys(A)

n=100;
xk=[1;0;0];

lamk=0;
c=(xk'*xk)^(1/2);
uk=xk/c;
epsL=1;
epsR=1;
nr_iteracji=1;

while(epsL>10^-6)&(epsR>10^-6)

xk1=A*uk;
lamk1=uk'*xk1;
c=[xk1'*xk1]^(1/2);
uk1=xk1/c;
epsL=abs((lamk1-lamk)/lamk1);
epsR=((uk1-uk)'*(uk1-uk))^(1/2);
nr_iteracji=nr_iteracji+1;

if (nr_iteracji>n)
printf("proces rozbiezny");
return;
end

lamk=lamk1;
uk=uk1;

end
iter=nr_iteracji;
lambda=lamk1;
x=uk1;

endfunction

[iter,lambda,x]=pys(A)

Prosze Very Happy