三角分解营销.pptVIP

3.3 矩阵的三角分解法 我们知道对矩阵进行一次初等变换，就相当于用相应的初等矩阵去左乘原来的矩阵。因此我们这个观点来考察Gauss消元法用矩阵乘法来表示，即可得到求解线性方程组的另一种直接法：矩阵的三角分解。 Gauss消元法的矩阵形式 Doolittle分解 Doolittle分解 若矩阵A有分解：A=LU，其中L为单位下三角阵，U为上三角阵，则称该分解为Doolittle分解，可以证明，当A的各阶顺序主子式均不为零时，Doolittle分解可以实现并且唯一。 A的各阶顺序主子式均不为零，即 Doolittle分解 Doolittle分解 Doolittle分解 Doolittle分解 Doolittle分解 Doolittle分解 例题 例题 例题 例题 例题 Doolittle分解 对称矩阵的Cholesky分解 在应用数学中，线性方程组大多数的系数矩阵为对称正定这一性质，因此利用对称正定矩阵的三角分解式求解对称正定方程组的一种有效方法，且分解过程无需选主元，有良好的数值稳定性。 对称矩阵的Cholesky分解 A对称：AT=A A正定：A的各阶顺序主子式均大于零。即 对称矩阵的Cholesky分解 由Doolittle分解，A有唯一分解 对称矩阵的Cholesky分解 定理3.2.4 设A为对称正定矩阵，则存在唯一分解A=LDLT，其中L为单位下三角阵，D=diag(d1,d2,…,dn)且di0(i=1,…,n) 对称矩阵的Cholesky分解 证明： 对称矩阵的Cholesky分解 对称矩阵的Cholesky分解 对称矩阵的Cholesky分解 推论：设A为对称正定矩阵，则存在唯一分解 其中L为具有主对角元素为正数的下三角矩阵。 证明： Cholesky分解的求法 Cholesky分解的求法 Cholesky分解的求法 Cholesky分解法 Cholesky分解法缺点及优点 优点：可以减少存储单元。 缺点：存在开方运算，可能会出现根号下负数。 改进Cholesky分解法 改进的cholesky分解A=LDLT 改进Cholesky分解 改进Cholesky分解 改进的cholesky分解算法 改进的cholesky分解算法 例题 例题 例题 例题 A=LDLT分解，既适合于解对称正定方程组，也适合求解A为对称，而各阶顺序主子式不为零的方程组 而对A=LLT只适合于对称正定方程组 追赶法 追赶法 追赶法 追赶法 例题 例题 3.5平方根法 在应用数学中，线性方程组大多数的系数矩阵为对称正定这一性质，因此利用对称正定矩阵的三角分解式求解对称正定方程组的一种有效方法，且分解过程无需选主元，有良好的数值稳定性。 平方根法 A对称：AT=A A正定：A的各阶顺序主子式均大于零。即 平方根法 由Doolittle分解，A有唯一分解 平方根法 定理3.2.4 设A为对称正定矩阵，则存在唯一分解A=LDLT，其中L为单位下三角阵，D=diag(d1,d2,…,dn)且di0(i=1,…,n) 平方根法 证明： 平方根法 平方根法 平方根法 推论：设A为对称正定矩阵，则存在唯一分解 其中L为具有主对角元素为正数的下三角矩阵。 证明 平方根法 平方根法 平方根法 平方根法 Cholesky分解法缺点及优点 优点：可以减少存储单元。 缺点：存在开方运算，可能会出现根号下负数。 改进平方根法 改进的平方根法分解A=LDLT 改进平方根法 改进平方根法 改进的平方根法 改进平方根法 例题 例题 例题 例题 A=LDLT分解，既适合于解对称正定方程组，也适合求解A为对称，而各阶顺序主子式不为零的方程组 而对A=LLT只适合于对称正定方程组 2.6范数与误差估计 用直接方法解n阶线性方程组Ax=b，由于原始数据A、b的误差及计算过程中的舍入误差，一般得不到方程的精确解 ，往往得到 它的近似解x，为了讨论解的精度，即误差向量x- ，的大小，也为了讨论用迭代法解线性方程组的收敛性问题，需要引入向量及矩阵的范数。  向量的范数  定义3.1 设是n维向量空间，如果 ，实值函数‖x‖满足条件 (1) 正定性：‖x‖≥0，当且仅当x=0时，‖x‖＝０； (2) 齐次性：λ∈R,‖λx‖=|λ|‖x‖； (3) 三角不等式：‖x+y‖≤‖x‖+‖y‖， 则称‖x‖为 上的向量范数(或向量的模)。 在数值计算中，常用的向量范数有三种。设