1不等于1？ MATLAB 中的矩阵条目比较

Question

我正在尝试通过乘以 1 来测试矩阵的求逆。我检查它是否等于 1，如果不等于，则打印一条消息和矩阵条目的值。

问题是打印的矩阵条目是 1。我什至使用我知道的舍入函数来确保浮点不准确没有问题。

这是我的代码：

 function [D1,f]=stiff(D,y)
 n = size(D,2);
 D1 = inv(D);
 f = D1 * y;
 for i = 1:n
        for j = 1:n
            if j == i % checking for diagonal entry
                yTempVar = D * D1 (1:n, i);
                if truncateFunction (yTempVar(j), 1) ~= 1
                    fprintf ("something's wrong")
                    truncateFunction (yTempVar(j), 1)
                 return
             end
         end
     end
 end

function B=truncateFunction(A,pow)
A(abs(A)<10^-pow)=0;
B=A;

对应的输出是这样的：

something's wrong
ans = 1.0000

不带舍入功能我也试过，我也试过对1.0000

这是我在实时脚本中用来测试的内容

D=rand(4)*0.1;
y=randi([10 20], 4,1);
stiff(D,y)

Answer 1

truncateFunction 并不像您想象的那样。它只会将低于给定的 10 次方的值四舍五入为零 - 它根本不会四舍五入到最接近的整数，除非该整数为零。

一般来说，对于浮点问题（这是），最好检查值是否在一个小值内而不是相等。

epsilon = .00001
if abs(yTempVar(j)-1) > epsilon
    fprintf ("something's wrong")
    truncateFunction (yTempVar(j), 1)
    return

1不等于1？ MATLAB 中的矩阵条目比较

1 doesn't equal 1? Matrix entry comparison in MATLAB

floating-point

matlab

equality