主元分析法 matlab 数值分析实习作业之不选主元法高斯分解(Matlab)

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