实验一用Matlab数据拟合..docx

用Matlab进行数据拟合 已知观测数据点如表所示-0.447 I.97X 3.2? 6.16 7.08 7.34 7.66 9.569.48分别用3次和6 已知观测数据点如表所示 -0.447 I.97X 3.2? 6.16 7.08 7.34 7.66 9.56 9.48 分别用3次和6次多项式曲线拟合这些数据点. p=polyfit(x,y,m) 其中,x, y为已知数据点向量，分别表示横，纵坐 标,m为拟合多项式的次数，结果返回m次拟合 多项式系数，从高次到低次存放在向量p中. y()=polyval(p,x()) 可求得多项式在xO处的值yO? 编写Matlab程序如下: x=0:0.1:l y=[-0.447,l .9783.28,6.16,7.07.34,7.66,9.56,9.49?3,11 ?2] plot(x,y/k.\f markersize \25) axis([O 1.3 -2 16]) p3=polyfit(x,y3) p6=polyfit(xty,6) 池u t20 02 04 06 08O:OJ:1 池 u t2 0 02 04 06 08 -0.447,1.9783.28,6.16,7.08,7.34,7.66,956,9.48,93,11.2] pRJt(x,y/k//markersize\25) axis([0 1.3 -2 16]) p3=polyfit(x,y3) p6=polyfit(x,y\6) t=0:0.1:1.2 s=polyval(p34) sl=polyval(p64) hold on piot(t^/r-\*linewid plot(t^/h--\linewi grid 官側2用切削机床进行金属品加工时，为了适当地调整 ,需要测定刀具的磨损速度?在一定的时间测量刀 N厚度?得数据如表所示： 切削时间t/h 0 1 2 3 4 5 6 7 8 刀具厚度y/cm 30.0 29」 28.4 28.1 28.0 27.7 27.5 27.2 27.0 切削时间t/h 9 10 11 12 13 14 15 16 刀具厚度y/cm 26.? 26.5 263 26」 25.7 253 24.X 24.0 描出散点图?在命令窗口输入: t=[0:l:16] y=[300 29.1 28.4 28.1 28.0 27.7 27.5 27.2 27.0 26.8 26.5 26.3 26A 25.7 253 24.8 24.0] 描出散点图，在命令窗口输入: 久 y=[30.0 29.1 28.4 21 28.0 2 % 26.5 26.3 26.1 25.7 plot 仏y/材) a=polyfit(t,yj) -03012 293804 hold on yl =-03() 12*1+293804 plot仏 yl), hold off 例2用切削机床进行金属品加工时，为了适当地调整 蠻.需要测定刀具的磨损速度?在一定的时间测量刀 天鮒厚度?得数据如表所示： 切削时间t/h 0 1 2 3 4 5 6 7 8 刀具厚度y/cm 30.0 29.1 2」 28」 28.0 27.7 27.5 27.2 27.0 切削时间t/h 9 10 11 12 13 14 15 16 刀具厚度y/cm 26.S 26.5 26.3 26」 25.7 25.3 24.8 24.0 拟合曲线为：y=?0?3012/+29?3804 刁M3 —个15.4cm X 30.48cm的混凝土柱在加压实验中的 ?应变关系测试点的数据如表所示 1.55 2.47 2.93 3.03 2.89 ￡ 500x10^ 1000x10^ 150()x10^ 2(MM)xl0* 2375x10* o7g/(N/m：) 3.103xl(f 2.465xl(T 1953x1 炉 1317x10 1.219x10 已知应力?应变关系可以用一条指数曲线来描述,即假设 式中，b表示应力,单位是N/nP; ￡表示应变. 1 式中，b表示应力,单位是N/n?; ￡表示应变. 解选取指数函数作拟合时?在拟合前需作变量代换. 化为k^k2的线性函数. 于是.In —= 111^! -k^c ￡ 令 z = In —a0 =-k“ a{ =ln[ € 即 Z =4声+ 4[ 令窗口输入:plot(xw*xyl/2) 令窗口输入: plot(xw*xyl/2) *1.0e-6 1000*1.0e-6 1500?1.0e-6 2000*1.0e-6 2375*1.0e-6] 5\y=|3.1()3*1.0e+3 2.465* 1.0e+3 1.953*1.0e+3 1.517*1.0e+3 *4.219* 1.0e+3] z=log(y) a