如何快速分析音频中的各种频率成分
从视频中提取音频
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)提供了一个高度可靠的存储系统…...
铭豹扩展坞 USB转网口 突然无法识别解决方法
当 USB 转网口扩展坞在一台笔记本上无法识别,但在其他电脑上正常工作时,问题通常出在笔记本自身或其与扩展坞的兼容性上。以下是系统化的定位思路和排查步骤,帮助你快速找到故障原因: 背景: 一个M-pard(铭豹)扩展坞的网卡突然无法识别了,扩展出来的三个USB接口正常。…...
使用VSCode开发Django指南
使用VSCode开发Django指南 一、概述 Django 是一个高级 Python 框架,专为快速、安全和可扩展的 Web 开发而设计。Django 包含对 URL 路由、页面模板和数据处理的丰富支持。 本文将创建一个简单的 Django 应用,其中包含三个使用通用基本模板的页面。在此…...
Zustand 状态管理库:极简而强大的解决方案
Zustand 是一个轻量级、快速和可扩展的状态管理库,特别适合 React 应用。它以简洁的 API 和高效的性能解决了 Redux 等状态管理方案中的繁琐问题。 核心优势对比 基本使用指南 1. 创建 Store // store.js import create from zustandconst useStore create((set)…...
逻辑回归:给不确定性划界的分类大师
想象你是一名医生。面对患者的检查报告(肿瘤大小、血液指标),你需要做出一个**决定性判断**:恶性还是良性?这种“非黑即白”的抉择,正是**逻辑回归(Logistic Regression)** 的战场&a…...
多场景 OkHttpClient 管理器 - Android 网络通信解决方案
下面是一个完整的 Android 实现,展示如何创建和管理多个 OkHttpClient 实例,分别用于长连接、普通 HTTP 请求和文件下载场景。 <?xml version"1.0" encoding"utf-8"?> <LinearLayout xmlns:android"http://schemas…...
AtCoder 第409场初级竞赛 A~E题解
A Conflict 【题目链接】 原题链接:A - Conflict 【考点】 枚举 【题目大意】 找到是否有两人都想要的物品。 【解析】 遍历两端字符串,只有在同时为 o 时输出 Yes 并结束程序,否则输出 No。 【难度】 GESP三级 【代码参考】 #i…...
【解密LSTM、GRU如何解决传统RNN梯度消失问题】
解密LSTM与GRU:如何让RNN变得更聪明? 在深度学习的世界里,循环神经网络(RNN)以其卓越的序列数据处理能力广泛应用于自然语言处理、时间序列预测等领域。然而,传统RNN存在的一个严重问题——梯度消失&#…...
全面解析各类VPN技术:GRE、IPsec、L2TP、SSL与MPLS VPN对比
目录 引言 VPN技术概述 GRE VPN 3.1 GRE封装结构 3.2 GRE的应用场景 GRE over IPsec 4.1 GRE over IPsec封装结构 4.2 为什么使用GRE over IPsec? IPsec VPN 5.1 IPsec传输模式(Transport Mode) 5.2 IPsec隧道模式(Tunne…...
听写流程自动化实践,轻量级教育辅助
随着智能教育工具的发展,越来越多的传统学习方式正在被数字化、自动化所优化。听写作为语文、英语等学科中重要的基础训练形式,也迎来了更高效的解决方案。 这是一款轻量但功能强大的听写辅助工具。它是基于本地词库与可选在线语音引擎构建,…...
springboot整合VUE之在线教育管理系统简介
可以学习到的技能 学会常用技术栈的使用 独立开发项目 学会前端的开发流程 学会后端的开发流程 学会数据库的设计 学会前后端接口调用方式 学会多模块之间的关联 学会数据的处理 适用人群 在校学生,小白用户,想学习知识的 有点基础,想要通过项…...
