用matlab写Guass-Jordan消元代码

好的，以下是一个简单的 Matlab Guass-Jordan 消元代码：

function [augMat] = gaussJordan(augMat)
% Gauss-Jordan elimination with partial pivoting
% Input: augmented matrix of the system [A b]
% Output: row-reduced augmented matrix [R d]

[m,n] = size(augMat);
% for each column in augmented matrix
for j = 1:n-1
% partial pivoting
[maxVal, maxRow] = max(abs(augMat(j:m,j)));
maxRow = maxRow + j-1;
if maxRow ~= j
augMat([j maxRow],:) = augMat([maxRow j],:);
end
% elimination
for i = j+1:m
factor = augMat(i,j) / augMat(j,j);
augMat(i,j:n) = augMat(i,j:n) - factor * augMat(j,j:n);
end
end
% back substitution
for i = m:-1:1
augMat(i,:) = augMat(i,:) / augMat(i,i);
for j = i-1:-1:1
factor = augMat(j,i) / augMat(i,i);
augMat(j,:) = augMat(j,:) - factor * augMat(i,:);
end
end

end

请注意这只是一个简单的示例，只适用于未知数数量少的情况。对于更多复杂的情况，您可能需要改进算法或使用其他库。