数学建模～插值与拟合概要.ppt

已知n个节点 其中 互不相同， 构造一个二元函数 通过全部已知节点,即 再用 计算插值，即 返回 要求x0,y0单调；x，y可取为矩阵，或x取行向量，y取为列向量，x,y的值分别不能超出x0,y0的范围． z=interp2(x0,y0,z0,x,y,’method’) 被插值点 插值方法 用MATLAB作网格节点数据的插值 插值节点 被插值点的函数值 ‘nearest’ 最邻近插值； ‘linear’ 双线性插值； ‘cubic’ 双三次插值； 缺省时 双线性插值. 例：测得平板表面3×5网格点处的温度分别为： 82 81 80 82 84 79 63 61 65 81 84 84 82 85 86 试作出平板表面的温度分布曲面z=f(x,y)的图形． 输入以下命令： x=1:5; y=1:3; temps=[82 81 80 82 84;79 63 61 65 81;84 84 82 85 86]; mesh(x,y,temps) 1.先在三维坐标画出原始数据，画出粗糙的温度分布曲线图. 2．以平滑数据,在 x、y方向上每隔0.2个单位的地方进行插值. 再输入以下命令: xi=1:0.2:5; yi=1:0.2:3; zi=interp2(x,y,temps,xi,yi,cubic); mesh(xi,yi,zi) 画出插值后的温度分布曲面图. 例 山地高程绘图 要在山区修一条公路， 首先测得一些地点的高程，数据见表. 表中的数据为在平面区域 0x2000, 0y2000 的坐标点上的海拔高程(单位：米） 试给出这个平面区域内地形的模型以便选择公路修建的位置。 0 400 800 1200 1600 2000 0 400 800 1200 1600 2000 370 470 550 600 670 690 510 620 730 800 850 870 650 760 880 970 1020 1050 740 880 1080 1130 1250 1280 830 980 1180 1320 1450 1420 880 1060 1230 1390 1500 1500 二维插值利用Matlab插值程序。 输入已知信息： (xi, yj, zij) x=0:4:20; %给出X轴的坐标 y=0:4:20; %给出Y轴的坐标 z=[37 51 65 74 83 88; 47 62 76 88 98 106; … ; 69 87 105 128 142 150]; %给出（xi，yj）点的高程 zij： [X，Y]=meshgrid(0:1:20,0:1:20); % 给出加密的插值坐标网格 Z=interp2(x,y,z,X,Y,’spline’); %在坐标上进行样条插值 画图： clf;%清空图形坐标系中的内容 mesh(X,Y,Z) %在网格上画出插值的结果 hold on %打开在同一坐标系中画图的功能 contour(X,Y,Z) %画平面等高线 con3=contour3(X,Y,Z) %画三维等高线 clabel(con3) %标高程 hold off %结束作图 插值函数griddata格式为: cz =griddata（x，y，z，cx，cy，‘method’） 用MATLAB作散点数据的插值计算 要求cx取行向量，cy取为列向量． 被插值点 插值方法 插值节点 被插值点的函数值 ‘nearest’最邻近插值 ‘linear’ 双线性插值 ‘