本菜鸡绝不放弃 2020-04-10 10:14 采纳率: 0%
浏览 665
已结题

时间序列做fft变换后计算出的欧氏距离为什么和用原数据计算出的欧氏距离不一样?

代码如下:

import numpy as np
from scipy.fftpack import fft
from scipy.spatial.distance import pdist

# 序列1
qy1 = [1,2,3,4,5,6,7,8]
# 序列2
qy2 = [2,3,4,5,6,7,8,9]

# 对序列1做傅里叶变换
fqy1 = fft(qy1)
qy1_5 = fqy1[:5]  # 取前五个系数
# 对序列2做傅里叶变换
fqy2 = fft(qy2)
qy2_5 = fqy2[:5]  # 取前五个系数

# 计算距离
X1 = np.vstack([qy1_5, qy2_5])  # 计算前五个傅里叶系数的距离
fd = pdist(X1)
print("傅里叶变换距离:", fd)
X2 = np.vstack([qy1, qy2])  # 计算两个序列的真实欧氏距离
td = pdist(X2)
print("真实距离:", td)

# 运行结果
# 傅里叶变换距离: [8.]
# 真实距离: [2.82842712]

有没有大佬告诉我哪里弄错了?
因为93年《Efficient Similarity Search In Sequence Databases 》这篇论文里提到“Another important observation is Parseval's theorem, which specifies that the Fourier transform preserves the Euclidean distance in the time or frequency domain.”应该是时域距离和频域距离是相等的(我的理解是这样的)。
求大佬指点!!!救命!!!

  • 写回答

1条回答 默认 最新

  • 关注
    评论

报告相同问题?

悬赏问题

  • ¥20 测距传感器数据手册i2c
  • ¥15 RPA正常跑,cmd输入cookies跑不出来
  • ¥15 求帮我调试一下freefem代码
  • ¥15 matlab代码解决,怎么运行
  • ¥15 R语言Rstudio突然无法启动
  • ¥15 关于#matlab#的问题:提取2个图像的变量作为另外一个图像像元的移动量,计算新的位置创建新的图像并提取第二个图像的变量到新的图像
  • ¥15 改算法,照着压缩包里边,参考其他代码封装的格式 写到main函数里
  • ¥15 用windows做服务的同志有吗
  • ¥60 求一个简单的网页(标签-安全|关键词-上传)
  • ¥35 lstm时间序列共享单车预测,loss值优化,参数优化算法