Python 怎么能ı；在M轴上找到cp和cp_p_g的交点？ 将numpy导入为np 将matplotlib.pyplot作为plt导入 Cp_0=-0.464 M=np.arange（0.3,1,0.1）#范围内的马赫数 Cp=（2/（1.4*M**2））*（（1+0.2*M**2）/（1.2））**3.5-1）#等熵流的气动关系 Cp_p_g=Cp_0/（（1-M**2）**0.5）#普兰特-格劳特规则 plt.绘图（M，Cp） plt.绘图（M，Cp\u p\g） 公共电网（“主要”） plt.title（“NACA**-***”）#我们工作的NACA配置文件中的标题各不相同 plt.xlabel（“M”） plt.ylabel（“Cp”）

---

我已经添加了两条曲线的绘图，cp和，cp_p_g，为了更好地理解问题

您可以使用

linspace

优化x网格，然后查找两条曲线的差异

import numpy as np
import matplotlib.pyplot as plt 

Cp_0=-0.464   
M=np.arange(0.3,1,0.1) #Mach number for range<0.3-1>
Cp=(2/(1.4*M**2))*(((1+0.2*M**2)/(1.2))**3.5-1) #Aerodynamic relation for isentropic flow
Cp_p_g=Cp_0/((1-M**2)**0.5) #prandtl-Glauert rule

plt.plot(M,Cp)
plt.plot(M,Cp_p_g)
plt.grid("major")
plt.title('NACA**-***') #Title varies on the NACA profiles that we work
plt.xlabel("M")
plt.ylabel("Cp")

M = np.linspace(0.3, 1, 200)

# Remaining code

idx = np.argmin(np.abs(Cp_p_g-Cp))

plt.axvline(M[idx], c='k', ls='--')
plt.axhline(Cp_p_g[idx], c='k', ls='--')