1-2 线性方程组求解.pptVIP

方程组求解 高斯消元法及算法实现 初等变分原理 最速下降法 共轭梯度法 参考文献 [1]李庆扬 关治 白峰杉, 数值计算原理(清华) [2]蔡大用 白峰杉, 现代科学计算 [3] 李庆扬 等, 数值分析 [4]Numerical Analysis (Seventh Edition) 数值分析 （第七版 影印版） [5]David Kincaid,数值分析(第三版) [6] John H. Mathews,数值方法(MATLAB版) 线性方程组的矩阵形式 a11x1+ a12x2+····+ a1nxn = b1 a21x1+ a22x2+····+ a2nxn = b2 ········································· an1x1+ an2x2+····+ annxn = bn A X = b ( i=1,2,···,n ) 线性方程组求解: 1. 直接方法; 2. 基本迭代法; 3. 子空间方法 X ? ? b 解线性方程组的克莱姆方法 1. 输入矩阵 A 和右端向量 b; 高斯消元法 第一步: 将方程组化简为三角形方程组; 第二步: 解三角形方程组,获方程组的解。 4. 计算并输出 x1 = D1 / D,····, xn=Dn/D, 结束。 3. 对 k=1,2,···,n 用 b 替换 A 的第 k 列数据,并计算替换后矩阵的行列式值 Dk; 2. 计算 A 的行列式 D,如果 D=0,则输出错信息结束,否则进行 3 ; 增广矩阵 计算: [m21 m31 m41]T = [a21 a31 a41]T / a11 用–m21乘矩阵第一行加到矩阵第二行; 用–m31乘矩阵第一行加到矩阵第三行; 用–m41乘矩阵第一行加到矩阵第四行; 实现第一轮消元 计算: [m32 m42]T = 用–m32乘矩阵第二行加到矩阵第三行; 用–m42乘矩阵第二行加到矩阵第四行; 实现第二轮消元、第三轮消元········· n阶方程组消元过程乘法次数: (n-1)n+(n-2)(n-1)+…+1×2=(n3-n)/3 除法次数: (n-1)+(n-2)+…+1=n(n-1)/2 上三角方程组 解上三角方程组 计算：xn = bn /ann (a11…ann≠0) xk=[bk－(ak , k+1xk+1+ … + ak n)] / ak k ( k =n－1,···,1 ) 除法: n次; 乘法: n(n-1)/2次, 乘、除法运算共 n(n+1)/2 次, 简记为 O( n2 ) n 2 3 4 5 6 高斯 6 17 36 65 106 克莱姆 8 51 364 2885 25206 定义1 设 Rn是n维向量空间,如果对任意x∈Rn,都有一个实数与之对应,且满足如下三个条件: (1)正定性: ||x||≥0,且||x||=0 = x = 0 ; 则称||x||为向量x的范数 . 一、向量的范数 例3. 范数意义下的单位向量: X=[x1, x2]T 二、矩阵的范数 例 5 Frobenius范数 极小化方法 一、与线性方程组等价的变分问题 三、共轭梯度法(共轭斜量法) 四、预条件共轭梯度法 二、最速下降法 设x, y∈R n, 记 ( x , y) = xT y ( x, y ) = ( y, x ); ( tx, y ) = t ( x, y); ( x+ y, z ) = ( x, z ) + ( y, z ); ( x, x) ≥ 0, 且( x, x) = 0 ? x = 0; 设A是n阶对称正定阵 ( Ax, y ) = ( x, Ay ) ; ( Ax，x ) ≥0, 且( Ax, x) = 0 ? x = 0 一、与线性方程组等价的变分问题 定理1 设A =( aij )n×n为实对称正定矩阵, b , x∈R n, 则 x使二次函数 取极小值 ? x 是线性方程组Ax = b 的解。 二、最速下降法 三、共轭梯度法 (CG) (共轭斜量法) 四、预条件共轭梯度法(PCG) 预条件共轭梯度法 实际计算，可通过变换，转化成用原方程组的量来计算。 预条件共轭梯度法MATLAB的三种调用格式: 1.不用预优矩阵的共轭梯度法 x=pcg(a,b,tol,kmax) 2.用预优矩阵的共轭梯度法 (1)x=pcg(a,b,tol,kmax,m) (2) r=chol(m) x=pcg(a,b,tol,kmax,r’,r,x0) 3.未给定预优矩阵的共轭梯度法 r=cholinc(sa,