For loop in a function (not getting right results)
4 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
I created the function Antoine to not bother solving the equation over and over. However, when I call the function in a for loop I dont get the right answers compared to the way i would do it in excel. I would appreciate any help.
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/932219/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/932209/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/932214/image.png)
p1_pure (using matlab) =
0 0 0 0 0.9800
using excel:
1.903628223
2.09045263
2.283986537
2.483798717
2.689463014
0 commentaires
Réponse acceptée
Voss
le 18 Mar 2022
Several things. First:
for i = length(T_K)
should be
for i = 1:length(T_K)
The first way just iterates once, using the last element of T_K, which is why the first four elements of p1_pure are zero.
Second: T_K in the code is in Kelvin, but T in the formula is in Celsius, so you would have to use T_K-273.15 in your MATLAB calculation to convert it.
Third: exp() is natural (base e) exponentiation, but the formula says log10 (not natural log), so your MATLAB code would have 10^() rather than exp().
Making those changes still doesn't give a result that matches what Excel gives you, so I would check the units of a, b and c, and check your Excel formula.
a = 4.00266;
b = 1171.53;
c = -48.784;
T_K = 300:10:340;
for i = 1:length(T_K)
p1_pure(i) = 10^(a-b/(T_K(i)-273.15+c)); % make this a function if you want
end
p1_pure
2 commentaires
Voss
le 18 Mar 2022
Modifié(e) : Voss
le 18 Mar 2022
You're welcome! And nope, log10() and exp() (or log()) are not interchangeable:
- exp() is the inverse function of log(), both of which have to do with the base e (i.e., natural) logarithm, also known as ln() in math (but not in MATLAB).
- 10^() is the inverse function of log10(), which have to do with the base-10 logarithm.
- There is also log2(), whose inverse function would be 2^() (or pow2()).
Plus de réponses (0)
Voir également
Catégories
En savoir plus sur Time Series Collections 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!