(数值分析 第二章线性方程组的直接解法 §1引 很多工程问题的计算最终都归结为解线性方程组.虽 然有些数学模型中不直接表现为解线性方程组但其数 值解法中将问题“离散化”或“线性化”将成为线性方程组 线性方程组的求解是数值分析课程中最基本的、最重 要的内容之 82Gaus消去法 Gaus消元法的基本思想 Gauss消去法的基本思想是对方程组所对应的增广 矩阵进行一系列的初等行变换,将其转化为上三角矩 阵,通过回代求得线性方程组的解 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
数值分析 1111 第二章 线性方程组的直接解法 §1 引 言 很多工程问题的计算最终都归结为解线性方程组. 虽 然有些数学模型中不直接表现为解线性方程组,但其数 值解法中将问题“离散化”或“线性化”将成为线性方程组 . 线性方程组的求解是数值分析课程中最基本的、最重 要的内容之一 . 一 、Gauss消元法的基本思想 §2 Gauss消去法 Gauss消去法的基本思想是对方程组所对应的增广 矩阵进行一系列的初等行变换, 将其转化为上三角矩 阵, 通过回代求得线性方程组的解. PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
(数值分析 二、 Gauss消元法公式 =k+1,…,n.(乘数因子) (k) j=k+1,k+2,…,n (k+1) k) i=k+1,k+2,……,n Gaus消去法的条件 定理 Gauss消去法消元过程能进行到底的充要条件是 系数阵A的1到n-1阶顺序主子式不为零;Ax=b能用 Gauss消元法解的充要条件是A的各阶顺序主子式不 为零 12 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
数值分析 1212 二、 Gauss消元法公式 ( ) ( ) ( 1) ( ) ( ) ( 1) ( ) ( ) , 1, , , , 1, 2, , . , 1, 2, , . k ik ik k kk k k k ij ij ik kj k k k i i ik k a l i k n a a a l a i j k k n b b l b i k k n + + ì = = + ï ï ï í = - = + + ï = - = + + ï ï î L L L .（乘数因子) 三、 Gauss消去法的条件 定理 Gauss消去法消元过程能进行到底的充要条件是 系数阵A的1到n-1阶顺序主子式不为零；Ax=b能用 Gauss消元法解的充要条件是A的各阶顺序主子式不 为零. PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
(数值分析 §3 Gauss主元素法 列主元消元法 算法( Gaussi列主元素法) 1选主元,对k=12,n-1确定i使k=maxl k: 若a1k=0则A为奇异矩阵,停机 2.若i=k,转3.否则换行 b,(≤j≤m) 3消元 a1,<a,-a,a 6.<6 k+1,k+2 13 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
数值分析 1313 §3 Gauss主元素法 一 、 列主元消元法 算法(Gauss列主元素法) max k i k ik k i n a a £ £ 1. 选主元, 对k=1,2,··· ,n-1确定i = k , 使 若aikk=0则A为奇异矩阵, 停机. , ,( ) k k kj i j k i a « a b « b k £ £j n 2. 若ik =k, 转3. 否则换行 ik ik kk a a a ¬ i = k + + 1, k n 2, , L , ij ij ik kj a ¬ - a a a i i ik k b ¬ - b a b 3.消元 PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
(数值分析 4若an=0,停机 5回代求解 ∑a b.< 丿=1+1 6输出(b1,b2…,bn) 二、 Gauss全主元消去法 在Gs消去法中,若每次选主元不局限在方框列 中,而在整个子矩阵 (k) kk (k) (k) PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
数值分析 4.若ann =0, 停机． 1 , , 1, 2, , 3, 2,1 n i ij j n j i n i nn ii b a b b b b i n n a a = + - ¬ ¬ = - - å L 5.回代求解 6.输出(b1 ,b2 , ··· ,bn ) T 二、 Gauss全主元消去法 在Gauss消去法中, 若每次选主元不局限在方框列 中, 而在整个子矩阵 ú ú ú û ù ê ê ê ë é ( ) ( ) ( ) ( ) k nn k nk k kn k kk a a a a L M M L PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
(数值分析 中选取便称为全主元Gas消去法,此时增加的步骤为 ①a|=maxq,确定i,若a1=0给出=0的信息 k 停止计算 ②作如下行交换和列交换 行交换:a)4>aW,b4b6,(≤j≤m) 列交换:ak)4a),(≤i≤n) 需要注意的是,在全主元的消去法中,由于进行了列 交换,x各分量的顺序已被打乱。因此必须在每次列交 换的同时,让机器存储列交换的信息,在回代得出解后再 将x的各分量换回到原来相应的位置处 15N PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
数值分析 中选取,便称为全主元Gauss消去法,此时增加的步骤为: ( ) ( ) , max k k k k i j ij k i j n a a £ £ = ( ) ( ) ( ) ( ) , , ( ) k k k k k k kj i j k i a « a b « b k £ £j n ( ) ( ) , ( ) k k k ik ij a « a k £ £i n 1515 ②作如下行交换和列交换 行交换: 列交换: ① , 确定ik , jk ,若aikjk =0 ,给出|A|=0的信息 停止计算. 需要注意的是, 在全主元的消去法中, 由于进行了列 交换, x各分量的顺序已被打乱。因此必须在每次列交 换的同时,让机器存储列交换的信息,在回代得出解后再 将x的各分量换回到原来相应的位置处. PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com