缺少部分参数的微分方程拟合之二

问题来源：非线性微分方程组多参数拟合-测试-CSDN问答

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

分析：S,E,R缺少初值数据，故作为拟合参数与α,p,β,ε,q,θ一起进行拟合。

Lu代码：

!!!using["luopt","math"]; //使用命名空间
f(t,S,E,I,R,dS,dE,dI,dR, params :: α,p,β,ε,q,θ)=
{