MATLAB在非线性曲线拟合中的应用研究.docVIP

MATLAB在非线性曲线拟合中的应用小结 摘要:归纳总结了非线性曲线拟合的方法、求解步骤和上机操作过程 关键词:曲线拟合非线性MATLAB 正文： 1．曲线拟合的基本原理 已知一组测定的数据（例如N个点（xi,yi）x和因变量y的一个近似解析表达式y=φ（x）δi=φ（xi）-yi，i=1，2，…N， 为最小，这就是常用的最小二乘法原理。 2 .MATLAB曲线拟合的相关方法 2.1.函数形式： （1）多项式拟合函数polyfit，调用格式为： p=polyfit（x,y,n） x，yn为拟合多项式的次数，函数返回值为拟合多项式的系数（按降幂排列）。n=1时，就为线性拟合。 例1：给出表1数据，试用最小二乘法求一次和二次拟合多项式。 表1 数据 x -1.00 -0.75 -0.50 -0.25 0 0.25 0.50 0.75 1.00 y -0.2209 0.3295 0.8826 1.4392 2.0003 2.5645 3.1334 3.7061 4.2836 在MATLAB命令窗口中输入： clear; close; x=-1:0.25:1; y=[-0.2209,0.3295,0.8826,1.4392,2.0003,2.5645,3.1334,3.7061,4.2836] p1=polyfit(x,y,1) p2=polyfit(x,y,2) y1=polyval(p1,x); y2=polyval(p2,x); plot(x,y,+,x,y1,r:,x,y2,k-.) 运行结果： 拟合多项式为：y*=2.0516+2.0131和y*=0.0313x2+2.2516x+2.20001 （2）非线性数据拟合函数lsqcurvefit调用格式为： c=lsqcurvefi（tfun,x0,xdata,ydata） fun为拟合函数的M-函数文件名，x0为初始向量，xdata,ydata为参与曲线拟合的实验数据。函数返回值c为非线性函数fun的拟合系数。 例2：2004年全国大学生数学建模竞赛C题（酒后驾车）中给出某人在短时间内喝下两瓶啤酒后，间隔一定的时间测量他的血液中酒精含量y（毫克/百毫升），得到数据如表2。 表2 酒精含量与饮酒时间的实验数据 时间（小时） 0.25 0.5 0.75 1 1.5 2 2.5 3 3.5 4 4.5 5 酒精含量 30 68 75 82 82 77 68 68 58 51 50 41 时间（小时） 6 7 8 9 10 11 12 13 14 15 16 酒精含量 38 35 28 25 18 15 12 10 7 7 4 （2） 根据实验数据，利用非线性拟合函数lsqcurvefit，确定模型（2）式中的参数c1，c2，c3M-函数文件Example2_1： function f=Example2_1(c,tdata) f=c(1)*(exp(-c(2)*tdata)-exp(-c(3)*tdata)); clear tdata=[0.25 0.5 0.75 1 1.5 2 2.5 3 3.5 4 4.5 5 6 7 8 9 10 11 12 13 14 15 16]; ydata=[30 68 75 82 82 77 68 68 58 51 50 41 38 35 28 25 8 15 12 10 7 7 4]; c0=[1 1 1]; for i=1:50; c=lsqcurvefit(Example2_1,c0,tdata,ydata); c0=c; end 得到最优解为：c= 117.05,0.1930,1.9546 从而得出拟合曲线： 2.2.图形窗口形式 （1）利用多项式拟合的交互图命令（GUI）polytool，调用格式为：polytoo（lx,y） 其中x,y分别为实验数据构成的向量，例如利用polytool求解例1的MATLAB命令如下： x=-1:0.25:1 y=[-0.2209,0.3295,0.8826,1.4392,2.0003,2.5645,3.1334,3.7061,4.2836]; polytool（x,y） Degree中输入1，点击导出数据Export,出现保存对话框ExporttoWorkspace,选中Parameters（参数）,Residuals（残差）后点击OK，在MATLAB的Workspace窗口中可以看到参数为：2.2516和2.0131，即拟合函数为 Y*=2.2516x+2.0131。 同样如果拟合的函数为二次函数，则只要在Degree中输入2，其它步骤相同，可得拟合函数为： Y*=0.0313x2+2.2516x+2.0001 通过查看Residuals（残差）值，可以发