R语言出现百思不得其解事件

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我在R里面得到了一个正定核矩阵K0，Q是一个对角矩阵，K=Q%*%K0%*%Q，奇怪的是，K竟然不对称，而且非正定。。。
为什么啊？？？首先K肯定对称啊，而且K肯定是正定的啊！！！求解！！！！！！！！！！！1

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：百思不得其解 R语言 不对称 非正定 毕业论文 R语言

黑呢 发表于 2016-4-7 13:44
我在R里面得到了一个正定核矩阵K0，Q是一个对角矩阵，K=Q%*%K0%*%Q，奇怪的是，K竟然不对称，而且非正定。。 ...

是不是有一个应该转置啊

对角矩阵的转至是本身

winnie03125 发表于 2016-4-7 13:51
是不是有一个应该转置啊

对角矩阵，不需要的

可以把K展示出来么？

另外，考虑矩阵运算，四舍五入的小数情况

jiangbeilu 发表于 2016-4-7 14:57
可以把K展示出来么？

另外，考虑矩阵运算，四舍五入的小数情况

维数有点大，不过可以肯定的是K0是正定的。
而且我觉得真的是矩阵运算的问题。我试了下K[1,2]=1.2e-220,K[2,1]=1.2e-220
但是K【2,1】-k[1,2]并不等于零（在r里面）
怎么解决r的矩阵运算问题啊。。。理论上K的特征值应该都是大于零的，但r做出来有负值！

jiangbeilu 发表于 2016-4-7 14:57
可以把K展示出来么？

另外，考虑矩阵运算，四舍五入的小数情况

理论上真没啥问题，是不是一开始要设置digits=2

黑呢 发表于 2016-4-7 17:25
理论上真没啥问题，是不是一开始要设置digits=2

R默认的是7位小数

jiangbeilu 发表于 2016-4-7 18:17
R默认的是7位小数

K仍然从出现了负的特征值。。。

K[1,2]=1.2e-220,K[2,1]=1.2e-220，这两个数早就可以约等于0了，小数点后面219个0，早就超出了R的精确范围了。
如果是11.02，11.03这样的数，相差0.01，你可以认为两个数不相等。