Vergleichen zweier Matrizen in Matlab

Question

Ich habe zwei Matrizen x und y, beide sind Ergebnisse von verschiedenen Algorithmen/Routinen, die das gleiche Ergebnis berechnen sollen. Während ich weiß, dass das isequal() überprüfen würde, ob x und y die gleiche Matrix sind, wären die Einträge in diesen Matrizen nicht genau die gleichen (d. H. Einige Einträge können im Worst-Case-Szenario mit 5% aus sein). Was wäre in diesem Szenario die beste Methode, sie zu vergleichen, um zu sehen, ob sie nah genug sind, um als dasselbe Ergebnis betrachtet zu werden? Vielen Dank im Voraus für die Ratschläge.

Answer 1

Ändern Edric's Lösung:

absTol = 1e-3; % You choose this value to be what you want! 
relTol = 0.05; % This one too! 
absError = x(:)-y(:); 
relError = absError./x(:); 
relError(~isfinite(relError)) = 0; % Sets Inf and NaN to 0 
same = all((abs(absError) < absTol) & (abs(relError) < relTol)); 

Die Variable gleiche wird falsch, wenn entweder die absolute oder der relative Fehler eines Elements ist größer als alle Toleranzen, die Sie wählen. Wenn Elemente von x zufällig genau 0 sind, dann könnten einige der Elemente von relError entweder unendlich oder nicht-eine-Zahl sein, also habe ich die ISFINITE Funktion verwendet, um diese Werte durch Setzen zu ignorieren zu 0.

Ich würde nicht vorschlagen IMAGESC Plots zu verwenden, da 1) die Daten skaliert wird, wenn es angezeigt wird, 2) die Colormap für die Anzeige hat eine diskrete Anzahl von Farbwerten (was ich denke, ist 256 standardmäßig, daher viel Rundung), und 3) feine Unterschiede in der Farbe sind möglicherweise nicht so offensichtlich, wie es beim visuellen Vergleich zweier Diagramme der Fall ist.

Answer 2

Diese Versuchen:

tf = abs((A-B)./B)<0.05 

Dies wird eine logische Matrix zurück, die für jedes Element der wahr ist, wenn der relative Unterschied zwischen A und B in Bezug auf B weniger als 5 Prozent ist.

Wenn Sie zu fragen, ob alle diese wahr sind (sie alle erfüllen die obige Bedingung):

all(tf(:)) 

Answer 3

Ich würde so etwas wie dies mit einer absoluten Toleranz sowie eine relative Toleranz zu tun:

function same = tol(x, y) 
absTol = 1e-3; 
relTol = 0.05; 
errVec = abs(x(:) - y(:)); 
same = all((errVec < absTol) | (errVec./x(:) < relTol)); 

Answer 4

Wenn Sie sehr kleine Wertepaare in x und y haben, würde das Ergebnis 0 ergeben, obwohl die Werte selbst ignoriert werden können. Also, ein Zusatz zur akzeptierten Lösung

relError(x < absTol) = 0;

kann verwendet werden, um sehr kleine Fehler zu verwerfen. Daher wird der relative Fehler für diese Werte nicht berücksichtigt.

Answer 5

nutzen ‚IsEqual (a, b) wobei a und b sind zwei Matrizen, wenn 1 es wahr ist

Answer 6

für Matrizen x und y enthält, Gleitkommazahlen, kann geprüft werden, ob Arrayelemente sind in einem gegeben Toleranz zueinander. Beispielcode:

tol = 0.05; 

result = abs(x - y) <= tol;