Wie erhält man eine Diagonalmatrix aus einem Vektor?

Question

Ich habe eine Kolumne:

0.0677
0.0584
0.0487
0.0453
0.0394

Welche Anweisung würde die folgende Ausgabe ergeben

0.0677   0          0         0          0
0        0.0584     0         0          0
0        0          0.0487    0          0
0        0          0         0.0453     0
0        0          0         0          0.0394

Answer 1

Diag ist die normale MATLAB-Lösung (wie von posdef hervorgehoben).

D = diag(vec);

ergibt eine Matrix mit Diagonalelementen nach Bedarf.

In manchen Anwendungen ist es vielleicht besser, eine dünnbesetzte Matrix zu erstellen, da eine Diagonalmatrix recht dünnbesetzt ist. Bei der Multiplikation von Matrizen hilft dies, die Zahl der unnötigen Operationen zu verringern.

n = length(vec);
D = spdiags(vec(:),0,n,n);

Wenn Sie die Zuordnung wirklich in expliziter Form vornehmen wollen, verwenden Sie einen einzelnen linearen Index wie diesen:

n = length(vec);
D = zeros(n);
D(cumsum([1,repmat(n+1,1,n-1)])) = vec;

Oder Sie können die Funktion sub2ind verwenden, um eine Reihe von Indizes in einen einzigen Index umzuwandeln.

Answer 2

Wenn ich mich richtig erinnere, gibt es einen Befehl, der in etwa so heißt diag(A)

Edit : Bitte sehr, eine Dokumentation über das diag http://www.mathworks.com/help/techdoc/ref/diag.html

Achten Sie besonders auf das Zitat:

X = diag(v) setzt v auf die Hauptdiagonale, wie oben mit k = 0.

Answer 3

Es ergibt sich die folgende Diagonalmatrix D dessen Diagonale der Vektor vec . Sie ist in MATLAB in vektorisierter Form geschrieben.

D      = zeros(numel(vec));
[I,J]  = ind2sub(size(D),1:numel(D));
ind    = [I(:) J(:)];
ind    = find(ind(:,1)==ind(:,2));
D(ind) = vec;

Answer 4

Nun, offensichtlich kann man das auf eine C-ähnliche Weise tun. Im Moment kann ich nicht herausfinden, eleganter Lösung.

vector;%Your vector
vec_length = length(vector);
A = zeros(vec_length);
for i=1:vec_length
   A(i,i)  = vector(i);
end;