作者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