查看: 2595|回复: 2

求拟合代码

字体大小: 正常放大

1 主题	6 听众	167 积分

升级 33.5%

TA的每日心情

	开心 2013-10-21 00:14

签到天数: 19 天

[LV.4]偶尔看看III

2013挑战赛参赛者

群组: 2013认证赛B题讨论群组

群组: 2013认证赛A题讨论群组

电梯直达

1^#

发表于 2013-5-18 13:34 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

t=[1,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,31,32,35,37,39,47];
x=[3,7,8,12,15,16,19,22,26,31,36,41,47,58,61,75,80,85,89,90,94,102,102,108,112,114,116,121,124,126,128,128,130,131];
plot(t,x,'*')
b0=[0.46,5.3,-0.45];
fun=inline('b(1)/(1+b(2)*exp(b(3)*t))','b','t');
[b,r,j]=nlinfit(t,x,fun,b0);
b
。。
请问大神们，为什么老是出现错误呢，麻烦帮我检查fun=inline('b(1)/(1+b(2)*exp(b(3)*t))','b','t');
这句哪儿出错了，急求。。

zan

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

使用道具举报

madio

3万主题	1312 听众	5万积分

TA的每日心情

	奋斗 2024-7-1 22:21

签到天数: 2014 天

[LV.Master]伴坛终老

自我介绍: 数学中国站长

群组: 数学建模培训课堂1

群组: 数学中国美赛辅助报名

群组: Matlab讨论组

群组: 2013认证赛A题讨论群组

群组: 2013认证赛C题讨论群组

2^#

发表于 2013-8-20 17:49 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

t=[1,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,31,32,35,37,39,47];
1 e0 N U% o\" _3 o( M6 p; f. H
x=[3,7,8,12,15,16,19,22,26,31,36,41,47,58,61,75,80,85,89,90,94,102,102,108,112,114,116,121,124,126,128,128,130,131];
3 g; ~6 w' {7 s3 b j$ W8 l2 b1 ~
plot(t,x,'*');3 I+ h5 f/ U; ?0 O, E
b0=[0.46,5.3,-0.45];9 u$ [\" B2 {7 B, I$ [/ P
fun=inline('b(1)./(1+b(2)*exp(b(3)*t))','b','t');. `. ]1 m2 b$ b0 q
[b,r,j]=nlinfit(t,x,fun,b0);
) p# W7 x\" L) X. e9 V
b