Politechnika Wrocławska

Wydział Elektryczny

Wojciech Calów

Rok studiów : III
Semestr : VI
Rok akad. : 2012/13

 Laboratorium Metod Numerycznych

Data: 

Temat : 

Metody rozwiązywania układu równań 

macierzowo

Ocena: 

1. Program

%  metoda LU

clear 

all

;

A=[5 -2 1 -1; 1 6 -2 3 ; 1 -3 4 -2 ; 2 1 -5 4 ];
b=[4 3 12 10];

x=gauss(A,b);

n = size(A,1);
    

for

 k=1:n

        L(k,k)=1;
    

end

 sum2=0;
 sum1=0;

 for

 k=1:n 

% ================================= 

      

    L=0;
    U=0;

    

for

 i=k:n 

% liczenie U =========

        

for

 m=1:(k-1)

            sum1 = sum1 + L(k,m) * U(m,i);            
        

end

 

        U(k,i)=A(k,i) - sum1;
    

end

   
    

for

 j=(k+1):n 

%liczenie L ========

        

for

 m=1:(k-1)

            sum2=sum2 + (L(j,m) * U(m,k));

        

end

        L(j,k) = (1/u(k,k))* (A(j,k)- sum2);        

    

end

 

% ===================================       

end

;

2. Wynik

x =

    0.5663
   -2.7831
    5.9398
   10.3373