I am trying to write my own rref and have the following code. It appears to have worked for the first column but not for all other columns for a test augmented matrix I entered. Can you please assist? Thanks!

Question

Jenel Darland 2017-4-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/337439-i-am-trying-to-write-my-own-rref-and-have-the-following-code-it-appears-to-have-worked-for-the-firs

回答： Sanjana Ramakrishnan 2017-5-8

%Loop through columns, then rows

%j=1:n-1 because matrix is augmented and I only want to check the square matrix values; however, row operations should be performed on the entire row of the matrix

for j=1:n-1

for i=1:m

%if i=j, the code already performed an operation to transform the value to 1, so I want to skip matrix(i,j) if i=j

if i~=j

        if matrix(i,j)~=0
            if matrix(i,j)>0
            matrix(i,:)=matrix(i,:)-matrix(i,j)*matrix(j,:);
            elseif matrix(i,j)<0
                matrix(i,:)=matrix(i,:)+(abs(matrix(i,j)))*matrix(j,:);
            end
        end
    end
    end
end