2010-11-19 周期境界条件の cspline (interpolation of periodic function by cspline) メモ Matlab Octave 標準の 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);