查看: 2635|回复: 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];
# C- I& O' C% K. ~, D% n+ I
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];( l( K1 Z! S [# U3 i+ F
plot(t,x,'*');% q/ A- @8 R, M& h% J3 c
b0=[0.46,5.3,-0.45];
6 }: x2 K2 q, r0 `6 C- s& ?
fun=inline('b(1)./(1+b(2)*exp(b(3)*t))','b','t');
$ q$ G2 s \+ D9 ~7 c
[b,r,j]=nlinfit(t,x,fun,b0);
6 h/ O9 ^0 j, c2 J
b