最小二乘样条逼近
返回样条的b形式f的订单样条
= spap2 (<一个href="#mw_80def079-53a5-40a2-8224-26f5258f536e" class="intrnllnk">结
,<一个href="#mw_539a8fa6-8aa7-4b82-bcc1-de1d24a6fff6" class="intrnllnk">k
,<一个href="#mw_6552bd4d-e96d-4487-a16d-a7598ef5c9f7" class="intrnllnk">x
,<一个href="#mw_43e7d63a-02ad-40c7-9417-0d0047fc2fab" class="intrnllnk">y
)k
用给定的结序列结
的
(*) y(:,j) = f(x(j)),所有j
在加权均方意义上,也就是总和
最小,默认权重为1。数据值y (:, j)
可以是标量、向量、矩阵或ND-arrays,以及|z|2是所有元素的平方和吗z。同一地点的数据点被它们的平均值所取代。
如果网站x
满足勋伯格-惠特尼条件
然后有一个给定顺序的唯一样条,并且结点序列完全满足(*)。的子序列不返回(**),否则不会返回样条x
。
spap2 (<一个href="#mw_5693a07c-ec6b-44ee-af07-3ff72043f012" class="intrnllnk">
,l
,<一个href="#mw_539a8fa6-8aa7-4b82-bcc1-de1d24a6fff6" class="intrnllnk">k
,<一个href="#mw_6552bd4d-e96d-4487-a16d-a7598ef5c9f7" class="intrnllnk">x
,<一个href="#mw_43e7d63a-02ad-40c7-9417-0d0047fc2fab" class="intrnllnk">y
)l
一个正整数,返回最小二乘样条近似值的b形式,但带有为您选择的结序列。通过应用求得结点序列<一个href="//www.tatmou.com/help/curvefit/aptknt.html">aptknt
到…的一个适当的子序列x
。得到的分段多项式包括l
多项式的分段k-2
连续的衍生品。如果你觉得内部结的不同分布可能会做得更好,跟着这个
sp1 = spap2 (newknt(花键),k, x, y));
允许指定权重样条
= spap2(…<一个href="#mw_6552bd4d-e96d-4487-a16d-a7598ef5c9f7" class="intrnllnk">x
,<一个href="#mw_43e7d63a-02ad-40c7-9417-0d0047fc2fab" class="intrnllnk">y
,<一个href="#mw_7d3e4106-535f-4c25-81f5-0ebc8276f53b" class="intrnllnk">w
)w
在误差测量中(如上所述)。w
必须是一个与x
,带有非负元素。当这些数据点被它们的平均值所取代时,与同一地点的数据点对应的所有权重将被求和。
spap2 (k {knorl1,…,knorlm}, {x1,…,xm}, y)
提供一个最小二乘样条近似网格数据。在这里,每一个knorli
要么是一个结序列,要么是一个正整数。此外,k
必须是一个米
向量,y
必须是(r + m
维数组,y (:, i1,…,im)
将基准面安装在网站
[x {1} (i1)…,x {m} (im)]
,所有i1
、……即时通讯
。然而,如果样条是标量值,那么,与单变量情况相反,y
是允许成为一个米
-维数组,在这种情况下y (i1,…,im)
基准面是否要安装在该处网站
[x {1} (i1)…,x {m} (im)]
,所有i1
、……即时通讯
。
spap2 (k {knorl1,…,knorlm}, {x1,…,xm}, y, w)
还允许指定权重。在这个米
变量的情况下,w
必须是一个单元格阵列米
条目,w{我}
大小相同的非负向量西
,否则w{我}
必须为空,在这种情况下我
变量。
spcol
调用来提供几乎块对角的搭配矩阵(Bj,k(x我)),<一个href="//www.tatmou.com/help/curvefit/slvblk.html">slvblk
在(加权)最小二乘意义上解决线性系统(*),使用块QR分解。
网格数据以张量积的方式进行拟合,一次一个变量,利用单变量加权最小二乘拟合与拟合值线性相关这一事实。