输入方程组的系数矩阵和常数向量A = 223; 477; -245;b = 3; 1; -7;使用matlab编写程序:使用高斯-若尔当消去法求解三阶方程组提示答案为2、-2、1:1先将原方程的增广矩阵进行输出2将经过高斯-若尔当消去法得到的对角矩阵和该对角矩阵对应的常数向量结合在一起进行输出3在输出结果4都要求格式化输出
% 输入方程组的系数矩阵和常数向量 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 著作权归作者所有。请勿转载和采集!