我正在进行以下计算:
N = 2**15
dx = 0.1
x = np.arange(-N/2,N/2)
u0 = np.zeros([N, 1])
L = N * dx
x0 = x[1] + 2 * delta
delta = 15
while x0 < L - delta:
l1 = 1.267;
x0 = x0 + delta
r = 1/(l1*np.cosh(x)**2)
u0 = r + u0
基本上,x0<L-delta这个循环将运行2^15分。
这可以很好地翻译成MATLAB,但python给了我以下错误:
RuntimeWarning:在cosh中遇到溢出
r=1/(l1*np.cosh(x)**2)
