| 查看: 651 | 回复: 8 | ||
| 【奖励】 本帖被评价3次，作者woshilsh增加金币 2.5 个 | ||
| 当前主题已经存档。 | ||
[资源]
【原创分享】曲线拟合简单操作！
|
||
|
%只考虑线性拟合 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %原始数据 t = [0 .3 .8 1.1 1.6 2.3]'; y = [0.5 0.82 1.14 1.25 1.35 1.40]'; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %多项式拟合 p=polyfit(t,y,2) %利用左除 X = [ones(size(t)) t t.^2]; a = X\y %regress函数 X = [ones(size(t)) t t.^2]; b=regress(y,X) %lsqcurvefit函数 fun=inline('x(1)*t.^2+x(2)*t+x(3)','x','t'); x=lsqcurvefit(fun,[0,0,0],t,y) %Curve Fitting Toolbox fit1= fit(t,y,'poly2') %Curve Fitting Toolbox(自定义多项式) mymodel = fittype('a*t^2+b*t+c','independent','t'); %mymodel = fittype('a*x^2+b*x+c'); fit1= fit(t,y,mymodel,'start',[0,0,0]) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %指数形式的拟合 X = [ones(size(t)) exp(-t) t.*exp(-t)]; a = X\y %lsqcurvefit函数 fun=inline('x(1)+x(2)*exp(-t)+x(3).*t.*exp(-t)','x','t'); x=lsqcurvefit(fun,[0,0,0],t,y) %Curve Fitting Toolbox mymodel = fittype('a+b*exp(-t)+c*t*exp(-t)','independent','t'); %mymodel = fittype('a+b*exp(-x)+c*x*exp(-x)'); fit1= fit(t,y,mymodel,'start',[0,0,0]) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %不含常数项的多项式拟合 %利用左除 X = [t t.^2]; a = X\y %regress函数 X = [t t.^2]; b=regress(y,X) %lsqcurvefit函数 fun=inline('x(1)*t.^2+x(2)*t','x','t'); x=lsqcurvefit(fun,[0,0],t,y) %Curve Fitting Toolbox mymodel = fittype('a*t^2+b*t','independent','t'); %mymodel = fittype('a*x^2+b*x'); fit1= fit(t,y,mymodel,'start',[0,0]) |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
2楼2008-11-27 09:20:20
3楼2008-11-27 10:25:48
4楼2008-11-28 10:01:38
5楼2008-11-28 10:09:18
6楼2009-02-17 06:42:15