moresplines

moresplines - ylabel('f(x)');

This preview shows page 1. Sign up to view the full content.

function moresplines(x,y); %Problem 16.5 %Fit data with a cubic spline with not-a-knot end conditions %v = ppval(pp,xx) returns the value of the piecewise polynomial f, %contained in pp, at the entries of xx. %pp = spline(x,Y) returns the piecewise polynomial form of the cubic spline %interpolant for later use with ppval and the spline utility unmkpp. x must %be a vector. gg=linspace(min(x),max(x),101); yy=interp1(x,y,gg,'spline'); plot(x,y,'o',gg,yy); title('(a) Piecewise cubic hermite interpolation'); xlabel('x');
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: ylabel('f(x)'); %-------------------------------------------------------%Cubic spline with zero-slope clamped end conditions yy=interp1(x,y,gg,'spline'); plot(x,y,'o',gg,yy); title('(a) Piecewise cubic hermite interpolation'); xlabel('x'); ylabel('f(x)'); %-------------------------------------------------------%Piecewise cubic hermite interpolation yy=interp1(x,y,gg,'pchip'); plot(x,y,'o',gg,yy); title('(c) Piecewise cubic hermite interpolation'); xlabel('x'); ylabel('f(x)');...
View Full Document

This note was uploaded on 09/27/2011 for the course EGM 3344 taught by Professor Raphaelhaftka during the Spring '09 term at University of Florida.

Ask a homework question - tutors are online