RError.com

RError.com Logo RError.com Logo

RError.com Navigation

  • 主页

Mobile menu

Close
  • 主页
  • 系统&网络
    • 热门问题
    • 最新问题
    • 标签
  • Ubuntu
    • 热门问题
    • 最新问题
    • 标签
  • 帮助
主页 / 问题 / 1602408
Accepted
krekninao
krekninao
Asked:2024-12-14 15:30:35 +0000 UTC2024-12-14 15:30:35 +0000 UTC 2024-12-14 15:30:35 +0000 UTC

通过模态频率确定有用信号的开始

  • 772

我有一个信号,例如: 在此输入图像描述

我将信号分成 1 秒的片段。在每个段中,我将信号展开为傅里叶级数并确定模态(即具有最大幅度)频率。我得到了一组分段和这些分段上的模态频率。例如,像这样: 在此输入图像描述

可以看出,当有用信号到达时(火车、电动火车等),5-6 Hz 的低频占主导地位。其余时间,高频噪声占主导地位。我需要自动确定有用信号的开始和结束。我该怎么做?有一个滑动窗口的想法,但我无法理解具有有用信号的部分的开始和结束的具体标准是什么。请告诉我可以完成此操作的算法。

фурье
  • 2 2 个回答
  • 25 Views

2 个回答

  • Voted
  1. MBo
    2024-12-14T16:23:40Z2024-12-14T16:23:40Z

    如果您的第一个图表与真实信号相关,那么您应该使用信号功率(信号^2在间隔内的积分或总和),并使用频率作为可靠性的证据。

    如果您依赖于频率,请构建一个频谱图 - 一组连续(或重叠)间隔的傅里叶变换幅度,并估计所需频率间隔(例如 5-15 Hz)中每个频谱的功率与总功率的比率功率(零,并且可能一定要从计算中删除第一个计数)。如果该比率高于阈值(取决于您的具体数据),则低频占主导地位。

    请注意,时尚不太适合您的任务。

    • 0
  2. Best Answer
    CrazyElf
    2024-12-14T16:43:32Z2024-12-14T16:43:32Z

    嗯,这是给您的一个例子。需要选择窗口大小和截止值。窗口越大,越能平滑异常值,但检测延迟也越大。那么,截止值的选择也需要大致相同的考虑因素。

    import pandas as pd
    import numpy as np
    import matplotlib.pylab as plt
    
    data = np.array([50]*40+[5]*15+[50]*30+[5]*10+[50]*40)
    data += np.random.randint(2, size=(len(data)))
    data[45] = 30
    data[50] = 40
    data[90] = 25
    
    series = pd.Series(data)
    roll = series.rolling(3).mean()
    
    found = roll < 20
    result = (found != found.shift(-1))[:-1]
    startstop = result[result].index
    
    plt.figure()
    plt.bar(range(len(data)), data)
    plt.scatter(startstop, data[startstop], color='red', marker='*')
    plt.savefig('data.png')
    
    plt.figure()
    roll.plot()
    plt.scatter(startstop, roll[startstop], color='red', marker='*')
    plt.savefig('dataroll.png')
    
    for a, b in zip(startstop[::2], startstop[1::2]):
        print(f'{a}:{b}')
    
    # 41:54
    # 86:94
    

    在此输入图像描述 在此输入图像描述

    • 0

相关问题

Sidebar

Stats

  • 问题 10021
  • Answers 30001
  • 最佳答案 8000
  • 用户 6900
  • 常问
  • 回答
  • Marko Smith

    我看不懂措辞

    • 1 个回答
  • Marko Smith

    请求的模块“del”不提供名为“default”的导出

    • 3 个回答
  • Marko Smith

    "!+tab" 在 HTML 的 vs 代码中不起作用

    • 5 个回答
  • Marko Smith

    我正在尝试解决“猜词”的问题。Python

    • 2 个回答
  • Marko Smith

    可以使用哪些命令将当前指针移动到指定的提交而不更改工作目录中的文件?

    • 1 个回答
  • Marko Smith

    Python解析野莓

    • 1 个回答
  • Marko Smith

    问题:“警告:检查最新版本的 pip 时出错。”

    • 2 个回答
  • Marko Smith

    帮助编写一个用值填充变量的循环。解决这个问题

    • 2 个回答
  • Marko Smith

    尽管依赖数组为空,但在渲染上调用了 2 次 useEffect

    • 2 个回答
  • Marko Smith

    数据不通过 Telegram.WebApp.sendData 发送

    • 1 个回答
  • Martin Hope
    Alexandr_TT 2020年新年大赛! 2020-12-20 18:20:21 +0000 UTC
  • Martin Hope
    Alexandr_TT 圣诞树动画 2020-12-23 00:38:08 +0000 UTC
  • Martin Hope
    Air 究竟是什么标识了网站访问者? 2020-11-03 15:49:20 +0000 UTC
  • Martin Hope
    Qwertiy 号码显示 9223372036854775807 2020-07-11 18:16:49 +0000 UTC
  • Martin Hope
    user216109 如何为黑客设下陷阱,或充分击退攻击? 2020-05-10 02:22:52 +0000 UTC
  • Martin Hope
    Qwertiy 并变成3个无穷大 2020-11-06 07:15:57 +0000 UTC
  • Martin Hope
    koks_rs 什么是样板代码? 2020-10-27 15:43:19 +0000 UTC
  • Martin Hope
    Sirop4ik 向 git 提交发布的正确方法是什么? 2020-10-05 00:02:00 +0000 UTC
  • Martin Hope
    faoxis 为什么在这么多示例中函数都称为 foo? 2020-08-15 04:42:49 +0000 UTC
  • Martin Hope
    Pavel Mayorov 如何从事件或回调函数中返回值?或者至少等他们完成。 2020-08-11 16:49:28 +0000 UTC

热门标签

javascript python java php c# c++ html android jquery mysql

Explore

  • 主页
  • 问题
    • 热门问题
    • 最新问题
  • 标签
  • 帮助

Footer

RError.com

关于我们

  • 关于我们
  • 联系我们

Legal Stuff

  • Privacy Policy

帮助

© 2023 RError.com All Rights Reserve   沪ICP备12040472号-5