24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wanggongming

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 186.6
散金: 1
帖子: 171
在线: 21.8小时
虫号: 276422
注册: 2006-09-02
性别: GG
专业: 计算机应用技术

[交流] 【求助】如何求解高精度double型矩阵的秩

求矩阵秩常用的方法是全选主元高斯消去法，我在徐士良的《C常用算法程序集（第二版）》中看到过C语言源代码，但是这个代码无法求解Double型高精度
矩阵（即每个元素的精度都在小数点6位之后）的秩！
有一个一个矩阵，其第一行和第二行元素完全相同（如下所示，由于两行完全相同，只列出第一行），从原理上说秩应该是1。但是求出的秩却是2。
0.288488
-0.010170
-0.007254
0.033046
-0.111335
0.308060
-0.424541
0.317609
-0.159026
0.042474
0.005092
0.000000
.005092
-0.042474
-0.159026
-0.317609
-0.424541
-0.308060
-0.111335
-0.033046
-0.007254
0.010170
0.288488
究其原因，应该是Double型矩阵进行高斯消元时，不好消成某行全为0的形式，所以导致求秩错误！
对于这种情况下的矩阵，不知谁有办法求出它的秩，望高手给予帮助！

[ Last edited by wenzhenzhong on 2008-12-4 at 23:38 ]

回复此楼