Matlab三次样条法插值

Posted Z.L.G

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了Matlab三次样条法插值相关的知识,希望对你有一定的参考价值。

以N方向为例:

 

1、将N方向数据导入Matlab,将十进制年转化为年积日

2、重新排序,将缺失数据的天数以NaN补齐

3、寻找出NaN所在的天数

nxx = find( isnan(n) );

4、直接进行三次样条插值,2192代表了6年的数据

x = 1:2192;
nyy = spline(x,n,nxx);

5、nyy向量里即是需要插补的值,将其值插入到n中

 j = 1;
for i = 1 : 2192
    if ( isnan(n(i)) == 1 )
        n(i) = nyy(j);
        j = j + 1;
    end
end

 

问题:

数据连续缺失多于5个的时候,采用此方法插值可能会出现插入的数据失真的情况,要考虑使用所谓的多点三次样条法插值。

详见参考文献:

 

以上是关于Matlab三次样条法插值的主要内容,如果未能解决你的问题,请参考以下文章

三次样条插值

三次样条插值法

如何利用matlab实现多种插值

三次样条插值介绍

三次样条插值

如何使用 CuPy 在 python 上进行三次样条插值?