我在时域中有两个波形，我需要在 MATLAB 中测量其中的cross-correlation 系数。我试过max(abs(xcorr(m,n,'coeff')))，但它似乎不能正常工作。
我还需要测量波形不同部分的互相关系数，例如每隔 1 分钟测量一次互相关系数。如果可能的话，将这些值输出到矩阵或其他东西。
我知道有很多问题要问，但我是 MATLAB 新手，发现这项任务令人生畏！
您可以就这个问题的任何部分给予我任何帮助，我们将不胜感激。
编辑： 这是我用来测试关联代码的代码：
x = rand(1,14400);
y = rand(1,14400);
r = max( abs(xcorr(x,y,'coeff')) )
【问题讨论】：
标签： matlab signal-processing correlation
根据this article中的方程可以这样算出互相关系数：
% Assuming: m and n are your signals organized as row vectors
r = cov([m;n]) / (std(m)*std(n));
如果您只想计算部分信号的系数，只需使用：
r = cov([m(1:100);n(1:100)]) / (std(m(1:100))*std(n(1:100)));
您是否也尝试过corrcoef 功能？
编辑 好的，我检查了corrcoef函数，它似乎工作正常，看看：
>> x = 100*randn(1000,1);
>> y=34*randn(1000,1);
>> corrcoef(x,y)
ans =
1.0000 -0.0543
-0.0543 1.0000
因此相关系数等于 -0.0543 - 相似度小（如预期）。
为了检查这一点，让我们计算相同信号的系数：
>> y=x;
>> corrcoef(x,y)
ans =
1 1
1 1
正如预期的那样，它等于 1。
编辑。 如您所见，corrcoef 的结果是这两个信号之间所有可能的相关系数的矩阵：
x y
x 1.0000 -0.0543
y -0.0543 1.0000
因此，对于互相关，您需要选择主对角线之外的元素之一（存在定位的自相关系数，在这种情况下始终等于 1）。
选择 ans(2,1) 或 ans(1,2) 没有区别 - 如果您计算 x 和 y 的相关性，或 y 和 x 的相关性，则没有区别。
所以最终的代码应该是这样的：
R = corrcoef(x,y); % Correlation matrix
r = R(2,1); % this is your Cross-Correlation coefficient
【讨论】：
尝试改用交叉协方差
x = rand(1,14400);
y = rand(1,14400);
r = max( abs(xcov(x,y,'coeff')) )
cross-covariance 序列是均值去除后的互相关
序列。就像 Joonas 提到的那样，rand 的 DC 偏移量为 0.5，会给你一个“不正确”的结果。
【讨论】：
我试过 max(abs(xcorr(m,n,'coeff'))) 但它似乎不能正常工作。
你这是什么意思？它输出什么，你期望什么？
互相关中的一个可能问题是波形中的直流偏置会破坏结果。据我所知，没有通用的方法可以解决这个问题。您必须以某种方式确保您的波形不包含任何直流偏置。
【讨论】：
这是我用来测试相关代码的代码：
x = rand(1,14400); y = rand(1,14400); r = max( abs(xcorr(x,y,'coeff')) )
问题在于rand 返回的数字均匀分布在区间 (0,1) 中。 换句话说，您的 DC 偏差（平均值）为 0.5！这就是为什么您会为看似随机的信号获得高相关系数：它们并不是完全随机的，因为每个信号都有一个相似的常数分量，会显示在相关系数中。
因此，请尝试改用 randn：它返回随机数，其元素正态分布，均值为 0，这正是您想要的。
【讨论】：