Question

MATLAB 差分方程求解

想求解下面这个差分方程...其中a1 a2 a3 a4 a5是参量 可以这么设
a1=20,a4=5,a2=a3=a5=0.1
s(t+1)=s(t)-a1*s(t)*i(t)-a2*a3*s(t)*i(t)/(s(t)+e(t)+i(t));
e(t+1)=e(t)+a1*s(t)*i(t)-e(t-a4)-a2*a3*(e(t)-i(t))/(s(t)+e(t)+i(t));
i(t+1)=i(t)+e(t-a4)-a2*i(t)^2/(s(t)+e(t)+i(t));
j(t+1)=j(t)+a2*i(t);
r(t+1)=r(t)+a5j(t);
关于初值，可以这么设
s=ones(1,50);
e=zeros(1,50);
i=zeros(1,50);
j=zeros(1,50);
r=zeros(1,50);

主要想问问 这种多元的差分方程想写编程怎么个思路...
解决了问题追加100...谢谢

Answer 1

这个问题要过细想，a1=20,a4=5,a2=a3=a5=0.1
s(t+1)=s(t)-a1*s(t)*i(t)-a2*a3*s(t)*i(t)/(s(t)+e(t)+i(t));
e(t+1)=e(t)+a1*s(t)*i(t)-e(t-a4)-a2*a3*(e(t)-i(t))/(s(t)+e(t)+i(t));
i(t+1)=i(t)+e(t-a4)-a2*i(t)^2/(s(t)+e(t)+i(t));
j(t+1)=j(t)+a2*i(t);
r(t+1)=r(t)+a5j(t);
关于初值，可以这么设
s=ones(1,50);
e=zeros(1,50);
i=zeros(1,50);
j=zeros(1,50);
r=zeros(1,50);