数值代数实验题目.docVIP

数值代数实验 数值线性代数实验一 实验名称：矩阵的LU分解. 实验目的：用不选主元的LU分解和列主元LU分解求解线性方程组 Ax=b, 并比较这两种方法. 实验内容与要求 （1）用所熟悉的计算机语言将不选主元和列主元LU分解编成通用的子程序，然后用编写的程序求解下面的84阶方程组 将计算结果与方程组的精确解进行比较，并就此谈谈你对Gauss消去法的看法。 （2）写出追赶法求解三对角方程组的过程，并编写程序求该实验中的方程组 数值线性代数实验二 实验名称：实对称正定矩阵的的Cholesky分解. 实验目的：用平方根法和改进的平方根方法求解线性方程组 Ax=b. 实验内容与要求 用所熟悉的计算机语言将Cholesky分解和改进的Cholesky分解编成通用的子程序，然后用编写的程序求解对称正定方程组Ax=b，其中 b随机的选取，系数矩阵为100阶矩阵 系数矩阵为40阶Hilbert矩阵，即系数矩阵A的第i行第j列元素为，向量b的第i个分量为 用实验一的程序求解这两个方程组，并比较所有的计算结果，然后评价各个方法的优劣。 数值线性代数实验三 实验名称：矩阵的分解 实验目的：应用改进的Gram—Schmidt方法和Householder变换的方法计算矩阵A的分解. 其中 rank = 实验内容与要求 输入：的各列 输出：的各列元素（存放在的相应位置上）以及的元素 数值线性代数实验四 实验名称：用迭代法求解方程组及超松弛迭代和最佳松弛因子的确定. 实验目的：应用Jacobi迭代法、Gauss－Seidel迭代法和超松弛迭代方法求解线性方程组 Ax=b, 并选择不同的松弛因子，观察松弛因子对松弛迭代法计算效果的影响. 实验内容与要求 将常微分方程（）离散化得到差分方程Ax=b，取，应用Jacobi迭代法、Gauss－Seidel迭代法和超松弛迭代方法求解线性方程组，分别取，用SOR迭代法计算对应的数值解，并与精确解进行比较. 写出这三种迭代法求解线性方程组的步骤，并对计算结果进行分析. 实验原理 将[0,1]区间n等分，方程离散化得差分方程 差分方程对应的系数矩阵和右端项分别为 其中 SOR迭代法的迭代矩阵和常数项分别为 SOR迭代法的迭代公式： 实验过程 利用Matlab语言实现SOR迭代公式的计算，分别对于，松弛因子??（1，2）范围内，先按步长0.1搜索，在最少迭代步数附近按步长0.01再搜索. 算法终止准则： 或者或者迭代次数. 记录四种情况下的计算结果，并找出最佳松弛因子的值，并观察松弛因子与间的关系.