python 直线求交_python绘制直线与d的交点

可以使用scipy中的interp1d函数在每组点之间创建一个线性插值。interp1d接受x和y值，并返回一个接受x值并返回相应y值的函数。一旦你有了应力和y的函数，你需要找出它们在哪里相交。一种有效的方法是使用像二分法这样的方法，它可以找到零。在

代码如下：将numpy作为np导入import matplotlib.pyplot as plt

from scipy.interpolate import interp1d

from scipy.optimize import bisect

data = '''0 0

0.1 0.1

0.2 0.2

0.3 0.3

0.4 0.31

0.5 0.32

0.6 0.35'''

data = [line.split() for line in data.split('\n')]

data = np.array(data, dtype = 'float')

strain = data[:,0]

stress = data[:,1]

E = 1

y = [0, 0.5]

x = 0.2, ((y[1] - y[0])/E+0.2)

#use interp1d to get interpolated points

y = interp1d(x, y)

stress = interp1d(strain, stress)

#set starting points

x1 = max(x[0], strain[0])

x2 = min(x[-1], strain[-1])

max_err = .01

#create function

f = lambda x : stress(x) - y(x)

#find x1 where f(x1) = 0

x1 = bisect(f, x1, x2, xtol = .001)

y1 = stress(x1)

plt.figure(num=None, figsize=(10, 6), dpi=100, facecolor='w', edgecolor='k')

plt.plot(strain, stress(strain))

plt.plot(x, y(x))

plt.scatter(x1, y1)

plt.show()

以及生成的图像：