双3次スプライン (2d) 評価関数
[zp[,dzpdx,dzpdy[,d2zpdxx,d2zpdxy,d2zpdyy]]]=interp2d(xp,yp,x,y,C [,out_mode])
doubleの mx行my列行列, スプラインを評価する点のx
座標.
doubleの mx行my列行列, スプラインを評価する点のy
座標.
doubleの1行nx列行列, 補間点のx
座標値.
i=1,2,...,nx-1について x(i)<x(i+1)が必要.
doubleの1行nx列行列, 補間点のy
座標値.
i=1,2,...,nx-1について y(i)<y(i+1)が必要.
双3次スプラインの係数. interp2d関数の入力引数は splin2d関数の出力引数です.
1行1列の文字列行列, [x(1),x(nx)]x[y(1),y(ny)]の外にある値の評価
mx行my列のdouble行列, スプラインのz
座標の評価値,
すなわち,i=1,2,...,mx および j = 1,2,...,my について,
zp(i,j)=s(xp(i,j),yp(i,j))
mx行my列のdouble行列, スプラインのx
に関する1階微分.
mx行my列のdouble行列, スプラインのy
に関する1階微分.
mx行my列のdouble行列, スプラインのx
に関する2階微分.
mx行my列のdouble行列, スプラインのx
とy
に関する2階微分.
mx行my列のdouble行列, スプラインのy
に関する2階微分.
双3次スプラインまたはサブスプライン関数を定義する
3つのベクトル (x,y,C)
(splin2d参照)を指定すると,
この関数は, 次のように(xp(i),yp(i))における
s (要すれば ds/dx,
ds/dy, d2s/dxx, d2s/dxy, d2s/dyy
も)を評価します:
out_mode
パラメータは,
捕外,すなわち(xp(i),yp(i))が
[x(1),x(nx)]x[y(1),y(ny)]にない場合,
の計算手法を定義します:
ゼロによる捕外が行われます
Nanによる捕外
捕外は以下のように定義されます :
捕外は(x,y)から最も近い双三次パッチにより行われます.
s
は周期的に拡張されます.
n = 7; // a n x n interpolation grid x = linspace(0,2*%pi,n); y = x; z = cos(x')*cos(y); C = splin2d(x, y, z, "periodic"); // now evaluate on a bigger domain than [0,2pi]x [0,2pi] m = 80; // discretization parameter of the evaluation grid xx = linspace(-0.5*%pi,2.5*%pi,m); yy = xx; [XX,YY] = ndgrid(xx,yy); zz1 = interp2d(XX,YY, x, y, C, "C0"); plot3d(xx, yy, zz1, flag=[2 6 4]) xtitle("extrapolation with the C0 outmode") | ![]() | ![]() |
n = 7; // a n x n interpolation grid x = linspace(0,2*%pi,n); y = x; z = cos(x')*cos(y); C = splin2d(x, y, z, "periodic"); // now evaluate on a bigger domain than [0,2pi]x [0, m = 80; // discretization parameter of the evaluation grid xx = linspace(-0.5*%pi,2.5*%pi,m); yy = xx; [XX,YY] = ndgrid(xx,yy); zz2 = interp2d(XX,YY, x, y, C, "by_zero"); plot3d(xx, yy, zz2, flag=[2 6 4]) xtitle("extrapolation with the by_zero outmode") | ![]() | ![]() |
n = 7; // a n x n interpolation grid x = linspace(0,2*%pi,n); y = x; z = cos(x')*cos(y); C = splin2d(x, y, z, "periodic"); // now evaluate on a bigger domain than [0,2pi]x [0,2pi] m = 80; // discretization parameter of the evaluation grid xx = linspace(-0.5*%pi,2.5*%pi,m); yy = xx; [XX,YY] = ndgrid(xx,yy); zz3 = interp2d(XX,YY, x, y, C, "periodic"); plot3d(xx, yy, zz3, flag=[2 6 4]) xtitle("extrapolation with the periodic outmode") | ![]() | ![]() |
n = 7; // a n x n interpolation grid x = linspace(0,2*%pi,n); y = x; z = cos(x')*cos(y); C = splin2d(x, y, z, "periodic"); // now evaluate on a bigger domain than [0,2pi]x [0,2pi] m = 80; // discretization parameter of the evaluation grid xx = linspace(-0.5*%pi,2.5*%pi,m); yy = xx; [XX,YY] = ndgrid(xx,yy); zz4 = interp2d(XX,YY, x, y, C, "natural"); plot3d(xx, yy, zz4, flag=[2 6 4]) xtitle("extrapolation with the natural outmode") | ![]() | ![]() |
Version | Description |
5.4.0 | 以前では, 入力引数の虚部は暗黙のうちに無視されていました. |