三维曲面插值

Halcom

三维曲面插值：

1. clc,clear,close all
   
2. [x,y] = meshgrid(-4:4);                           %产生已知数据栅格点
   
3. z = peaks(x,y);                                   %计算已知点上的函数值
   
4. [x1,y1] = meshgrid(-4:0.25:4);                %产生更精细的插值栅格点
   
5. z1 = interp2(x,y,z,x1,y1,'nearest');
   
6. subplot(131),surf(x1,y1,z1)                %画基于最邻近法插值的三维表面图
   
7. xlabel('最邻近法插值')
   
8. z2 = interp2(x,y,z,x1,y1,'linear');
   
9. subplot(132),surf(x1,y1,z2)                    %画基于二维分段线性插值的三维表面图
   
10. xlabel('分段线性插值')
    
11. z3 = interp2(x,y,z,x1,y1,'cubic');
    
12. subplot(133),surf(x1,y1,z3)                    %画基于二维三次多项式插值的三维表面图
    
13. xlabel('三次多项式插值')

复制代码