Fitting experimental data to solve for rate constants
Show older comments
I have y (concentration) and x (time) data for a reaction that can be modeled by the following function:
F(x) = k*L0*((m/(g*h))+((h-m)/(g*(g-h))*exp(-h*x))+((m-g)/(g*(g-h)))*exp(-g*x))+2*k*l*L0*(1/(g*h)+(exp(-g*x)/(g*(g-h)))-(exp(-h*x)/(h*(g-h))))
where g*h=k*(l+m) and g+h=k+l+m.
My only variables are x and y and all other parameters are constants, all I want to do is fit my data to this function to solve for k, l, and m.
Could anybody please help me with writing this code or suggest a resource that might help me?
Answers (1)
Roshni Garnayak
on 7 Nov 2019
0 votes
You can use Symbolic Toolbox to solve your problem. Refer to the following link for the User’s Guide:
For a list of functions available in the toolbox, refer to the following link:https://www.mathworks.com/help/symbolic/referencelist.html?type=function&s_cid=doc_ftr
Categories
Find more on Chemistry in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!