matlab第六章 插值问题.ppt

6. 1 引例 82 81 80 82 84 79 63 61 65 81 84 84 82 85 86 试作出平板表面温度分布曲面。 平板表面温度分布曲面图 T: 68, 68, 87, 87 P: 9.798, 13.324, 9.090, 13.355 K: 0.085, 0.090, 0.076, 0.081 6. 2 一维插值 MATLAB常用的一维插值方法 1．最近邻点插值：nearest 2．线性插值：linear 3．三次样条插值：spline 4．三次插值：cubic 最近邻点插值：nearest 函数y= ?(x) 满足条件： ?(xk)= yk =f(xk), k=0,1,…,n 对任意xk ≤ ? ≤ xk+1 ?(?) 线性插值：linear 函数y= ?(x) 满足条件： ?(xk)= yk =f(xk), k=0,1,…,n 对任意xk ≤ ? ≤ xk+1 ?(?) = MATLAB 一维插值命令 interp1 格式：cy=interp1(x,y,cx,method) 其中x, y为已知数据点的横、纵坐标向量， x必须单调； cx为待求插值点的横坐标向量，cx不能超出x的范围； 输出cy为与cx同维数的向量，即插值点处的函数值； method为可选参数，可选值为nearest，linear(默认值)，spline，cubic 求例1中五点半、十一点和十九点的温度 t=0:2:24; T=[12,9,9,10,18,24,28,27,25,20,18,15,13]; x=[5.5,11,19]; y=interp1(t,T,x) y = 9.7500 26.0000 19.0000 例1中24小时环境温度的变化曲线 t=0:2:24; T=[12,9,9,10,18,24,28,27,25,20,18,15,13]; cx=0:1/1440:24; % 间隔一分钟一个插值点 cy1=interp1(t,T,cx); cy2=interp1(t,T,cx,spline); plot (t,T, o,cx,cy1, r ,cx,cy2, k) 24小时环境温度的变化曲线 6. 3 二维插值 二维插值一般可分为两种 网格节点插值 散乱节点插值 1. 插值基点为网格节点时的二维插值 已知二元函数z=f(x,y)在互异的m×n个点(xk,yl)(k=1,2,…,m, l=1,2,…,n)处的函数值zkl，构造函数z=G(x,y)，使其满足条件： G(xk,yl)=zkl k=1,2,…,m, l=1,2, …,n 用G(x,y)作为f(x,y)的近似值，此方法称为二维插值法。 MATLAB二维网格节点插值命令interp2 格式: cz=interp2(x,y,z,cx,cy,method) 其中x, y为已知数据网格点的横、纵坐标向量，其分量必须单调递增; z是m×n维矩阵(m,n分别为向量x, y的维数)，其元素是已知数据网格点相应的函数值; cx，cy是需要插值处网格点的横、纵坐标向量，cx不能超出x的范围，cy不能超出y的范围。cx与cy的方向必须不同(一个行向量,一个列向量)。 输出cz是(cx的维数)×(cy的维数)矩阵 method的取值同interp1 例2 平板表面温度分布曲面的作图程序 x=1:5; y=1:3; z=[82, 81, 80, 82, 84; 79, 63, 61, 65, 81; 84, 84, 82, 85, 86]; cx=1:0.2:5; cy=1:0.2:3; cz=interp2(x,y,z,cx,cy,cubic); mesh(cx,cy,cz) 例2 平板表面温度分布曲面图 2. 插值基点为散乱节点时的二维插值 例3 T=99,P=10.3时导热系数K求解程序 T=[68, 68, 87, 87, 106, 106, 140, 140]; P=[9.798, 13.324, 9.090, 13.355, 9.798, 14.277, 9.656, 12.463]; K=[0.085, 0.090, 0.076, 0.081, 0.070, 0.0750, 0.061, 0.065]; V=griddata(T,P,K,99,10.3) V = 0.0729 已知空间互异的n个点(xk,yk,zk) (k=1,…,n), 求点(x,y)处