SEIRS模型参数估计与趋势分析.pdfVIP

附录

(1)估计参数

函数

SEIRS_估计

全局tdataxdatax0

tdata=1:15;

xdata=[20,191,220,215,206,204,203,202,204,201,203,201,202,201,201];

x0=[10,20,3];//给出初值

b=0.01;

afminsearch(@fun_SEIRS,b)

end

函数S=fun_SEIRS(b)全局变量tdataxdatax0函数

dy=f_SEIRS(t,y)K=1000;d=0.03;e=0.4;delta=0.

7;g=1;dy=zeros(3,1);dy(1)

b*(K-y(1)-y(2)-y(3))*y(2)-(d+e)*y(1);dy(2)

e*y(1)-(d+g)*y(2);dy(3)=g*y(2)-(d+

delta)*y(3);

end

sol=ode23(@f_SEIRS,[1,15],x0);

x_est=deval(sol,tdata,2);

S=norm(x_est-xdata)//计算误差

end

_

SEIRS估计

S=