런타임 경고: 분할에서 잘못된 값이 발견되었습니다.

런타임 경고: 분할에서 잘못된 값이 발견되었습니다.

저는 "스프링 속의 공" 모델에 대한 오일러의 방법을 사용하여 프로그램을 만들어야 합니다.

from pylab import*
from math import*
m=0.1
Lo=1
tt=30
k=200
t=20
g=9.81
dt=0.01
n=int((ceil(t/dt)))
km=k/m
r0=[-5,5*sqrt(3)]
v0=[-5,5*sqrt(3)]
a=zeros((n,2))
r=zeros((n,2))
v=zeros((n,2))
t=zeros((n,2))
r[1,:]=r0
v[1,:]=v0
for i in range(n-1):
    rr=dot(r[i,:],r[i,:])**0.5
    a=-g+km*cos(tt)*(rr-L0)*r[i,:]/rr
    v[i+1,:]=v[i,:]+a*dt
    r[i+1,:]=r[i,:]+v[i+1,:]*dt
    t[i+1]=t[i]+dt

    #print norm(r[i,:])

plot(r[:,0],r[:,1])
xlim(-100,100)
ylim(-100,100)
xlabel('x [m]')
ylabel('y [m]')

show()

다음 오류가 계속 발생합니다.

a=-g+km*cos(tt)*(rr-L0)*r[i,:]/rr
RuntimeWarning: invalid value encountered in divide

알 수가 없네요, 코드가 뭐가 잘못된 건가요?

당신의 코드가 "0으로 나누기" 또는 "NaN으로 나누기"를 시도하고 있다고 생각합니다.이를 알고 있으면서도 문제가 발생하지 않도록 하려면 다음을 시도할 수 있습니다.

import numpy as np
np.seterr(divide='ignore', invalid='ignore')

자세한 내용은 다음을 참조하십시오.

• http://docs.scipy.org/doc/numpy/reference/generated/numpy.seterr.html

Python 인덱싱은 (1이 아닌) 0에서 시작되므로 할당 "r[1,:] = r0"은 r의 두 번째(즉, 인덱스 1) 요소를 정의하고 첫 번째(인덱스 0) 요소를 0 쌍으로 둡니다.for 루프에서 i의 첫 번째 값은 0이고, 첫 번째 항목의 점 곱의 제곱근을 r로 r을 구하면(즉 0) 후속 라인에서 r을 기준으로 나눗셈하면 오류가 발생합니다.

0으로 나눗셈하는 것을 방지하기 위해 div0 오류가 발생하는 출력 'out'을 사전 초기화할 수 있습니다.np.where완전한 선은 조건에 관계없이 평가되므로 절단하지 않습니다.

사전 초기화를 사용한 예:

a = np.arange(10).reshape(2,5)
a[1,3] = 0
print(a)    #[[0 1 2 3 4], [5 6 7 0 9]]
a[0]/a[1]   # errors at 3/0
out = np.ones( (5) )  #preinit
np.divide(a[0],a[1], out=out, where=a[1]!=0) #only divide nonzeros else 1

다음으로 나눕니다.rr0.0일 수 있습니다.확인 여부rr0이고 분모에 사용하는 것 이외에 합리적인 방법을 사용합니다.

언급URL : https://stackoverflow.com/questions/14861891/runtimewarning-invalid-value-encountered-in-divide