数学建模竞赛常用方法之数据处理.ppt

例：测得平板表面3×5网格点处的温度分别为： 82 81 80 82 84 79 63 61 65 81 84 84 82 85 86 试作出平板表面的温度分布曲面z=f(x,y)的图形． 1.先在三维坐标画出原始数据，画出粗糙的温度分布曲线图. 2．以平滑数据,在 x、y方向上每隔0.2个单位的地方进行插值. [x,y] = meshgrid(1:5,1:3); z=[82 81 80 82 84;79 63 61 65 81;84 84 82 85 86]; mesh(x,y,z) xi = 1:0.2:5; yi = 1:0.2:3; zi = interp2(x,y,z,xi',yi) mesh(xi,yi,zi) 通过此例对最近邻点插值、双线性插值方法和双三次插值方法的插值效果进行比较． * 数据的极差归一化变换 * 一、 极差归一化变换公式 * 二、极差归一化变换的MATLAB实现 调用格式： R = rscore(X) [R,xmin,xrange] = rscore(X) […] = rscore(X,dim) 1. 自编的rscore函数 * % 调用rand函数产生一个10行，4列的随机矩阵，每列服从不同的均匀分布 >> x = [rand(10,1), 5*rand(10,1), 10*rand(10,1), 500*rand(10,1)] % 调用rscore函数对x按列进行极差归一化变换， % 返回变换后矩阵R，以及矩阵x各列的最小值构成的向量xmin，各列的极差构成的向量xrange >> [R,xmin,xrange] = rscore(x) 【例3.3-1】调用rand函数产生一个随机矩阵，然后调用rscore函数对其按列进行极差归一化变换 * 调用格式： [Y,PS] = mapminmax(X,YMIN,YMAX) [Y,PS] = mapminmax(X,FP) Y = mapminmax('apply',X,PS) X = mapminmax('reverse',Y,PS) dy_dx = mapminmax('dy_dx',X,Y,PS) dx_dy = mapminmax('dx_dy',X,Y,PS) 2. 自带的mapminmax函数 * % 调用rand函数产生一个10行，4列的随机矩阵，每列服从不同的均匀分布 >> x = [rand(10,1), 5*rand(10,1), 10*rand(10,1), 500*rand(10,1)] % 调用rscore函数对x按列进行极差归一化变换， % 返回变换后矩阵R，以及矩阵x各列的最小值构成的向量xmin，各列的极差构成的向量xrange >> FP.ymin = 0; >> FP.ymax = 1; >> y = mapminmax(x',FP)' 【例3.3-2】调用rand函数产生一个随机矩阵，然后调用mapminmax函数对其按列进行极差归一化变换 数 据 插 值 基本定义： 插值问题：给定一批数据点，需确定满足特定要求的曲线或曲面，如果要求所求曲线（面）通过所有给定的数据点，这就是插值问题。 数据拟合：如果不要求曲线（面）通过所有给定的数据点，而是要求它反映对象整体的变化趋势，可得到更简单实用的近似函数，这就是数据拟合，又称曲线拟合或曲面拟合。 一维插值的定义 已知 n+1个节点 其中 互不相同，不妨设 求任一插值点 处的插值 ? ? ? ? ? 节点可视为由 产生, 表达式复杂, 或无封闭形式, 　或未知. ? 构造一个(相对简单的)函数 通过全部节点, 即 再用 计算插值，即 ? ? ? ? ? ? 返回 1 .拉格朗日(Lagrange)插值 两点一次(线性)插值多项式(一次拉格朗日