Problem:
mam obliczyć pozycję peryhelium orbity, tz. r(t_min) = (x,y,z),
mając dane tylko dwa punkty: r0, v0, oraz r1 i v1, i wiem że minimum jest pomiędzy tymi punktami.
|r| = r_min.
Zatem biorę interpolację kubiczną, bo taka akurat tu podchodzi:
mamy tylko dwa punkty i dwie pochodne, czyli prędkości, w tych samych punktach.
y(t=0), y'(0), oraz y(1), y'(1); gdzie y jak i v to wektory: x,y,z;
Zatem obliczam ten wielomian 3-go stopnia, i tam wyjdzie coś takiego:
gdzie: a, b, c i d należy wyznaczyć, co jest dość proste, np. od razu można zgadnąć, że:
y(0) = d, oraz: y'(0) = c = v(0);
Pozostają tylko a i b do wyznaczenia, co też można sobie łatwo wyliczyć:
y(1) = a+b+c+d
y'(1) = 3a + 2b + d;
...
OK. Mamy ten wielomian 3-stopnia... i cóż z tego?
Gdzie jest to minimum:
|r(t_min)| = r_min, jak wyznaczyć: t = t_min, żeby wyliczyć pozycję: r(t_min) = ?