Matlab作常微分方程组相图
用streamslice函数
eg.作出下列系统的相图 x′=x2−2x−xy y′=y2−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.的意义是点乘
效果:
从图中可以很明显的看到4个奇点
第一行改变图像区间[x,y]=meshgrid(linspace(0,5));
可以得到(3,1)处焦点的大图:
探索过程中一些有用的操作:
- clear:清除workspace
- clc:清除command window
- clf:清空figure(就是图像显示框)
|
请发表评论