音频干扰检测(时域方法)

2024-03-27 13:12

本文主要是介绍音频干扰检测(时域方法),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

请注意注释掉的代码:逐个包络比对就不能加窗了。

import librosa

import numpy as np

from scipy.signal import windows

import matplotlib.pyplot as plt

# 读取音频文件

audio_file = 'sine.wav'

signal, sample_rate = librosa.load(audio_file, sr=None, mono=False)

# 检查通道数并处理信号

if signal.ndim > 1:

    num_channels = signal.shape[0]

    print(f"音频文件有 {num_channels} 个通道")

    # 如果是4通道,取第X个通道进行处理,这里示例取第4个通道(索引为3)

    if num_channels == 2:

        signal = signal[0, :]

else:

    # 如果信号是单通道,直接使用

    print("音频文件是单通道")

# 计算每个周期的采样点数

cycle_samples = int(sample_rate / 1000)

# # 创建汉宁窗

# window_length = cycle_samples * 1  # 窗长度为10个周期

# window = windows.hann(window_length)

# # 对信号的开头和结尾分别应用汉宁窗

# windowed_signal = signal.copy()

# windowed_signal[:window_length//2] *= window[:window_length//2]

# windowed_signal[-window_length//2:] *= window[window_length//2:]

# 计算周期数

num_cycles = len(signal) // cycle_samples

# 存储异常周期的时间点和幅值

anomaly_times = []

anomaly_amplitudes = []

# 逐个周期比较包络

for i in range(num_cycles - 1):

    start = i * cycle_samples

    end = (i + 1) * cycle_samples

    current_cycle = signal[start:end]

    next_cycle = signal[end:end+cycle_samples]

   

    # 计算当前周期和下一个周期的包络差异

    diff = np.abs(current_cycle - next_cycle)

   

    # 如果差异大于阈值,则认为是异常周期

    if np.max(diff) > 0.1:

        anomaly_time = start / sample_rate

        anomaly_times.append(anomaly_time)

        anomaly_amplitudes.append(np.max(np.abs(current_cycle)))

# 打印异常周期的时间点和幅值

for time, amplitude in zip(anomaly_times, anomaly_amplitudes):

    print(f"异常周期时间点: {time:.3f}s, 幅值: {amplitude:.3f}")

# 绘制时域波形图

time = np.arange(len(signal)) / sample_rate

plt.figure(figsize=(8,4))

plt.plot(time, signal, label='Signal')

# 标注异常周期

for t in anomaly_times:

    plt.axvline(x=t, color='r', linestyle='--', label='Anomaly Detected')

plt.xlabel('Time(s)')

plt.ylabel('Amplitude')

plt.title('Windowed Waveform with Anomalies Highlighted')

plt.legend()

plt.show()

这篇关于音频干扰检测(时域方法)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/852204

相关文章

SpringBoot日志配置SLF4J和Logback的方法实现

《SpringBoot日志配置SLF4J和Logback的方法实现》日志记录是不可或缺的一部分,本文主要介绍了SpringBoot日志配置SLF4J和Logback的方法实现,文中通过示例代码介绍的非... 目录一、前言二、案例一:初识日志三、案例二:使用Lombok输出日志四、案例三:配置Logback一

Python实现无痛修改第三方库源码的方法详解

《Python实现无痛修改第三方库源码的方法详解》很多时候,我们下载的第三方库是不会有需求不满足的情况,但也有极少的情况,第三方库没有兼顾到需求,本文将介绍几个修改源码的操作,大家可以根据需求进行选择... 目录需求不符合模拟示例 1. 修改源文件2. 继承修改3. 猴子补丁4. 追踪局部变量需求不符合很

mysql出现ERROR 2003 (HY000): Can‘t connect to MySQL server on ‘localhost‘ (10061)的解决方法

《mysql出现ERROR2003(HY000):Can‘tconnecttoMySQLserveron‘localhost‘(10061)的解决方法》本文主要介绍了mysql出现... 目录前言:第一步:第二步:第三步:总结:前言:当你想通过命令窗口想打开mysql时候发现提http://www.cpp

Mysql删除几亿条数据表中的部分数据的方法实现

《Mysql删除几亿条数据表中的部分数据的方法实现》在MySQL中删除一个大表中的数据时,需要特别注意操作的性能和对系统的影响,本文主要介绍了Mysql删除几亿条数据表中的部分数据的方法实现,具有一定... 目录1、需求2、方案1. 使用 DELETE 语句分批删除2. 使用 INPLACE ALTER T

MySQL INSERT语句实现当记录不存在时插入的几种方法

《MySQLINSERT语句实现当记录不存在时插入的几种方法》MySQL的INSERT语句是用于向数据库表中插入新记录的关键命令,下面:本文主要介绍MySQLINSERT语句实现当记录不存在时... 目录使用 INSERT IGNORE使用 ON DUPLICATE KEY UPDATE使用 REPLACE

CentOS 7部署主域名服务器 DNS的方法

《CentOS7部署主域名服务器DNS的方法》文章详细介绍了在CentOS7上部署主域名服务器DNS的步骤,包括安装BIND服务、配置DNS服务、添加域名区域、创建区域文件、配置反向解析、检查配置... 目录1. 安装 BIND 服务和工具2.  配置 BIND 服务3 . 添加你的域名区域配置4.创建区域

mss32.dll文件丢失怎么办? 电脑提示mss32.dll丢失的多种修复方法

《mss32.dll文件丢失怎么办?电脑提示mss32.dll丢失的多种修复方法》最近,很多电脑用户可能遇到了mss32.dll文件丢失的问题,导致一些应用程序无法正常启动,那么,如何修复这个问题呢... 在电脑常年累月的使用过程中,偶尔会遇到一些问题令人头疼。像是某个程序尝试运行时,系统突然弹出一个错误提

电脑提示找不到openal32.dll文件怎么办? openal32.dll丢失完美修复方法

《电脑提示找不到openal32.dll文件怎么办?openal32.dll丢失完美修复方法》openal32.dll是一种重要的系统文件,当它丢失时,会给我们的电脑带来很大的困扰,很多人都曾经遇到... 在使用电脑过程中,我们常常会遇到一些.dll文件丢失的问题,而openal32.dll的丢失是其中比较

使用Python实现文本转语音(TTS)并播放音频

《使用Python实现文本转语音(TTS)并播放音频》在开发涉及语音交互或需要语音提示的应用时,文本转语音(TTS)技术是一个非常实用的工具,下面我们来看看如何使用gTTS和playsound库将文本... 目录什么是 gTTS 和 playsound安装依赖库实现步骤 1. 导入库2. 定义文本和语言 3

python中字符串拼接的几种方法及优缺点对比详解

《python中字符串拼接的几种方法及优缺点对比详解》在Python中,字符串拼接是常见的操作,Python提供了多种方法来拼接字符串,每种方法有其优缺点和适用场景,以下是几种常见的字符串拼接方法,需... 目录1. 使用 + 运算符示例:优缺点:2. 使用&nbsjsp;join() 方法示例:优缺点:3