0线性方程的数值解法试卷.docx

姓名评分实验报告 实验报告 课程名称： 数值分析 课题名称： 线性方程的数值解法 专 业： 姓 名： 班 级： 《数值分析》实验报告 实验报告 一、实验名称 线性方程的数值解法 二、 实验目的 掌握解线性方程的直接法，并编程实现LU分解、三角分解和追赶法； 掌握解线性方程的迭代法，分析比较雅可比迭代和赛德尔迭代的特点； 理解不同ω对松弛迭代法的影响; 学会比较解线性方程的各种算法. 三、 实验要求 编程实现LU分解、三角分解和追赶法并比较其特点； 利用不同实验体验各种线性方程解法的特性. 四、 实验原理 矩阵的LU分解法： 设方程组Ax=b的系数矩阵非奇异，且各阶顺序主子Dk≠0,则矩阵A可以分解为：A=LU； L=10…0l211…0…………ln1ln2…1; U=u11u12…u1n0u22…u2n…………00…unn Ly=b; Ux=b; 三角分解法： 若A为对称正定阵，则可以把A分解为：A=GGT或A=LDLT; G=g110…0g21g22…0…………gn1gn2…gnn ; L=10…0l211…0…………ln1ln2…1 ; D=u11?00u220???0?unn ; 追赶法： A=b1c1…0a1b2…0………cn-