| 查看: 751 | 回复: 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号出结果
已经有4人回复
-
基金中了
已经有15人回复
-
帮忙看看fileCode
已经有8人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有3人回复
-
FileCode能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
静等基金结果
已经有16人回复
4楼2008-11-28 10:01:38
2楼2008-11-27 09:20:20
3楼2008-11-27 10:25:48
5楼2008-11-28 10:09:18
100