| 查看: 756 | 回复: 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]) |
回复此楼» 猜你喜欢
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
时间戳又变了8-15
已经有15人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有58人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有9人回复
6楼2009-02-17 06:42:15
2楼2008-11-27 09:20:20
3楼2008-11-27 10:25:48
4楼2008-11-28 10:01:38
