作者ghost17612 (就是要ROCK)
看板MATLAB
標題[討論] 請問關於fitting
時間Thu May 23 17:47:56 2013
請問有辦法 一次直接fitting 三組數據嗎?
例如說
Y1 = [..... ] (對應到的分率)
Y2 = [..... ]
Y3 = [..... ]
X = [..... ] (時間)
T = [T1 T2 T3] (三個溫度下)
然後Eqution 像是
d(Y1)/d(x) = a*exp(Ef/8.3145/T1)*f(Y)-b*exp((Ef+5330)/8.3145/T1)*g(Y)
d(Y2)/d(x) = a*exp(Ef/8.3145/T2)*f(Y)-b*exp((Ef+5330)/8.3145/T2)*g(Y)
d(Y3)/d(x) = a*exp(Ef/8.3145/T3)*f(Y)-b*exp((Ef+5330)/8.3145/T3)*g(Y)
這樣一次Fitting出 a,b和Ef
目前我只會用"nlinfit"分別fitting三個equation的數據在三個溫度下
詳細的code怎麼寫如果需要 我可以再附上
現在的想法是
想寫出一個code可以一次fitting三組data (Y1,X) (Y2,X) (Y3,X)
然後是可以得到Ef a b 在最小SSE的情況下
非常感謝!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.22.253
※ 編輯: ghost17612 來自: 140.112.22.253 (05/23 17:48)
1F：→ aia:會用一條方程式nlinfit出a,b,Ef的話，令一個誤差平方和函數為 05/25 00:41
2F：→ aia:(d(Y1)/d(x)-Y1).^2+(d(Y2)/d(x)-Y2).^2+(d(Y3)/d(x)-Y3).^2 05/25 00:43
3F：→ aia:將此目標函數最小化即可 05/25 00:44
4F：→ aia:上面d(Y1)/d(x)等就是指你寫的那一串方程式 05/25 00:44