DOSEIの日記

技術メモ+日常ログ

周期境界条件の cspline (interpolation of periodic function by cspline)

標準の spline は、周期境界条件では使えない。
Spline Toolbox の csape で可能。何の略かは不明。 Octave-forge にも互換パッケージがある。

例:

n=10; % 標本点数
m=1000; % リサンプリング点数

% サンプリング点
x=(0:(n-1))./(n-1);
xx=(0:(m-1))./(m-1);

% 適当な標本データ
y=rand(n,1);
y(1)=y(n); % 両端一致

% 周期スプライン補間
s=csape(x,y,"periodic");
yy=fnval(s,xx);

plot(x,y,"-o",xx,yy);