2020-10-04 Matlab中相图的画法

Matlab中相图的画法

最近在Youtube中看到一段很好的相图代码，把它拿出来分享给大家。
这是对如下微分方程所代表的动态系统绘制相图：
d y d x = t + sin ⁡ ( y ) \frac{dy}{dx}=t+\sin (y) dxdy​=t+sin(y)
Matlab程序如下：

[t,y]=meshgrid(0:.2:4,0:.2:4);
m=t+sin(y);
L=sqrt(1+m.^2);
quiver(t,y,1./L,m./L)
axis tight
hold on
[e-mail](u,v) u+sin(v)
[u,v]=ode45(f,[0 5],1:3);
plot(u,v)
hold off
axis([0 4 0 4])
axis tight equal
axis([0 4 0 4])

演示结果