level 1
古棋香茗
楼主
ti2 = Table[i, {i, 0, 1100, 50}]
ci2 = {0.0000069659, 0.0000097818, 0.0000260849, 0.0000671390, 0.0001474686, 0.0002614418, 0.0003763042, 0.0004837562, 0.0005680875, 0.0006190716, 0.0006496028, 0.0006620524, 0.0006712414, 0.0006709450, 0.0006651648, 0.0006533080, 0.0006407102, 0.0006155146, 0.0006072148, 0.0005931349, 0.0005778693, 0.0005614181, 0.0005465971}
data = Table[{ti2[[i]], ci2[[i]]}, {i, 1, Length[ci2]}]
eq = c'[t] == k1*(0.001575 - 2*c[t])^m*(0.05 - 5*c[t])^n - k2*c[t]^l
pfun = ParametricNDSolveValue[{eq, c[ti2[[1]]] == ci2[[1]]}, c, {t, 0, 1100}, {k1, n, m, k2, l}]
sol = NonlinearModelFit[data,pfun[k1, n, m, k2, l][t], {k1, n, m, k2, l}, t]
请大佬们赐教
2020年08月06日 08点08分
1
ci2 = {0.0000069659, 0.0000097818, 0.0000260849, 0.0000671390, 0.0001474686, 0.0002614418, 0.0003763042, 0.0004837562, 0.0005680875, 0.0006190716, 0.0006496028, 0.0006620524, 0.0006712414, 0.0006709450, 0.0006651648, 0.0006533080, 0.0006407102, 0.0006155146, 0.0006072148, 0.0005931349, 0.0005778693, 0.0005614181, 0.0005465971}
data = Table[{ti2[[i]], ci2[[i]]}, {i, 1, Length[ci2]}]
eq = c'[t] == k1*(0.001575 - 2*c[t])^m*(0.05 - 5*c[t])^n - k2*c[t]^l
pfun = ParametricNDSolveValue[{eq, c[ti2[[1]]] == ci2[[1]]}, c, {t, 0, 1100}, {k1, n, m, k2, l}]
sol = NonlinearModelFit[data,pfun[k1, n, m, k2, l][t], {k1, n, m, k2, l}, t]
请大佬们赐教

