0001 function zx=LSM_fx_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 M1= 1/2*[ -1 0 1 0;
0022 4 -10 8 -2;
0023 -3 9 -9 3;
0024 0 0 0 0];
0025 Fij = [f(i-1,j-1), f(i-1,j), f(i-1,j+1), f(i-1,j+2);...
0026 f(i ,j-1), f(i ,j), f(i ,j+1), f(i ,j+2);...
0027 f(i+1,j-1), f(i+1,j), f(i+1,j+1), f(i+1,j+2);...
0028 f(i+2,j-1), f(i+2,j), f(i+2,j+1), f(i+2,j+2);...
0029 ];
0030 zx=p'*M1*Fij*M'*q;
0031 end