数学建模Matlab数据拟合详解.pptVIP

用Matlab進行数据拟合;例1已知观测数据點如表所示;t=0:0.1:1.2

s=polyval(p3,t)

s1=polyval(p6,t)

holdon

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.94x3.;练习:;插值問題;1.分段线性插值;输入:节點x0,y0,插值點x(均為数组,長度自定义);

输出:插值y(与x同長度数组).;例5对在[-1,1]上,用n=20的等距分点进行分段线性插值,绘制f(x)及插值函数的图形.;例6对在[-5,5]上,用n=11个等距分点作分段线性插值和三次样条插值,用m=21个插值点作图,比较结果.;01.00001.00001.0000

0.50000.80000.75000.8205

1.00000.50000.50000.5000

1.50000.30770.35000.2973

2.00000.0.0.

2.50000.13790.15000.1401

3.00000.10000.10000.1000

3.50000.07550.07940.0745

4.00000.05880.05880.0588

4.50000.04710.04860.0484

5.00000.03850.03850.0385;解在命令窗口输入:;例8在飞机的机翼加工時,由于机翼尺寸很大,壹般在图紙上只能標出部分要點的数据.某型号飞机的机翼上缘轮廓线的部分数据如下:;例8在飞机的机翼加工時,由于机翼尺寸很大,壹般在图紙上只能標出部分要點的数据.某型号飞机的机翼上缘轮廓线的部分数据如下:;例9天文學家在198月份的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.結论