這個方程要解出來還要初始條件吧
不妨設(shè)(x(0),y(0))=(A,B);
還有用ode45只能得到數(shù)值解,不能給出解析式表達
那么先建立一個函數(shù)文件
test_fun.m
function dx=tsst_fun(x,t)
dx(1)= c-b*x(1)-w*x(1)*x(2);
dx(2)=w*x(1)*x(2)-(b+m)*x(2);
那么調(diào)用ode45
x0=(0,0);
t0=0:0.01：10；
；%根據(jù)你所想要的求得值設(shè)定t0,間隔是任意的,與求解所用的步長無關(guān),
[x,t]=ode45(@test_fun,t0,x0);
得到了