Jacobi Verfahren < Numerik < Hochschule < Mathe < Vorhilfe
|
Aufgabe | Umsetzung des Jacobi Verfahrens in Matlab |
geg:Ax=b
A (2250,2250);
b (2250,1);
D = diag(diag(A));
[L,U]=lu(A);
Tg = -inv(D-L)*U;
cg = inv(D-L)*b;
A=L+D+U das eingesetzt in Ax=b
x =( Tg*x + cg);
k = 1;
x = zeros(2250,1); startwert
Tg = -inv(D-L)*U;
cg = inv(D-L)*b;
Ist das so richtig? bzw. bei mir verändert sich das nur im 2. iterationschritt dann bleibts gleich ich speicher x zwischen und setz es dann wieder ein ist das der Sinn oder muss man da noch mehr machen?
while k <= N
x =( Tg*x + cg);
k = k+1;
[X_rec]=plot_array(x,data);
figure(2)
subplot(1,2,1)
imagesc(X); colorbar
subplot(1,2,2)
imagesc(X_rec); colorbar
end
Ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt.
|
|
|
|
Hallo melchior2K,
lu berechnet eine Faktorisierung der Matrix also A=L*U. Das ist nicht das was du brauchst. Dann könnte man ja auch direkt lösen - ganz ohne iteratives Verfahren.
GRüße
mathemaduenn
|
|
|
|