- 7
- 0
- 约6.25千字
- 约 92页
- 2019-02-28 发布于北京
- 举报
相似矩阵,就是同一个线性变换的不同的描述矩阵. 或者说相似矩阵都是同一个线性变换的描述 . 定理 设 非奇异,则存在正交矩阵P和Q, 使得 (2) 其中 证 因为A非奇异,所以 为实对称正定矩阵,于是存 在正交矩阵Q使得, 其中 为 特征值 令 , 引理: 1.设 则 是对称矩阵, 且其特征值是非负实数. 2. 3. 设 则 的充要条件是 谢谢大家, 加强联系与合作! 因为第三列是前两列的 和,所以 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) 计算结果为: D = 2.421457493421318e+000 3.406534035359026e-001 1.875146052457622e-016 因为有“三”个非零奇异值,所以A的秩 为“3”. 然而,注意到在IEEE双精度的标准下, 其中一个奇异值是微小的. 也许应该将它看作 零.因为这个原因,引人数值秩的概念. 如果矩阵 有 个“大”的奇异值,而其它都很“微小”,则称 的数值秩为 . 为了确定哪个奇异值是“微小”的,需要引人阈值或容忍度 .就MATLAB而言,可以把 设为阈值,大于这个阈值的奇异值的数目就是 A的数值秩,把小于这个阈值的奇异值看作零. 利用MATLAB的命令rank计算 的秩,它的结果是2,就是这个道理. 求矩阵 的奇异值分解 解:MATLAB程序为: A=[0,-1.6,0.6;0 ,1.2,0.8;0,0,0;0,0,0] [U,S,V]=svd(A) 计算结果 A = 0 -1.6000 0.6000 0 1.2000 0.8000 0 0 0 0 0 0 U = 0.8000 0.6000 0 0 -0.6000 0.8000 0 0 0 0 1.0000 0 0 0 0 1.0000 S = 2.0000 0 0 0 1.0000 0 0 0 0 0 0 0 V = 0 0 1.0000 -1.0000 0.0000 0 0.0000 1.0000 0 奇异值分解的几何意义 研究将一个空间映射到不同空间,特别是 不同维数的空间时,例如超定或欠定方程组所 表示的情况,就需要用矩阵的奇异值来描述算 子对空间的作用了. 考察二维平面上的单位圆 在映射A下的变换过程,其中 MATLAB程序为: A=[sqrt(3)\sqrt(2),sqrt(3)\sqrt(2);-3\sqrt(2),3\sqrt(2); 1\sqrt(2),1\sqrt(2)] [U,S,V]=svd(A) 记: 基
原创力文档

文档评论(0)