数学建模社区-数学中国

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

---

作者: 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];
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];
* }5 _4 f; n4 p3 K# i: kfun=inline('b(1)/(1+b(2)*exp(b(3)*t))','b','t');
: a4 I" x5 W. b5 |5 `0 f$ I [b,r,j]=nlinfit(t,x,fun,b0);
b
4 h/ t* u- R' h6 i$ p8 C。。
% ?2 b# d' d- s% s* y% d请问大神们，为什么老是出现错误呢，麻烦帮我检查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];
   ( C" z1 J" c4 R2 H0 j: P! \
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,'*');
   3 u+ h& m( |! k( I
4. b0=[0.46,5.3,-0.45];
   
5. fun=inline('b(1)./(1+b(2)*exp(b(3)*t))','b','t');
   + j/ _/ X: i& k# ^1 |3 ?. G" V
6. [b,r,j]=nlinfit(t,x,fun,b0);
   
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