MATLAB插值实验报告--数学实验.doc

文档介绍：新乡学院数学与信息科学系
实验报告
实验项目名称 插值实验
所属课程名称 数学实验
实 验 类 型 综
1200 1100 1550 1600 1550 1380 1070
1200 1100 1350 1450 1200 1150 1010
1390 1500 1500 1400 900 1100 1060 950
1320 1450 1420 1400 1300 700 900 850
1130 1250 1280 1230 1040 900 500 700
y/x
1200 1600 2000 2400 2800 3200 3600 4000
解：原始图程序：
x=1200:400:4000;
y=1200:400:3600;
z=[1480 1500 1550 1510 1430 1300 1200 980;
1500 1550 1600 1550 1600 1600 1600 1550;
1500 1200 1100 1550 1600 1550 1380 1070;
1500 1200 1100 1350 1450 1200 1150 1010;
1390 1500 1500 1400 900 1100 1060 950;
1320 1450 1420 1400 1300 700 900 850;
1130 1250 1280 1230 1040 900 500 700];
mesh(x,y,z);
title('原始图'),xlabel('x'),ylabel('y'),zlabel('z')
　　　　
最邻近插值图：[xi,yi]=meshgrid(1200:10:4000,1200:10:3600);
zi=interp2(x,y,z,xi,yi,'nearest');
mesh(xi,yi,zi)
title('最邻近插值图'),xlabel('x'),ylabel('y'),zlabel('z')

线性插值：[xi,yi]=meshgrid(1200:10:4000,1200:10:3600);
zi=interp2(x,y,z,xi,yi,'linear');
mesh(xi,yi,zi)
title('线性插值'),xlabel('x'),ylabel('y'),zlabel('z')

三次样条插值：[xi,yi]=meshgrid(1200:10:4000,1200:10:3600);
zi=interp2(x,y,z,xi,yi,'spline');
mesh(xi,yi,zi)
title('三次样条插值'),xlabel('x'),ylabel('y'),zlabel('z')

立方差值：[xi,yi]=meshgrid(1200:10:4000,1200:10:3600);
zi=interp2(x,y,z,xi,yi,'cubic');
mesh(xi,yi,zi)
title('立方插值'),xlabel('x'),ylabel('y'),zlabel('z')

等高线图：x=1200:400:4000;
y=1200:400:3600;
z=[1480 1500 1550 1510 1430 1300 1200 980;
1500 1550 1600 1550 1600 1600 1600 1550;
1500 1200 1100 1550 1600 1550 1380 1070;
1500 1200 1100 1350 1450 1200 1150 1010;
1390 1500 1500 1400 900 1100 1060 950;
1320 1450 1420 1400 1300 700 900 850;
1130 1250