Using "find" for finding decimal values

Question

MiauMiau 2012-12-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/56518-using-find-for-finding-decimal-values

采纳的回答： Matt Fig

在 MATLAB Online 中打开

Hi

I use the matlab command importdata:

X = importdata('filename.csv');

to read in a csv file with three columns.

Now, for finding a specific values in the matrix X, I simple use the find command as follows:

idx = find(X(:,1 ) == 17)

But, the same seems not to be possible for decimal numbers. This for instance would not work:

idx = find(X(:,1 ) == 17.9203)

even though 17.9203 is to be found in the original csv file. What is the problem and what can I do?

Thanks

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Matt Fig 2012-12-13

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/56518-using-find-for-finding-decimal-values#answer_68448

编辑：Matt Fig 2012-12-13

在 MATLAB Online 中打开

MiauMiau, the numbers you are looking for are obviously different from the short format you see. Do this:

[~,idx] = min(abs(X(:,1) - 6.0018));
Y = X(idx,1);
[~,idx] = min(abs(X(:,1) - 17.9203));
fprintf('%15.15f   %15.15f\n',Y,X(idx,1))

And show us the ouput in a comment on this answer, don't add another answer!

I have a feeling you will need to set your tolerance much higher, like 10^-4.

4 个评论
显示 2更早的评论隐藏 2更早的评论

MiauMiau 2012-12-13

Thanks - it works perfectly, wow! I thought there is a number like 17.9203 because I used two input files which were supposed to be identical - obviously there were not. May I ask how you found out the tolerance has to be set to 10^-3? I don't understand that yet. Thank you very much!

Matt Fig 2012-12-14

I found out by looking at how close your guess value was to the actual value. You may have to do some tweaking but you get the idea now.

请先登录，再进行评论。

Answer 2

Azzi Abdelmalek 2012-12-13

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/56518-using-find-for-finding-decimal-values#answer_68415

在 MATLAB Online 中打开

Try

idx = find(abs(X(:,1 )-17.9203)<eps)

2 个评论
显示无隐藏无

Jan 2012-12-13

编辑：Jan 2012-12-13

在 MATLAB Online 中打开

Substracting two numbers in the magnitude of 20 can have a roundoff error of eps(20). Then this is better:

idx = find(abs(X(:,1 ) - 17.9203) <= eps(17.9203)) % EDITED: was "<"

Even 10*eps(17.9203) would be a reliable limit.

[EDITED] Thanks, Kye. I cannot test this currently. But I assume that even "<=" can fail, when the values are at the limits of 2^n. Does eps(16 + eps(16)) reply eps(32)? Then a factor of 2 would be obligatory.

Kye Taylor 2012-12-13

should be <=

请先登录，再进行评论。

Answer 3

Jan 2012-12-13

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/56518-using-find-for-finding-decimal-values#answer_68422

编辑：Jan 2012-12-13

在 MATLAB Online 中打开

There is no exact representation of decimal floating point numbers in binary format for all values. You find a lot of corresponding discussion in this forum:

0.1 + 0.2 - 0.3 == 0
>> false

This is no bug, but the expected behaviour, when floating point numbers are represented with a limited precision.

A consequence is, that you cannot compare numbers like 17.9203 and 17.92029999999999999999 sufficiently and even the display in the command window can be rather confusing.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 4

MiauMiau 2012-12-13

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/56518-using-find-for-finding-decimal-values#answer_68427

在 MATLAB Online 中打开

I was not aware of that, that is pretty horrible! Anyway, both of your suggestions did not work, so also for:

idx = find(abs(X(:,1 ) - 17.9203) < eps(17.9203))

I did get an empty matrix returned - although the value is to be found in X. Also, when I changed the number a bit, I then got too many answers returned. But what strikes me most, is that in the first case, I just did not get any answer at all!..?!

7 个评论
显示 5更早的评论隐藏 5更早的评论

Azzi Abdelmalek 2012-12-13

Yes, I said Matlab because we are working with Matlab. I think, representing real numbers for numerical programming is almost similar to what an ADC (analogic-digital converter) do. There is a quantification of a real number, which means there is an error of quantification which depends on the number of used bits and method.

Jan 2012-12-13

@Azzi: I've stressed this detail, because the OP seems to be confused about this topic already. I did not assume, that you struggle with floating point arithmetics.

请先登录，再进行评论。

Answer 5

MiauMiau 2012-12-13

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/56518-using-find-for-finding-decimal-values#answer_68446

在 MATLAB Online 中打开

I am really sorry for my use of the word "horrible".

If you meant the following command though ( I have used another nummerical value):

idx = find(abs(X(:,1 ) - 6.0018) <= eps(6.0018))

I - again - got an empty 0,1-matrix.

Also, I was not criticising the obviously necessary restriction on measurement and computation precision but I was amazed that there is not a simpler command to handle that. Anyway. What shall I try next?