Transform 2d spline function f(t) into f(x)

后端 未结 3 1394
刺人心
刺人心 2020-12-11 04:55

So I\'ve got a special case set of cubic splines, whose 2d control points will always result in a curve that will never cross itself in the x axis. That is, the cu

相关标签:
3条回答
  • 2020-12-11 05:06

    Often people will use a root finding technique (like Newton's Method) if an numerical approximation is good enough.

    0 讨论(0)
  • 2020-12-11 05:23

    Well, you could solve your fx(t)=x for t. That would be a cubic equation; ugly but still possible to solve explicitely. If your spline is as you describe it, then two of the solutions will be conjugate complex, so the only remaining one is the one to take. Use that to compute y=fy(t). I doubt you can accompish anything easier if you want exact solutions.

    You can use the general formula from Wikipedia to compute the solution of the cubic equation.

    0 讨论(0)
  • 2020-12-11 05:26

    Your question states that you want evenly spaces x coordinates, and approximate solutions are all right. So I propose the following algorithm:

    • Decide on the grid points you want, e.g. one every 0.1 x units.
    • Start with l = 0 and r = 1.
    • Compute fx(l) and fx(r) and consider the interval denoted by these endpoints.
      • If the interval is sufficiently small and contains exactly one grid point, use the central parameter t=(l+r)/2 as a good approximation for this grid point, and return that as a one-element list.
      • If there is at least one grid point in that interval, split it in two using (l+r)/2 as the splitting point, and concatenate the resulting lists from both computations.
      • If there is no grid point in the interval, then skip the current branch of the computation, returning an empty list.

    This will zoom in on the grid points, bisecting the parameter space in each step, and will come up with suitable parameters for all your grid points.

    0 讨论(0)
提交回复
热议问题