Treppennormalform in matlab < Matlab < Mathe-Software < Mathe < Vorhilfe
|
Hallo!
Ich habe folgenden Code in matlab geschrieben, um eine Matrix in die Treppennormalform zu bringen.
function [B,Q] = tnf(A) % B Treppennormalform, Q Tansformationsmatrix
[Z,T] = stufen(A); % Zeilenstufenform von A, T.A=Z
[n,m] = size(Z);
T1 = eye(n);
for j = 2:m
for k = 2:n % unterste 1 in der Spalte j suchen
if abs(Z(n-k+1,j)-1) < eps
break; % welche Zeile
end
end
for l = 2:j
if abs(Z(k,l-1)) > eps %
break;
end
for i = 1:k-1
if abs(Z(i,j)) < eps
break;
else
T1 = (add(n,i,k,-Z(i,j)))'*T1;
Z = T1*Z
end
end
end
end
Q = T1*T
B = Q*A
return;
Das Programm rechnet mit der Eingabematrix nach Erreichen der TNF weiter.
Kann mir jemand da helfen?
gruss,
logarithmus
|
|
|
|
Status: |
(Mitteilung) Reaktion unnötig | Datum: | 01:04 Sa 21.05.2005 | Autor: | Stefan |
Hallo logarithmus!
Es tut mir leid, dass dir keiner im von dir vorgesehenen Zeitraum bei der Beantwortung deiner Frage weiterhelfen konnte.
Vielleicht ja beim nächsten Mal!
Viele Grüße
Stefan
|
|
|
|