Diese berechnete Geschwindigkeit aus dem gegebenen Arrays x und t, mit zentralen Unterschieden:
v = (x(3:$) - x(1:$-2)) ./ (t(3:$) - t(1:$-2))
Um zu sehen, was das ist, in jedem Bereich auf dem ersten Index konzentrieren:
(x(3) - x(1)) ./ (t(3) - t(1))
ist klar, Dies ist die Geschwindigkeit im 2. Moment der Zeit. Die Formel führt diese Berechnung für alle Zeiten durch, wenn es möglich ist; Die zentrierte Differenzformel gilt nicht im ersten und letzten Moment. Man will kann verkürztes Zeitbereich einzuführen, dies zu reflektieren:
tr = t(2:$-1)
Ebenso für die Beschleunigung:
a = (x(3:$) - 2*x(2:$-1) + x(1:$-2)) ./ (t(3:$) - t(1:$-2)).^2
Dies kann nun mit plot(tr,v)
oder plot(tr,a)
aufgetragen werden. Und ihre Werte zu sehen, wenn die Zeit 10 ist, verwenden
v(tr==10)
und
a(tr==10)