计算机仿真技术(MATLAB数据处理(插值、拟合)).pptxVIP

第九章 MATLAB数据处理;第八章 MATLAB数据处理;一、数据拟合;一、数据拟合;一、数据拟合;一、数据拟合;一、数据拟合;一、数据拟合;一、数据拟合;一、数据拟合;一、数据拟合;已知应力-应变关系可以用一条指数曲线来描述, 即假设;在实际应用中常见的拟合曲线有:;2. 非线性曲线拟合: lsqcurvefit.;例4 已知观测数据点如表所示;例4 已知观测数据点如表所示;编写下面的程序调用拟合函数.;例4 已知观测数据点如表所示;f(x)= 3ex+ 4.03x2 + 0.94 x3.;练习:;二、数据插值;已知 n+1个节点;1.分段线性插值;输入: 节点x0, y0, 插值点x (均为数组,长度自定义); 输出: 插值y (与x同长度数组).;例 5 对 在[-1, 1]上, 用n=20的等距分点进行分段线性插值, 绘制 f(x)及插值函数的图形.;例 6 对 在[-5, 5]上, 用n=11个等距分点作分段线性插值和三次样条插值, 用m=21个插值点作图,比较结果.; 0 1.0000 1.0000 1.0000 0.5000 0.8000 0.7500 0.8205 1.0000 0.5000 0.5000 0.5000 1.5000 0.3077 0.3500 0.2973 2.0000 0.2000 0.2000 0.2000 2.5000 0.1379 0.1500 0.1401 3.0000 0.1000 0.1000 0.1000 3.5000 0.0755 0.0794 0.0745 4.0000 0.0588 0.0588 0.0588 4.5000 0.0471 0.0486 0.0484 5.0000 0.0385 0.0385 0.0385;解 在命令窗口输入:;例 8 在飞机的机翼加工时, 由于机翼尺寸很大, 通常在图纸上只能标出部分关键点的数据. 某型号飞机的机翼上缘轮廓线的部分数据如下:;例 8 在飞机的机翼加工时, 由于机翼尺寸很大, 通常在图纸上只能标出部分关键点的数据. 某型号飞机的机翼上缘轮廓线的部分数据如下:;例9 天文学家在1914年8月份的7次观测中, 测得地球与金星之间距离(单位: m), 并取其常用对数值与日期的一组历史数据如下所示, 试推断何时金星与地球的距离(单位: m)的对数值为 9.9352.;解 由于对数值 9.9352 位于 24 和 26 两天所对应的对数值之间, 所以对上述数据用三次样条插值加细为步长为1的数据:;练习:;练习:;（1）打开曲线拟合工具界面 通过cftool命令打开曲线拟合工具界面 ;Data按钮：可输出、查看和平滑数据； Fitting按钮：可拟合数据、比较拟合曲线和数据集； Exclude按钮：可以从拟合曲线中排除特殊的数据点； Ploting按钮：在选定区间后，单击按钮，可以显示拟合曲线和数据集； Analysis按钮：可以做内插法、外推法、微分或积分拟合。; 在输入数据之前，数据变量必须存在于 matlab的工作区间。可以通过load命令输 入变量。单击曲线拟合工具界面中的Data 按钮，打开Data对话框，在对话框中进行设 置，可以输入数据。 ;包括两个选项卡：Data Sets 和 Smooth. Data Sets选项卡： .Import workspace vectors 把向量输 入工作区，要注意的是变量必须具有相同的 维数，无穷大的值和不定值被忽略。 X data 用于选择观测数据 Y data 用于选择X的响应数据 Weight 用于选择权重，与响应数据相联系的向量，如果没选择，默认值为1.;三、拟合工具箱;census 有两个变量：cdate和pop。 cdate是一个年向量，包括1790-1990年， pop是对应年份的美国人口。 whos -file census Name Size Bytes Class Attributes cdate 21x1 168 double pop 21x1 168 double