| 查看: 705 | 回复: 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]) |
回复此楼» 猜你喜欢
面上本子正文33页，违规吗？会被低分嘛？
已经有12人回复
-
Sci. Bull. 悲剧经验
已经有7人回复
-
提交了我也来说说感想
已经有5人回复
-
今年审到国自然15份，谈谈感受
已经有20人回复
-
国自然评分
已经有4人回复
-
评审有感
已经有22人回复
-
反应很差，大量原料没有反应
已经有6人回复
-
如果工作了想读博，可以边工作边读全日制嘛？
已经有5人回复
-
如何实现卤原子转化
已经有8人回复
-
河北省自然科学基金
已经有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
30