【Matlab的三次样条插值程序】

%示例开始---------------------------------------
%[x]=[2.1,12.5,33.5,77.5,98];
%[y]=[50,100,150,200];
%[z]=[
%85.99992857,91.24228571,91.38814286,83.01521429;
%80.71485714,82.33292857,83.56392857,81.49857143;
%88.30728571,87.38005921,82.24071429,80.05985714;
%99.81578571,98.06721429,100.9337143,93.43307143;
%102.6605356,100.7432857,108.8566429,95.85807143;
%];
%
%[xi]=[3:1:98];
%[yi]=[50:1:200];
%示例结束----------------------------------------

%初始化说明
%调用前先输入一维数组x、y，xi，yi和二维数组z
%x，y为插值节点，z为节点函数值，xi，yi为插值点
%要求，x长度=z的行数，y的长度=z的列数

%插值参数调整
%(1)最近邻插值     method=nearest
%(2)线性插值       method='linear'
%(3)三次样条插值   method='spline'
%(4)三次插值       method='pchip'或'cubic'

%程序开始，三次样条插值---------------------------
[rz,cz]=size(z);[rx,cx]=size(x);[ry,cy]=size(y);
[rxi,cxi]=size(xi);[ryi,cyi]=size(yi);
if cx == rz && cy == cz