数值分析上机课资料.doc

请先看一下数模教材 306页及317页内容 一维插值 基本格式为： interpl (x9y,cx/method，) 其中X, y分别表示为数据点的横、纵坐标向量，x必须单调。 CX为需要插值的横坐标数据(或数组)，CX不能超过X的范围。 method有以下四个值中情况：▲ nearesf——最邻近点插值 4inear, 线性插值 spline, 三次样 条插值 ? cubic——三次插值 例1：在12h内，每隔lh测量一次温度，温度依次为：5, 8, 9, 15, 25, 29, 3b 30, 22, 25, 27, 24o 试估计在 3.2h, 6.5h, 7.1h时的温度值。 解：输入命令： hours=l:12; temps=[5 8 9 15 25 29 31 30 22 25 27 24]; q=[3.2 6.5 7.1 11.7]; t=interp 1 (hours,temps,q) T= interp 1 (hours,temps,q, spline) 例2：在例1的条件下，每隔O.lh估计一次温度。 解：输入命令： hours=l:12; temps=[5 8 9 15 25 29 31 30 22 25 27 24]; h=l:0.1:12; t=interp 1 (hours,temps,h, spline1); plot(hours,temps, +\h,t,hours,temps, T:) xlabelChouF),ylabel(degi*ees Celsius) 二维插值 1?插值基点为网格节点 问题：l_A 知 mxn 个节点：(石,儿，命)，(i = j = 1,2,...,〃)且 x2 ...xm;yl y2 ... yn。求点(%*,/)(工(兀,牙))处的插值/。 Matlab提供了插值命令为: cz=interp2(x,y,z,cx,cy, Method*) 其中向量兀和y是自变量，八y的分量值必须是单调递增的。 *是肌维向量，指明所给数据网格点的横坐标；y是“维向量，指 明所给数据网格点的纵坐标；乙是加xn阶矩阵，标明相应于所给 数据网格点的函数值。 向量gey是给定的网格点的横坐标和纵坐标，指明函数 cz=interp2(x,y,z,cx,cy,cz, method)返回在网格(cx.cy)处的函数值, 欽与◎应是方向不同的向量，即一个是行向量，另一个是列向量。 Method可选参数如下： inearest, 最邻近点插值 ninear9……线性插值(缺省值) Mpline，三次样条插值 cubic，——三次插值 例3：测得平板表面3x5网格点处的温度分别为: 82 81 80 82 84 79 63 61 65 81 84 84 82 85 86 做出平板表面温度分布曲面。 解：(1)先在三维坐标画出原始数据，看一下该数据的粗糙 程度。输入命令如下： x=l:5; y=i：3； temps=[82 81 80 82 84;79 63 61 65 81;84 84 82 85 86]; mesh(x,y,temps) (2)在八y方向上每隔0.2个单位的地方进行插值，以平滑数 据。再输入以下命令： xi=l:0.01:5; yi=l:0.01:3; zi=interp2(x,y,temps,xi,yiVcubicf); !注：书中可能有误，去掉一个。 mesh(xi,yi,zi) TOC \o 1-5 \h \z 2.插值基点为散乱节点 4 问题：已知72个节点：(兀,)2),(心1,2, ???,〃)， O O 求点(门)「)0(无,兀))处的插值八 O ? ° ° Matlab提供了插值函数为： °??° cz=griddata(x,y,z,cx9cy9 method ) 英中*、y、z均是维向量，指明所给数据点的横坐标、纵 坐标和竖坐标。向量cx,cy是给定的网格点的横坐标和纵坐标，指 明函数 cz= griddata(x9y,z,cx9cy/method,)返回在网格(o)处的 函数值，与刁应是方向不同的向量，即一个是行向量，另一个 是列向量。 Method可选参数如下： nearest9 最邻近点插值 6linear, 线性插值(缺省值) cubic 三次插值 matlab中提供的插值方法 例5：在某海域测得一些点gy)处的水深込由下表给出，在矩形 区域(75,200)x(-50,150)内画出海底曲面的图形。 X 129 140 103.5 88 185.5 195 105 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 -3