R 具有不规则网格的Julia三次样条插值
我正在将一些代码从R移植到Julia,并努力翻译R的R 具有不规则网格的Julia三次样条插值,r,julia,R,Julia,我正在将一些代码从R移植到Julia,并努力翻译R的spline函数。我需要一个Julia函数: function spline_j(x,y,xout) spline(x,y,,"fmm",,,xout) 这将产生与调用R函数相同的返回: function spline_j(x,y,xout) spline(x,y,,"fmm",,,xout) i、 e.使用Forsyth、Malcolm和Moler的方法,这是R中的默认方法 Myx和y始终是一维的,但是x的点之间没有规则的间隔。这种非
spline
函数。我需要一个Julia函数:
function spline_j(x,y,xout)
spline(x,y,,"fmm",,,xout)
这将产生与调用R函数相同的返回:
function spline_j(x,y,xout)
spline(x,y,,"fmm",,,xout)
i、 e.使用Forsyth、Malcolm和Moler的方法,这是R中的默认方法
Myx
和y
始终是一维的,但是x
的点之间没有规则的间隔。这种非规则性似乎排除了使用纯Julia插值
包的可能性,因为“目前仅线性插值
支持不规则网格”
该软件包支持不规则的x
,因此spline_j
的候选者是:
using Dierckx
function spline_j(x, y, xout)
spl = Dierckx.Spline1D(x, y)
spl(xout)
end
如果方法
是“自然的”
,则匹配R的样条曲线
函数
是否可以在Julia中复制R的
“fmm”
方法?您可以查看这个包含三次样条曲线的软件包:谢谢Michael。不幸的是,datainterpolation
没有我想要的(至少在我第一次阅读文档时)。但你的评论让我觉得我需要学习如何最好地搜索Julia软件包。这个想法很快就把我带到了一个关于朱莉娅的演讲中,作者是。。。你。哈哈——幸好从那以后情况有所改善。现在可以搜索。我只是在那里搜索“三次样条曲线”,它让我找到了