Matlab - Función spline

Buenas tardes,

Estoy liado con la función spline, y no puedo resolver mi problema.

Para un caso más sencillo, a partir de la siguiente gráfica



Que se trata de la rigidez de un resorte, de la recta obtengo la rigidez del resorte: K = 95 N/mm

Esta rigidez hay que "quitársela", de modo que uso:

Kr=95;
FzaR=(Fzacelula)+(Kr*zRmm)+PreQ;

Donde zRmm es la posición.


En otro caso, que se tiene resortes de paso variable, se obtiene:




Por lo tanto, ahora no es una recta F=kx. Para ello utilizo una spline con los valores de la curva. Mi problema es que no sé como hacer en este caso para "quitar" esa curva cómo hice antes. He intentado:

x=[1.181, 19.65, 34.85, 52.73, 64.16, 72.94, 79.98, 86.32, 90.31, 92.74, 93.66];
y=[-57.85, -184.2, -296.1, -454.5, -580.3, -710.5, -866, -1088, -1318, -1544, -1742];

FzaR2=(Fzacelula2)+spline(x,y,zRmm2)+PreQ2;

Pero no es correcto.

Muchas gracias, Saludos!

no veo errores en lo que comentas salvo ira un menos en vez de un mas:


Saludos

Daniel

Pues algo hago mal, porque, para el primer caso, en el que es una recta:

Original:



Aplicando:


Kr=95;
FzaR=(Fzacelula)+(Kr*zRmm)+PreQ;

Donde Kr es la pendiente y zRmm datos de posición en mm del archivo, obtengo:



Que es lo correcto.

En cambio, si aplico con spline, para varios puntos de la recta con pendiente Kr:

x=[2.427, 7.735, 15.92, 25.48, 35.54];
y=[-370.7, -382.8, -399.9, -415.5, -436.9];

FzaR=(Fzacelula)-spline(x,y,zRmm)+PreQ;

El resultado es:



El resultado no es el mismo.