Gauss列主元消去法程序设计.docVIP

Gauss列主元消去法程序设计

设，第i个方程减去第k个方程的倍，，得到.

经过n-1次消元，可把方程组化为上三角方程组.

(2)回代过程：

以解如下线性方程组为例测试结果.

2.列主元消去法

由高斯消去法可知，在消元过程中可能出现的情况，这是消去法将无法进行，即使主元素但很小时，用其作除数，会导致其他元素数量级的严重增长和舍入误差的扩散，最后也使得计算解不可靠.这时就需要选取主元素，假定线性方程组的系数矩阵A是菲奇异的.

(1)消元过程：

对于，进行如下步骤：

按列选主元，记

交换增广阵A的p,k两行的元素。

A(k,j)=A(p,j)(j=k,…,n+1)

3)交换常数项b的p,k两行的元素。

b(k)=b(p)

4)计算消元

(2)回代过程

(3)以解如下线性方程组为例测试结果.

三、实验环境

MATLABR2014a

四、实验步骤

1.高斯列主元消去法流程图：

开始

开始

输入系数阵a和常数项b

按列选主元

交换元素

计算消元

回代

输出线性方程组的解

结束

2.程序设计：

(一)高斯消去法：

a=input(请输入系数阵：);

b=input(请输入常数项：);

n=length(b);

A=[a,b];

x=zeros(n,1);%初始值

fork=1:n-1

fori=k+1:n%第k次消元

m(i,k)=A(i,k)/A(k,k);

forj=k+1:n

A(i,j)=A(i,j)-A(k,j)*m(i,k);

end

b(i)=b(i)-m(i,k)*b(k);

end

end

x(n)=b(n)/A(n,n);%回代

fori=n-1:-1:1;

s=0;

forj=i+1:n;

s=s+A(i,j)*x(j);

end

x(i)=(b(i)-s)/A(i,i)

end

(二)高斯列主元消去法：

a=input(请输入系数阵：);

b=input(请输入常数项：);

n=length(b);

A=[a,b];

x=zeros(n,1);%初始值

fork=1:n-1

ifabs(A(k,k))10^(-4);%判断是否选主元

y=1

else

y=0;

end

ify;%选主元

fori=k+1:n;

ifabs(A(i,k))abs(A(k,k))

p=i;

elsep=k;

end

end

ifp~=k;

forj=k:n+1;

s=A(k,j);

A(k,j)=A(p,j);%交换系数

A(p,j)=s;

end

t=b(k);

b(k)=b(p);%交换常数项

b(p)=t;

end

end

fori=k+1:n

m(i,k)=A(i,k)/A(k,k);%第k次消元

forj=k+1:n

A(i,j)=A(i,j)-A(k,j)*m(i,k);

end

b(i)=b(i)-m(i,k)*b(k);

end

end

x(n)=b(n)/A(n,n);%回代

fori=n-1:-1:1;

s=0;

forj=i+1:n;

s=s+A(i,j)*x(j);

end

x(i)=(b(i)-s)/A(i,i)

end

五、实验结果

Gauss1

请输入系数阵：[10,-7,0;-3,2,6;5,-1,5]

请输入常数项：[7;4;6]

x=

0

-1.0000

1.0000

x=

-0.0000

-1.0000

1.0000

X=（0，-1，1）

Gauss2

请输入系数阵：[10^(-5),1;2,1]

请输入常数项：[1;2]

y=

1

x=

0.5000

1.0000

X=（0.5，1）

六、实验讨论、结论

本实验通过matlab程序编程实现了高斯消去法及高斯列主元消去法的求解，能加深对高斯消去法基本思路与计算步骤的理解。当主元素特别小时，需要选取主元，否则会影响结果，这时就需要采用高斯列主元消