南京邮电大学数值代数实验.pdfVIP

.\ 数值代数实验 数值线性代数实验一 一、实验名称：矩阵的 LU 分解 . 二、实验目的：用不选主元的 LU 分解和列主元 LU 分解求解线性方程组 Ax=b, 并比较这 两种方法 . 三、实验内容与要求 （1）用所熟悉的计算机语言将不选主元和列主元 LU 分解编成通用的子程序，然后用编写 的程序求解下面的 84 阶方程组 将计算结果与方程组的精确解进行比较，并就此谈谈你对 Gauss 消去法的看法 . （2 ）写出追赶法求解三对角方程组的过程，并编写程序求该实验中的方程组 Gauss消去法： 用消去法解方程组的基本思想是用逐次消去未知数的方法把原来方程组 Ax=b 化为与其等价的三角方程组，而求解三角方程组就容易了。换句话说，上 述过程就是用行的初等变换将原方程组系数矩阵化为简单形式， 从而将求解原方 程组的问题转化为求解简单方程组的问题。 利用 Gauss消去法对线性方程组 Ax=b 进行求解。 用 MATLAB 建立 m 文件 DelGauss.m,程序如下： function x=DelGauss(a,b) [n,m]=size(a); nb=length(b); det=1; x=zeros(n,1); for k=1:n-1 for i=k+1:n if a(k,k)==0 return end .\ m=a(i,k)/a(k,k); for j=k+1:n a(i,j)=a(i,j)-m*a(k,j); end b(i)=b(i)-m*b(k); end det=det*a(k,k); end det=det*a(n,n); for k=n:-1:1 for j=k+1:n b(k)=b(k)-a(k,j)*x(j); end x(k)=b(k)/a(k,k); End 在 matlab 中输入如下： 结果如下： .\ .\ 方程组的精确解为 x 1 2 84 =x =…=x =1.0000，与 Gauss消去法求得的解差距很大， 所 得结果不够准确， 计算简单但其消元过程有时不能进行到底而使求解出现解失真 的情况。 数值线性代数实验二 一、实验名称：实对称正定矩阵的 A 的 Cholesky 分解 . 二、实验目的：用平方根法和改进的平方根方法求解线性方程组 Ax=b . 三、实验内容与要求 用所熟悉的计算机语言将 Cholesky 分解和改进的 Cholesky 分解编成通用的子程序，然后用 编写的程序求解对称正定方程组 Ax=b ，其中 （1） b 随机的选取，系数矩阵为 100 阶矩阵 （2） 系数矩阵为 40 阶 Hilbert 矩阵，即系数矩阵 A 的第 i 行第 j 列元素为 ，向量 b 的第 i 个分量为 （3） 用实验一的程序求解这两个方程组，并比较所有的计