Gauss列主元消去法程序设计.pdfVIP

《Gauss 列主元消去法》实验报告 实验名称： Gauss 列主元消去法程序设计 ??? 成绩： ___________ 专业班级：数学与应用数学 1202 班 ? 姓名：王晓阳 ???学号： 实?验? 日?期 ：?2014?年 11 月 10 日 实验报告日期： ?2014 年?11 月 10 日 一．实验目的 1. 学习 Gauss消去法的基本思路和迭代步骤 . 2. 学会运用 matlab 编写高斯消去法和列主元消去法程序，求解线性方程组 . 3. 当 akk k 绝对值较小时，采用高斯列主元消去法 . 4. 培养编程与上机调试能力 . 二、实验内容 用消去法解线性方程组的基本思想是用逐次消去未知数的方法把原线性方程组 Ax b 化为与其等价的三角形线性方程组，而求解三角形线性方程组可用回代的方法求 解 . 1. 求解一般线性方程组的高斯消去法 . (1) 消元过程： k k k 设 a 0 ，第 i 个方程减去第 k 个方程的 m a / a 倍， (i k 1, , n) ，得 kk ik ik kk 到 A k 1 x b k 1 . 经过 n-1 次消元，可把方程组 A 1 x b 1 化为上三角方程组 A n x b n . (2) 回代过程： 以解如下线性方程组为例测试结果 . 2. 列主元消去法 由高斯消去法可知，在消元过程中可能出现 akk k 0 的情况，这是消去法将无法进行， 即使主元素 akk k 0 但很小时，用其作除数， 会导致其他元素数量级的严重增长和舍入误 差的扩散，最后也使得计算解不可靠 . 这时就需要选取主元素，假定线性方程组的系数矩 阵 A 是菲奇异的 . (1) 消元过程： 对于 k 1,2, , n 1 ，进行如下步骤： 1) 按列选主元，记 2) 交换增广阵 A 的 p,k 两行的元素。 A(k,j)=A(p,j) ( j=k, …,n +1) 3) 交换常数项 b 的 p,k 两行的元素。 b(k)=b(p) 4) 计算消元 (2) 回代过程 (3) 以解如下线性方程组为例测试结果 . 三、实验环境 MATLAB R2014a 四、实验步骤 1. 高斯列主元消去法流程图： 开始 输 入 系 数 按列 选主 交 换 元 计 算 消 回代