MATLAB程序 插值

例1 用不同插值方法对一维数据进行插值。

程序代码

%interp1_example.m

%用不同插值方法对一维数据进行插值，并比较其不同

x=0:1.2:10;

y=sin(x);

xi=0:0.1:10;

yi_nearest=interp1(x,y,xi,'nearst');%最邻近插值

yi_linear=interp1(x,y,xi);% 默认插值方法线性插值

yi_spline=interp1(x,y,xi,'spline');% 三次样条插值

yi_cubic=interp1(x,y,xi,'cubic');% 三次多项式插值

yi_v5cubic=interp1(x,y,xi,'v5cubic'); %MATLAB5中使用的三次多项式插值hold on;

subplot(2,3,1);

plot(x,y,'ro',xi,yi_nearest,'b-');

title('最邻近插值');

subplot(2,3,2);

plot(x,y,'ro',xi,yi_linear,'b-');

title('线性插值');

subplot(2,3,3);

plot(x,y,'ro',xi,yi_spline,'b-');

title('三次样条插值');

subplot(2,3,4);

plot(x,y,'ro',xi,yi_cubic,'b-');

title('三次多项式插值');

subplot(2,3,5);

plot(x,y,'ro',xi,yi_v5cubic,'b-');

title('三次多项式插值（MATLAB5）');

运行结果

例2 利用一维快速傅里叶插值实现数据增采样程序代码

%interpt_example.m

%一维快速傅里叶插值实现数据增采样

x=0:1.2:10;

y=sin(x);

n=2*length(x);%增采样1倍

yi=interpft(y,n);%一维快速傅里叶插值

xi=0:0.6:10.4;

hold on;

plot(x,y,'ro');

plot(xi,yi,'b.-');

title('一维快速傅里叶插值');

legend('原始数据','插值结果');

运行结果

例3 举例比较各种二维插值插值算法的不同

程序代码

%interp2_example2.m

%采用二次插值对三维高斯型分布函数进行插值

[x,y]=meshgrid(-3:0.8:3)%原始数据

z=peaks(x,y);

[xi,yi]=meshgrid(-3:0.25:3);

zi_nearest=interp2(x,y,z,xi,yi,'nearst');%最邻近插值zi_linear=interp2(x,y,z,xi,yi);%默认插值方法线性插值zi_spline=interp2(x,y,z,xi,yi,'spline');%三次样条插值zi_cubic=interp2(x,y,z,xi,yi,'cubic');%三次多项式插值hold on;

subplot(2,3,1);

surf(x,y,z);

title('原始数据');

subplot(2,3,2);

surf(xi,yi,zi_nearest);

title('最邻近插值');

subplot(2,3,3);

surf(xi,yi,zi_linear);

title('线性插值');

subplot(2,3,4);

surf(xi,yi,zi_spline);

title('三次样条插值');

subplot(2,3,5);

surf(xi,yi,zi_cubic);

title('三次多项式插值');

figure;

subplot(2,2,1);

contour(xi,yi,zi_nearest);

title('最邻近插值');

subplot(2,2,2);

contour(xi,yi,zi_linear);

title('线性插值');

subplot(2,2,3);

contour(xi,yi,zi_spline);

title('三次样条插值');

subplot(2,2,4);

contour(xi,yi,zi_cubic);

title('三次多项式插值');

运行结果

二维插值

二维插值后的等高线