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

已知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,θ这些参数的值

新手小白,请各位大佬赐教!!!

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

https://stackoverrun.com/cn/q/3378893