a=interp1(x,y,z,'spline');
interp1是matlab里的插值函数,用法:
Yi=interp1(x,y,xi,'方法参数')
(1)x是原始数据自变量单增向量y是原始数据函数向量或矩阵,若为矩阵则按列插值。
(2)xi是要插值的自变量矩阵,yi是由插值方法计算出的函数矩阵。
(3)方法参数有:Nearset(近点插值)linear(线性)spline(三次样条)cublic(内插)
把NaN去掉就好了呀,不过要注意坐标向量也要把相应的位置去掉,
以保证坐标和数据一一对应,假矩阵的一列是y,坐标是x=(1:length(y))'
这样做就可以了
m=isnan(y);
x(m)=[];
y(m)=[];
yi=interp1(x,y,xi);