Question

在matlab中用ode45命令求解微分方程组，画随时间变化的曲线。

前两张图是我要求解的微分方程和想得到的类似图形；后两张图是我自己编的，但没运行出来。请帮我看看应该怎么改或怎么写程序。感激不尽！虽然我没有积分了，但是被采纳者我发个qq红包之类的。诚意求解答，急！！！最好在6月20号之前。谢谢！





Answer 1

function test()

p0=0.6;q0=0.7;

a1=0.07; b1=0.04; s1=0.3; k1=0.5; m1=1800; w1=0.5; c1=3;

a2=0.04; b2=0.02; s2=0.4; k2=0.4; m2=1500; w2=0.4; c2=2;

u1 = (a1 - b1) * s2 * k2 * m2 + w1;

u2 = (a2 - b2) * s1 * k1 * m1 + w2;

[t, y] = ode45(@func, [0, 5], [p0; q0], [], u1, -c1, u2, -c2);

plot(t, y);

function f=func(t, y, ua, ub, va, vb)

f = [y(1) * ( 1 - y(1)) * (y(2) * ua + ub); y(2) * (1 - y(2)) * (y(1) * va + vb)];

Answer 2

里面的α等希腊字母是输入法打上的吧，在matlab中貌似不支持