求问用matlab做三次样条拟合

qinkong

x = 0:10;
y = sin(x);
xx = 0:.25:10;
yy = spline(x,y,xx);
plot(x,y,'o',xx,yy)

这是matlab帮助里面的一个example。我想做过去几十年中某个气象要素比如年降水的三次样条拟合；我不知道怎样用matlab的spline函数实现，上面的x=0:10相当于选择好的分割点么？50年的年降水序列如何合理的选择分割点的数目？分割点之间的数据信息如何利用？

总之我不明白怎样用这个函数实现年降水序列的三次样条拟合。

恳请了解的朋友答疑解惑！

谢谢！

沙颖凯

x, y为已知的自变量，xx为需要插值的自变量，yy为插值得到的因变量。

你问的分割点什么的不知道是想说什么。不过我觉得很奇怪，插值一般是生成序列或者确定特定点取值用的，既然已经有现成的序列为什么要插值，年降水序列按理说可以直接进行分析。

qinkong

沙颖凯 发表于 2014-8-3 13:26
x, y为已知的自变量，xx为需要插值的自变量，yy为插值得到的因变量。

你问的分割点什么的不知道是想说什 ...

我做的不是年降水序列，是另一种不过说起来比价复杂。另外我不是想做插值，而是想分析一下该要素的年变化趋势。因此想用三次样条拟合一下生成一个相对平滑的拟合曲线，来判断变化趋势。

滑动平均应该也可以做，但是我担心滑动平均不能比较准确的看出突变的时点在哪里。

沙颖凯

qinkong 发表于 2014-8-3 14:01
我做的不是年降水序列，是另一种不过说起来比价复杂。另外我不是想做插值，而是想分析一下该要素的年变化 ...

年变化趋势应该包括很多方面吧，线性和非线性的增减，周期的变化，方差或者说振幅的变化之类，样条感觉只能平滑一下啊，不一定能看出很多问题，不过也看具体数据啦。

突变点可以考虑看看MK检验有没有明显的突变点。

小萌萌

我就是想问问为什么要用三次样条。。而且这是插值不是拟合呀。。