溢出exp，python

Question

无法确定此错误RuntimeWarning: overflow encountered in exp出现的原因。我试图实现的功能是：

Id = lambda t_u, yp: Is * (np.exp((Vin(t_u) - L*yp)/(n_Ut*Ut)) - 1.0)

值：

Vin = lambda t: Vo * np.sin(2*np.pi*w*t)

L = 50e-3  # 50 mH
Vo = 5  # 5 V
w = 50  # 50 Hz
Is = 1e-9   # 1 nA
Ut = 25e-3  # 25 mV
n_Ut = 1.0

函数Id是我试图用Runge-Kutta-Method解决的ODE的一部分。

编辑：

Methoden\Serien\6>python circuit.py
Traceback (most recent call last):
  File "circuit.py", line 148, in <module>
    perform_experiment(exrk_o5(), "exrk_o5")
  File "circuit.py", line 48, in perform_experiment
    t, y = run_method(method, label)
  File "circuit.py", line 69, in run_method
    t, y = method.__call__(circuit_rhs, y[0, :], t_end, n_steps)
  File "C:\Users\-\numerical_methodes\6\rk.py", line 49, in __call__
    t[k+1],y[k+1, :] = self.step(f, y[k, :], t[k], dt)
  File "C:\Users\-\numerical_methodes\6\rk.py", line 102, in step
    dydt[:, i] = rhs(t + c[i] * dt, y0 + dt * np.dot(A[i, :].T, dydt.T))
  File "circuit.py", line 38, in circuit_rhs
    dydt=  np.array ([y[1] , Id(t, y [1]) /(C*L) - y [1]/( R*C) - y [0]/( C*L)])
  File "circuit.py", line 27, in <lambda>
    Id = lambda t_u, yp: Is * (np.exp((Vin(t_u) - L*yp)/(n_Ut*Ut)) - 1.0)
FloatingPointError: overflow encountered in exp

其中y [1] = 0且y [0] = 0，作为第一个值。 rk.py文件只是Runge-Kutta-Method的实现。

Answer 1

Vin可以返回高达1.0的值;然后你减去一些东西，导致一个完全未知的数量。你将它乘以40，减去1，并将epsilon提高到该功率。根据yp的值，在此过程中有足够的溢出空间。 yp的足够负值将导致问题而不参考中间计算。

numpy.exp（n）溢出了700-800范围内的某个值。这表明yp可能<= -350。