• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    公众号

numpy和matlab计算协方差矩阵的不同(matlab是标准的,numpy相当于转置后计算) ...

原作者: [db:作者] 来自: [db:来源] 收藏 邀请

matlab是标准的,numpy相当于转置后计算

>> x = [2,0,-1.4;2.2,0.2,-1.5;2.4,0.1,-1;1.9,0,-1.2]

x =

2.0000    0            -1.4000
2.2000    0.2000    -1.5000
2.4000    0.1000    -1.0000
1.9000    0            -1.2000

>> cov(x)

ans =

0.0492        0.0142      0.0192
0.0142        0.0092      -0.0058
0.0192        -0.0058     0.0492

 

>> xt=x\'

xt =

2.0000     2.2000         2.4000       1.9000
0             0.2000         0.1000       0
-1.4000    -1.5000        -1.0000     -1.2000

>> cov(xt)

ans =

2.9200   3.1600   2.9500   2.6700
3.1600   3.4300   3.1750   2.8850
2.9500   3.1750   3.0100   2.7050
2.6700   2.8850   2.7050    2.4433

 

np.array(x).T
array([[ 2. , 2.2, 2.4, 1.9],
         [ 0. , 0.2, 0.1, 0. ],
         [-1.4, -1.5, -1. , -1.2]])

>>> print(np.cov(np.array(x).T))
[[ 0.04916667 0.01416667 0.01916667]
[ 0.01416667 0.00916667 -0.00583333]
[ 0.01916667 -0.00583333 0.04916667]]

 

np.array(x)
array([[ 2. , 0. , -1.4],
          [ 2.2, 0.2, -1.5],
          [ 2.4, 0.1, -1. ],
          [ 1.9, 0. , -1.2]])

>>> print(np.cov(np.array(x)))
[[ 2.92 3.16 2.95 2.67 ]
[ 3.16 3.43 3.175 2.885 ]
[ 2.95 3.175 3.01 2.705 ]
[ 2.67 2.885 2.705 2.44333333]]


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
Rust: 流程控制发布时间:2022-07-18
下一篇:
【Rust日报】2020-11-06 - Rust Reference中文版发布时间:2022-07-18
热门推荐
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap