interp2innongrid.asv

基于matlab编程

function interp2z = Interp2InNonGrid(a, b)

%step 1
db_yoldmatrix = [

   -2.9295   -2.6113   -2.7456   -2.6449   -2.9099   -2.7162   -2.7606
   -2.1809   -0.9963   -1.6443   -0.7211   -2.6984   -0.3680   -1.0798
   -1.8052   -0.4590   -0.4026   -0.2334   -1.1700   -0.0101   -0.5283
   -1.2936   -0.1847    0.1789    0.4070   -0.5078    0.8609   -0.3604
   -0.2077    0.0931    0.4788    1.0851    0.6172    1.7332    1.3598
    0.9687    0.4968    0.8432    1.7000    1.6077    2.8251    1.4674
    2.8634    2.9300    2.9622    2.7653    2.8462    2.9405    2.3923 ] ;


%step 2
db_newcoly =  [-2.5, -2.0, -1.5, -1.0, -0.5, 0.0, 0.5, 1.0, 1.5, 2.0, 2.2 ] ;
db_oldvectorx= [ -2.19736, -1.75720, -1.24647, 0.77271, 1.09999, 2.03543, 2.60028 ] ;

%step 3
db_newmatrixz = zeros(11, 7) ;

for i = 1:7
db_newmatrixz(:,i) = interp1( db_yoldmatrix(:,i), db_zoldmatrix(:,i), db_newcoly ) ;
db_newmatrixy(:,i) = db_newcoly ;
end

for i= 1:11 ;
db_newmatrixx(i ,:) = db_oldvectorx ;
end

interp2z = interp2(db_newmatrixx, db_newmatrixy, db_newmatrixz, a, b, 'cubic') ;