aptknt
可接受的结序
语法
结= aptknt(tau,k)
[knot,k] = aptknt(tau,k)
描述
返回适合在数据点进行插值的结序列结= aptknt(tau,k)
τ
通过k阶的样条,提供这个结序列τ
至少有k
,是非递减的,满足τ(我)
<τ(i + k - 1)
对所有我
.在这种情况下,只有一条样条是有序的k
带结序列结
与这些地点的给定值相匹配。这是因为顺序结
返回满足勋伯格-惠特尼条件
节(i) <节(i) <节(i+k), i=1:长度(tau)
只有在极端结点处才相等,每一个结点都有精确的多重性k
.
如果τ
小于k
条目,然后k
减为值长度(τ)
.如果出现以下情况,则会产生错误τ
非递减和/或τ(我)
=τ(i + k - 1)
对于一些我
.
也返回实际的[knot,k] = aptknt(tau,k)
k
使用(等于输入中较小者k
而且长度(τ)
).
例子
如果τ
是等距的,例如,等于linspace (a、b、n)
对于一些n
> = 4,y
序列的大小是否与τ
,然后Sp = spapi(aptknt(tau,4),tau,y)
给出带有非结结束条件的三次样条插值。这与命令生成的三次样条相同样条(τ,y)
,但是是B-form而不是ppform。
不过值得注意的是
如果τ
是非常不均匀的,那么使用结果结序列插值到数据的站点τ
可能导致不满意的结果。
算法
的(k - 1)
分平均sum(τ(i + 1:我+ k - 1)) / (k - 1)
序列的τ
,由aveknt(τ,k)
,加ak
倍τ(1)
和一个k
倍τ(结束)
.换句话说,该命令给出的结果与augknt([τ(1)aveknt(τ,k),τ(结束),k)
,只要τ
至少有k
条目和k
大于1。
这个话题有用吗?<我nput type="button" value="是的" onClick="openFeedbackWindow('//www.tatmou.com/programs/bounce_hub_help.html?s_cid=Help_Topic_Survey&surveyParams='+document.location.href+'-YES-R2017a');">