悬赏 30 个论坛币 未解决
我现在有一个微分方程组
Fx =-x*(x - 1)*(cs - k*x - 3*k*y + 2*k*y^2 + 2*k*x*y)
Fy =-y*(y - 1)*(cs - 3*k*x - k*y + 2*k*x^2 + 2*k*x*y)
其中x,y是自变量,cs,k是常数(下面的例子中,cs=0.5, k=5)。
现在想绘制x,y在[0,1]上的相图,我自己写的matlab代码如下,不知道对不对,希望有经验的朋友指点一下,有积分感谢!
syms x y cb cs k;
[x,y] = meshgrid(0:0.05:1, 0:.05:1);
Fx=-x*(x - 1)*(cs - k*x - 3*k*y + 2*k*y^2 + 2*k*x*y);
Fx=subs(Fx,{cs,k},{0.5,5});
Fy=-y*(y - 1)*(cs - 3*k*x - k*y + 2*k*x^2 + 2*k*x*y);
Fy=subs(Fy,{cs,k},{0.5,5});
streamslice(x,y,Fx,Fy)