数学实验-二维插值.pptVIP

* 下一页 上一页 主 页 Gongqu 数学实验之四 ——（2）二维插值 * 实验目的 [1] 了解二维插值的基本原理，了解三种网格节点数据的插值方法的基本思想； [2] 掌握用MATLAB计算二维插值的方法，并对结果作初步分析； [3] 通过实例学习如何用插值方法解决实际问题。 * 二维插值主要内容 范例2：船在该海域会搁浅吗？ 布置实验 引例1,引例2 二维插值的基本原理 二维插值方法 用MATLAB作二维插值计算 范例1：山区的山形地貌 返回 * 某山区的地形图 * * 引例1：如何绘制山区地貌图 要在某山区方圆大约27平方公里范围内修建一条公路，从山脚出发经过一个居民区，再到达一个矿区。横向纵向分别每隔400米测量一次，得到一些地点的高程：(平面区域0=x=5600,0=y=4800)，首先需作出该山区的地貌图和等高线图。 * 引例1：如何绘制山区地貌图 * * ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? x y 引例1：如何绘制山区地貌图 * 引例2：船在该海域会搁浅吗？ 在某海域测得一些点(x,y)处的水深z由下表给出，船的吃水深度为5英尺，在矩形区域（75，200）*（-50，150）里的哪些地方船要避免进入。 * ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? y x 引例2：船在该海域会搁浅吗？ 返 回 * 二维插值的提法 已知 m?n个节点 (xi, xj, zij) ( i=1, 2, …,m; j=1, 2, …, n ) 其中xi, yj互不相同，不妨设 a=x1x2…xm=b c=y1y2…yn=d 求任一插值点(x*, y*) ( ≠(xi, yj) )处的插值Z* 第一种（网格节点） 插值的基本原理 * 二维插值的提法 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? x y x1 x2 x3 x4 x5 y1 y2 y3 ? (x*, y*) * 第二种（散乱节点） 已知n个节点 其中 互不相同， 求任一插值点 处的插值 二维插值的提法 * 二维插值方法的基本思路 构造一个二元函数 通过全部已知节点,即 再用 计算插值，即 或 返 回 * 二维插值方法 网格节点插值法： 最邻近插值； 分片线性插值; 双线性插值； 双三次插值。 散点数据插值法： 修正Shephard法 返 回 * 注意：最邻近插值一般不连续。具有连续性的最简单的插值是分片线性插值。 1．最邻近插值 x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? y (x1, y1) (x1, y2) (x2, y1) (x2, y2) 返 回 * x y ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? (xi, yj) (xi, yj+1) (xi+1, yj) (xi+1, yj+1) 2．分片线性插值 * 四个插值点（矩形的四个顶点）处的函数值： f (xi, yj)=f1， f (xi+1, yj)=f2， f (xi+1, yj+1)=f3， f (xi, yj+1)=f4 第一片(下三角形区域)： 插值函数为： 2．分片线性插值 * 第二片(上三角形区域)： 插值函数为： 注意：(x, y)当然应该是在插值节点所形成的矩形区域内。显然，分片线性插值函数是连续的； 2．分片线性插值 返 回 * 双线性插值是一片一片的空间二次曲面构成。 插值函数的形式如下： 3．双线性插值 x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? y (x1, y1) (x1, y2) (x2, y1) (x2, y2) 返 回 * 用MATLAB作插值计算 网格节点的插值计算； 散点数据的插值计算； 用MATLAB作插值计算小结 返 回 * 例：测得平板表面3*5网格点处的温度分别为： 82 81 80 82 84