Matlab增量錯誤？

Question

我不知道發生了什麼，但是當我在代碼中使用i = 0.1:0.1:7時，我得到整數1和2，然後跳過3，但得到4,5,6和7沒有問題。

x=zeros((t_f/h)+1,1); 
x(1,1)=0; 

table=zeros(t_f,1); 

for i=0.1:0.1:7; 
     x(round(i/h)+1,1)=i; 
     if ~mod(i,1) 
      table(i,1)=i; 
     end 
end 

然後測試我回到這些

table 
a=[x(1) x(11) x(21) x(41) x(51) x(61) x(71)] 
a=[x(1) x(11) x(21) x(31) x(41) x(51) x(61) x(71)] 

它沒有找到3爲一個整數，因爲我從來是3，這是3.000 ...但1，2，4，5，6和7是整數。

Answer 1

它是一個浮點數表示問題。以這種方式構建數字時，您期望得到的數字恰好爲3（4.440892e-16）。這是預期的行爲，而不是bug，請參閱：http://docs.oracle.com/cd/E19957-01/806-3568/ncg_goldberg.html