12 / 2 页下一页

查看: 14660|回复: 15

logistic人口模型的求解问题

[复制链接]

字体大小: 正常放大

106347084

2 主题	4 听众	117 积分

升级 8.5%

TA的每日心情

	开心 2016-11-10 11:01

签到天数: 1 天

[LV.1]初来乍到

电梯直达

1^#

发表于 2009-3-25 16:31 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

logistic 方程的系数求解

已知logistic模型为 x(t)=xm/(1+(xm/x0-r)*exp(-r*t))其中XM,X0,R是参数。
已经1978-2006年的人口数据，x=[1:29]
y=[3750.96
3792.33
3826.58
3871.51
3924.32
3963.1
3993.09
4029.56
4070.07
4121.19
4169.85
4208.88
4234.91
4261.37
4285.91
4313.3
4341.2
4369.63
4400.09
4422.28
4446.86
4467.46
4501.22
4519.84
4535.98
4551.58
4577.22
4602.11
4629.43]'希望能拟合出那3个参数？
请教高手求解~~~~

本人试过CFTOOL，lsqcurvefit都拟合不来？(或者说完全失真)
请高手实际能做出后得到答案？最好是能解释下为什么那2个工具不行～～～～

zan

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

使用道具举报

g853513345

0 主题	4 听众	225 积分

升级 62.5%

该用户从未签到

群组: 数学建模

群组: LINGO

2^#

发表于 2009-7-28 20:06 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

高手怎么没有回答呢

我也说一句

发表

使用道具举报

g853513345

0 主题	4 听众	225 积分

升级 62.5%

该用户从未签到

群组: 数学建模

群组: LINGO

3^#

发表于 2009-7-28 20:06 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

我也急需此问题的解决办法

我也说一句

发表

使用道具举报

fc013

7 主题	2 听众	819 积分

升级 54.75%

TA的每日心情

	开心 2015-6-12 19:46

签到天数: 6 天

[LV.2]偶尔看看I

群组: 数学建模

群组: 数学趣味、游戏、IQ等

群组: 学术交流A

4^#

发表于 2009-7-29 20:36 |只看该作者 |招呼Ta 关注Ta

:):victory:

我也说一句

发表

使用道具举报

shuqg02

0 主题	5 听众	62 积分

升级 60%

该用户从未签到

5^#

发表于 2009-8-1 19:09 |只看该作者 |招呼Ta 关注Ta

我也想帮你，可我也不会。

我也说一句

发表

使用道具举报

水木年华zzu

21 主题	7 听众	3435 积分

升级 47.83%

TA的每日心情

	怒 2014-5-25 20:58

签到天数: 20 天

[LV.4]偶尔看看III

群组: Matlab讨论组

群组: 小草的客厅

群组: 数学趣味、游戏、IQ等

群组: C 语言讨论组

群组: 我行我数

6^#

发表于 2009-8-11 16:07 |只看该作者 |招呼Ta 关注Ta

先去对数然后再拟合

我也说一句

发表

使用道具举报

Kadyniost

0 主题	4 听众	187 积分

升级 43.5%

该用户从未签到

7^#

发表于 2009-8-11 20:55 |只看该作者 |招呼Ta 关注Ta

.................。。。

我也说一句

发表

使用道具举报

dcxyfcm

0 主题	0 听众	2 积分

升级 40%

该用户从未签到

8^#

发表于 2010-12-21 21:23 |只看该作者 |招呼Ta 关注Ta

clc;
clear;clf
t=1971:1990;
N=[8.523,8.718,8.921,9.086,9.242,9.372,9.497,9.626,9.754,9.871...
10.007,10.165,10.301,10.436,10.585,10.751,10.930,11.103,11.270,11.433];
M=1./N-1/20;
plot(t,N,'k.','markersize',20);
axis([1971,2010,0,20]);
grid;
hold on
pause(0.5)
n=20;
a=sum(t(1:n));
b=sum(t(1:n).*t(1:n));
c=sum(log(M(1:n)));
d=sum(t(1:n).*log(M(1:n)));
A=[n a;a ,b];
B=[c;d];
p=inv(A)*B;
x=1971:2010;
y=1./(0.05+exp(p(1)+p(2)*x));
plot(x,y,'r-','linewidth',2)
p,y

我也说一句

发表

使用道具举报

zsy312

2 主题	4 听众	64 积分

升级 62.11%

TA的每日心情

	开心 2013-9-23 18:22

签到天数: 4 天

[LV.2]偶尔看看I

9^#

发表于 2010-12-22 23:22 |只看该作者 |招呼Ta 关注Ta

大家以后再发布代码的时候千万用高级模式点代码，比如：

x(t)=xm/(1+(xm/x0-r)*exp(-r*t))

复制代码

不然拷过去就有一堆花乱码，修改很麻烦。

我也说一句

发表

使用道具举报

zsy312

2 主题	4 听众	64 积分

升级 62.11%

TA的每日心情

	开心 2013-9-23 18:22

签到天数: 4 天

[LV.2]偶尔看看I

10^#

发表于 2010-12-22 23:31 |只看该作者 |招呼Ta 关注Ta

修改数据又折腾了一阵子：
y=[3750.96
3792.33
3826.58
3871.51
3924.32
3963.1
3993.09
4029.56
4070.07
4121.19
4169.85
4208.88
4234.91
4261.37
4285.91
4313.3
4341.2
4369.632
4400.09
4422.28
4446.86
4467.462
4501.22
4519.843
4535.98
4551.58
4577.22
4602.114
4629.43]'

我也说一句

发表

使用道具举报

12 / 2 页下一页

返回列表

帐号		密码		只需要一步，快速开始		注册地址	找回密码

logistic人口模型的求解问题

新人进步奖

优秀斑竹奖

QQ

电话咨询

关于我们| 联系我们| 诚征英才| 对外合作| 产品服务|