How to use intersect command with a tolerance value?
显示 更早的评论
I have two vector with different lenght and I want to compare the values of the vectors one by one with a certain tolerance because the numbers are similar but not equal, e.g. 222.456 and 222.389. I was using the command intersect but it only gives me the numbers that are exactly the same.
回答（1 个）
While ismembertol might be useful, I'm still puzzled by the way it defines the tolerance.
If the vectors are not too large (some thousand elements should be fine), a simple loop approach should be sufficient:
function [AI, BI] = CommonElemTol(A, B, Tol)
A = A(:);
B = B(:);
nA = numel(A);
M = zeros(1, nA);
% Collect the index of the first occurrence in B for every A:
for iA = 1:nA
dist = abs(A(iA) - B); % EDITED: Of course abs() is needed
Ind = find(dist < Tol, 1); % Absolute tolerance
% Ind = find(dist ./ A(iA) < Tol, 1); % Relative tolerance
if ~isempty(Ind)
M(iA) = Ind;
end
end
AI = find(M); % If any occurrence was found, this A exists
if isempty(AI) % prevent: Empty matrix: 0-by-1
AI = [];
end
BI = M(AI); % at this index in B
end
3 个评论
Vishal Dwivedi
2019-2-15
Stephen23
2019-2-15
I suspect that
dist = (A(iA) - B);
should be
dist = abs(A(iA) - B);
otherwise there is no absolute involved.
Jan
2019-2-15
@Stephen: Of course. Thanks.
类别
在 帮助中心 和 File Exchange 中查找有关 Operators and Elementary Operations 的更多信息
2019-2-12
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
