How to make a spline with end point limits?

3 vues (au cours des 30 derniers jours)
Abdul Mannan
Abdul Mannan le 15 Sep 2015
Commenté : Abdul Mannan le 15 Sep 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 commentaire
Abdul Mannan
Abdul Mannan le 15 Sep 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

Connectez-vous pour commenter.

Réponses (0)

Catégories

En savoir plus sur Splines dans Help Center et File Exchange

Produits

Community Treasure Hunt

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

Start Hunting!

Translated by