在Python中通过坡度字段绘制曲线
我正在使用以下Python代码(摘自另一篇我记不起来的SO帖子)来绘制微分方程dy/dx=1/(x+y)的斜率场: 现在我想画一条曲线(比如y=-x-1)穿过这个相同的斜率场在Python中通过坡度字段绘制曲线,python,matplotlib,graph-visualization,Python,Matplotlib,Graph Visualization,我正在使用以下Python代码(摘自另一篇我记不起来的SO帖子)来绘制微分方程dy/dx=1/(x+y)的斜率场: 现在我想画一条曲线(比如y=-x-1)穿过这个相同的斜率场 我该怎么做呢?plt.plot(x,-x-1)在plt.show()之前。@gboffi明白了!非常感谢你的帮助。 import numpy as np from matplotlib import pyplot as plt # Differential equation # diff = y'= y/x (or sa
我该怎么做呢?
plt.plot(x,-x-1)
在plt.show()之前。
@gboffi明白了!非常感谢你的帮助。
import numpy as np
from matplotlib import pyplot as plt
# Differential equation
# diff = y'= y/x (or say x+y)
def diff(x,y):
if (x+y == 0):
return 0
else:
return 1/(x+y) # try also x+y
x = np.linspace(-3,3,25)
y = np.linspace(-3,3,25)
# use x,y
for j in x:
for k in y:
slope = diff(j,k)
domain = np.linspace(j-0.07,j+0.07,2)
def fun(x1,y1):
z = slope*(domain-x1)+y1
return z
plt.plot(domain,fun(j,k),solid_capstyle='projecting',solid_joinstyle='bevel')
plt.title("Slope field y'")
plt.grid(True)
plt.show()
print("End of the program")