如何快速分析音频中的各种频率成分
从视频中提取音频
from moviepy.editor import VideoFileClip# Load the video file and extract audio
video_path = "/mnt/data/WeChat_20241026235630.mp4"
video_clip = VideoFileClip(video_path)# Extract audio and save as a temporary file for further analysis
audio_path = "/mnt/data/extracted_audio.wav"
video_clip.audio.write_audiofile(audio_path)audio_path
分析分离出的音频
import numpy as np
import matplotlib.pyplot as plt
from scipy.io import wavfile
from scipy.fft import fftaudio_path="/mnt/data/extracted_audio.wav"
# Load the audio file
sample_rate, audio_data = wavfile.read(audio_path)# If the audio has two channels (stereo), convert it to mono by averaging the channels
if audio_data.ndim > 1:audio_data = audio_data.mean(axis=1)# Select a segment of the audio for analysis (e.g., the first 5 seconds)
duration_seconds = 5
num_samples = duration_seconds * sample_rate
audio_segment = audio_data[:num_samples]# Perform FFT on the audio segment to get the frequency spectrum
audio_fft = fft(audio_segment)
frequencies = np.fft.fftfreq(len(audio_segment), 1 / sample_rate)
magnitude = np.abs(audio_fft)# Filter positive frequencies for plotting
positive_frequencies = frequencies[:len(frequencies) // 2]
positive_magnitude = magnitude[:len(magnitude) // 2]# Plot the frequency spectrum
plt.figure(figsize=(12, 6))
plt.plot(positive_frequencies, positive_magnitude)
plt.title("Frequency Spectrum of Audio")
plt.xlabel("Frequency (Hz)")
plt.ylabel("Magnitude")
plt.xlim(0, 5000) # Focus on lower frequencies (common range for mechanical noises)
plt.grid(True)
plt.show()
这是音频的频谱图,展示了在0到5000 Hz范围内的频率成分。通过观察图中的频率峰值,我们可以得出以下初步分析:
主要频率峰值:
- 如果在高频范围(1000 Hz以上)有明显峰值,可能来源于电机的旋转噪音,尤其是电机的无刷运行特性。
- 如果在中低频范围(200-1000Hz)出现峰值,可能是齿轮啮合、传动系统或轮胎摩擦噪声。
- 低频段(100 Hz以下)的峰值一般是由振动噪声或结构传递噪音引起。
噪声来源可能性:
- 电机噪声:如果有高频噪声,电机或其驱动电路可能是主要噪声源。
- 传动系统噪声:中频峰值可能与齿轮啮合或皮带传动有关,可能由于齿轮不对中、松动或缺乏润滑等原因引起。
- 轮胎和路面摩擦:在较低频率的情况中,轮胎与路面接触的摩擦噪声较为常见。
常见频率峰值的音频噪声源
常见的频率峰值可以帮助识别音频噪声的来源。以下是一些常见的噪声源及其对应的频率范围和特点:
- 电机噪声
无刷直流电机(BLDC):主要噪声峰值在高频段(1000 Hz以上),因其高速旋转和换向特性,可能会在电机转速的倍频位置产生峰值。
有刷电机:通常会产生中频至高频噪声(300-2000 Hz),由于电刷和换向器的摩擦,带来高频成分。
步进电机:通常在200-600 Hz之间产生噪声峰值,其噪声源于步进脉冲和扭矩振动,尤其在低速时更明显。 - 齿轮传动噪声
齿轮啮合频率:齿轮传动系统中最明显的频率峰值通常是齿轮啮合频率(Gear Mesh Frequency, GMF),典型频率范围在200-1000 Hz之间。
倍频噪声:齿轮系统可能会在啮合频率的倍频位置产生额外的噪声峰值,尤其在齿轮不对中或磨损情况下。 - 轴承噪声
滚珠轴承缺陷频率:轴承噪声通常位于中高频(500 Hz以上),主要来源于内外圈、滚珠或保持架的缺陷。
轴承故障频率: - 外圈缺陷:产生的频率取决于滚珠通过外圈的频率,通常较低频段。
- 内圈缺陷:通常位于较高频段,频率与轴承的转速和滚珠通过内圈的频率有关。
保持架缺陷:产生频率较低,因为保持架运动速度较慢。 - 振动和结构共振
机械共振频率:机械结构和支架可能会在某些频率上产生共振峰值,通常是低频段(20-200 Hz)。共振频率依赖于结构的材料、形状和支撑方式。
固有频率:当外界激励频率接近结构的固有频率时,可能会引起大幅振动,并导致明显的低频噪声。 - 风扇和冷却系统噪声
风扇叶片通过频率(Blade Pass Frequency, BPF):当风扇运行时,叶片通过的频率会产生一个显著的频率峰值,通常在100-500 Hz之间,频率峰值取决于叶片数和风扇速度。
空气涡流噪声:高速空气流动会产生随机噪声,频率分布较宽,但一般集中在中高频段(500 Hz以上)。 - 皮带和链条传动噪声
皮带传动噪声:频率通常在中低频(100-500 Hz),主要由皮带振动和摩擦引起。
链条传动噪声:与皮带相似,链条传动也会在100-500 Hz范围内产生噪声,特别是链条松紧度不当或磨损时。 - 轮胎和地面摩擦噪声
低频摩擦噪声:轮胎与地面的接触摩擦通常在低频段(20-200 Hz),特别在硬地面或粗糙地面上行驶时更明显。
振动噪声:轮胎不平衡或路面不平整会在较低频率产生明显的振动噪声(通常低于100 Hz)。 - 电磁干扰噪声
电源开关噪声:常见于开关电源,通常在数千赫兹到十几千赫兹的频率范围,可能通过电磁辐射方式产生噪声。
变频器噪声:变频器控制的电机可能产生明显的高频噪声(1000 Hz以上),尤其是PWM(脉宽调制)频率。
通过识别频谱图中的这些常见频率峰值,可以帮助快速定位噪声来源。
相关文章:

如何快速分析音频中的各种频率成分
从视频中提取音频 from moviepy.editor import VideoFileClip# Load the video file and extract audio video_path "/mnt/data/WeChat_20241026235630.mp4" video_clip VideoFileClip(video_path)# Extract audio and save as a temporary file for further anal…...
MongoDB 6.0 主从复制配置
以下是 MongoDB 6.0 版本配置主从的详细安装步骤: 1. 安装 MongoDB:可以从官网下载 MongoDB 6.0 的安装包并进行安装,或者使用相应的包管理工具进行安装。 2. 配置主节点:在主节点的 MongoDB 配置文件(默认路径为 …...

NPU 神经网络处理单元
Ⅰ 什么是 NPU? 当前正处于神经网络和机器学习处理需求爆发的初期。传统的 CPU(中央处理器)/GPU(图形处理器)可以执行类似任务,但专门为神经网络优化的 NPU(神经处理单元)比 CPU/GP…...

安宝特分享 | AR技术引领:跨国工业远程协作创新模式
在当今高度互联的工业环境中,跨国合作与沟通变得日益重要。然而,语言障碍常常成为高效协作的绊脚石。安宝特AR眼镜凭借其强大的多语言自动翻译和播报功能,正在改变这一局面,让远程协作变得更加顺畅。 01 多语言翻译优势 安宝特A…...

Vulkan 开发(五):Vulkan 逻辑设备
图片来自《Vulkan 应用开发指南》 Vulkan 开发系列文章: 1. 开篇,Vulkan 概述 2. Vulkan 实例 3. Vulkan 物理设备 4. Vulkan 设备队列 在 Vulkan 中,逻辑设备(Logical Device)是与物理设备(Physical D…...
Kafka 解决消息丢失、乱序与重复消费
一、引言 在分布式系统中,Apache Kafka 作为一种高吞吐量的分布式发布订阅消息系统,被广泛应用于日志收集、流式处理、消息队列等场景。然而,在实际使用过程中,可能会遇到消息丢失、乱序、重复消费等问题,这些问题可能…...

计算机专业毕业生面试工具推荐:白瓜面试
随着毕业季的临近,计算机专业的毕业生们即将步入职场,面试成为了他们必须面对的挑战。在这个过程中,选择合适的面试工具可以大大提高求职成功率。今天,我要向大家推荐一款专为计算机专业毕业生设计的面试工具——白瓜面试。 为什…...

数字IC开发:布局布线
数字IC开发:布局布线 前端经过DFT,综合后输出网表文件给后端,由后端通过布局布线,将网表转换为GDSII文件;网表文件只包含单元器件及其连接等信息,GDS文件则包含其物理位置,具体的走线࿱…...

高空作业未系安全带监测系统 安全带穿戴识别预警系统
在各类高空作业场景中,安全带是保障作业人员生命安全的关键防线。然而,由于人为疏忽或其他原因,作业人员未正确系挂安全带的情况时有发生,这给高空作业带来了巨大的安全隐患。为有效解决这一问题,高空作业未系安全带监…...

k8s的配置和存储(ConfigMap、Secret、Hostpath、EmptyDir以及NFS的服务使用)
ConfigMap 简介 在 Kubernetes 中,ConfigMap 是一种用于存储非敏感信息的 Kubernetes 对象。它用于存储配置数据,如键值对、整个配置文件或 JSON 数据等。ConfigMap 通常用于容器镜像中的配置文件、命令行参数和环境变量等。 ConfigMap 可以通过三种方…...

JS轮播图实现自动轮播、悬浮停止轮播、点击切换,下方指示器与图片联动效果
代码: <!DOCTYPE html> <html lang"zh-CN"> <head><meta charset"UTF-8"><meta name"viewport" content"widthdevice-width, initial-scale1.0"><title>Document</title><s…...

使用 Kafka 和 MinIO 实现人工智能数据工作流
MinIO Enterprise Object Store 是用于创建和执行复杂数据工作流的基础组件。此事件驱动功能的核心是使用 Kafka 的 MinIO 存储桶通知。MinIO Enterprise Object Store 为所有 HTTP 请求(如 PUT、POST、COPY、DELETE、GET、HEAD 和 CompleteMultipartUpload…...

力扣题86~90
题86(中等): python代码 # Definition for singly-linked list. # class ListNode: # def __init__(self, val0, nextNone): # self.val val # self.next next class Solution:def partition(self, head: Optional[Li…...

【JavaEE】【多线程】定时器
目录 一、定时器简介1.1 Timer类1.2 使用案例 二、实现简易定时器2.1 MyTimerTask类2.2 实现schedule方法2.3 构造方法2.4 总代码2.5 测试 一、定时器简介 定时器:就相当于一个闹钟,当我们定的时间到了,那么就执行一些逻辑。 1.1 Timer类 …...

CI/CD 的原理
一、CI/CD 的概念 CI/CD是一种软件开发流程,旨在通过自动化和持续的集成、测试和交付实现高质量的软件产品。 CI(Continuous Integration)持续集成 目前主流的开发方式是协同开发,即多位开发人员同事处理同意应用不同模块或功能。 如果企业在同一时间将…...

进一步认识ICMP协议
在日常工作中,我们经常需要判断网络是否连通,相信大家使用较多的命令就是 ping啦。ping命令是基于 ICMP 协议来实现的,那么什么是 ICMP 协议呢?ping命令又是如何基于 ICMP 实现的呢? 今天这篇文章,我们就来…...

NUUO网络视频录像机upload.php任意文件上传漏洞复现
文章目录 免责声明漏洞描述搜索语法漏洞复现nuclei修复建议 免责声明 本文章仅供学习与交流,请勿用于非法用途,均由使用者本人负责,文章作者不为此承担任何责任 漏洞描述 NUUO网络视频录像机(Network Video Recorder࿰…...
WebGL 3D基础
1. 归一化函数 对一个向量进行归一化处理,即调整向量的模长(长度)为1,同时保持其方向不变。 // 归一化函数 function normalized(arr) {let sum 0;for (let i 0; i < arr.length; i) {sum arr[i] * arr[i];}const middle …...

Docker 部署MongoDb
1. 编写docker-compose.conf 文件 version: 3 services:mongo:image: mongo:latest # 指定 MongoDB 版本,确保 > 3.6container_name: mongo-replicarestart: alwayscommand: ["mongod", "--replSet", "rs0", "--oplogSize&…...
【Hadoop】hadoop的路径分不清?HDFS路径与本地文件系统路径的区别
/usr/local/hadoop /user/hadoop /home/hadoop/ 这里有些路径名很相似,帮我区分? 在Hadoop生态系统中,理解文件存储的位置对于有效管理数据至关重要。Hadoop分布式文件系统(HDFS)提供了一个高度可靠的存储系统…...

Mybatis逆向工程,动态创建实体类、条件扩展类、Mapper接口、Mapper.xml映射文件
今天呢,博主的学习进度也是步入了Java Mybatis 框架,目前正在逐步杨帆旗航。 那么接下来就给大家出一期有关 Mybatis 逆向工程的教学,希望能对大家有所帮助,也特别欢迎大家指点不足之处,小生很乐意接受正确的建议&…...

YSYX学习记录(八)
C语言,练习0: 先创建一个文件夹,我用的是物理机: 安装build-essential 练习1: 我注释掉了 #include <stdio.h> 出现下面错误 在你的文本编辑器中打开ex1文件,随机修改或删除一部分,之后…...

Opencv中的addweighted函数
一.addweighted函数作用 addweighted()是OpenCV库中用于图像处理的函数,主要功能是将两个输入图像(尺寸和类型相同)按照指定的权重进行加权叠加(图像融合),并添加一个标量值&#x…...

学校招生小程序源码介绍
基于ThinkPHPFastAdminUniApp开发的学校招生小程序源码,专为学校招生场景量身打造,功能实用且操作便捷。 从技术架构来看,ThinkPHP提供稳定可靠的后台服务,FastAdmin加速开发流程,UniApp则保障小程序在多端有良好的兼…...

什么是库存周转?如何用进销存系统提高库存周转率?
你可能听说过这样一句话: “利润不是赚出来的,是管出来的。” 尤其是在制造业、批发零售、电商这类“货堆成山”的行业,很多企业看着销售不错,账上却没钱、利润也不见了,一翻库存才发现: 一堆卖不动的旧货…...
ffmpeg(四):滤镜命令
FFmpeg 的滤镜命令是用于音视频处理中的强大工具,可以完成剪裁、缩放、加水印、调色、合成、旋转、模糊、叠加字幕等复杂的操作。其核心语法格式一般如下: ffmpeg -i input.mp4 -vf "滤镜参数" output.mp4或者带音频滤镜: ffmpeg…...

现代密码学 | 椭圆曲线密码学—附py代码
Elliptic Curve Cryptography 椭圆曲线密码学(ECC)是一种基于有限域上椭圆曲线数学特性的公钥加密技术。其核心原理涉及椭圆曲线的代数性质、离散对数问题以及有限域上的运算。 椭圆曲线密码学是多种数字签名算法的基础,例如椭圆曲线数字签…...
实现弹窗随键盘上移居中
实现弹窗随键盘上移的核心思路 在Android中,可以通过监听键盘的显示和隐藏事件,动态调整弹窗的位置。关键点在于获取键盘高度,并计算剩余屏幕空间以重新定位弹窗。 // 在Activity或Fragment中设置键盘监听 val rootView findViewById<V…...

pikachu靶场通关笔记22-1 SQL注入05-1-insert注入(报错法)
目录 一、SQL注入 二、insert注入 三、报错型注入 四、updatexml函数 五、源码审计 六、insert渗透实战 1、渗透准备 2、获取数据库名database 3、获取表名table 4、获取列名column 5、获取字段 本系列为通过《pikachu靶场通关笔记》的SQL注入关卡(共10关࿰…...
Device Mapper 机制
Device Mapper 机制详解 Device Mapper(简称 DM)是 Linux 内核中的一套通用块设备映射框架,为 LVM、加密磁盘、RAID 等提供底层支持。本文将详细介绍 Device Mapper 的原理、实现、内核配置、常用工具、操作测试流程,并配以详细的…...