CN105139863A - 一种音频频域连续性图谱计算方法 - Google Patents
一种音频频域连续性图谱计算方法 Download PDFInfo
- Publication number
- CN105139863A CN105139863A CN201510362518.4A CN201510362518A CN105139863A CN 105139863 A CN105139863 A CN 105139863A CN 201510362518 A CN201510362518 A CN 201510362518A CN 105139863 A CN105139863 A CN 105139863A
- Authority
- CN
- China
- Prior art keywords
- audio frequency
- frequency domain
- audio
- fourier transform
- continuity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种音频频域连续性图谱计算方法，其中，包括如下步骤：步骤a、将一时长为T′的音频文件分割为时长为T的片段，所有音频片段标记为Ai,i＝1,…,k，k为片段总数；步骤b、对每个Ai以窗长2L进行傅里叶变换计算，提取前半部分长度为L的傅里叶变换因子，所有音频片段的长度为L的傅里叶变换因子标记为Fi,j，i＝1，...，k，j＝1，...，L；步骤c、计算Fi，j的行平均向量，标记为mF；步骤d、计算Fi，j中的所有行向量Fi，：与向量mF的向量相关性，得到所有的相关系数标记为Ci，i＝1，...，k；步骤e、在二维图谱中绘制Ci。本发明的有益效果是：提出一种音频频域连续性分析的计算方法，以及呈现音频频率域连续性的图谱表征方法。
Description
技术领域
本发明涉及的是一种计算音频频域连续性图谱表征的方法。
背景技术
随着视听资料成为我国《刑事诉讼法》和《民事诉讼法》等规定的法定证据形式，录音资料作为视听资料的重要组成部分，其真实性司法鉴定新技术新方法的研究具有重要的理论研究意义和实际应用价值。录音真实性鉴定主要涵盖录音资料的客观性、原始性、连续性和完整性等属性方面的技术鉴定。现有的针对录音资料的检验方法和技术主要包括听觉检验以及频谱检验等技术；针对频谱检验方法，目前的技术仍停留在录音频谱图的计算和频谱图表征初级层面上，缺乏对录音频谱连续性的有效分析计算方法和针对录音频谱连续性分析的直观表征图谱呈现方法。
发明内容
针对音频取证技术在音频频域分析所存在的上述问题，本专利提出一种音频频域连续性分析的计算方法，以及呈现音频频率域连续性的图谱表征方法。
本发明解决技术问题所采用的技术手段为：
一种音频频域连续性图谱计算方法，其中，包括如下步骤：
步骤a、将一总时长为T′的音频文件分为每段时长为T的小段，最后时长不足T的部分也作为一段，所述音频文件包含的所述小段数为k个，所有音频小段标记为Ai,i＝1,…,k；
步骤b、对每个Ai以窗口长度为2L进行傅里叶变换计算，提取前半部分长度为L的傅里叶变换因子，所有音频小段计算的长度为L的傅里叶变换因子标记为Fi，j，i＝1，...，k,j＝1，...，L；
步骤c、计算Fi，j的行平均向量，标记为mF；
步骤d、计算Fi，j中的所有行向量Fi，：与向量mF的向量相关性，得到所有的相关系数向量标记为Ci，i＝1，...，k；
步骤e、在二维图谱中绘制Ci。
上述音频频域连续性图谱计算方法，其中，所述步骤e中的二维图谱的横轴为时间，纵轴为相关系数。
本发明的有益效果是：
1、本发明提出了一种分析音频在频率域上的连续性计算方法。
2、本发明提出了一种呈现音频频率域连续性的图谱表征方法。
3、本发明提出的呈现音频频率域连续性的图谱表征方法，其图谱直观、易懂，特征表现明确、清晰，将成为分析录音资料的有效检验图谱。
附图说明
图1是本发明一种音频频域连续性图谱计算方法的流程图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明，但不作为本发明的限定。
本实施例提供了一种音频频域连续性图谱计算方法，只需输入一段音频文件便可迅速的得出反映其频域连续性的图谱。整个工作流程如图1所示，本实施例包括如下步骤：
步骤a、首先选择一个合适的常数T对音频文件进行没有重叠的分割，在此以1秒作为一个分割单元进行说明。须注意，选择1秒作为分割单元仅仅是为了可以更直观的说明，而不是将本发明的实施方式限制于此，以下所举出的具体数字均可作此理解。
不妨设给定的音频材料的播放时间T′＝120秒，采样率为96KHz，则我们可以把该音频文件分成k＝120个播放时间都为1秒的小片段，标记为Ai,i＝1,…,120。
步骤b、对Ai中的每个小片段以窗口长度为2*2048进行傅里叶变换计算，提取前半部分长度为2048的傅里叶变换因子，所有音频小段计算的长度为2048的傅里叶变换因子标记为Fi，j，i＝1，...，120,j＝1，...，2048；
步骤c、计算Fi，j的行平均向量，标记为mF；
步骤d、计算Fi，j中的所有行向量Fi，：与向量mF的向量相关性，这里可以采用皮尔逊相关系数计算方法得到各个行向量Fi，：与向量mF的相关系数，标记为Ci，i＝1，...，120；
步骤e、以时间片段[0,120]范围为横轴，以相关系数[0,1]范围为纵轴在二维图谱中绘制Ci。
以上所述仅为本发明较佳的实施例，并非因此限制本发明的申请专利范围，所以凡运用本发明说明书及图示内容所做出的等效结构变化，均包含在本发明的保护范围内。
Claims (2)
1.一种音频频域连续性图谱计算方法，其特征在于，包括如下步骤：
步骤a、将一总时长为T′的音频文件分为每段时长为T的小段，最后时长不足T的部分也作为一段，所述音频文件包含的所述小段数为k个，所有音频小段标记为Ai,i＝1,…,k；
步骤b、对每个Ai以窗口长度为2L进行傅里叶变换计算，提取前半部分长度为L的傅里叶变换因子，所有音频小段计算的长度为L的傅里叶变换因子标记为Fi，j，i＝1，...，k,j＝1，...，L；
步骤c、计算Fi，j的行平均向量，标记为mF；
步骤d、计算Fi，j中的所有行向量Fi，：与向量mF的向量相关性，得到所有的相关系数向量标记为Ci，i＝1，...，k；
步骤e、在二维图谱中绘制Gi。
2.如权利要求1所述音频频域连续性图谱计算方法，其特征在于，所述步骤e中二维图谱的横轴为时间，纵轴为相关系数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510362518.4A CN105139863B (zh) | 2015-06-26 | 2015-06-26 | 一种音频频域连续性图谱计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510362518.4A CN105139863B (zh) | 2015-06-26 | 2015-06-26 | 一种音频频域连续性图谱计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105139863A true CN105139863A (zh) | 2015-12-09 |
| CN105139863B CN105139863B (zh) | 2020-07-21 |
Family
ID=54725184
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510362518.4A Active CN105139863B (zh) | 2015-06-26 | 2015-06-26 | 一种音频频域连续性图谱计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105139863B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106652998A (zh) * | 2017-01-03 | 2017-05-10 | 中国农业大学 | 一种基于fft短时傅里叶算法的语音综合电路结构及其控制方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6253026A (ja) * | 1985-09-02 | 1987-03-07 | Nec Corp | 適応直交変換符号化方式とその装置 |
| KR20000023379A (ko) * | 1998-09-24 | 2000-04-25 | 이데이 노부유끼 | 정보 처리 장치 및 방법, 정보 기록 장치 및 방법, 기록매체 및 제공 매체 |
| CN101425291A (zh) * | 2007-10-31 | 2009-05-06 | 株式会社东芝 | 语音处理装置及语音处理方法 |
| CN103137138A (zh) * | 2011-11-25 | 2013-06-05 | 司法部司法鉴定科学技术研究所 | 一种音频重复插入的检测方法 |
-
2015
- 2015-06-26 CN CN201510362518.4A patent/CN105139863B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6253026A (ja) * | 1985-09-02 | 1987-03-07 | Nec Corp | 適応直交変換符号化方式とその装置 |
| KR20000023379A (ko) * | 1998-09-24 | 2000-04-25 | 이데이 노부유끼 | 정보 처리 장치 및 방법, 정보 기록 장치 및 방법, 기록매체 및 제공 매체 |
| CN101425291A (zh) * | 2007-10-31 | 2009-05-06 | 株式会社东芝 | 语音处理装置及语音处理方法 |
| CN103137138A (zh) * | 2011-11-25 | 2013-06-05 | 司法部司法鉴定科学技术研究所 | 一种音频重复插入的检测方法 |
Non-Patent Citations (1)
| Title |
|---|
| 曾锦华: "录音设备识别司法鉴定技术研究", 《鉴定科学》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106652998A (zh) * | 2017-01-03 | 2017-05-10 | 中国农业大学 | 一种基于fft短时傅里叶算法的语音综合电路结构及其控制方法 |
| CN106652998B (zh) * | 2017-01-03 | 2021-02-02 | 中国农业大学 | 一种基于fft短时傅里叶算法的语音综合电路结构及其控制方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105139863B (zh) | 2020-07-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2458335A3 (en) | Road estimation device and method for estimating road | |
| US20220366213A1 (en) | Convolution acceleration operation method and apparatus, storage medium and terminal device | |
| JP2012095309A5 (zh) | ||
| CN106571146B (zh) | 噪音信号确定方法、语音去噪方法及装置 | |
| CN104484410A (zh) | 应用于大数据系统的数据融合方法及系统 | |
| EP2458332A3 (en) | Road estimation device and method for estimating road | |
| CN106303524B (zh) | 基于预测残差异常模式的视频双重压缩检测方法 | |
| CN104484412A (zh) | 基于多形式处理的大数据分析系统 | |
| CN116486834A (zh) | 基于特征融合和改进卷积神经网络的轧制声音分类方法 | |
| CN105139863A (zh) | 一种音频频域连续性图谱计算方法 | |
| CN110489810B (zh) | 一种基于数据块的趋势自动提取方法 | |
| CN104484409A (zh) | 用于大数据处理的数据挖掘方法 | |
| CN107037262B (zh) | 一种大数据频谱分析系统及其方法 | |
| CN112817049B (zh) | 一种声波时差的计算方法 | |
| CN103117063A (zh) | 一种基于软件实现的音乐内容截幅检测方法 | |
| Xu et al. | IPSO-VMD based signal feature extraction and internal defect detection of hardwood logs through acoustic impact test | |
| CN115508335A (zh) | 基于傅里叶变换的拉曼光谱曲线数据增强方法 | |
| Cantzos et al. | Fractal analysis, information-theoretic similarities and svm classification for multichannel, multi-frequency pre-seismic electromagnetic measurements | |
| US10303847B2 (en) | Single molecule identification using intensity time sequencing, line charting and run-length coding | |
| CN105030231A (zh) | 一种胎心监护加速特征参数的提取方法 | |
| CN105631029A (zh) | 基于用户主观行为分析的网络监控系统 | |
| CN111027599A (zh) | 基于随机抽样的聚类可视化方法及装置 | |
| Milillo et al. | Spatial statistics and interpolation methods for TOF SIMS imaging | |
| Wu et al. | Synchro-reassigning generalized Chirplet basis transform for analyzing signals with crossover frequencies | |
| US20160277732A1 (en) | Method of motion segmentation in video using randomized voting and system having the same |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information |
Address after: 200063 No. 1347 Guangfu West Road, Shanghai, Putuo District Applicant after: ACADEMY OF FORENSIC SCIENCE Address before: 200063 No. 1347 Guangfu West Road, Shanghai, Putuo District Applicant before: INSTITUTE OF FORENSIC SCIENCE, MINISTRY OF JUSTICE PRC |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |