Mercurial > repos > public > sbplib
diff +parametrization/dataSpline.m @ 1106:00203fcc962f
Merged in feature/dataspline (pull request #14)
Feature/dataspline
Approved-by: Jonatan Werpers <jonatan.werpers@it.uu.se>
Approved-by: Martin Almquist <malmquist@stanford.edu>
| author | Martin Almquist <malmquist@stanford.edu> |
|---|---|
| date | Tue, 09 Apr 2019 20:24:21 +0000 |
| parents | 36d092a00040 |
| children | 60c875c18de3 |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/+parametrization/dataSpline.m Tue Apr 09 20:24:21 2019 +0000 @@ -0,0 +1,19 @@ +% dataSpline calculates a Curve through the points f_i using cubic spline interpolation. +% The spline curve is parametrized with the arc length parametrization +% to facilitate better grids. +% +% f - m x D matrix of m points in D dimensions +function C = dataSpline(f) + m = size(f, 1); + + t = linspace(0,1,m); + + pp_g = spapi(4, t, f'); + pp_gp = fnder(pp_g); + + g = @(t) fnval(pp_g, t); + gp = @(t) fnval(pp_gp, t); + + C = parametrization.Curve(g, gp); + C = C.arcLengthParametrization(); +end