Ответы даны все правильно. Я просто хотел рассказать о замечании gnovice's о тестировании с плавающей точкой.
При сравнении чисел с плавающей запятой для равенства необходимо использовать значение допуска. Обычно используются два типа сравнения толерантности: абсолютный допуск и относительная толерантность. (source)
Абсолютное сравнение толерантность a
и b
выглядит следующим образом:
|a-b| < tol
Сравнение относительной толерантности выглядит следующим образом:
|a-b| < tol*max(|a|,|b|) + tol_floor
Вы можете осуществить вышеуказанные два, как анонимные функции:
%# absolute tolerance equality
isequalAbs = @(x,y,tol) (abs(x-y) <= tol);
%# relative tolerance equality
isequalRel = @(x,y,tol) (abs(x-y) <= (tol*max(abs(x),abs(y)) + eps));
Тогда вы можете использовать их как:
%# let x and y be scalars/vectors/matrices of same size
x == y
isequalAbs(x, y, 1e-6)
isequalRel(x, y, 1e-6)
Всегда, всегда векторизовать код MATLAB, если это возможно. – Doresoom
Должно быть перемещено на http://math.stackexchange.com/? – aaronsnoswell