Hilfe bei MATLAB-Schnippsel

[Dennis S. (eltio)]

Guten Morgen,
ich habe Verständnisprobleme bei folgendem MATLAB-Code:

   % find local minimums
   dif = sign(diff(f));
   pos = find(dif > [dif(1) dif(1: end - 1)]);

In "dif" befinden sich 1 oder 0 abhängig vom Vorzeichen, ist klar. Kann 
mir jemand erklären was in der zweiten Zeile passiert? Das Argument des 
find-Befehls will mir einfach nicht einleuchten...

Dank & Gruß,
Dennis

[Walter T. (nicolas)]

In "dif" befinden sich 1en, 0en und -1en.

[Dennis S. (eltio)]

Richtig... aber für die Frage auch zweitrangig, oder?

[Marcel B. (marcel_b)]

Das "[dif(1) dif(1: end - 1)]" erzeugt doch nur einen neuen Vektor, mit 
dem dann dif verglichen wird. Warum das in diesem Falle gemacht wird, 
geht aus dem Code nicht hervor!

[Justus S. (jussa)]

Dennis S. schrieb:
> Kann
> mir jemand erklären was in der zweiten Zeile passiert?

das was darüber steht, es werden lokale Minima gesucht, also wo
(f(n) - f(n+1)) > 0 und (f(n+1) - f(n+2)) < 0

[Dennis S. (eltio)]

Hmm, okay ich glaube jetzt habe ich es: der Vektor "dif" wird ohne das 
letzte Element kopiert und das erste Elemtent des alten Vektors an den 
Anfang gestellt.. also
dif_alt = [1 2 3 4 5 6]
dif_neu =  1 1 2 3 4 5

Richtig?

Gruß Dennis

[Justus S. (jussa)]

jup, die Vektoren müssen halt gleich lang sein...man könnte auch
dif(1: end-1) und dif(2:end) vergleichen, kommt auf das selbe heraus, 
lediglich die Indizes sind verschoben...

[Dennis S. (eltio)]

Du meinst man könnte da auch so schreiben:

   % find local minimums
   dif = sign(diff(f));
   pos = find(dif(1: end-1) > dif(2:end));

Das wäre für meine Aufgabe recht praktisch... :-)

[Justus S. (jussa)]

ja, im Vergleich zur ersten Variante sind die Indizes in pos halt um 1 
geringer

[Dennis S. (eltio)]

Damit könnte ich leben. Vielen Dank für die Hilfe! :-)

[Dennis S. (eltio)]

Justus Skorps schrieb:
> ja, im Vergleich zur ersten Variante sind die Indizes in pos halt um 1
> geringer

Ohne dir auf die Füße treten zu wollen, aber kann es sein, dass du dich 
irrst? Ich hatte gerade die Gelegenheit das auzuprobieren und es kommen 
zwei unterschiedliche Ergebnisse heraus (73 gegenüber 38 und 500).

Gruß Dennis

[stebbl (Gast)]

Der Ursprünliche Code vergleicht folgende Elemente:

dif(1) > dif(1), dif(2) > dif(1), dif(3) > dif(2),... , dif(n) > 
dif(n-1)

die zweite variante vergleicht:

dif(1) > dif(2), dif(2) > dif(3), ... , dif(n-1) > dif(n)

=> umgekehrt

pos = find(dif(2:end) > dif(1: end-1));

oder

pos = find(dif(1:end-1) < dif(2: end));

Schöne Grüße
Stephan