So sánh hai ma trận trong Matlab

Question

Tôi có hai ma trận x và y, cả hai đều là kết quả từ các thuật toán/thường trình khác nhau được cho là để tính cùng một kết quả. Trong khi tôi biết rằng isequal() sẽ kiểm tra nếu x và y là cùng một ma trận, các mục trong các ma trận đó sẽ không chính xác giống nhau (nghĩa là một số mục có thể giảm 5% trong trường hợp xấu nhất). Trong kịch bản này, phương pháp nào tốt nhất để so sánh chúng để xem liệu chúng có đủ gần để được xem là kết quả tương tự không? Cảm ơn trước cho lời khuyên.

Answer 1

Sửa đổi Edric's giải pháp:

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)); 

Các biến cùng sẽ là false nếu hoặc tuyệt đối hoặc lỗi tương đối của bất kỳ phần tử nào lớn hơn bất kỳ dung sai nào bạn chọn. Ngoài ra, nếu bất kỳ phần tử nào của x xảy ra chính xác 0, thì một số thành phần của relError có thể kết thúc hoặc là vô hạn hoặc không phải số, vì vậy tôi đã sử dụng hàm ISFINITE để bỏ qua các giá trị đó bằng cách đặt chúng đến 0.

Tôi sẽ không đề xuất sử dụng IMAGESC để so sánh ô, vì 1) dữ liệu được thu nhỏ khi hiển thị, 2) bản đồ màu cho màn hình có số lượng màu riêng biệt (mà tôi nghĩ là 256 theo mặc định, do đó rất nhiều làm tròn) và 3) các biến thể tinh tế về màu sắc có thể không rõ ràng so với sự so sánh trực quan của hai ô.

Answer 2

Hãy thử điều này:

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

này sẽ trả về một ma trận logic đó sẽ đúng đối với mỗi phần tử nếu sự chênh lệch tương đối giữa A và B liên quan đến B với ít hơn 5 phần trăm.

Nếu bạn muốn hỏi nếu tất cả trong số này là đúng (tất cả họ đều đáp ứng các điều kiện nêu trên):

all(tf(:)) 

Answer 3

tôi sẽ xem xét làm một cái gì đó như thế này với một sự khoan dung tuyệt đối cũng như khả năng chịu đựng tương đối:

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

Answer 4

Khi bạn có các cặp giá trị rất nhỏ trong x và y, kết quả sẽ trả về 0 mặc dù các giá trị có thể tự bỏ qua. Do đó, việc thêm vào giải pháp được chấp nhận

relError(x < absTol) = 0;

có thể được sử dụng để loại bỏ các lỗi rất nhỏ. Do đó, lỗi tương đối không được xem xét cho các giá trị này.

Answer 5

làm sử dụng 'isequal (a, b), nơi a và b là hai ma trận, nếu 1 đó là sự thật

Answer 6

Đối với ma trận x và y có chứa các giá trị dấu chấm động, bạn có thể kiểm tra xem phần tử mảng nằm trong một chịu được sự khoan dung lẫn nhau. Mã mẫu:

tol = 0.05; 

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