数值计算方法线代方程组直接法分析.pptVIP

1 n阶线性代数方程组的一般形式 Gauss消去与矩阵LU分解 一 Gauss消去 1 直接法的关键思想 如果方程组是“上三角方程组”或“下三角方程组”就可以很容易地求出方程组的解。因此，直接法的关键思想就是如何把一般方程组约化为上/下三角方程。 2 上三角方程组与回代过程 3 下三角方程组与前推过程 4 Gauss消去过程 则Gauss消去过程如下： 二 LU分解 1 矩阵的三角分解 对方程组Ax=b的顺序Gauss消去过程的结果，就是把矩阵A分解成两个三角距阵L和U的乘积：A=LU。 利用这个特点可以进行线性方程组的直接三角分解法。 解方程组的直接三角分解法有3种形式： （1）A=LU 2 解方程组的直接三角分解法 function [L,U,flag]=LU_decom(A) [n,m]=size(A); if n~=m error(The rows and columns of matrix A must be equal!); return; end L=eye(n); U=zeros(n); flag=OK; for k=1:n for j=k:n z=0; for q=1:k-1 z=z+L(k,q)*U(q,j); end U(k,j)=A(k,j)-z; end if abs(U(k,k))eps flag=failure; return; end for i=k+1:n z=0; for q=1:k-1 z=z+L(i,q)*U(q,k); end L(i,k)=(A(i,k)-z)/U(k,k); end end 课堂练习 答案 Cholesky分解 function [L,flag]=chol_factor(A) n=length(A); L=zeros(n); flag=OK; for k=1:n delta=A(k,k); for j=1:k-1 delta=delta-L(k,j)^2; end if delta1e-10 flag=failure; return; end L(k,k)=sqrt(delta); for i=k+1:n L(i,k)=A(i,k); for j=1:k-1 L(i,k)=L(i,k)-L(i,j)*L(k,j); end L(i,k)=L(i,k)/L(k,k); end end 向量范数与矩阵范数 一 向量范数 1 定义 设向量x∈Rn，若与x对应的 一个实值函数（并记为）||x||满足： ||x||0，其中||x||=0当且仅当x=0,称为正定性； ||kx||=|k|||x||,k∈R,称为齐次性； ||x+y||≤||x||+||y||,且x,y∈Rn,称为三角不等式。 线性方程组的条件 定义3.2 设A为非奇异矩阵，称 cond(A)=||A-1||||A|| 为矩阵A的条件数。 定理3.6 对任意的非奇异矩阵 A，cond(A)是由任意的矩阵范数定义的条件数，则 for for for 计算 的k列 计算 的k+1行 矩阵 分解的实际计算公式： 该方法也称为平方根法。 思想 ?Cholesky分解的计算公式 设 由 对应元素相等得 Cholesky分解公式 因对称性无需存储 Step1 Step2 Step3 Stepn 的计算过程: 逐 列 计 算 元素 仍然存放在矩阵 的相应位置上 例：用Cholesky分 解法求解下列方程组 解： 系数矩阵为 Step1 Step2 Step3 求解方程组 求解方程组 ? Cholesky分解法求解方程组中需说明的几个问题 ?工作量：约为 分解的一半； ?不必选主元： 的正定性和稳定性 ?稳定性：是数值稳定的； ?缺陷：存在开平方运算。 例：设x=(2,-4,3)T,求||x||1,||x||2和||x||∞ 称为矩阵ATA的谱半径 例：给定矩阵 求矩阵 的1、2、 范数。 若 是实对称矩阵，则 矩阵 的特征值为 。 ，求 例：设 2 1 || || , || || , || || 4 3 2 1 A A A A ￥ ú ? ù ê ? é - - = 对方程组的右端b作一