线性插值用c语言程序源代码,实现快速线性插值的实现

Fast Linear Interpolation

F = lininterp1f(X,Y,XI,Ydefault) returns the value of the 1-D function Y at the points XI using linear interpolation. Length(F)=length(XI). The vector X specifies the coordinates of the underlying interval. Ydefault is returned for values of XI outside the coordinates in X.

For lininterp1f to work properly:

X must be a monotonically increasing array;

Y must be an array with length(Y)=length(X);

XI must be an array.

Ydefault must be a scalar value or an empty matrix [].

Example:

>> x = [1:1:1000];

>> y =log(sqrt(x+1.001)-1.001);

>> xv =[5:.001:100];

>> yinterp =lininterp1f(x,y,xv,[]);

Warning: not much in the way of error checking, since this slows things down, so pay attention to the argument passed to the function!!!

Installation (Windows users):