I have two equal vectors a = b. When I check for equality:
y = eq(a, b) find(y==0)
I get most of values, meaning they are NOT equal.
What would be the reason for this?
N.
No products are associated with this question.
To compare two vectors use:
r = isequal(a, b)
or to find the difference between vectors (or other arrays) of the same size:
any(a - b)
Note, that the later replies FALSE even for ones(1,10) and ones(1,1).
If you only assume, that the vectors have equal values, check a - b and note, that the limited precision of the type double leads to effects like:
03. - 0.2 ~= 0.1
Then check "almost equal" by:
limit = eps(max(max(abs(a(:)), abs(b(:))))); % or: eps(max(abs(a(:)), abs(b(:)))); % or: the same times 10 a_eq_b = isequal(size(a), size(b)) & all(abs(a - b) < limit);
The choise of the limit depends on the physical meaning of the variables. There is no unique definition to distinguish measurement noise from noise caused by numerical artifacts for a tiny values.
5 Comments
Direct link to this comment:
http://mathworks.com/matlabcentral/answers/44021#comment_90476
What does the "find(a==0)" do in your code?
Direct link to this comment:
http://mathworks.com/matlabcentral/answers/44021#comment_90890
Sorry, it was wrong put, "a" should be the output of eq(a,b).
Direct link to this comment:
http://mathworks.com/matlabcentral/answers/44021#comment_90908
Then please edit the question instead of adding this important information as a comment.
Direct link to this comment:
http://mathworks.com/matlabcentral/answers/44021#comment_90983
Thanks, I did.
Direct link to this comment:
http://mathworks.com/matlabcentral/answers/44021#comment_90985
@Nuchto: Is your problem solved?