Comparing two matrix in Matlab

Question

john vattic 2014-6-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/137776-comparing-two-matrix-in-matlab

评论： Anuj Kumar 2020-12-24

I have two matrices x and y, both are results from different algorithms/routines that are supposed to calculate the same result. While I know that the isequal() would check if x and y are the same matrix, the entries in those matrices would not be exactly the same (i.e. some entries may be with 5% off in worst case scenario). In this scenario, what would be the best method of comparing them to see if they are close enough to be considered the same result? Thanks in advance for the advices.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Anuj Kumar 2020-12-24

you can try this,assumming x and y has same dimensions

len=lenght(x);

tolerance=x-y;

check=tolerance<0.001;

sum_Check=sum(check)

if sum_Check==len %then results can be accepted

请先登录，再进行评论。

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

Answer 1

Titus Edelhofer 2014-6-26

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/137776-comparing-two-matrix-in-matlab#answer_142591

在 MATLAB Online 中打开

Hi,

this question is rather tricky to answer in general. A simple approach is to compute the relative error. That works fine as long as there are no zeros in it, i.e.,

absError = abs(A-B);
relError = max(absError(:) ./ abs(A(:)));

This also assumes, that A and B are at least comparable in size, so it does not make a difference if you devide by A or B.

Taking care of all eventualities (NaN, inf, zeros ...) makes it more difficult to answer.

Titus

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

john vattic 2014-6-26

i got inf as result of relError

请先登录，再进行评论。

Answer 2

Anastasios 2014-6-26

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/137776-comparing-two-matrix-in-matlab#answer_142593

在 MATLAB Online 中打开

I assume you have some floating point variables, in order to have this 5% error. Try to define a function to check if they are relatively equal, such as

function test = isequalRel(x,y,tol);
test = ( abs(x-y) <= ( tol*max(abs(x),abs(y)) + eps) );
if (!  test)
    printf("fail x=%e y=%e tol=%e\n", x, y, tol);
end

So basically you do the following :

|a-b| < tol*max(|a|,|b|) + tol_floor