Using splines to find the Area Under the Curve

Using splines to find the Area Under the Curve - In these two worksheet the Area Under the Curve is found for a biomass batch reactor and a pulse - Messaggi

#1 Pubblicato: 15/05/2023 20:34:26
Javier Ochoa

Javier Ochoa

3 apprezzamenti in 26 post.

Gruppo: User

Hi,

These two examples solve the problem calculations for the total biomass formed in a batch bioreactor and the pattern flow of a reactor with a pulse of some tracer. The splines method and the Simpson integration formula found in the forum let me solve both problems.


Kind regards,


splines_integral_data.sm (30,28 KiB) scaricato 1673 volta(e).

splines_integral_data1.sm (31,52 KiB) scaricato 1649 volta(e).
1 utenti hanno apprezzato questo post
Oscar Campo 26/05/2023 07:10:00
#2 Pubblicato: 22/05/2023 15:02:07
CBG

CBG

78 apprezzamenti in 312 post.

Gruppo: User

Hello

Fitting Javier's data to the Gauss function by the Levenberg-Marquardt Algorithm

Fit_Javier_LMA.png

Javier_fit_CBG.sm (37,14 KiB) scaricato 1639 volta(e).


Best Regards

Carlos
1 utenti hanno apprezzato questo post
Javier Ochoa 22/05/2023 15:20:00
#3 Pubblicato: 22/05/2023 15:56:44
Javier Ochoa

Javier Ochoa

3 apprezzamenti in 26 post.

Gruppo: User

Hi Jean and Carlos,
These worksheets help me to better understand the problem which I tried to solve.

Thanks in advance,

#4 Pubblicato: 23/05/2023 16:29:20
Javier Ochoa

Javier Ochoa

3 apprezzamenti in 26 post.

Gruppo: User

Hi Jean and Carlos,

I am curious about your thoughts for the result I found using the al_nlqsolve routine from ALGLIB.

May be I was lucky to find a similar solution.

Kind regards,


LevenbergMarquardt_solver.sm (27,68 KiB) scaricato 1687 volta(e).
#5 Pubblicato: 24/05/2023 00:30:27
CBG

CBG

78 apprezzamenti in 312 post.

Gruppo: User

Wrote

Hi Jean and Carlos,

I am curious about your thoughts for the result I found using the al_nlqsolve routine from ALGLIB.

May be I was lucky to find a similar solution.

Kind regards,


LevenbergMarquardt_solver.sm (27,68 KiB) scaricato 1687 volta(e).



Hi Javier
As you can see in the graph, the green line is that of your solution,
and it is not correct, since it does not fit the data.

javier_fit2.png

It is always important to draw the curve

Javier fit.sm (24,31 KiB) scaricato 1625 volta(e).

Best Regards

Carlos

1 utenti hanno apprezzato questo post
Javier Ochoa 24/05/2023 09:31:00
#6 Pubblicato: 24/05/2023 03:41:54
Alvaro

Alvaro

1.139 apprezzamenti in 1.805 post.

Gruppo: User

Hola Javier. Your Levenberg Marquardt algorithm works pretty well, but you choose a bad guess starting point, with a very high tolerance. BTW, I made some few modifications and use also other target function.

LevenbergMarquardt_solver.sm (32,85 KiB) scaricato 1667 volta(e).
LevenbergMarquardt_solver.pdf (233,35 KiB) scaricato 1004 volta(e).

Best regards.
Alvaro.
2 utenti hanno apprezzato questo post
Javier Ochoa 24/05/2023 09:31:00, Вячеслав Мезенцев 24/05/2023 04:21:00
#7 Pubblicato: 24/05/2023 15:52:25
Alvaro

Alvaro

1.139 apprezzamenti in 1.805 post.

Gruppo: User

Wrote

Hi,

These two examples solve the problem calculations for the total biomass formed in a batch bioreactor and the pattern flow of a reactor with a pulse of some tracer. The splines method and the Simpson integration formula found in the forum let me solve both problems.


Kind regards,


splines_integral_data.sm (30,28 KiB) scaricato 1673 volta(e).

splines_integral_data1.sm (31,52 KiB) scaricato 1649 volta(e).



Hola Javier. You can also use this code for integrate and differentiate tabulated data.

Click to enlarge

splines_integral_data.sm (85,66 KiB) scaricato 1666 volta(e).

splines_integral_data.pdf (180,67 KiB) scaricato 938 volta(e).

Best regards.
Alvaro.
1 utenti hanno apprezzato questo post
Javier Ochoa 29/05/2023 09:25:00
#8 Pubblicato: 25/05/2023 13:16:26
Jean Giraud

Jean Giraud

983 apprezzamenti in 6.866 post.

Gruppo: User

.. refreshed/augmented.
Cheers ... Jean.

Inst_Frechet GAUSSIAN.sm (92,39 KiB) scaricato 1724 volta(e).
1 utenti hanno apprezzato questo post
Javier Ochoa 29/05/2023 09:25:00
  • Nuovi post
  • Nessun nuovo post