Fit a spline function to a set of data points

49 visualizaciones (últimos 30 días)
Deepa Maheshvare
Deepa Maheshvare el 3 de Jul. de 2022
Comentada: Deepa Maheshvare el 4 de Jul. de 2022
Hi All,
I have a set of coordinate points
x= [0,4,6,10,15,20]
y = [18,17.5,13,12,8,10]
I want to fit a spline function.
Could someone please suggest the the curve fit function that can be used?
  7 comentarios
Deepa Maheshvare
Deepa Maheshvare el 3 de Jul. de 2022
x = [0,4,6,10,15,20]
y = [18,17.5,13,12,8,10]
vq = interp1(x,y,'linear')
returns NaN
Deepa Maheshvare
Deepa Maheshvare el 3 de Jul. de 2022
I tried
x = [0,4,6,10,15,20]
v = [18,17.5,13,12,8,10]
vq = interp1(x,v, 'linear', 'pp')
xq = x
figure
vq1 = interp1(x,v,xq);
plot(x,v,'o',xq,vq1,':.');
title('(Default) Linear Interpolation');
returns
vq =
struct with fields:
form: 'pp'
breaks: [0 4 6 10 15 20]
coefs: [5×2 double]
pieces: 5
order: 2
dim: 1
orient: 'first'
Here I am not sure why the order (which I infer as the polynomial degree) is 2 for linear interpolation. Could you please explain this a bit?

Iniciar sesión para comentar.

Respuesta aceptada

Sam Chak
Sam Chak el 3 de Jul. de 2022
Thought you want the 'spline'? Is there a reason to use the 'linear' method?
x = [0, 4, 6, 10, 15, 20];
y = [18, 17.5, 13, 12, 8, 10];
xq = 0:0.01:20;
subplot(2,1,1)
yq1 = interp1(x, y, xq, 'linear');
plot(x, y, 'o', xq, yq1, '--');
ylim([5 22]), grid on, title('Linear Interpolation');
subplot(2,1,2)
yq2 = interp1(x, y, xq, 'spline');
plot(x, y, 'o', xq, yq2, '--');
ylim([5 22]), grid on, title('Spline Interpolation');
  3 comentarios
Torsten
Torsten el 3 de Jul. de 2022
Seems that order = 2 means that two coefficients are needed to describe the (linear) spline function on each subinterval and is set to 1 + degree of the interpolating polynomial. A bit confusing in my opinion.
Deepa Maheshvare
Deepa Maheshvare el 4 de Jul. de 2022
Thanks a lot for the clarification.

Iniciar sesión para comentar.

Más respuestas (1)

Simon Chan
Simon Chan el 3 de Jul. de 2022
Use function spline

Categorías

Más información sobre Interpolation en Help Center y File Exchange.

Productos


Versión

R2021a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by