% 输入方程组的系数矩阵和常数向量 A = [2,2,3; 4,7,7; -2,4,5]; b = [3; 1; -7];

% 将增广矩阵进行输出 augmented_matrix = [A,b]; disp("增广矩阵为:"); disp(augmented_matrix);

% 高斯-若尔当消去法 n = length(b); for k = 1:n-1 for i = k+1:n factor = A(i,k) / A(k,k); A(i,k+1:n) = A(i,k+1:n) - factor * A(k,k+1:n); b(i) = b(i) - factor * b(k); end end disp("经过高斯-若尔当消去法后得到的对角矩阵为:"); disp(A); disp("对应的常数向量为:"); disp(b);

% 回代求解 x = zeros(n,1); x(n) = b(n) / A(n,n); for i = n-1:-1:1 x(i) = (b(i) - A(i,i+1:n)*x(i+1:n)) / A(i,i); end

% 输出结果 fprintf("方程组的解为:x1=%.2f, x2=%.2f, x3=%.2f\n",x(1),x(2),x(3))


原文地址: https://www.cveoy.top/t/topic/cGoW 著作权归作者所有。请勿转载和采集!

免费AI点我,无需注册和登录