数值分析上机作业2.docVIP

实验一

实验名称：矩阵的LU分解.

实验目的：用不选主元的LU分解和列主元LU分解求解线性方程组Ax=b,并比拟这两种方法.

实验内容与要求

〔1〕用所熟悉的计算机语言将不选主元和列主元LU分解编成通用的子程序，然后用编写的程序求解下面的84阶方程组

将计算结果与方程组的精确解进行比拟，并就此谈谈你对Gauss消去法的看法。

〔2〕写出追赶法求解三对角方程组的过程，并编写程序求该实验中的方程组

实验步骤：

程序：

不选主元：

functionx=GAuss(A,b)

n=length(b);

x=b;

fori=1:n-1

forj=i+1:n

mi=A(j,i)/A(i,i);

b(j)=b(j)-mi*b(i);

fork=i:n

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

end

AB=[A,b];

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

选主元：

functionx=GAussZhu(A,b)

n=length(b);

x