MATLAB在处理数据中的应用.pdfVIP

冶金仿真软件应用 MATLAB在处理数据中的应用 概 述 在实际的试验和工程测量中，通过不同的方式测得 的一些离散数据点，称为采样点。对这些数据进行利用 之前，还要对其进行分析与处理，如剔除误差较大的或 明显不正确的点，以提高数据的准确性。 有时，由于条件限制不能通过现有的测量手段得到 希望的数据点，则可以通过测量其他的量，并对这一量 的测量数据进行运算，便可间接得到所希望的数据量。 冶金仿真软件应用 2007年 秋 本章主要内容 4.1 数据插值* 4.2 曲线拟合* 4.3多项式求根 4.4 线性规划 冶金仿真软件应用 2007年 秋 4.1 数据插值 在流动、传热及燃烧等问题中运用数值计算进行求解 时通常也是得到的一些不连续的点，如果要得到这些分散 点以外的其他地方的数值，就可以运用这些已知的点进行 插值。插值可以分为一维插值、二维插值和多维插值。 冶金仿真软件应用 2007年 秋 一维数据插值 在MATLAB中，实现这些插值的函数是interp1，其调用 格式为： Y1=interp1(X,Y,X1,method) 函数根据X,Y的值，计算函数在X1处的值。X,Y是两个等 长的已知向量，分别描述采样点和样本值，X1是一个向量 或标量，描述欲插值的点，Y1是一个与X1等长的插值结果。 method是插值方法，允许的取值有‘linear’、‘nearest’、 ‘cubic’、‘spline’ 。 冶金仿真软件应用 2007年 秋 例4-1 某观测站测得某日6:00时至18:00时之间每隔2小时的 室内外温度(℃)， h 6 8 10 12 14 16 18 t 内 18 20 22 25 30 28 24 t外 15 19 24 28 34 32 30 （1）求该日室内外11：30时的近似温度(℃) 。 （2 ）求该日室内外6:30至17:30时之间每隔2小时各点 的近似温度(℃) 。 冶金仿真软件应用 2007年 秋 解：设时间变量h为一行向量，温度变量t为一个两列 矩阵，其中第一列存放室内温度，第二列储存室外温度。 命令如下： （1） h =6:2:18; t=[18,20,22,25,30,28,24;15,19,24,28,34,32,30]; X1 =11.5 Y1=interp1(h,t,X1,‘spline’) % ‘linear’、‘nearest’、 %‘cubic’、‘spline’ （2 ） X2 =6.5:2:15.5 Y2=interp1(h,t,X2,spline) %用3次样条插值计算 冶金仿真软件应用 2007年 秋