协方差、协方差矩阵的数学概念及算法计算

在讲解协方差之前,我们先一起回忆一下样本的均值、方差、标准差的定义。

方差,协方差和协方差矩阵

1、概念

方差(Variance)是度量一组数据的分散程度。方差是各个样本与样本均值的差的平方和的均值: 


协方差(Covariance)是度量两个变量的变动的同步程度,也就是度量两个变量线性相关性程度。如果两个变量的协方差为0,则统计学上认为二者线性无关。注意两个无关的变量并非完全独立,只是没有线性相关性而已。计算公式如下: 

如果协方差大于0表示一个变量增大是另一个变量也会增大,即正相关,协方差小于0表示一个变量增大是另一个变量会减小,即负相关。 
协方差矩阵(Covariance matrix)由数据集中两两变量的协方差组成。矩阵的第(i,j)(i,j)个元素是数据集中第ii和第jj个元素的协方差。例如,三维数据的协方差矩阵如下所示: 

2、练习

计算下表数据的协方差矩阵:

 

Python代码如下:

可以由python中的numpy包计算均值和协方差:

import numpy as np
X = [[2, 0, -1.4],
    [2.2, 0.2, -1.5],
    [2.4, 0.1, -1],
    [1.9, 0, -1.2]]
print(np.mean(X,axis=0))
print(np.cov(np.array(X).T))

计算结果如下:

  • 15
    点赞
  • 57
    收藏
    觉得还不错? 一键收藏
  • 3
    评论
MUSIC (MUltiple SIgnal Classification) 算法是一种用于估计信号源位置的高分辨率频谱分析算法。MUSIC 算法的核心是计算信号源在空间域中的协方差矩阵,然后对其进行特征分解。在这个过程中,协方差矩阵计算是非常关键的。 假设有 $M$ 个传感器,每个传感器接收到的信号为 $x_i(n)$,其中 $n$ 表示时间步,$i$ 表示传感器编号。则协方差矩阵 $R$ 的定义为: $$ R = E\{x(n)x^H(n)\} $$ 其中 $E\{\cdot\}$ 表示取期望,$^H$ 表示共轭转置。在实际应用中,我们通常使用样本协方差矩阵 $\hat{R}$ 来替代真实协方差矩阵 $R$,其定义为: $$ \hat{R} = \frac{1}{N}\sum_{n=1}^{N}x(n)x^H(n) $$ 其中 $N$ 表示采样点数。由于样本协方差矩阵的估计误差会影响 MUSIC 算法的性能,因此需要采用一些方法来提高协方差矩阵的估计精度。 计算协方差矩阵的过程可以分成两步。首先,我们需要将接收到的信号进行预处理,以消除信号中的噪声和干扰。常见的预处理方法包括:空间平滑、信号滤波、时域平滑等。然后,我们可以使用预处理后的信号计算协方差矩阵。常见的计算方法包括:样本协方差矩阵法、修正样本协方差矩阵法等。 在 MUSIC 算法中,我们需要对协方差矩阵进行特征分解。具体来说,我们需要对 $\hat{R}$ 进行特征值分解,得到矩阵 $E$ 和向量 $\Lambda$,使得 $\hat{R} = E\Lambda E^H$。然后,我们将 $E$ 的前 $M$ 列作为信号子空间,$E$ 的后 $M-N$ 列作为噪声子空间。最后,我们可以使用信号子空间来估计信号源的位置。 总之,计算协方差矩阵是 MUSIC 算法的核心步骤之一。通过对接收到的信号进行预处理,并采用合适的协方差矩阵估计方法,可以提高 MUSIC 算法的性能。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值