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

2024-04-27 05:10:51 发布

您现在位置:Python中文网/ 问答频道 /正文

在对以前工作的一些回答中,我无法访问有关样条曲线和从def创建的多条直线之间的交点的文档的某些链接。我想找出系统需求曲线和泵曲线相交的点。我尝试过几种方法;但是,numpy数组是空列表。我有8个正确的列表,只是没有在他们

有关于样条曲线和曲线相交的文档,其中x值必须增加(在我的例子中是这样),但我无法阅读文档。你知道吗

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

我创建了一条样条曲线,并希望找到我创建的曲线的交点。我试着用np.1d段在我的for循环中,创建了我希望是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()

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


Tags: 文档列表ratedefnpplt数组flow