Python绘制斜率场

网上找了一圈,没找到Python里现成的绘制斜率场(方向场)的函数,利用Matplotlib.pyplot里的quiver() DIY了一个斜率场的绘图函数,通过sympy加了一个寻找经过指定点的解并绘图的功能。

博主用Python学ODE的笔记见:http://ryancheunggit.gitbooks.io/study-notes-for-differential-equations/

例如: 给定微分方程frac{dy}{dt}=y-t要绘制-3leq t leq 3,-3leq y leq 3范围内的斜率场,并分别找出经过如下几点:(2,4),(1,3),(0,2),(2,0),(2,1),(1,1),(0,0)的解的函数图。可以如此:

效果图:

04-03LinesInField

函数代码如下:

 

Leave a Reply

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