Matlab作常微分方程组相图

Matlab作常微分方程组相图

用streamslice函数

eg.作出下列系统的相图
 x=x22xxy y=y24y+xy \ x'=x^2-2x-xy\\ \ y'=y^2-4y+xy
Matlab代码:
[x,y]=meshgrid(linspace(-10,10));
streamslice(x,y,x.(-2+x-y),y.(-4+y+x));
xlabel(‘x’);ylabel(‘y’);

*这里x.和y.不可改成x,y,x.的意义是点乘

效果:
Matlab作常微分方程组相图
从图中可以很明显的看到4个奇点
第一行改变图像区间[x,y]=meshgrid(linspace(0,5));
可以得到(3,1)处焦点的大图:
Matlab作常微分方程组相图

探索过程中一些有用的操作:
  • clear:清除workspace
  • clc:清除command window
  • clf:清空figure(就是图像显示框)