Python绘制相线

写了个函数,用Python绘制自治微分方程frac{dy}{dt}=f(y(t))的相线。

调用时,只要给一个微分方程的右边部分就行了:

如果微分方程右边是周期函数,因为sympy默认只给出很少几个正根,如果需要指定更多的平衡解位置,可以手工限定点的位置。

例如:画出frac{dy}{dt}=y^2cosy的相线,用默认参数,则只会显示0,frac{pi}{2}frac{3pi}{2}

figure_1

如果需要更多位置,可以手工指定所需位置。

Leave a Reply

Your email address will not be published. Required fields are marked *