Intento graficar una funcion de una curva par vs velocidad de un motor de inducción, pero me salta el error del título, estoy tratando de generar la gráfica para poder ver su curva. ¿Cómo puedo solucionar este error ó qué estoy haciendo mal? el error me salta en el bucle for
, en la ecuacion de torque(ii)
Existe un código en MATLAB
que es por el que me estoy guiando, donde ya muestra esa caracteristica pero debo hacerlo en python
.
import math
import numpy as np
import matplotlib.pyplot as plt
#valores requeridos del programa
r1 = 0.641 #resistencia del estator.
r2 = 0.332 #resistencia del rotor.
x1 = 1.106 #reactancia del estator.
x2 = 0.464 #reactancia del rotor.
xm = 26.3 #reactancia del circuito tanque.
voltaje_fase = 460/math.sqrt(3)
n_sincronismo = 1800
n_sincronismo_rad = 188.5
#Calcuamos el voltaje e impedancia de Thevenin.
vth = voltaje_fase * (xm /math.sqrt(r1**2 + (x1 + xm)**2))
zth = ((xm*1j)*(r1 + x1*1j)) / (r1 + (x1 + xm)*1j)
rth = zth.real
xth = zth.imag
#Calculamos la caracteristica par vs velocidad.
s = np.arange(0, 51, 0.001)
nm = (1 - s)*n_sincronismo
for ii in range(1,51):
torque(ii) = ((3*vth**2*r2)/s(ii)) / (n_sincronismo_rad*(rth +
(r2*2)/s(ii))**2 + (xth + x2)**2)
plt.plot(nm,torque)
plt.grid(True)
plt.show()
[]
, no estorque(ii)
os(ii)
, estorque[ii]
ys[ii]
..append()
o una comprensión de lista