Python 样条曲线与函数中多条直线的交点

Python 样条曲线与函数中多条直线的交点,python,numpy,matplotlib,spline,Python,Numpy,Matplotlib,Spline,在对以前工作的一些回答中,我无法理解有关样条曲线和从def创建的多条直线之间的交点的文档的一些链接。我想抓住系统需求曲线和泵曲线相交的点。我尝试了一些方法;但是,numpy数组是空列表。我正确地列出了8个列表,但没有任何内容 有关于样条曲线和曲线之间的交点的文档,其中x值必须增加(在我的例子中是这样),但我无法阅读文档 对于我的代码…这是我尝试过的 我创建了一条样条曲线,希望找到我创建的曲线的交点。我尝试在for循环中使用np.intersect1d,并创建了我希望长度为1的8数组;但是,它们是

在对以前工作的一些回答中,我无法理解有关样条曲线和从def创建的多条直线之间的交点的文档的一些链接。我想抓住系统需求曲线和泵曲线相交的点。我尝试了一些方法;但是,numpy数组是空列表。我正确地列出了8个列表,但没有任何内容

有关于样条曲线和曲线之间的交点的文档,其中x值必须增加(在我的例子中是这样),但我无法阅读文档

对于我的代码…这是我尝试过的

我创建了一条样条曲线,希望找到我创建的曲线的交点。我尝试在for循环中使用np.intersect1d,并创建了我希望长度为1的8数组;但是,它们是空数组…请提供帮助

plt.figure()
xnew = np.linspace(Flow_Rate.min(),Flow_Rate.max(),300)
power_smooth = spline(Flow_Rate,Pump_Head,xnew)
plt.plot(xnew,power_smooth)
plt.show()
Changing_diameter = [0.06985, .0762, .08255, .0889, .09525, .1016, .10795, .1143]

def H_l(x):    
    reyn = ((density*(x*(1.6667*10**-5))/Area)*d_diameter)/(v_vis)
    return Length/d_diameter*(x/Area)**2/(2*g)*64./reyn

Head = H_l(Q_si)
Changing_head = []


for i in range(len(Changing_diameter)):
    d_diameter = Changing_diameter[i]
    Head = H_l(Q_si)
    Changing_head = []
    plt.xlim(-5,100)
    plt.ylim(-5,50)
    plt.plot(Q_imp,Head,label=d_diameter)
    np.array(Q_imp)
    np.array(Head)
    print(np.intersect1d(Q_imp, Head, assume_unique = False))

plt.legend()
plt.show()

我的预期结果将是一个列表,该列表的长度为一个项目,只有样条线和创建的曲线之间的交点

欢迎使用StackOverflow。请尽量提供一个清晰的答案。没有mcve是很难帮助的。但是,您的代码的一个问题可能是
np.intersect1d
返回两个数组之间的交集,您需要具有相同的元素,并且您的问题可能不是这样(但我无法验证此假设…)