Hello,
I have the following system of odes:
ds/dt = k1*C-(k2*S*(0.1-C)) ;
dc/dt = (k1*S*(0.1-C))-((k2+k3)*C);
with initial conditions s0 = 5, c0 = 0;
I have two data sets ( re run) for the same components (s and c) at slightly different times in the data shown below
first set second set
time(s1) s1 time(c1) c1 time(s2) s2 time(c2) c2
0.0000 5.0000 0.0000 0.0000 0.0000 4.9662 0.0000 0.0000
1.0000 5.4574 0.9586 0.0961 1.0195 5.5637 0.9586 0.0921
2.0000 5.9366 1.9733 0.0964 2.0516 6.0899 1.9733 0.0951
3.0000 6.4172 2.9608 0.0967 2.9798 6.3740 2.9608 0.0954
4.0000 6.8990 4.0740 0.0969 4.0622 7.0063 4.0740 0.0987
5.0000 7.3818 5.0490 0.0971 4.8049 7.0937 5.0490 0.0981
6.0000 7.8653 6.1736 0.0973 6.2603 8.2064 6.1736 0.1001
7.0000 8.3497 6.9074 0.0974 6.7812 8.0887 6.9074 0.0961
8.0000 8.8348 7.7648 0.0976 7.8129 8.6282 7.7648 0.0947
9.0000 9.3205 8.6280 0.0977 9.2680 9.5981 8.6280 0.0936
10.0000 9.8069 10.2719 0.0978 9.9876 9.7947 10.2719 0.1005
11.0000 10.2939 10.6762 0.0979 11.2959 10.5707 10.6762 0.0950
12.0000 10.7812 11.8659 0.0980 11.8752 10.6691 11.8659 0.0969
13.0000 11.2692 13.0673 0.0981 12.7048 11.0133 13.0673 0.0986
14.0000 11.7575 13.6205 0.0982 13.3521 11.2134 13.6205 0.0955
15.0000 12.2461 15.2129 0.0982 15.2599 12.4583 15.2129 0.0996
16.0000 12.7352 15.9752 0.0983 15.8873 12.6455 15.9752 0.0981
17.0000 13.2246 16.4081 0.0984 16.9180 13.1608 16.4081 0.0949
18.0000 13.7143 18.5075 0.0984 18.1977 13.8650 18.5075 0.1012
19.0000 14.2043 18.2412 0.0985 18.1629 13.5785 18.2412 0.0945
20.0000 14.6946 19.5881 0.0985 19.6316 14.4240 19.5881 0.0965
21.0000 15.1852 20.4485 0.0986 21.5727 15.5993 20.4485 0.0960
22.0000 15.6760 22.0679 0.0986 22.4322 15.9839 22.0679 0.0989
23.0000 16.1670 22.0604 0.0987 22.1383 15.5613 22.0604 0.0946
24.0000 16.6582 23.7728 0.0987 23.1124 16.0421 23.7728 0.0978
25.0000 17.1497 24.0121 0.0987 23.9809 16.4505 24.0121 0.0948
26.0000 17.6412 24.9919 0.0988 24.7203 16.7730 24.9919 0.0949
27.0000 18.1330 27.7680 0.0988 26.7924 17.9936 27.7680 0.1016
28.0000 18.6249 27.4164 0.0988 28.4356 18.9147 27.4164 0.0968
29.0000 19.1170 29.3002 0.0989 29.6465 19.5432 29.3002 0.0999
30.0000 19.6092 31.3933 0.0989 30.6600 20.0406 31.3933 0.1035
I am trying to use lsqcurvefit to estimate the values of the unknown parameters (k1, k2 and k3) using the two sets of the experimental data. I am somewhat new to MATLAB coding as I only use it once in a while.
Could anyone help me set my data up and set up the code. If someone will write the code too that will be super helpful. Thanks.
I have the following system of odes:
ds/dt = k1*C-(k2*S*(0.1-C)) ;
dc/dt = (k1*S*(0.1-C))-((k2+k3)*C);
with initial conditions s0 = 5, c0 = 0;
I have two data sets ( re run) for the same components (s and c) at slightly different times in the data shown below
first set second set
time(s1) s1 time(c1) c1 time(s2) s2 time(c2) c2
0.0000 5.0000 0.0000 0.0000 0.0000 4.9662 0.0000 0.0000
1.0000 5.4574 0.9586 0.0961 1.0195 5.5637 0.9586 0.0921
2.0000 5.9366 1.9733 0.0964 2.0516 6.0899 1.9733 0.0951
3.0000 6.4172 2.9608 0.0967 2.9798 6.3740 2.9608 0.0954
4.0000 6.8990 4.0740 0.0969 4.0622 7.0063 4.0740 0.0987
5.0000 7.3818 5.0490 0.0971 4.8049 7.0937 5.0490 0.0981
6.0000 7.8653 6.1736 0.0973 6.2603 8.2064 6.1736 0.1001
7.0000 8.3497 6.9074 0.0974 6.7812 8.0887 6.9074 0.0961
8.0000 8.8348 7.7648 0.0976 7.8129 8.6282 7.7648 0.0947
9.0000 9.3205 8.6280 0.0977 9.2680 9.5981 8.6280 0.0936
10.0000 9.8069 10.2719 0.0978 9.9876 9.7947 10.2719 0.1005
11.0000 10.2939 10.6762 0.0979 11.2959 10.5707 10.6762 0.0950
12.0000 10.7812 11.8659 0.0980 11.8752 10.6691 11.8659 0.0969
13.0000 11.2692 13.0673 0.0981 12.7048 11.0133 13.0673 0.0986
14.0000 11.7575 13.6205 0.0982 13.3521 11.2134 13.6205 0.0955
15.0000 12.2461 15.2129 0.0982 15.2599 12.4583 15.2129 0.0996
16.0000 12.7352 15.9752 0.0983 15.8873 12.6455 15.9752 0.0981
17.0000 13.2246 16.4081 0.0984 16.9180 13.1608 16.4081 0.0949
18.0000 13.7143 18.5075 0.0984 18.1977 13.8650 18.5075 0.1012
19.0000 14.2043 18.2412 0.0985 18.1629 13.5785 18.2412 0.0945
20.0000 14.6946 19.5881 0.0985 19.6316 14.4240 19.5881 0.0965
21.0000 15.1852 20.4485 0.0986 21.5727 15.5993 20.4485 0.0960
22.0000 15.6760 22.0679 0.0986 22.4322 15.9839 22.0679 0.0989
23.0000 16.1670 22.0604 0.0987 22.1383 15.5613 22.0604 0.0946
24.0000 16.6582 23.7728 0.0987 23.1124 16.0421 23.7728 0.0978
25.0000 17.1497 24.0121 0.0987 23.9809 16.4505 24.0121 0.0948
26.0000 17.6412 24.9919 0.0988 24.7203 16.7730 24.9919 0.0949
27.0000 18.1330 27.7680 0.0988 26.7924 17.9936 27.7680 0.1016
28.0000 18.6249 27.4164 0.0988 28.4356 18.9147 27.4164 0.0968
29.0000 19.1170 29.3002 0.0989 29.6465 19.5432 29.3002 0.0999
30.0000 19.6092 31.3933 0.0989 30.6600 20.0406 31.3933 0.1035
I am trying to use lsqcurvefit to estimate the values of the unknown parameters (k1, k2 and k3) using the two sets of the experimental data. I am somewhat new to MATLAB coding as I only use it once in a while.
Could anyone help me set my data up and set up the code. If someone will write the code too that will be super helpful. Thanks.