分段线性插值matlab程序

一个很经典的分段线性插值的matlab程序，程序短小精悍，其中应用了向量思想，还有数组的逻辑坐标。

% function v=piecelin（xyu）
%PIECELIN Piecewise linear interpolation
% v=piecelin（xyu） finds the piecewise linear L（x）
% with L（x（j））=y（j） and returns v（k）=L（u（k））
% First divided difference
clear
clc
x=1:6;
y=[16 18 21 17 15 12];
u=[2.5 2 4.5 3.7 5 5.8];
delta=diff（y）./diff（x）;

n=length（x）;
k=ones（size（u））;
for j=2:n-1

    k（x（j）<=u）=j
    
end

s=u-x（k）;
v=y（k）+s.*delta（k）;
figure 
plot（xy‘*-‘）
hold on
plot（uv‘o‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        455  2010-01-05 15:41  piecelin.m

----------- ---------  ---------- -----  ----

                  455                    1