使用ode系列的函数求解微分方程组.
d = @(t,x)[x(2); -x(1).^2-2*x(1)-x(2)/2];
[t x] = ode45(d,[0 5.1],[0; 4]); % 这个tfinal自己调出来的一个近似值
plot(x(:,1),x(:,2))
用数值法ode45的函数计算。
fun = @(t,x)[x(1).*(1-x(1)).*(11-5*x(2)) ; x(2).*(1-x(2)).*(0.5+9.5*x(2))];
[t x] = ode45(fun, [0,5], [0; 0]); %初值和区间可修改
plot(x(:,1), x(:,2))