主要内容gydF4y2Ba

newkntgydF4y2Ba

新断裂分布gydF4y2Ba

语法gydF4y2Ba

newknots = newknt (f, newl)gydF4y2Ba
newknt (f)gydF4y2Ba
[…,distfn] = newknt(...)

描述gydF4y2Ba

newknots = newknt (f, newl)gydF4y2Ba返回的结序列内部结的基本时间间隔gydF4y2BafgydF4y2Ba成gydF4y2BanewlgydF4y2Ba块,这样一个分段线性单调函数的高阶导数gydF4y2BafgydF4y2Baequidistributed。gydF4y2Ba

目的是选择一个结序列适用于近似函数的罚款gydF4y2BaggydF4y2Ba的粗糙近似gydF4y2BafgydF4y2Ba被认为包含足够的信息吗gydF4y2BaggydF4y2Ba让这个可行。gydF4y2Ba

newknt (f)gydF4y2Ba用途gydF4y2BanewlgydF4y2Ba它的默认值,即多项式块的数量gydF4y2BafgydF4y2Ba。gydF4y2Ba

[…,distfn] = newknt(...)也返回,gydF4y2BadistfngydF4y2Ba的ppform equidistributed分段线性单调函数。gydF4y2Ba

例子gydF4y2Ba

如果误差的最小二乘近似gydF4y2BaspgydF4y2Ba一些数据gydF4y2Bax, ygydF4y2Ba花键的秩序gydF4y2BakgydF4y2Ba看起来不均匀,你可能会尝试一个更公平的分配节通过使用gydF4y2Ba

spap2 (newknt (sp), k, x, y);gydF4y2Ba

另一个例子,看到gydF4y2Ba解决非线性ODE边界层的搭配gydF4y2Ba。gydF4y2Ba

算法gydF4y2Ba

这是Fortran程序gydF4y2BaNEWNOTgydF4y2Ba在gydF4y2Ba动力gydF4y2Ba。与gydF4y2BakgydF4y2Ba的顺序分段多项式函数gydF4y2BafgydF4y2Ba在gydF4y2Ba页gydF4y2Ba,函数|gydF4y2BaDgydF4y2BakgydF4y2BafgydF4y2Ba|由分段常数函数近似得到的地方、gydF4y2Ba离散,gydF4y2Ba分化的变化gydF4y2BaDgydF4y2BakgydF4y2Ba1gydF4y2BafgydF4y2Ba。新的打破序列选择细分gydF4y2Ba分段多项式的基本时间间隔gydF4y2BafgydF4y2Ba以这样一种方式gydF4y2Ba

∫gydF4y2Ba newknotsgydF4y2Ba (gydF4y2Ba 我gydF4y2Ba )gydF4y2Ba newknotsgydF4y2Ba (gydF4y2Ba 我gydF4y2Ba +gydF4y2Ba 1gydF4y2Ba )gydF4y2Ba |gydF4y2Ba DgydF4y2Ba kgydF4y2Ba fgydF4y2Ba |gydF4y2Ba 1gydF4y2Ba /gydF4y2Ba kgydF4y2Ba =gydF4y2Ba 常量gydF4y2Ba ,gydF4y2Ba 我gydF4y2Ba =gydF4y2Ba kgydF4y2Ba :gydF4y2Ba kgydF4y2Ba +gydF4y2Ba newlgydF4y2Ba −gydF4y2Ba 1gydF4y2Ba