文档

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');">