流行病模型的参数估计

Question

我有一个流行病模型：

function ypsirid = ypsirid(t,y)
global q m b r
ypsirid(1) = q - m*y(1) - b*y(1)*y(2);
ypsirid(2) = (b*y(1)*y(2)) - (m+r)*y(2);
ypsirid(3) = (r*y(2)) - m*y(3);
ypsirid = ypsirid(:)

任何人都可以建议使用matlab进行参数估计的最简单方法。谢谢！