线性方程组的求解.doc

南昌航空大学数学与信息科学学院实验报告

第PAGE6页

一、实验目的

了解列主元消去法、LU分解法、Jacobi迭代法、Gauss-Seidel迭代法,并用以上方法解决线性方程组。

二、实验原理、程序框图、程序代码等

1．实验原理：

（1）列主元高斯消去法：

设已完成k-1步消元，即有，利用，确定第k不列主元所在的位置，在交换行和k行后，再进行消元。根据矩阵理论，交换和k两个方程的位置，相当于再系数矩阵和有段顶上都左乘初等排列矩阵，因此列主元素的消去过程相当于对交换后的新矩阵进行消元，即

同时，右端向量变化为

列主元gauss消去法的矩阵解释为

（2）LU三角分解法

先直接将矩阵A分解为A=LU，则求解的问题就等价于解两个三角形方程组：

=1\*GB3①Ly=b,求y

=2\*GB3②Ux=y，求x

求矩阵A的LU分解分选主元和不选主元两种情况。

(i)不选主元的三角分解法

直接利用矩阵乘法，得到矩阵的三角分解计算公式为：

由上两式得到矩阵A的LU分解后，求解Ux=y的计算公式为

(ii)选主元的三角分解法

对不选主元的三角分解法或很小时，该分解方法将会失效。与Gauss列主元消去法相对应，可以通过行交换实现矩阵A的LU分解，从而构造出线性方程组的选主元分解法。

设第k-1步分解已经完成，这时有

A→

第K步分解需要用到做除数，可以引进变量则有。比较的大小，如有，此时可以交换A的第K行和行元素，然后再进行第k步分解计算。

（3）jacobi迭代法

考虑n阶线性代数方程组

其矩阵形式为

设该方程组的系数矩阵A非奇异且，可将A分解为：

其中

，U=－

然后化为如下等价形式

简记为

选取初始向量,通过上式可得到线性代数方程组的迭代格式

其分量形式为

即

以上计算过程称迭代法，矩阵陈为jacobi迭代法的迭代矩阵。

(4)Gauss-Seidel迭代法

将jacobi迭代格式改为如下形式

其矩阵形式为

于是得Gauss-Seidel迭代公式为

称为Gauss-Seidel迭代法的迭代矩阵。

2.实验题目：

一．用列主元消去法解方程组：

二．用LU分解法界方程组，其中

三．分别用Jacobi迭代法和：Gauss-Seidel迭代法求解方程组

3.程序代码：

MATLAB程序如下：

1.列主元消去法解

functionx=hanshu1(a,b)

n=length(b);

x=ling(n,1);

c=ling(1,n);

d=0

fori=1:n-1

max=abs(a(i,i));

m=i;

forj=i++:n

ifmaxabs(a(j,i))

max=abs(a(j,i));

m=j;

end

end

if(m≠i)

fork=i:n

c(k)=a(i,k);

a(i,k)=a(m,k);

a(m,k)=c(k);

end

d=b(i);

b(i)=b(m);

b(m)=d;

end

fork=i+1:n

forj=i+1:n

a(k,j)=a(k,j)-a(i,j)*a(k,i)/a(i,i);

end

b(k)=b(k)-b(i)*a(k,i)/a(i,i);

a(k,i)=0;

end

end

x(n)=b(n)/a(n,n);

fori=n-1:-1:1

sum=0;

forj=i++:n

sum=sum+a(i,j)*x(j)