非线性微分方程组多参数拟合

Question

已知4个微分方程如下：

dS/dt= − αSI − pS
dE/dt= αSI − βE − εE + qR
dI/dt=βE − θI
dR/dt= pS +εE+ θI − qR

已知的参数是：

t=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15]
I=[11714 21352 21789 22403 16308 17017 17146 33229 60907 38875 33437 47282 51749 40530 31617]'

怎么才可以求出方程中的其他数据：S,E,R ,α,p,β,ε,q,θ这些参数的值

新手小白，请各位大佬赐教！！！

Answer 1

参考这里，已解决。
https://blog.csdn.net/forcal/article/details/119281810?