| 首先，矩阵开方如何定义：M=A*A，还是M=A*AT（AT表示A的转秩）？？ 其次，当M=A*A时，已知M，怎么求A？？ 我在运用一个算法的时候发现这样一个问题，求矩阵的开方，利用Malab试一些方法，归集如下： 1）利用求未知数的方法，解N元二次方程组，理论上是可解得，但是我昨天用Matlab试过了，结果没有出来，电脑差不多快崩溃了； 2）利用Matlab中的函数Sqrtm（），这个解法比较方便，但是其算法究竟是什么，不得而知； 3）利用奇异值开方求矩阵的开放，这边要利用奇异值的性质，同样有三种方法求奇异值，经Matlab仿真表明均可行： a.利用eig（）函数，但是这样的函数必须是实矩阵，而且是正定的，当然eig求出来的是特征值，我认为特征值是奇异值的一种特殊情况，所以条件比较苛刻； b.利用Svd（）函数，这是标准的求奇异值的函数，可以在我以前的博文中看到定义； c.用乔来斯基分解Chol（）函数，先将原矩阵乔来斯基分解，然后再求上对角阵的奇异值，最后同上面几种解法，这个方法在老外的一本矩阵分析中出现，但是有人持怀疑态度，我也没有看到原文，但是我用Matlab试过，是可行的。 可能还有其他的解法，我个人认为，只有想办法把矩阵的奇异值求出来就解决了问题，因而问题的焦点是究竟由多少种求解的方法？ 这个问题可能真是个小问题，但是到现在我还没有找到答案。 故，特来求助，望哪位大虾能给予理论上的指导！！ |
| 欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) | Powered by Discuz! X2.5 |