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