%x为已知二维矩阵13*22, 采用线性插值。for j=1:size(x,2)x(isnan(x(:,j)),j)=interp1(find(~isnan(x(:,j))),x(~isnan(x(:,j)),j),find(isnan(x(:,j))));end