How to make a spline with end point limits?

3 visualizaciones (últimos 30 días)
Abdul Mannan
Abdul Mannan el 15 de Sept. de 2015
Comentada: Abdul Mannan el 15 de Sept. de 2015
I want to make a spline of some data x (input) and y (output), with constrained slope (first derivative) and zero initial and final 2nd derivative. For example
if true
x = [ 0, 0.125, 2*0.125, 3*0.125, 4*0.125, 5*0.125, 6*0.125, 7*0.125, 8*0.125 ];
y = [0, sin(2*pi*f/8), 1, sin(2*pi*3*f/8), 0, sin(2*pi*5*f/8),-1,sin(2*pi*7*f/8), 0];
end
I know that I can find fifth order spline with this like
if true
k = 8;
sp = spapi( optknt(x,k), x, y );
xx = 0:0.0001:1;
yy = fnval(xx,sp); % Evaluating data for xx using given spline
dsp = fnder(sp);
dyy = fnval(xx,dsp); % First Derivative
ddsp = fnder(sp,2); % Second Derivative
ddyy = fnval(xx,ddsp);
plot(xx, yy, xx, dyy, xx, ddyy), shg % Plotting all data.
end
Although, it is good, but I want to put constrained on second derivative at end points. How can I do that?
  1 comentario
Abdul Mannan
Abdul Mannan el 15 de Sept. de 2015
OK. This is explained in help in "spapi" MATLAB Command. You can do it as follows.
if true
x = [ 0, 0.125, 2*0.125, 3*0.125, 4*0.125, 5*0.125, 6*0.125, 7*0.125, 8*0.125 ];
y = [0, sin(2*pi*f/8), 1, sin(2*pi*3*f/8), 0, sin(2*pi*5*f/8),-1,sin(2*pi*7*f/8), 0];
k = 6;
sp = spapi(augknt(x,6,2),[x,x,min(x),max(x)],[y,[0 2*pi*cos(2*pi*f/8), 2*pi*cos(2*pi*f/4),2*pi*cos(2*pi*3*f/8),2*pi*cos(2*pi*4*f/8),2*pi*cos(2*pi*5*f/8),2*pi*cos(2*pi*6*f/8),2*pi*cos(2*pi*7*f/8),0],0,0]);
xx = 0:0.001:1;
yy = fnval(xx,sp);
dsp = fnder(sp);
dyy = fnval(xx,dsp);
ddsp = fnder(sp,2);
ddyy = fnval(xx,ddsp);
plot(xx, yy, xx, dyy, xx, ddyy), shg
end

Iniciar sesión para comentar.

Respuestas (0)

Categorías

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

Etiquetas

Productos

Community Treasure Hunt

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

Start Hunting!

Translated by