Assume we are given a set of $n$ points from $\mathbb R^2$, $(x_1,y_1),(x_2,y_2)\dots(x_n,y_n)$.
We want to construct a path connecting all these points using a pair of parametric equations \begin{cases} x = f(t) \\ y = g(t) \\ \end{cases} such that \begin{cases} f(0) = f(1) = x_1 \\ g(0) = g(1) = y_1 \\ \end{cases} and find a set of points $(t_2,t_3\dots t_n),$ for which \begin{cases} f(t_k) = x_k \\ g(t_k) = y_k \\ \end{cases} and the length of the path from $(x_1,y_1)$ to itself via all the other points is minimized.
Now, finding A path is easy. Divide the unit line into $n$ parts, take the starting point of each part to be $t_k$ and construct two interpolation polynomials of degree $n$, $f_p$ and $g_p$, such that $f_p(t_k)=x_k$ and $g_p(t_k)=y_k$ and set \begin{cases} x = f_p(t) \\ y = g_p(t) \\ \end{cases}
Example: Let the original points be $(0,0),(1,1),(2,0)$ and select points of time to be $(1/3,2/3)$. Now the parametric interpolation function is \begin{cases} x = -\frac{27}{2} (-1 + t) t^2 \\ y = \frac{9}{2} (-1 + t) t (-2 + 3 t) \\ \end{cases}
The length of this path is acquired via formula $L=\int_0^1 \sqrt{{\frac{dx}{dt}}^2+{\frac{dy}{dt}}^2} dt$ and it is approximately $5.422$.
However, if we select times to be $(1/4,3/4),$ we get function \begin{cases} x = -\frac{8}{3} (-1 + t) t (1 + 4 t) \\ y = \frac{8}{3} (-1 + t) t (-3 + 4 t) \\ \end{cases} and the length is now approximately $5.411$.
Thus optimizing the length of the path requires two components, selecting the times $t_k$ and also selecting the original interpolation functions. The times can be optimized via normal differentiation methods, but how can we optimize over the interpolation functions? The method for solving this type of problem is calculus of variations. Can it be applied into this type of problem with interpolation and parametric equations?
This is some kind of "continuous" variant of the travelling salesman problem and I was hoping maybe I could get some faster, more robust methods for solving TSP in general from this, but looking at the simplest of examples, I notice that both the interpolation functions and the required integrations become exceedingly complex very fast. So is there some methods for simplifying the process?
Edit: Is there a quick way to get good selection of $t_k$s?