第二章线性方程组的数值解法.pptVIP

数值计算方法 第二章 线性方程组的数值解法 第二章 线性方程组的数值解法 §2.0 引 言 §2.1 Gauss消去法 §2.2 矩阵的三角分解 §2.3 QR分解和奇异值分解 §2.0 引 言 在自然科学和工程技术中很多问题的解决常常归结为解线性代数方程组。 例如：电学中的网络问题 用最小二乘法求实验数据的曲线拟合问题 解非线性方程组问题 用差分法或者有限元方法解常微分方程 偏微分方程边值问题等 都导致求解线性代数方程组。 §2.0 引 言 这些方程组的系数矩阵大致分为两种 一种是低阶稠密矩阵(例如，阶数大约为≤150) 另一种是大型稀疏矩阵(即矩阵阶数高且零元素较多) §2.0 引 言 设有线性方程组Ax = b，其中 为非奇异阵， ， 关于线性方程组的数值解法一般有两类： 直接法与迭代法。 §2.0 引 言 1. 直接法 就是经过有限步算术运算，可求得方程组精确解的方法(若计算过程中没有舍入误差)。 但实际计算中由于舍入误差的存在和影响，这种方法也只能求得线性方程组的近似解。 本章将阐述这类算法中最基本的高斯消去法及其某些变形。 这类方法是解低阶稠密矩阵方程组的有效方法，近十几年来直接法在求解具有较大型稀疏矩阵方程组方面取得了较大进展。 §2.0 引 言 2. 迭代法 就是用某种极限过程去逐步逼近线性方程组精确解的方法。 迭代法具有需要计算机的存贮单元较少、程序设计简单、原始系数矩阵在计算过程中始终不变等优点，但存在收敛性及收敛速度问题。 迭代法是解大型稀疏矩阵方程组(尤其是由微分方程离散后得到的大型方程组)的重要方法。 第6章介绍迭代法解线性方程组。 §2.1 Gauss消去法 高斯(Gauss)消去法是解线性方程组最常用的方法之一 它的基本思想是通过逐步消元（行的初等变换），把方程组化为系数矩阵为三角形矩阵的同解方程组，然后用回代法解此三角形方程组（简单形式）得原方程组的解。 例如： §2.1 Gauss消去法 下面讨论一般的解n阶方程组的高斯消去法。 1. 消去过程 将原方程组记为 A(1)x =b(1) 其中A(1)=(aij(1))n?n=(aij)n?n ，b(1)=b (1) 第一次消元。 §2.1 Gauss消去法 1. 消去过程 (1) 第一次消元。 其中 §2.1 Gauss消去法 (2) 第k次消元。 §2.1 Gauss消去法 1. 消去过程 (2) 第k次消元。 注：为减少计算量，令 ， 则 §2.1 Gauss消去法 1. 消去过程 (3) 当k = n – 1时得 完成第n – 1次消元后得到与原方程组等价的三角形方程组 A(n)x = b(n) 注：当det(A) ≠ 0时，显然有aii(i) ≠ 0，(i = 1，…，n) 称aii(i)为主元素。 §2.1 Gauss消去法 2. 回代过程 求解三角形方程组A(n)x = b(n)， 得到求解公式 注：求解过程称为回代过程。 §2.1 Gauss消去法 3. 算法设计 令 i = k+1,k+2,…,n §2.1 Gauss消去法 3. 算法设计 function x = gauss(a) m=size(a);n=m(1); for k=1:n-1 for i=k+1:n a(i,:)=a(i,:)-a(k,:)*a(i,k)/a(k,k); end end for j=n:-1:1 a(j,n+1)=(a(j,n+1)-a(j,j+1:n)*a(j+1:n,n+1))/a(j,j) end x=a(1:n,n+1); §2.1 Gauss消去法 3. 算法设计 i = k+1,k+2,…,n §2.1 Gauss消去法 4. Gauss消去法的计算量 以乘除法的次数为主 (1) 消元过程： 第k步时(n – k) + (n – k) (n – k + 1) = (n – k) (n – k + 2) 共有 §2.1 Gauss消去法 4. Gauss消去法的计算量 (1) 消元过程： (2) 回代过程： 求xi中, 乘n–i次,除1次,共n–i+1次(i = 1,…,n–1) 共有 §2.1 Gauss消去法 4. Gauss消去法的计算量 (1) 消元过程： (2) 回代过程： (3) 总次数为 注：当n = 20时约为2670次，比克莱姆法则9.7?