第4讲插值(简化)精选.ppt

* * 例 在某海域测得一些点(x,y)处的水深z由下表给出，船的吃水深度为5英尺，在矩形区域（75，200）*（-50，150）里的哪些地方船要避免进入。 * 4.作出水深小于5的海域范围,即z=5的等高线. * x=[129 140 103.5 88 185.5 195 105.5 157.5 107.5 77 81 162 162 117.5]; y=[7.5 141.5 23 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.5 84 -33.5]; z=[-4 -8 -6 -8 -6 -8 -8 -9 -9 -8 -8 -9 -4 -9]; cx=75:0.5:200; cy=-70:0.5:150; cz=griddate(x,y,z,cx,cy,cubic); meshz(cx,cy.cz),rotate3d xlabel(x),ylabel(y),zlabel(z); %pause figure(2),contour(ca,cy,cz,[-5 -5]);grid hold on plot(x,y,+) xlabel(x),ylabel(y) * 实验作业 山区地貌：在某山区测得一些地点的高程如下表：(平面区域1200=x=4000,1200=y=3600)，试作出该山区的地貌图和等高线图，并对几种插值方法进行比较。 * * * * * * * * * * * * * * * * xi=linspace(0,5,50);%加密横坐标数据到50个 yi=linspace(0,6,80);%加密纵坐标数据到80个[xii,yii]=meshgrid(xi,yi);%生成新网格数据zii=interp2(x,y,z,xii,yii,cubic);%立方插值 mesh(xii,yii,zii)%加密后的地貌 * hold on?%?保持图形 plot3(x,y,z+0.1,ob)??%原始数据用‘O’绘出 * 已知 m?n个节点 其中 互不相同，不妨设 构造一个二元函数 通过全部已知节点,即 再用 计算插值，即 * 例：用随机数据生成地貌图再进行插值 x=rand(100,1)*4-2; y=rand(100,1)*4-2; z=x.*exp(-x.^2-y.^2); ti=-2:.25:2; [xi,yi]=meshgrid(ti,ti);%加密数据 zi=griddata(x,y,z,xi,yi);%默认线性插值 mesh(xi,yi,zi) hold on plot3(x,y,z,o) ? * * 第二种（散乱节点）： ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? y x 0 * 已知n个节点 其中 互不相同， 构造一个二元函数 通过全部已知节点,即 再用 计算插值，即 * 注意：最邻近插值一般不连续。具有连续性的最简单的插值是分片线性插值。 最邻近插值 x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? y (x1, y1) (x1, y2) (x2, y1) (x2, y2) O 二维或高维情形的最邻近插值，与被插值点最邻近的节点的函数值即为所求。 * 将四个插值点（矩形的四个顶点）处的函数值依次简记为： 分片线性插值 x y ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? (xi, yj) (xi, yj+1) (xi+1, yj) (xi+1, yj+1) O f (xi, yj)=f1，f (xi+1, yj)=f2，f (xi+1, yj+1)=f3，f (xi, yj+1)=f4 * 插值函数为： 第二片(上三角形区域)：(x, y)满足 插值函数为： 注意：(x, y)当然应该是在插值节点所形成的矩形区域内。显然，分片线性插值函数是连续的； 分两片的函数表达式如下： 第一片(下三角形区域)： (x, y)满足 * 双线性插值是一片一片的空间二次曲面构成。 双线性插值函数的形式如下： 其中有四个待定系数，利用该函数在矩形的四个顶点（插值节点）的函数值，得到四个代数方程，正好确定四个系数。 双线性插值 x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? y (x1, y1) (x1, y2) (x2, y1) (x2, y2) O * 要求x0,y0单调；x，y可取为矩阵，或x取行向量，y取为列向量，x,y的值分别不能超出x0,y0的范围。 z=interp2(x0,y0,z0,x,y,’method’) 被插值点 插值方法 用MATLAB作网格节点数据的插值 插值节点 被插值点的函数值 ‘nearest’ 最邻近插值