主要内容gydF4y2Ba

newkntgydF4y2Ba

新断点分布gydF4y2Ba

语法gydF4y2Ba

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

描述gydF4y2Ba

newknot = newknt(f,newl)gydF4y2Ba返回其内部结点截断的基本间隔的结序列gydF4y2BafgydF4y2Ba成gydF4y2BanewlgydF4y2Ba块,这样就可以使某个分段线性单调函数与的高导数相关gydF4y2BafgydF4y2Baequidistributed。gydF4y2Ba

目的是选择一个适合于函数的精细近似的结序列gydF4y2BaggydF4y2Ba其在gydF4y2BafgydF4y2Ba是否包含足够的信息gydF4y2BaggydF4y2Ba让这个可行。gydF4y2Ba

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

[…,distfn] = newknt(...)也返回,在gydF4y2BadistfngydF4y2Ba,则该分段线性单调函数的ppform为等分布。gydF4y2Ba

例子gydF4y2Ba

如果误差是最小二乘近似gydF4y2BaspgydF4y2Ba一些数据gydF4y2Bax, ygydF4y2Ba通过一个样条的顺序gydF4y2BakgydF4y2Ba似乎不均匀,你可以尝试更公平的分布结使用gydF4y2Ba

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

另一个例子,请参见gydF4y2Ba利用搭配法求解带边界层的非线性ODEgydF4y2Ba.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