Python在振动信号处理中的应用（十）：三分之一倍频程谱的计算

最新推荐文章于 2024-04-27 17:30:34 发布
最新推荐文章于 2024-04-27 17:30:34 发布
阅读量5k 收藏 3
点赞数
分类专栏： Python在振动信号处理中的应用 文章标签： python 开发语言 后端
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/pengchen_wuhan/article/details/121346380
版权

文章目录

一、概述

  三分之一倍频程谱是一种频率分析方法，它具有谱线少频带宽的特点，常用于声学、人体振动、机械振动等测试分析以及频带范围较宽的随机振动测试分析等。

二、算法原理

  倍频程实际上是频域分析中频率的一种相对尺度。倍频程谱是由一系列频率点以及对应这些频率点附近的频带内信号的平均幅值（有效值）所构成。这些频率点称为中心频率 f c f_c f

了解本专栏 订阅专栏 解锁全文
白银时代_
关注
  • 0
    点赞
  • 踩
  • 3
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 10
    评论
专栏目录
订阅专栏
三分之一倍频程程序解读[借鉴].pdf
10-13
三分之一倍频程程序解读[借鉴].pdf
matlab实现三分之一倍频程,三分之一倍频程程序
weixin_30899023的博客
03-18 3879
方法一：%A计权声压级频谱分析clc;clear;close all;y=wavread('abc.wav');fs=51200;%采样频率p0=2e-5;%参考声压f=[1.00 1.25 1.600 2.00 2.50 3.15 4.00 5.00 6.30 8.0]; %基准中心频率f1=[20.00 25.0 31.5 40.0 50.0 63.0 80];fc=[f1,100*f,100...
基于MATLAB的频谱、能量谱、三分之一倍频程分析
weixin_41406486的博客
03-16 1万+
以信号为例，信号在时域下的图形可以显示信号如何随着时间变化，而信号在频域下的图形（一般称为频谱）可以显示信号分布在哪些频率及其比例。频域的表示法除了有各个频率下的大小外，也会有各个频率的相位，利用大小及相位的资讯可以将各频率的弦波给予不同的大小及相位，相加以后可以还原成原始的信号。时域：描述数学函数或物理信号对时间的关系；例如一个信号的时域波形可以表达信号随着时间的变化[1]。这符合现实世界中人们对信号的认识。频域：是指在对函数或信号进行分析时，分析其和频率有关部分，而不是和时间有关的部分，和时域。
matlab 三分之一倍频程,三分之一倍频程谱
weixin_30249953的博客
03-18 2507
三分之一倍频程谱是一种频率分析方法，它具有谱线少频带宽的特点。倍频程实际上是频域分析中频率的一种相对尺度。倍频程谱是由一系列频率点以及对应这些频率点附近的频带内信号的平均幅值(有效值)所构成。这些频率点称为中心频率fc，中心频率附近的频带出于下限频率fl和上限频率fu之间。三分之一倍频程谱是按逐级式频率进行分析的，它是由多个带通滤波器并联组成，为的是使这些带通滤波器的带宽覆盖整个分析频带。根据国籍...
频谱、能谱、功率谱、倍频程谱、1/3 倍频程谱
热门推荐
Ivan 的专栏
01-13 5万+
在做声学信号处理时经常会遇到下面几个概念：能谱、功率谱、倍频程谱、1/3 倍频程谱。这些概念有区别也有联系。大家谈论问题时经常将其中一些概念混用。最近有点时间，我将这几种术语做一个梳理，记录在这里。其中，准确的说明功率谱密度的概念需要用到一点平稳随机过程的知识，考虑到多数人对随机过程不太了解，我这里尽量只用最简单的傅立叶分析的基本概念来说明，这样虽然不太严谨，但是对于我们平常的应用来说也够用了。
一个MATLAB关于三分之一倍频程处理的程序.rar_matlab 左1/3倍频_matlab1/3倍频程_proudscj_三
07-14
使用matlab计算噪声处理：1/3倍频程
python实现1/n倍频程计算
jin666669的博客
05-15 4887
倍频程定义 在音频分析领域，经常要分析音频信号的频谱，这时最常用的是倍频程功率谱和1/3倍频程功率谱。 所谓倍频程功率谱，是将音频分为一个个的频段，然后分别计算每个频段内的功率谱。相邻频段的宽度为二比一的关系。 1/3倍频程是将倍频程再细分为三段。 1/n倍频程是将倍频程再细分为n段。 [倍频程介绍参考链接](https://blog.csdn.net/liyuanbhu/article/deta...
2倍频程、1倍频程、1/3倍频程介绍
学海无涯的专栏
12-22 4万+
什么是倍频程？有何用？ 倍频程，又称倍波程，指在滤波特性曲线上，频率或波长之比为2或1/2的两个频率或波长之间的间隔。常用于音响设备的音频信号分析。 声音信号的频率范围为20Hz~20kHz。分析音频设备的信号频率特性时，一般不可能分析到每一个频点，因此，将20Hz~20kHz频率分为若干个频段，每个频段的上限频率是下限频率的两倍。 音频信号的倍频程的中心频率为： 31.25、62.5、125、250、500、1K、2K、4K、8K、16KHz十个频率 声音信号的频谱，如图所示： 可知整个频谱的
三分之一倍频程谱
尘埃飞舞
10-26 1万+
三分之一倍频程谱是一种频率分析方法，它具有谱线少频带宽的特点。 倍频程实际上是频域分析中频率的一种相对尺度。倍频程谱是由一系列频率点以及对应这些频率点附近的频带内信号的平均幅值（有效值）所构成。这些频率点称为中心频率fc，中心频率附近的频带出于下限频率fl和上限频率fu之间。 三分之一倍频程谱是按逐级式频率进行分析的，它是由多个带通滤波器并联组成，为的是使这些带通滤波器的带宽覆盖整个分析频带。...
octave 倍频程、1/3倍频程分析
12-17
非常详细的octave 倍频程、1/3倍频程分析方法程序代码
wav音频文件三分之一倍频程.zip
11-22
wav文件三分之一倍频程分析 并保存结果到csv
pyfilterbank:为python实现分数倍频程滤波器组。 基于Numpy和CFFI
05-10
分数倍频程滤波器组（适用于IEC-61260：1995） 频谱加权滤波器（IEC 61672：2003） A加权 B和C加权 梅尔频率滤波器组（三角滤波器组） 作为转换矩阵 计划中的：使用STFT 二阶截面/双二阶滤波器和滤波器设计 ...
pyRiemann:Python软件包，用于协方差矩阵处理和生物信号分类，并在脑计算机接口中得到应用
04-29
这项工作仍在进行中，敬请期待。 此代码是BSD许可的（3子句）。 文献资料 该文档位于 安装 使用PyPI pip install pyriemann 或使用pip + git作为最新版本的代码： pip install git+...
[python3] 读取一个正在更新的日志文件
言之。
04-24 1149
等同于–follow=name --retry，根据文件名进行追踪，并保持重试，即该文件被删除或改名后，如果再次创建相同的文件名，会继续追踪。这两种方法都会持续监听日志文件的变化，并实时读取新增的日志内容。你可以根据实际需求选择其中一种方法。要读取一个正在更新的日志文件（即实时写入的日志文件），你可以使用 Python 的。的技巧来实现实时读取。函数打开文件，并使用。
Day25-Java基础之常用类1
最新发布
m0_46053885的博客
04-27 889
同时我们发现Math类中的方法都是静态的，因此在使用的时候我们可以直接通过类名去调用。对于计算机而言，其实是没有数据类型的概念的，都是0101010101，数据类型是编程语言自己规定的，所以在实际存储的时候，先把具体的数字变成二进制，每32个bit为一组，存储在数组中。比较内存地址值一般情况下是没有意义的，我们希望比较的是对象的属性，如果两个对象的属性相同，我们认为就是同一个对象；如果我们的数据是一个浮点类型的数据，有的时候计算机并不会将这个数据完全转换成一个二进制数据，而是将这个将其转换成一个无限的。
实现自定义注解、实现自定义幂等性注解
qq_44721738的博客
04-23 1172
添加 Spring AOP 依赖。创建自定义注解。创建一个新的 Java 注解类，通过@interface关键字来定义，并可以添加元注解以及属性。@Target(ElementType.METHOD) //表示作用于方法上@Retention(RetentionPolicy.RUNTIME) // 表示这个注解在运行时是可见的，这样 AOP 代理才能在运行时读取到这个注解编写 AOP 拦截（自定义注解）的逻辑代码。@Aspect@Component// 方法执行前的处理。
Pandas 2.2 中文官方教程和指南（二十二）
龙哥盟
04-24 1165
原文：pandas.pydata.org/docs/ 时间增量 原文：pandas.pydata.org/docs/user_guide/timedeltas.html 时间增量是时间之间的差异，以不同的单位表示，例如天、小时、分钟、秒。它们可以是正数也可以是负数。 Timedelta是datetime.timedelta的子类，并且行为类似，但也允许与np.timedelta64类型兼容，以及一系列自定义表示、解析和属性。 解析 您可以通过各种参数构造一个Timedelta标量，包括ISO 860
Pandas 2.2 中文官方教程和指南（十五）
龙哥盟
04-24 768
原文：pandas.pydata.org/docs/ 处理文本数据 原文：pandas.pydata.org/docs/user_guide/text.html 文本数据类型 在 pandas 中有两种存储文本数据的方式： object -dtype NumPy 数组。 StringDtype 扩展类型。 我们建议使用StringDtype来存储文本数据。 在 pandas 1.0 之前，object dtype 是唯一的选项。这在很多方面都是不幸的： 你可能会在object dty
三分之一倍频程 python
07-02
### 回答1： 三分之一倍频程是指频率范围的一种度量方式，它表示的是频率范围中的一个值，位于频率范围的三分之一处。 Python是一种编程语言，它具有广泛的应用领域，包括数据分析、人工智能、机器学习等。在编写Python代码时，有时需要处理频率范围的计算或处理，包括计算三分之一倍频程。 为了计算三分之一倍频程，我们首先需要确定频率范围的两个边界值。假设这两个值分别为f1和f2。接下来，我们可以使用以下的公式计算三分之一倍频程： 三分之一倍频程 = f1 + (f2 - f1) / 3 这个公式的含义是，将频率范围平均分成三个部分，然后取第一个部分的末尾作为三分之一倍频程的值。 在Python中，可以使用变量和数学操作符来计算三分之一倍频程。下面是一个示例代码： ```python f1 = 100 # 频率范围的起始值 f2 = 300 # 频率范围的结束值 third_octave_bandwidth = f1 + (f2 - f1) / 3 # 计算三分之一倍频程 print("三分之一倍频程为:", third_octave_bandwidth) ``` 上述代码中的third_octave_bandwidth变量用于保存计算得到的三分之一倍频程的值。最后，通过使用print函数打印该值，我们可以在控制台上看到计算结果。 总结起来，三分之一倍频程是频率范围的其中一种度量方式，而在Python中，我们可以使用变量和数学操作符来计算和处理三分之一倍频程的值。 ### 回答2： 三分之一倍频程是指音频信号中频率的范围，其大小是整个频率范围的三分之一。Python 是一种广泛应用于程序开发和数据科学领域的编程语言。 在音频信号处理中，计算三分之一倍频程可以通过以下步骤实现： 1. 首先，需要获取音频信号的采样数据。可以使用Python中的库如`numpy`或`scipy`来导入和处理音频数据。 2. 接下来，使用傅里叶变换将时域信号转换为频域信号。Python中的`numpy.fft`模块提供了进行快速傅里叶变换的函数。 3. 计算信号的幅度谱，可以使用`numpy.fft.fft`函数获取频域信号的实部和虚部，然后通过计算幅度谱得到信号在每个频率上的幅度大小。 4. 根据幅度谱的大小降序排序，找到频率值为三分之一处的位置。可以使用`numpy.argsort`函数得到排序后的索引，然后取第三分之一处对应的频率值。 5. 最后，计算该频率值对应的频率范围。将该频率值乘以2，得到整个频率范围的三分之一倍频程。 总结起来，通过Python中的库和函数，可以方便地计算音频信号的三分之一倍频程。实现步骤包括获取音频数据、进行傅里叶变换、计算幅度谱、排序、找到三分之一位置的频率值，并最终计算出频率范围。这样可以帮助我们分析和处理音频信号中特定频率范围的信息。 ### 回答3： 三分之一倍频程是指将某个频率范围划分为三等分，每个部分之间的频率间隔相等。在python中，要实现这个功能，可以通过使用numpy库来实现。 首先，需要导入numpy库： ```python import numpy as np ``` 然后，定义一个函数来计算三分之一倍频程： ```python def calculate_third_octave(freq_range): # 计算频率范围的起始频率和终止频率 start_freq = freq_range[0] end_freq = freq_range[1] # 计算三分之一倍频程频率间隔 interval = (end_freq - start_freq) / 3 # 生成三分之一倍频程的频率数组 freq_array = np.arange(start_freq, end_freq, interval) return freq_array ``` 接下来，可以调用这个函数来计算三分之一倍频程： ```python freq_range = (20, 20000) # 设置频率范围为20Hz到20000Hz third_octave_freq = calculate_third_octave(freq_range) # 计算三分之一倍频程 print(third_octave_freq) # 输出三分之一倍频程的频率数组 ``` 这样，就可以得到三分之一倍频程的频率数组。输出结果可能类似于以下内容： ``` [ 20. 100. 500. 2500. 12500.] ``` 以上就是使用python实现三分之一倍频程的方法。

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

白银时代_

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值