issue using the rref function ?!
1 次查看(过去 30 天)
显示 更早的评论
Hello! I would like to use rref function in order to remove the linearly dependent elements from my matrix , but when using rref I get the number of basic columns as 1 which is not correct! any suggestion of any potential error in my code will be very welcomed ! thank you in advance !
load('matFile.mat', 'MyMatrix')
>> [rows, columns]=size(MyMatrix);
>> reshapedMatrix = reshape(MyMatrix,(rows * columns),1);
>> Gram = reshapedMatrix * reshapedMatrix';
>> GramDeterm = det(Gram);
>> if (GramDeterm <0.1)
[E,basiccol] = rref(Gram);
MatrixOut = Gram(:,basiccol)';
end
2 个评论
Steven Lord
2017-7-31
Do not use det to try to determine if a matrix is singular or near-singular!
In this example, A is a nonzero scalar multiple of the identity and as such is most definitely NOT singular. But dA is 0. Why? Underflow. The condition number, on the other hand, correctly identifies this matrix as well conditioned.
A = 0.1*eye(350);
dA = det(A)
cA = cond(A)
In this example the first and second rows of B are almost the same, and dB reflects this fact. But the magnitude of the elements in C (a nonzero scalar multiple of B) makes the determinant reasonably large so C appears to be nonsingular. But the condition numbers cB and cC tell the tale.
B = [1 1; 1 1+eps];
C = flintmax * B;
dB = det(B)
dC = det(C)
cB = cond(B)
cC = cond(C)
回答(0 个)
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!