Matlab die Indizierung in 2D-array mit

Ich habe eine 2D matrix, sagen wir M=zeros(10,10);

Habe ich eine andere Spalte der matrix, V=[1;2;3;4;5;6;5;4;3;2];

Ich würde gerne in der Lage sein zu set M(i,j) = 1 für alle j >= V(i)

Ich weiß kann ich dies in eine Schleife

for i=1:10
   M(i,V(i):10) = 1;
end

aber es scheint, dass es möglich wäre die Verwendung einer form von Matlab-Indizierung zu vermeiden, die mit einer Schleife. Zum Beispiel so etwas wie :

M(:,V:10)=1;

oder 

M(:,V(:):10)=1;

aber weder von diesen produziert das erwartete Ergebnis.

Gibt es einige syntaktische Zucker ich verwenden können, um dies zu erreichen, oder sollte ich wieder Schleife?

InformationsquelleAutor Dave Durbin | 2012-06-13
  1. 2

    Da Sie bist suchen für syntaktischen Zucker, hier ist eine Art von esoterischen Weg, es zu tun. 🙂

    Vorausgesetzt, die Länge der V ist die Größe der beiden Dimensionen in die gewünschte matrix M erstellen Sie zuerst eine identity-matrix der gleichen Größe), dann ist index entsprechend und nehmen cumsum:

    V = [1;2;3;4;5;6;5;4;3;2]; #% 10x1 vector
E = eye(length(V), length(V)); #%10x10 identity matrix
M = cumsum(E(V,:),2)

M =

     1     1     1     1     1     1     1     1     1     1
     0     1     1     1     1     1     1     1     1     1
     0     0     1     1     1     1     1     1     1     1
     0     0     0     1     1     1     1     1     1     1
     0     0     0     0     1     1     1     1     1     1
     0     0     0     0     0     1     1     1     1     1
     0     0     0     0     1     1     1     1     1     1
     0     0     0     1     1     1     1     1     1     1
     0     0     1     1     1     1     1     1     1     1
     0     1     1     1     1     1     1     1     1     1

    Ok, jetzt: weniger Spaß, aber (auf meiner Maschine) schneller als alle anderen Optionen geprüft wurden bisher:

    n=10000;
V = randi(n-1, 1, n); #% as in @KevinRatelle's answer (but not transposed)

tic;
Vlinear = reshape(V + (0:n-1)*n, 1, []); #% find linear indices of first "ones"
M = zeros(n);
M(Vlinear)=1;
M=cumsum(M);
toc
    • Viele Punkte für Stil und Einfallsreichtum. Vielen Dank an alle für die Antworten - es war eine nützliche Lernerfahrung. Ich werde lernen, zu denken, die parallel...
    InformationsquelleAutor tmpearce
  2. 2

    Seine kaum subtiler und nicht wirklich besser als eine Schleife, ich glaube nicht, aber:

    [J,I] = meshgrid(1:10,1:10); 
V = [1;2;3;4;5;6;5;4;3;2];
M = J>V(I);

    Genießen.

    InformationsquelleAutor mwengler
  3. 1

    Habe ich versucht, das loop-Methode, und die 'meshgrid' - Methode. Ich Frage mich, über die Zeit zu berechnen, die für große Matrizen (seit das Problem mit Schleifen in matlab ist in der Regel die Zeit).

    Zuerst optimierte ich den code wie folgt Aussehen :

    V = V*ones(1,n);
N = ones(1,n)'*(1:n);
M = N>=V;

    Eigentlich N eine meshgrid, aber es scheint viel schneller zu tun, es auf diese Weise...

    Habe ich versucht, dieses :

    n = 10000;
V = randi(n-1,1,n)';

tic;
M = zeros(n);
for i=1:n
    M(i,V(i):n) = 1;
end
toc

tic;
[J,I] = meshgrid(1:n,1:n);
M = J>=V(I);
toc

tic;
V = V*ones(1,n);
N = ones(1,n)'*(1:n);
M = N>=V;
toc

    Und Ergebnisse waren :

    Elapsed time is 1.726872 seconds.
Elapsed time is 5.206657 seconds.
Elapsed time is 1.548600 seconds.

    Aber Methoden mit Hilfe von Matrizen anstelle von Schleifen sind Speicher für große n. Ich würde personnaly stick der Schleife.

    • Ich habe ein paar Optionen, darunter eine, die ist ziemlich schnell (schneller als diese) mit linearer Indizierung und cumsum. Wenn du Wert auf Geschwindigkeit testen Sie gegen diesen direkt, ich wäre daran interessiert zu wissen, wie es Tarife.
    • Auch in Ihrem nicht-loop-tests, beachten Sie, dass M wird eine logische matrix, es sei denn, Sie warf es zu verdoppeln, das wird Sie verlangsamen.
    • Ihre Methode ausprobiert habe (mit cumsum) und für N=10000, dauerte es 0.852 Sekunde, so dass Ihre Methode die Schnellste von allen, doppelt so schnell wie die Schleife ein
    InformationsquelleAutor Kevin Ratelle
  4. 0

    versuchen Sie dies:

    v = [1;2;3;4;5;6;5;4;3;2];
n = 10;
M = repmat((1:n)', 1, numel(v)) > repmat(v', n, 1);
    InformationsquelleAutor Jack Schwager
Schreibe einen Kommentar