def func_i_beta(r, l, lam):
ib_p = []
for i in range(len(parm.TIME)):
ib = parm.get_i_beta_at(i)
w = theta = parm.get_theta_at(i)
ib_p.append(r / l * ib + parm.P * lam / l * w * math.cos(theta) +
1 / l * parm.U_BETA)
ib_p = np.asarray(ib_p)
return ib_p * SCALER
def func_w(j, lam):
w = []
for i in range(len(parm.TIME)):
ia = parm.get_i_alpha_at(i)
ib = parm.get_i_beta_at(i)
theta = parm.get_theta_at(i)
ci = parm.get_ci_at(i)
w.append(-(lam / j * (-ia * math.sin(parm.P) * theta +
ib * math.cos(parm.P) * theta) -
parm.F / j * theta - ci / j))
w = np.asarray(w)
return w * SCALER