我正在用Matlab编写Python脚本。碰巧在Matlab中,代码运行得更快(比如10分钟更少),所以也许你可以帮我找出我做错了什么。我用Python运行了所有会话,发现最慢的是以下两个会话:
X = xi * math.cos(theta) + zi * (math.sin(theta))
Y = yi
Z = xi * ((-1) * math.sin(theta)) + zi * (math.cos(theta))
x2 = [a * b * c * X[i] / np.sqrt((a * c * Y[i])**2 + b**2 * (
c**2 * (X[i])**2 + a**2 * (Z[i])**2)) for i in range(len(X))]
y2 = [a * b * c * Y[i] / np.sqrt((a * c * Y[i])**2 + b**2 * (
c**2 * X[i]**2 + a**2 * Z[i]**2)) for i in range(len(X))]
zt = c**2 * (1 - (x2 / a)**2 - (y2 / b)**2)
i = 0
z2 = []
for i in range(len(Z)):
z2.append((-1 if Z[i] < 0 else 1) * np.sqrt(zt[i]))
< >席、彝、子是100万余个元素的浮标
谢谢
对于z2零件
相关问题 更多 >
编程相关推荐