How to solve non-linear langmuir type equations which have 4 variables?
2 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Hi
I want to solve the langmuir type equations which have 4 unknown constants.
The equation is
k*a*p*c*(aer)^2*(x-xeq)/(1+(a*p)^0.5+b*0.03)^3-y=0
The variables are k a b c p, aer, xeq are the parameters for the conditions.
p = [0.97 0.8 0.5 0.3]
aer = [0.6 0.5 0.4 0.36]
xeq = [2.2e-10 2.6e-10 4.2e-10 7.1e-10]
x and y are the experimental data at each conditions.
For example,
if [p, aer, xeq] = [0.97, 0.6, 2.2e-10], we have 2~3 pairs of [x, y].
So with this condition, I want to estimate the 4 unknown values, a, b, c, k.
I have start the study of MATLAB just yesterday, so the answer with detail like whole function definition and codes will be very helpful for me.
Thanks.
2 commentaires
Torsten
le 25 Août 2017
Modifié(e) : Torsten
le 25 Août 2017
Use MATLAB's "lsqcurvefit".
Note that k and c cannot be determined independently from each other.
Reformulate your equation as
kc*a*p*(aer)^2*(x-xeq)/(1+(a*p)^0.5+b*0.03)^3-y=0
with the three variables kc, a and b to be determined.
Best wishes
Torsten.
Réponses (1)
John D'Errico
le 25 Août 2017
Modifié(e) : John D'Errico
le 25 Août 2017
So, for each combination of p,aer,xeq, you wish to estimate a,b,c,k variables. To do this at each such combination you say that you have 2-3 pairs of data points, thus (x,y) pairs?
You cannot estimate 4 variables using 2 or 3 data points! That is impossible. Any solution that is found (assuming one is possible) will be only one of infinitely many solutions. So the solution will be not a useful solution.
While Torsten told you to use lsqcurvefit, that would be appropriate only if you have more data points that you have parameters to estimate.
The answer is you need to get more data. To have any realistic chance of being able to get good estimates for 4 variables, even 10 data pairs would be the minimum I would normally suggest, and more is always better. How many you really need depends on how much noise there is in your data. But at a bare minimum, you absolutely need at least 4 pairs of (x,y).
1 commentaire
Torsten
le 25 Août 2017
As I understood the OP, he/she has 2-3 pairs of [x,y] for the four combinations of [p aer xeq], thus 8-12 pairs of [x y] in total.
Best wishes
Torsten.
Voir également
Catégories
En savoir plus sur 非線形方程式系 dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!