1、引言
求解线性方程组在许多领域中都有重要应用,写成矩阵的形式:
求解
1、设置增广矩阵,利用高斯消元法,通过初等行列变换可以求
但这种方法不利于使用计算机计算。
2、利用矩阵对角化求
这种方法关键在于求解矩阵
2、经典Jacobi方法
1846年数学家Jacobi提出的经典Jacobi方法用于求解实对称矩阵的特征值。它的核心思想是采用一些列的Jacobi平面旋转矩阵将对称阵
把
Jacobi旋转矩阵定义如下:
一个Jacobi旋转矩阵,对角线上只有第i行i列,j行j列为c,其余为1,未标出的元素均为0。
其中
解得:
当一次变换结束后,
其中
3、复杂度分析
当
但由于每次做旋转变换后会影响矩阵
可能的改进方法
(1)循环Jacobi方法[3]
按照行顺序或者列顺序依次做Jacobi旋转变换,持续多轮,直到满足收敛条件。
![45539b5cff75ddcb5b0547f3f4e4b1e8.png](https://i-blog.csdnimg.cn/blog_migrate/5e367d0dd2a19535ac1089753647025a.png)
(2)过关Jacobi方法[4]
在顺序遍历的基础上,在加入门限值,大于某个门限值,才做旋转变换,其中门限值与
4、非对称矩阵的Jacobi方法
对于非对称矩阵
5、讨论
高效使用Jacobi方法的关键在于,如何使用尽量少的旋转角度完成对角化矩阵,贪心算法是否是最优方法还值得进一步探讨。是否在矩阵中存在一个固定的最优的
参考文献
[1] 郭强. 并行JACOBI方法求解矩阵奇异值的研究[D]. 苏州大学.
[2] C.F. Van Loan G. H. Golub. Matrix COmputations[M]. John Hopkins University Press, Baltimore and London, second edition, 1993.
[3] E. R. Hansen. On Cyclic Jacobi Methods[J]. J.Soc,Indust.Appl.Math, 1963, 11(2): 448–459.
[4] B. N. Parlett. The Symmetric Eigenvalue Problem[M]. Prentice-Hall, 1980.