0001 function z=LSM_f_cubic_interpolation(x,y,f)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011 i=floor(x);
0012 j=floor(y);
0013 u=x-i;
0014 v=y-j;
0015 p=[1 u u^2 u^3]';
0016 q=[1 v v^2 v^3]';
0017 M= 1/2*[ 0 2 0 0 ;
0018 -1 0 1 0;
0019 2 -5 4 -1;
0020 -1 3 -3 1];
0021 Fij = [f(i-1,j-1), f(i-1,j), f(i-1,j+1), f(i-1,j+2);...
0022 f(i ,j-1), f(i ,j), f(i ,j+1), f(i ,j+2);...
0023 f(i+1,j-1), f(i+1,j), f(i+1,j+1), f(i+1,j+2);...
0024 f(i+2,j-1), f(i+2,j), f(i+2,j+1), f(i+2,j+2);...
0025 ];
0026 z=p'*M*Fij*M'*q;
0027 end