数学建模社区-数学中国

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

---

作者: 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];
: ]  l! G0 W8 a! Ax=[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];
$ p0 ?8 e6 l5 I8 a/ u2 ?- zplot(t,x,'*')
# l( ?% u3 b% |: J- [ b0=[0.46,5.3,-0.45];
4 A% E: I: l  F! Q: u4 g, Ofun=inline('b(1)/(1+b(2)*exp(b(3)*t))','b','t');
0 t' \4 Z* e: [7 A0 O1 {2 ` [b,r,j]=nlinfit(t,x,fun,b0);
7 A) a& R. x( a$ y: T) pb
8 i0 e1 A+ u9 a* V6 y/ i) ^$ C。。
1 j  o/ n# X9 p, q* O1 I& z1 y请问大神们，为什么老是出现错误呢，麻烦帮我检查fun=inline('b(1)/(1+b(2)*exp(b(3)*t))','b','t');
" b! ~0 P8 ~" @1 ]这句哪儿出错了，急求。。

---

作者: 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];
   ) \" }% _1 I0 M& w' L7 n
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];
   , {" R7 I$ \1 [& q8 P' R) M1 l7 L
3. plot(t,x,'*');
   
4. b0=[0.46,5.3,-0.45];
   ! m$ l0 T' z) J8 }/ Q) Q) q& W: a
5. fun=inline('b(1)./(1+b(2)*exp(b(3)*t))','b','t');
   
6. [b,r,j]=nlinfit(t,x,fun,b0);
   + r+ p' u3 Z* N
7. b

复制代码

改了一下，你试试看
4 \3 E8 J! Z8 j9 `  H

---

作者: 净心、精心    时间: 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