输入方程组的系数矩阵和常数向量A = 223; 477; -245;b = 3; 1; -7;使用matlab编写程序:使用高斯-若尔当消去法求解三阶方程组提示答案为2、-2、1:1先将原方程的增广矩阵进行输出2将经过高斯-若尔当消去法得到的对角矩阵和该对角矩阵对应的常数向量结合在一起进行输出该对角矩阵应该除主对角元素外全为03在输出结果4都要求格式化输出
% 输入方程组的系数矩阵和常数向量 A = [2,2,3; 4,7,7; -2,4,5]; b = [3; 1; -7];
% 将增广矩阵输出 augmented_matrix = [A,b]; disp("The augmented matrix is:") disp(augmented_matrix)
% 高斯-若尔当消去法 n = size(A,1); for i = 1:n-1 for j = i+1:n m = A(j,i)/A(i,i); A(j,i:n) = A(j,i:n) - mA(i,i:n); b(j) = b(j) - mb(i); end end
% 将对角矩阵和常数向量输出 diagonal_matrix = [A,zeros(n,1)]; for i = 1:n diagonal_matrix(i,:) = diagonal_matrix(i,:)/diagonal_matrix(i,i); diagonal_matrix(i,i) = 1; end disp("The diagonal matrix and corresponding constant vector are:") disp([diagonal_matrix,b])
% 回代求解 x = zeros(n,1); for i = n:-1:1 x(i) = b(i) - A(i,i+1:n)*x(i+1:n); end
% 输出结果 disp("The solution is:") fprintf("x1 = %f, x2 = %f, x3 = %f", x(1), x(2), x(3)
原文地址: https://www.cveoy.top/t/topic/cGqt 著作权归作者所有。请勿转载和采集!