数学建模社区-数学中国

标题: 求拟合代码 [打印本页]

---

作者: woshiluyu666    时间: 2013-5-18 13:34
标题: 求拟合代码
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];
! ~. u. k' J. e7 r% _( k: d0 @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,'*')
+ a  X8 f% ~/ t: B b0=[0.46,5.3,-0.45];
& s) |) ^4 P! a0 ?fun=inline('b(1)/(1+b(2)*exp(b(3)*t))','b','t');
[b,r,j]=nlinfit(t,x,fun,b0);
9 n1 n9 I6 j& S1 sb
。。
' W' O" ]; E! S0 b请问大神们，为什么老是出现错误呢，麻烦帮我检查fun=inline('b(1)/(1+b(2)*exp(b(3)*t))','b','t');
这句哪儿出错了，急求。。

---

作者: madio    时间: 2013-8-20 17:49

1. 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];
   
2. 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. plot(t,x,'*');
   
4. b0=[0.46,5.3,-0.45];
   ( G# _* ^6 n- j3 _6 k0 D6 t7 S& R
5. fun=inline('b(1)./(1+b(2)*exp(b(3)*t))','b','t');
   ( d3 m% [! z( S2 w' ~
6. [b,r,j]=nlinfit(t,x,fun,b0);
   - t; q$ |) n+ J
7. b

复制代码

改了一下，你试试看

---

作者: 净心、精心    时间: 2013-8-20 20:17
fun=inline('b(1)./(1+b(2)*exp(b(3)*t))','b','t');换成点除

---

欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)

Powered by Discuz! X2.5