24小时热门版块排行榜    

查看: 1000  |  回复: 11
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

yanggis

木虫 (著名写手)

小木虫GIS家族族长

[交流] 【求助】求助关于矩阵求逆的诡异结果

我现在有个3*3的矩阵,如下:
20                         37700               71131000
37700                 71131000        1.34333E+11
71131000        1.34333E+11        2.53926E+14
就这个矩阵发生了比较诡异的事情:
用R软件计算这个矩阵的逆,居然算不出结果。用Matlab算出的结果,我自己用Fortran也计算出了结果,我用Excel也计算了一个结果。诡异的是这三个计算结果居然都不一样。而且相差很远。
有兴趣的可以试试。
有谁知道原因吗?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

magic7004

金虫 (职业作家)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
wangen994(金币+1,VIP+0):欢迎参与讨论,呵呵 1-19 10:40
最小的数只有20,最大的数2.53926E+14,如果不选择精度足够足够高的数据类型和算法的话,算出来一定不准的。

估计是FORTRAN、EXCEL和MATLAB用的算法不一样或者数据精度不一样,所以结果不一样。
流氓不可怕,可怕的是流氓有文化,有文化又BH的流氓无敌~~!
11楼2010-01-19 08:30:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

黄胖胖

木虫 (正式写手)


波不动(金币+1,VIP+0):感谢应助! 12-31 20:04
是不是数值的原因啊.
2楼2009-12-31 15:26:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yanggis

木虫 (著名写手)

小木虫GIS家族族长


波不动(金币+1,VIP+0):感谢应助! 12-31 20:04
这个数值应该是没有问题的。这个矩阵是我做时间序列分析时得到的一个矩阵。就矩阵自己而言不存在有没有问题的事情啊。
3楼2009-12-31 17:51:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

黄胖胖

木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
那你可以试着算算简单的整数  看看三种软件计算的结果是否一样?看看是否程序上有没有问题.
4楼2009-12-31 20:49:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见