1,用不选主元法今儿选主元法的高斯消去法求解下列方程组,并记下变换后的增广矩阵:
(3)
0.729x1+0.81x2+0.9x3 = 0.6867
X1+x2+x3 = 1
1.331x1+1.21x2+1.1x3 = 1
Matlab程序:
(1)不选主元的高斯消去法求解
创建主程序:
clear
clc
A =input('输入线性方程组 Ax = b 的A = ');
b = input('输入线性方程组 Ax = b 的b = ');
x = fguass(A,b);
disp('通过不选主元的高斯消去法的到 x = ')
disp(x)
fguass函数
function x = fguass(A,b)
tic
A_b = [A,b];
[m,n] = size(A_b);
for k = 1:1:m-1;
for i = k+1:1:m
c = A_b(i,k)/A_b(k,k);
A_b(i,:) = A_b(i,:) - c*A_b(k,:);
end
end
disp('不选主元的高斯消去后的增广矩阵 A_b = ')
disp(A_b)
x(m) = A_b(m,n)/A_b(m,m);
for k = m-1:-1:1
x(k) = (A_b(k,n) - A_b(k,k+1:m)*x(k+1:m)')/A_b(k,k);
end
toc