[討論] 請問關於fitting

看板MATLAB作者 (就是要ROCK)時間12年前 (2013/05/23 17:47), 編輯推噓0(004)
留言4則, 1人參與, 最新討論串1/1
請問有辦法 一次直接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)

05/25 00:41, , 1F
會用一條方程式nlinfit出a,b,Ef的話,令一個誤差平方和函數為
05/25 00:41, 1F

05/25 00:43, , 2F
(d(Y1)/d(x)-Y1).^2+(d(Y2)/d(x)-Y2).^2+(d(Y3)/d(x)-Y3).^2
05/25 00:43, 2F

05/25 00:44, , 3F
將此目標函數最小化即可
05/25 00:44, 3F

05/25 00:44, , 4F
上面d(Y1)/d(x)等就是指你寫的那一串方程式
05/25 00:44, 4F
文章代碼(AID): #1HdUNFW0 (MATLAB)