Sprawozdanie

Dawid Kasiński

Grupa T11A

28.11.2011r

Przedmiot: Technika obliczeniowa i symulacja (Labolatoria)

Temat: Temat 4

Zadanie 1. Algorytm do metody Jacobiego, program do rozbieżności.

licz=0

for k=1:100

N=4

a=0.5

mA=rand(N)

mA=mA*2-1

vb=rand(N,1)

vb=mA'*vb

mA=mA'*mA+a*eye(N)

diaA=diag(diag(mA))

invDiaA=diag(1./diag(mA))

mA_d=mA-diaA

x=vb

y=[ ]

for i=1:50

x=(invDiaA*(vb-mA_d*x));

y=[y,x]

if i>10,

dy=max(abs(y(:,end)-y(:,end-1)))

end;

end

if dy<0.00001 licz=licz+1; end

figure(1)

plot(y')

[V,D]=eig(mA)

end

licz

Zadanie 2. Algorytm do metody Gausa Saidla, program do rozbieżności.

licz=0

for k=1:100

N=4

mA=rand(N)

mA=mA*2-1

vb=rand(N,1)

vb=mA'*vb

mA=mA'*mA

diaA=diag(diag(mA))

invdiaA=diag(1./diag(mA))

mA_d=mA-diaA

x=vb

y=[ ]

for j=1:100

for i=1:N

x(i)=invdiaA(i,i)*(vb(i,1)-(mA_d(i,:))*x)

end;

if j>10

dy=max(abs(y(:,end)-y(:,end-1)))

end

dy

y=[y,x];

end

if dy<0.1 licz=licz+1; end

figure(2)

plot(y')

[V,D]=eig(mA)

end

licz