matlab插值曲线拟合.ppt

a = polyfit（xdata，ydata，n） 其中n表示多项式的最高阶数 xdata，ydata 为要拟合的数据，它是用向量的方式输入。 输出参数a为拟合多项式 y = a1xn + … + anx + an+1的系数a = [a1, …, an, an+1]。 多项式在x处的值y可用下面程序计算。 y = polyval (a, x) 例a 已知观测数据点如表所示 x y 0 -0.447 0.1 1.978 0.2 3.28 0.3 6.16 0.4 7.08 0.5 7.34 0.6 7.66 0.7 9.56 0.8 9.48 0.9 9.3 1 11.2 分别用3次和6次多项式曲线拟合这些数据点. x=0:0.1:1 y=[-0.447,1.978,3.28,6.16,7.08,7.34,7.66,9.56,9.48,9.3,11.2] plot(x,y,k.,markersize,25) axis([0 1.3 -2 16]) p3=polyfit(x,y,3) p6=polyfit(x,y,6) 编写Matlab程序如下: 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 x=0:0.1:1 y=[-0.447,1.978,3.28,6.16,7.08,7.34,7.66,9.56,9.48,9.3,11.2] plot(x,y,k.,markersize,25) axis([0 1.3 -2 16]) p3=polyfit(x,y,3) p6=polyfit(x,y,6) 例b 用切削机床进行金属品加工时, 为了适当地调整机床, 需要测定刀具的磨损速度. 在一定的时间测量刀具的厚度, 得数据如表所示: 切削时间 t/h 0 30.0 1 29.1 2 28.4 3 28.1 4 28.0 5 27.7 6 27.5 7 27.2 8 27.0 刀具厚度 y/cm 切削时间 t/h 9 26.8 10 26.5 11 26.3 12 26.1 13 25.7 14 25.3 15 24.8 16 24.0 刀具厚度 y/cm 解: 在命令窗口输入: t=[0:1:16] y=[30.0 29.1 28.4 28.1 28.0 27.7 27.5 27.2 27.0 26.8 26.5 26.3 26.1 25.7 25.3 24.8 24.0] plot(t,y,*) 解: 在命令窗口输入: t=[0:1:16] y=[30.0 29.1 28.4 28.1 28.0 27.7 27.5 27.2 27.0 26.8 26.5 26.3 26.1 25.7 25.3 24.8 24.0] plot(t,y,*) a = -0.3012 29.3804 hold on plot(t, y1), hold off a=polyfit(t,y,1) y1=-0.3012*t+29.3804 例b 用切削机床进行金属品加工时, 为了适当地调整机床, 需要测定刀具的磨损速度. 在一定的时间测量刀具的厚度, 得数据如表所示: 切削时间 t/h 0 30.0 1 29.1 2 28.4 3 28.1 4 28.0 5 27.7 6 27.5 7 27.2 8 27.0 刀具厚度 y/cm 切削时间 t/h 9 26.8 10 26.5 11 26.3 12 26.1 13 25.7 14 25.3 15 24.8 16 24.0 刀具厚度 y/cm 拟合曲线为: y=-0.3012t+29.3804 例c 一个15.4cm×30.48cm的混凝土柱在加压实验中的应力-应变关系测试点的数据如表所示 1.55 2.47 2. 93 3. 03 已知应力-应变关系可以用一条指数曲线来描述, 即假设 式中, 表示应力, 单位是 N/m2; 表示应变. 2.89 已知应力-应变关系可以用一条指数曲线来描述, 即假设 式中, 表示应力, 单位是 N/m2; 表示应变. 解 选取指数函数作拟合时, 在拟合前需作变量代换, 化为 k1, k2 的线性函数. 于是, 令 即 在命令窗口输入: x=[500*1.0e-6 1000*1.0e-6 1500*1.0e-6 2000*1.0e-6 2375*1.0e-6] y=[3.103*1.0e+3 2.465*1.0e+3 1.953*1.0e+3 1.517*1.0e+3 1.219*1.0e+3] z=log(y) a=polyfit(x,z,1) k1