ValueError：操作数无法与形状一起广播（51,51）（100，）

Question

我试图分配，但是我得到一个错误：ValueError：操作数无法与形状一起广播（51,51）（100，）

def get_velocity(strength, xs, ys, X, Y,):

u = strength /(2*np.pi)*(X-xs)/((X-xs)**2 +(Y-ys)**2)
v = strength /(2*np.pi)*(Y-ys)/((X-xs)**2 +(Y-ys)**2)

return u, v

我想对某些给定点使用该函数

x_naca0012 = np.loadtxt('NACA0012_x.txt')

y_naca0012 = np.loadtxt('NACA0012_y.txt')

strength_naca0012 = np.loadtxt('NACA0012_sigma.txt')

使用该功能：

u_naca0012, v_naca0012 = get_velocity(strength_naca0012, x_naca0012, y_naca0012, X, Y)

这里出错了，我搜索错误并找到一些关于不匹配数组的东西，但我检查了每个数组的长度并匹配（100）。我注意到NACA0012_y的第50个元素是正值，第51个元素是负值。我觉得它在那里搞砸了。