所以我的代码遇到了一个非常奇怪的问题。这是有问题的部分:
vrt = np.matrix.transpose(vr) #np.shape(vrt) = (1832,360)
term2 = C*np.cos((beta*np.pi/180)-azrad) #np.shape(term2) = (1832,360)
test = np.shape(vrt) == np.shape(term2) #test = True
ci = vrt - term2
ValueError: operands could not be broadcast together with shapes (360,1832) (1832,360)
两个数组的大小相同(test = True),但是当我尝试减去两个矩阵时,它现在说第一个矩阵不再是(1832,360)。当我尝试这个而不采用vr的转置时,显然test = False,但错误完全相同(ValueError:操作数不能与形状一起广播(360,1832)(1832,360))。