永发信息网

以下数据请帮忙用逻辑斯蒂拟合,需要MATLAB的代码。验证后给分。

答案:2  悬赏:20  手机版
解决时间 2021-04-04 18:00
数据如下:
x=[1900:10:1990]';
y=[76.0 92.0 106.5 123.2 131.7 150.7 179.3 204.0 226.5 251.4]';

我仿照别人的拟合方式写了如下代码:
function logistic
tdata=[1900:10:1990]';
cdata=[76.0 92.0 106.5 123.2 131.7 150.7 179.3 204.0 226.5 251.4]';
x0=[75;0.0033];
x=lsqcurvefit(@curvefun1,x0,tdata,cdata);
f=curvefun1(x,tdata);
plot(tdata,cdata,'r+',tdata,f);
function f = curvefun1(x,t)
y0=76;
f=x(1)./(1+(x(1)/y0-1)*exp(-t*x(2))); %其中x(1)=Xm;x(2)=r
但是拟合效果很差。其中x0对曲线的形状影响很大,请问x0应该怎么定。
命令窗口出现如下提示:
Local minimum found.

Optimization completed because the size of the gradient is less than
the default value of the function tolerance.


请问具体怎么改可以把拟合曲线改好。
搞定再加分!!!
最佳答案
你这个本来不合适用logistic
你自己画图看看:
x=[1900:10:1990]';
y=[76.0 92.0 106.5 123.2 131.7 150.7 179.3 204.0 226.5 251.4]';
plot(x,y)
这个用logistic还不如用直线拟合
全部回答
这样,在一个程序里: function hh tdata=[1 2 3 4 5 6 7 8 9 10 11 12]; cdata=[3.1 3.6 4.2 4.8 6.7 9.1 13 18.9 29 51.1 78.2 113]; x0=[300,0.6]; x=lsqcurvefit(@curvefun1,x0,tdata,cdata) f=curvefun1(x,tdata) plot(tdata,cdata,'o',tdata,f) function f = curvefun1(x,t) a=3.1 f=x(1)./(1+(x(1)/a-1)*exp(-t*x(2))) %其中x(1)=xm;x(2)=r 如果分开: function f = curvefun1(x,t) a=3.1 f=x(1)./(1+(x(1)/a-1)*exp(-t*x(2))) %其中x(1)=xm;x(2)=r 主程序: tdata=[1 2 3 4 5 6 7 8 9 10 11 12]; cdata=[3.1 3.6 4.2 4.8 6.7 9.1 13 18.9 29 51.1 78.2 113]; x0=[300,0.6]; x=lsqcurvefit('curvefun1',x0,tdata,cdata) f=curvefun1(x,tdata) plot(tdata,cdata,'o',tdata,f)
我要举报
如以上问答信息为低俗、色情、不良、暴力、侵权、涉及违法等信息,可以点下面链接进行举报!
大家都在看
从你的全世界路过诗词,贾宝玉诗有哪些
耐克羽毛球拍价钱
绿源简易款电动车有保险丝嘛?
三生有幸遇见你,一世热忱只为君!这句话是什
对于分式3x-1÷x-a,当x=a时,以下结论中正确
平遥哪里有卖散一口想牛肉冠云的
青岛哪里有卖装喜糖的袋子 在李沧区李村附近
梦幻西游144级DT三修20.飞升后修炼变成多少?
我要写一篇作文,关于追求成功的人心中充满阳
学习数控模具,以后可以做什么工作?以后的发展
来临之际是什么意思
无枭印是什么意思啊!
我现在生活的城市,租房很贵,吃饭也很贵,一
()修造长城也具有保卫人们安居乐业的目的,()
感悟生命作文评语,作文 生命的感悟 家长评语
推荐资讯
笑曰的曰是什么意思
近几年,我国企业海外“并购”、“收购”案例
阅读下面文字,完成1—4题。 困境也是
为什么添加或删除组件删不了maya 2010?
化学物品的读音问题
韩国人是不是都整容,为什么极限挑战张艺兴的
宝鸡石油公司原来是一家具有优良传统的国有企
z282列车三号车厢让坐吗?听说是餐厅
摸不透的“透”字是什么意思?急急急!
纳兰家的兴衰
开店必须请财神?郁闷中——
工程施工业当月工程结算收入如何确认?
正方形一边上任一点到这个正方形两条对角线的
阴历怎么看 ?