Re: [討論] 給你10筆資料(t,y) 求y=1-a*e^(b*t) 的a b已回收
※ 引述《poo123456 (poo)》之銘言：
: t = [0:1:9];
: y = [0 0.393 0.632 0.777 0.865 0.918 0.950 0.970 0.982 0.989 ];
: 求解 y=1-a*e^(b*t) 的 a b
: 怎麼用MATLAB打 用手算一看就知道a=1,b=-0.499(大概)
: 但是怎麼用MATLAB打呢 感謝@@
兩種做法，一種是先線性化
ln(1-y) = ln(a) + b*t
然後做線性迴歸
====================Matlab code==================
A = [ones(size(transpose(t))), transpose(t)];
ybar = log(1-transpose(y));
b = A\ybar;
=================================================
上面做出來的b(1)就是ln(a)，b(2)就是b。
另外一種法是用非線性迴歸，這個比較麻煩一點，可以參
考 #1BwOrYjT (推銷一下自己的文章)
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 220.136.52.186
推
01/13 16:57, , 1F
01/13 16:57, 1F
→
01/13 16:59, , 2F
01/13 16:59, 2F
→
01/13 17:00, , 3F
01/13 17:00, 3F
→
01/13 17:01, , 4F
01/13 17:01, 4F
推
01/13 17:01, , 5F
01/13 17:01, 5F
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
14
35