| 查看: 765 | 回复: 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]) |
回复此楼» 猜你喜欢
29号明天会评吗
已经有4人回复
-
国社科又开始会评了，不知道这次命运如何
已经有10人回复
-
梦想
已经有4人回复
-
中青基了要发朋友圈吗？
已经有6人回复
-
基金不中，共勉
已经有12人回复
-
能否申诉？
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上意见出来了
已经有19人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
5楼2008-11-28 10:09:18
2楼2008-11-27 09:20:20
3楼2008-11-27 10:25:48
4楼2008-11-28 10:01:38
