最近在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 f=@(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])