SVD矩阵的奇异值分降拟.ppt

称上式为矩阵A的奇异值分解. Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 在矩阵理论中，奇异值分解实际上是“对称矩阵正交相似于对角矩阵”的推广. 奇异值分解中 是 的特征向量，而 的列向量是 的特征向量，并且 与 的非零特征值完全相同. 但矩阵 的奇异值分解不惟一. 注意 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 数值秩 在没有误差时，奇异值分解可以确定矩阵的秩. 但是误差的存在使得确定变得非常困难. 例如，考虑矩阵 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 因为第三列是前两列的 和，所以 A 的秩是2. 如果不考虑到这个关系， 运用IEEE标准的双精度浮点计算模式， 用MATLAB命令SVD计算A 的奇异值： format long e A=[1/3,1/3,2/3;2/3,2/3,4/3;1/3,2/3,1;2/5,1/5,3/5;3/7,1/7,4/7]； D= svd(A) Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 计算结果为： D = 2.421457493421318e+000 3.406534035359026e-001 1.875146052457622e-016 因为有“三”个非零奇异值，所以A的秩 为“3”. 然而，注意到在IEEE双精度的标准下, 其中一个奇异值是微小的. 也许应该将它看作 零.因为这个原因，引人数值秩的概念. Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 如果矩阵 有 个“大”的奇异值，而其它都很“微小”，则称 的数值秩为 . 为了确定哪个奇异值是“微小”的，需要引人阈值或容忍度 .就MATLAB而言，可以把 设为阈值，大于这个阈值的奇异值的数目就是 A的数值秩，把小于这个阈值的奇异值看作零. 利用MATLAB的命令rank计算 的秩，它的结果是2，就是这个道理. Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 数无形时少直观, 形无数时难入微, 数形结合百般好, 隔离分家万事休. --------华罗庚 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 将抽象思维形象化 将理论知识实用化 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 二、矩阵的四个