Matlab数据拟合.ppt

;数据拟合;用Matlab进行数据拟合;例1 已知观测数据点如表所示;t=0:0.1:1.2 s=polyval(p3,t) s1=polyval(p6,t) hold on plot(t,s,r-,linewidth,2) plot(t,s,b--,linewidth,2) grid;例2 用切削机床进行金属品加工时, 为了适当地调整机床, 需要测定刀具的磨损速度. 在一定的时间测量刀具的厚度, 得数据如表所示:;解: 描出散点图, 在命令窗口输入:;解: 描出散点图, 在命令窗口输入:;例2 用切削机床进行金属品加工时, 为了适当地调整机床, 需要测定刀具的磨损速度. 在一定的时间测量刀具的厚度, 得数据如表所示:;例3 一个15.4cm×30.48cm的混凝土柱在加压实验中的应力-应变关系测试点的数据如表所示;已知应力-应变关系可以用一条指数曲线来描述, 即假设;在命令窗口输入:;已知应力-应变关系可以用一条指数曲线来描述, 即假设;在实际应用中常见的拟合曲线有:;2. 非线性曲线拟合: lsqcurvefit.;例4 已知观测数据点如表所示;例4 已知观测数据点如表所示;编写下面的程序调用拟合函数.;例4 已知观测数据点如表所示;f(x)= 3ex+ 4.03x2 + 0.94 x3.;练习:;插值问题;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. 设 在区间[-2, 2]上用10等分点作为节点, 分别用三种插值方法:;1. 设 在区间[-2, 2]上用10等分点作为节点, 分别用两种插值方法:;练习:;例5 给药方案;例5 给药方案;例5 给药方案;例5 给药方案;2. 给药方案设计;2. 给药方案设计;2. 由实验数据作曲线拟合以确定参数;问题化为由数据 ti , yi ( i=1,…,8 ) 拟合直线;拟合曲线为:;3. 结论