matlab mod() 함수가 잘못된 답을 제공합니다.

Question

다음 코드를 사용하여 dt * j가 정수 인 경우 결과를 저장하기 위해 Matlab을 사용하려고했습니다. .

j=1; 
dt = 1E-05; 
a=[]; 
while dt*j <=20 
    if mod(dt*j,1) ==0 
     a=[a;[dt*j,j]]; 
    end 
j=j+1; 
end 

그러나 Matlab은 잘못된 결과를 제공합니다. 아래에 표시된 것처럼 일부 정수 (3,6,7,11 등 ...)가 누락되었습니다.

a = 
1  100000 
2  200000 
4  400000 
5  500000 
8  800000 
9  900000 
10  1000000 
16  1600000 
17  1700000 
18  1800000 
19  1900000 
20  2000000 

dt = 1E-4를 사용하여 다시 시도했으며 올바른 결과를 제공합니다. 아무도 여기서 무슨 일이 일어나는지 말해 줄 수 없습니까? 고맙습니다.

Answer 1

직접 ==와 부동 소수점 숫자를 비교,하지만 그들은 정말 가까운 대신 확인하지 마십시오 내가 추측하고

j=1; 
dt = 1E-05; 
a=[]; 
while dt*j <=20 
    if abs(mod(dt*j,1)) < 1e-10 
     a=round([a;[dt*j,j]]); 
    end 
j=j+1; 
end 

Answer 2

는 이상한 반올림 오류가있다. 나는 당신의 코드를 시도했고 나는 같은 결과를 얻고있다. 하나의 접근법은 mod (dt * j, 1) == 0을 사용하는 것이 아니라 mod (dt * j, 1) < dt 또는 mod (dt * j, 1) < thresh를 사용하여 임계 값을 설정할 수 있습니다. 트릭을 수행 할 정도로 제로에 가깝습니다.

배율 부동 소수점 형식 페이지의 위키 피 디아 페이지에서 다음과 같이 말합니다. 1.0000000000000002, 가장 작은 수> 1. 두 배를 사용하여 추정하면 약 2E-16의 수치 해상도를 얻을 수 있습니다. 따라서 2E-15의 임계 값을 사용하여 안전 할 수 있습니다.