python - 从比较信号的角度规范互相关的基础知识
问题描述
我试图了解如何使用互相关来确定两个信号的相似性。本教程对基础知识提供了非常清晰的解释,但我仍然不明白如何有效地使用归一化来防止强信号在具有不同能级的信号时主导互相关测量。同一位导师 David Dorran在这里讨论了归一化的问题,并解释了如何使用点积对相关性进行归一化,但我仍然有一些疑问。
我编写了这个 python 例程来在一组信号中的每对信号之间进行互相关:
import numpy as np
import pandas as pd
def mycorrelate2d(df, normalized=False):
# initialize cross correlation matrix with zeros
ccm = np.zeros(shape=df.shape, dtype=list)
for i, row_dict1 in enumerate(
df.to_dict(orient='records')):
outer_row = list(row_dict1.values())
for j, row_dict2 in enumerate(
df.to_dict(orient='records')):
inner_row = list(row_dict2.values())
x = np.correlate(inner_row, outer_row)
if normalized:
n = np.dot(inner_row, outer_row)
x = x / n
ccm[i][j] = x
return ccm
假设我有 3 个幅度增加的信号:[1, 2, 3], [4, 5, 6] 和 [7, 8, 9]
我想交叉关联这三个信号以查看哪些信号对是相似的,但是当我将这三个信号传递到我编写的例程中时,我似乎没有得到相似度的度量。互相关值的大小只是能量信号的函数。时期。甚至一个信号与其自身的互相关产生的值也低于同一信号与另一个更高能量信号的互相关。
df_x3 = pd.DataFrame(
np.array([[1, 2, 3],
[4, 5, 6],
[7, 8, 9]]).reshape(3, -1))
mycorrelate2d(df_x3)
这产生:
array([[array([ 3, 8, 14, 8, 3]),
array([12, 23, 32, 17, 6]),
array([21, 38, 50, 26, 9])],
[array([ 6, 17, 32, 23, 12]),
array([24, 50, 77, 50, 24]),
array([ 42, 83, 122, 77, 36])],
[array([ 9, 26, 50, 38, 21]),
array([ 36, 77, 122, 83, 42]),
array([ 63, 128, 194, 128, 63])]], dtype=object)
现在,我传入相同的 3 个信号,但这次我表示我想要标准化结果:
mycorrelate2d(df_x3, normalized=True)
这产生:
array([[array([ 0.2142, 0.5714, 1., 0.5714, 0.2142]),
array([ 0.375, 0.71875, 1., 0.5312, 0.1875]),
array([ 0.42, 0.76, 1., 0.52, 0.18])],
[array([ 0.1875, 0.5312, 1., 0.7187, 0.375]),
array([ 0.3116, 0.6493, 1., 0.6493, 0.3116]),
array([ 0.3442, 0.6803, 1., 0.6311, 0.2950])],
[array([ 0.18, 0.52, 1., 0.76, 0.42]),
array([ 0.2950, 0.6311, 1., 0.6803, 0.3442]),
array([ 0.3247, 0.6597, 1., 0.6597, 0.3247])]],
dtype=object)
现在所有的最大值都是 1!!因此,我们从具有基于虚假差异的最大值变为完全没有最大值之间的差异!我欣然承认,我不明白如何使用互相关来检测信号之间的相似性。某人将信号与互相关进行比较的分析工作流程是什么?
解决方案
因此,您用于规范化的公式并不完全正确。标准化发生在我们在 NCC 中关联之前,然后我们将答案除以向量长度,如此维基百科公式https://en.wikipedia.org/wiki/Cross-correlation#Zero-normalized_cross-correlation_(ZNCC)所示
所以你需要类似的东西
import numpy as np
def mycorrelate2d(df, normalized=False):
# initialize cross correlation matrix with zeros
ccm = np.zeros((3,3))
for i in range(3):
outer_row = df[i][:]
for j in range(3):
inner_row = df[j][:]
if(not normalized):
x = np.correlate(inner_row, outer_row)
else:
a = (inner_row-np.mean(inner_row))/(np.std(inner_row)*len(inner_row))
#print(a)
b = (outer_row-np.mean(outer_row))/(np.std(outer_row))
#print(b)
x = np.correlate(a,b)
ccm[i][j] = x
return ccm
df_x3 =np.array([[1, 2, 3],
[4, 5, 6],
[7, 8, 9]]).reshape(3, -1)
print(mycorrelate2d(df_x3,True))
df_x3 =np.array([[1, 2, 3],
[9, 5, 6],
[74, 8, 9]]).reshape(3, -1)
print(mycorrelate2d(df_x3,True))
输出是
[[1. 1. 1.]
[1. 1. 1.]
[1. 1. 1.]]
[[ 1. -0.72057669 -0.85933941]
[-0.72057669 1. 0.97381599]
[-0.85933941 0.97381599 1. ]]
推荐阅读
- google-sheets - 在 Google 表格公式中为 Row() 函数添加前缀
- automated-tests - 试图让 testcafe、saucelabs 和 circleci 一起工作但没有运气
- java - 应用程序因 java.lang.IllegalStateException 而崩溃
- c# - 如何使 Unity 按钮更改 ScrollView 中的某些内容
- scala - Spark Structured Streaming 输出未显示在 inteliJ 控制台中
- python - 变量名称的长度会影响程序运行的速度吗?
- ruby-on-rails - Ruby on Rails 教程:route_path 和 route_url 之间的区别
- ios - 在 UITableViewCell 上显示和使用拖动/重新排序按钮而不处于“编辑”模式?
- assembly - 汇编除以零溢出
- grafana - 计算值在转换为 1 之前设置为 0 的平均时间