JuliaのMatlabスプライン関数

Question

Matlab（https://uk.mathworks.com/help/matlab/ref/spline.html）の同じアプローチを使用して、Juliaの2つの点を補間しようとしています。 Interpolations（https://github.com/tlycken/Interpolations.jl）ライブラリを試しましたが、作業用スクリプトの作成にいくつか問題があります。

私は2つの点（のは1.5と10.5をしましょう）とそれらの間の5 NAとDataArrays.DataArray {Float64,1}持っている：スプラインを実行するのに十分であるMATLABで

using DataFrames 
using Interpolations 

a = @data([1.5, NA, NA, NA, NA, NA, 10.5]); 

を関数。ジュリアでは、の補間関数は3次補間を可能にします。しかし、それはNAと一緒に働いていないようです。どうしたらいいですか？また、それらのポイントを補間するために同じ/類似のアルゴリズムを使用していますか？

Answer 1

これらが均等に配置されていると仮定すると、linspaceを使用できる線形補間を想定しています。あなただけの開始、終了、およびその間の値の数が必要です。どんなのNAの間の補間を行うために、より一般的に

linspace(a[1],a[end],sum(isna(a))) 

を、あなたは見つけることができる非NAは

idxs = find(~isna(a)) 

と値してください

for i in 1:length(idxs)-1 
    tmpidxs = idxs[i]:idxs[i+1] 
    a[idxs[i]+1:idxs[i+1]-1] = linspace(a[idxs[i]],a[idxs[i+1]],length(tmpidxs))[2:end-1] 
end 

あなたはそれをきれいにしたり、必要に応じて関数に入れることができます。それは、それぞれの既知の値の間にそれが線形に間隔を置いていると仮定するのは奇妙な仮定です。