查看: 2663|回复: 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];
8 j: J# o% L, t* y
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];2 \9 r. N' g0 U2 P
plot(t,x,'*');\" @, X$ P K5 z# A\" m\" K& Z2 B
b0=[0.46,5.3,-0.45];
5 `/ I2 g1 R8 s# j' o
fun=inline('b(1)./(1+b(2)*exp(b(3)*t))','b','t');
& S' a* w$ w7 g ^( V N# i
[b,r,j]=nlinfit(t,x,fun,b0);7 }( s2 I5 B0 u$ B4 C6 | K
b