资源描述：
《曲线拟和函数lsqcurvefit和nlinfit》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、曲线拟和函数lsqcurvefit&nlinfit琢磨了好久matlab自带的曲线拟和工具箱,发现这货只能解决从离散数据得到各种类型的拟和效果,但是反之貌似没法实现,google一下有这两个函数可以用:lsqcurvefit和nlinfitlsqcurvefit(非线性最小二乘法)help了一下,发现官方的文档过于详尽,节选部分吧.LSQCURVEFITsolvesnon-linearleastsquaresproblems.X=LSQCURVEFIT(FUN,X0,XDATA,YDATA)startsatX0an
2、dfindscoefficientsXtobestfitthenonlinearfunctionsinFUNtothedataYDATA(intheleast-squaressense).FUNacceptsinputsXandXDATAandreturnsavector(ormatrix)offunctionvaluesF,whereFisthesamesizeasYDATA,evaluatedatXandXDATA.NOTE:FUNshouldreturnFUN(X,XDATA)andnotthesum-of-
3、squaressum((FUN(X,XDATA)-YDATA).^2).((FUN(X,XDATA)-YDATA)issquaredandsummedimplicitlyinthealgorithm.)这里面对函数的各个参数的意义和用法都做了解释.[x,resnorm]=lsqcurvefit(fun,x0,xdata,ydata,...)fun是我们需要拟合的函数，这是重点x0是我们对函数中各参数的预测值，这也是重点xdata则是横轴坐标的值ydata是纵轴的值需要注意的是lsqcurvefit里面的fun是我们
4、需要拟和的函数,需要另外编写.下面用以前某年的CUMCM真题修改作为示例.--------------------------------------------------------------------------------------------tj/s
5、100
6、200
7、300
8、400
9、500
10、600
11、700
12、800
13、900
14、1000
15、-----------------------------------------------------------------
16、---------------------------Cj/(mg*cm-3)
17、4.54
18、4.99
19、5.35
20、5.65
21、5.90
22、6.10
23、6.26
24、6.39
25、6.50
26、 6.59
27、--------------------------------------------------------------------------------------------在已知部分参数的情况下,求函数求函数的最下值的点(K,a,b) 1.编写M文件(curvefun.m)functionf=curvefun(x
28、,tdata)f=x(1)+x(2)*exp(-0.02*x(3)*tdata);end2.编写程序(test1.m)tdata=linspace(100,1000,10);cdata=1e-05.*[454499535565590610626639650659];x0=[0.20.050.05];x=lsqcurvefit(@curvefun,x0,tdata,cdata);f=curvefun(x,tdata)plot(tdata,cdata,'*',tdata,f,'r-')3.输出结果x=0.0069-0.
29、00290.0809即表示k=0.2542,a=0.0063,b=-0.0034nlinfit从matlab给出的帮助文档来看,nlinfit与lsqcurvefit同属与非线性最小二乘拟和,一般来说都是能得到比较接近的结果.但是由于nlinfit使用的是牛顿方法,在使用是需要给出你和参数的假设初值,有些问题对初值比较敏感,不同的初值会导致差异比较大.下面示例nlinfit的用法:混凝土的抗压强度随着养护时间的延长而增加,现将一批混凝土作成12个试块,记录了养护日期x(日)以及抗压强度y(kg/cm2)的数据---
30、----------------------------------------------------------------------------------------------------养护时间x
31、 2
32、 3
33、 4
34、 5
35、 7
36、 9
37、 12
38、 14
39、 17
40、 21
41、 28
42、 56
43、-----
