22.10.2013 Views

Rešene naloge iz numerične matematike

Rešene naloge iz numerične matematike

Rešene naloge iz numerične matematike

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Rešitev.<br />

% Naloga1.m<br />

clear all<br />

a = diag([4,-2,2i,0]-9) + ones(4)*9<br />

b = [ones(3), zeros(3,2); zeros(2,3), 4*eye(2)]<br />

c = zeros(5);<br />

c(1,:) = ones(1,5) * 2;<br />

c(:,1) = ones(5,1) * 2;<br />

c(3:4,3:4) = 7*ones(2);<br />

c(5,5) = 10;<br />

c<br />

2. V Octave vnesite naslednje matrike:<br />

⎡<br />

⎤<br />

1 1 1 1<br />

⎢<br />

⎢2<br />

1 1 1 ⎥<br />

A1 = ⎢<br />

⎥<br />

⎣1<br />

1 1 3⎦<br />

1 1 1 1<br />

, A2<br />

⎡<br />

⎤<br />

1 + i 1 1 − i <br />

⎢<br />

⎥ A1 I<br />

= ⎣ 2i 1 − 3i 1 ⎦ , A3 = .<br />

0 2A1<br />

4 − i 1 + 2i i<br />

Izpišite tretjo vrstico in prvi stolpec matrike A1. V matriki A3 vse<br />

ničle v spodnjem levem kotu spremenite v 2. Izračunajte determinanto<br />

in lastne vrednosti nove matrike A3. Poiščite maksimalni element po<br />

absolutni vrednosti v matriki A2. Naj bo x drugi, y pa tretji stolpec<br />

matrike A2. Izračunajte skalarni produkt vektorjev x in y.<br />

Rešitev.<br />

% Naloga2.m<br />

clear all<br />

a1 = ones(4);<br />

a1(2,1) = 2;<br />

a1(3,4) = 3<br />

a2 = [1+i,1,1-i;2i,1-3i,1;4-i,1+2i,i]<br />

a3 = [a1, eye(4); zeros(4), 2*a1]<br />

a1(3,:) % tretja vrstica matrike a1<br />

a1(:,1) % prvi stolpec matrike a1<br />

a3(5:end,1:4) = ones(4)*2 % spremenimo matriko a3<br />

det(a3) % determinanta<br />

eig(a3) % lastne vrednosti<br />

max(max(abs(a2))) % maksimalni element po abs. vrednosti<br />

x = a2(:,2);<br />

14

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!