语音信号的自相关基音周期检测
（ｃｏｌｆＣ ｍ ｕ ｒ ｎ ｆｒ ａｉ ｎｉ ｅ ｎ ， ｈｎ ｈｉＵ ｉｒ ｙｏ ｌ ｔｃＰｗ ｒ Ｓ ａｇａ ２Ｏ９ ， ｈｎ ） Ｓｈｏ ｏ ｏ ｐｔ ｄＩ ｏ ｔ ｎＥ ｇｎ ｒ ｇ Ｓａ ｇａ ｎｖ ｓ ｅ ｒ ｏｅ， ｈｎ ｈｉ ＯＯＯ Ｃ ｉ ｅａ ｎ ｍ ｏ ｅｉ ｅ ￣ ｆＥ ｃ ｉ ａ
Ｅ ｃｅ ｃ ｆｔ ｅ ａ ｇ ｒ ｈ ｉ ｌｏ ｉ ｏ ｔｎ ｘ ｅ ｔ ｏ ｈ ｕ ｌ ｙ ｏ ｅ ｓｇ ａ ｒ ｃ ｓ ．Ｐ ｔ ｈ ｐ ｒ ｄ ｉ ｆ ｉ ｎ ｙ ｏ ｌｏ ｔ ｍ ｓａ ｓ ｈ ｉ ｍｐ ｒ ｔ ｃ ｐ ｒｔ ｅ ｑ ａｉ ｆ ｈ ｉｎ ｌｐ ｅ ｓ ｉ ｅ ｉ ａ ｅ ｆ ｔ ｔ ｏ ｃ ｏ ｏ ｅ ｄ ｌ ｓ ｕ ｄ ｉ ｏ ｔｉ ｅ ｒ ｕ ｈ ｔ ｅ ｓ ｏ — ｍｅ ａ ｔｃ ｒｅ ａ ｉｎ ｆ ｎ ｔ ｎ，ａ ｄ ｐ ｌ ｒ ｙ ｃ ｒｅａ ｉｎ ｆｔ ｕ ｌ ｏ ｎ ｓ ｂ ａｎ ｄ ｔ ｏ ｇ ｈ ｈ ｒ ｔ ｕ ｏ ｏ ｒ ｌｔ ｕ ｃ ｉ ｈ ｈ ｔ ｉ ｏ ｏ ｎ ｏ ａ ｔ ｏ ｒ ｌ ｔ ｉ ｏ ａ ｄ ｐ ａ ａｕ ｓｉ ｔ ｎ ａ ｅ ｕ ｅ ｏ ｒ ｉｅ ｔ ｅ ｅ ｃ ｅ ｃ ｎ ｔ ｅ ｐ ｏ ｅ ｓ ｏ ｅ ａ ｔｃ ｒ ｅａｉ ｎ ｎ ｅ ｋ ｖ ｌｅ ｅｔ ｍａ ｉ ｒ ｓ ｄ ｔ ａｓ ｆ ｉｎ ｙ ｉ ｒ ｃ ｓ ｆｔ ｕ ｏ ｏ ｒ ｌｔ ． ｏ ｈ ｉ ｈ ｈ ｏ Ｋｅ ｒ ｓ ａ ｔ ｃ ｒｅａ ｉｎ；ｐ ｔ ｈ ｐ ｒ ｄ；ｐ ｌ ｒ ｙ ｃ ｒｅａ ｉｎ ｙ ｗｏ ｄ ： ｕ ｏ ｏ ｌ ｔ ｏ ｉ ｅｏ ｃ ｉ ｏ ａ ｔ ｏ ｌ ｔ ｉ ｏ
基音周期的估计称为基音周期检测， 基音周 期检测的最终 目的是得到与声音振动频率吻合较
好 的基 音周 期 变 化 轨 迹 曲线 ． 语 音 信 号 的 处 理 在
如图 １ 所示 ．
中只有准确捕获语 音信号参 数 ， 才能高效地识别 语音． 而在这些语音信号参数中 ， 基音周期提取的 精确性和效率直接影响到合成语音能否真实快速
为在 ｌ ３ Ｏ～ ０ｍｓ内语 音 特 性 基 本 上 是 不 变 的 ， 或
算法 中采用布莱克曼窗 ， 具体实现代码如下 ： ｙ３＝ ｌ ｒｈ，，３ ； １ ｆｔ （ ｌ， ） ％滤 波处 理 ； ｉｅ ，
Ｌ Ｌ＝ｌｎｔ（ １ ） ２ ０ ｅｇ ｙ３ ／４ ； ｈ
的谐波分量 ， 基音频率最低可达 ８ ｚ最高可达 ０Ｈ ， ５０Ｈ ， ０ ｚ但基音频率大多数分布在 １０ ２０Ｈ ， ０ ～ ０ ｚ为 了提高检测的准确度 ， 算法 中引入了一个 ６ ～ ０ ０ ５０ Ｈ 的带通滤波模块滤除语音帧的高次谐波分量 ； ｚ
１ 语音信 号基 音周 期检测一般过程
１ ３ ｔ ｍ． ６ ．ｏ
上
海
电
力
学
院
学
报
（） ２ 取样模块 从频率为 ８００Ｈ 的语音信 ０ ｚ 号 中截取 Ｌ ｎｔ（ 点数 ） ｅｇ 样 ｈ 长的个样点值 ， 一般取
样 点数 设计 成 帧 长度 的整数 倍 ；
Ⅳ—— 窗 长 ； （ ）＝ ｎ ２ ｎ （ ）．
（） ３ 分帧模块 主要将语 音去样点分割为若 干个语音帧 ， 语音信号属于短时平稳信号 ， 一般认
Ａｂ ｔａ ｔ Ａｕｏ ｏｒｌｔ ｎ ｐｔｈ ｐ ｒ ｄ ｄ ｔｃｉｎ ｉ ｅ ｔｃ ｎｑ ｅ ｉ ｖ ｉｅ ｉｎ ｌ ｐ ｏ ｅｓ ｓｒ ｃ ： ｔｃ ｒｅａｉ ｉ ｅ ｏ ｅｅ ｔ ｓ ａ ｋ ｙ ｅ ｈ ｉｕ ｎ ｏｃ ｓ ａ ｒｃ ｓ ． ｏ ｃ ｉ ｏ ｇ
文章 编号 ：１０ —４ ２ （０ １ ０ Ｏ ９ ｏ ０ ６ ７ ９ ２ １ ）３一 ２ ７一 ４
语 音 信 号 的 自相 关基 音周 期 检 测
赵 萍
２０９ ） ０００
（ 上海 电力学院 计算机与信息工程 学院 ， 上海 摘
要 ：自相关基音周期 检测是语音信 号处理 的关键 技术 ， 了保 证 信号 处理 的质量外还 要注重 算 法 的效 除
图 １ 基 音周 期检测的处理过程
地再现原始语音信号． 本文在基音周期检测一般 方法 的基础上 ， 自相关运算过程加以修改 ， 对 并给
出使用该种方法对两段语音信号进行基音周期检 测 的结果 ．
图 ｌ中各部 分 的功 能及原 理 主要 有 ：
（） １ 带通滤波模块
语音信号包含非常丰富
率． 通过短 时 自相关函数得到浊音语 音的基音 周期 ， 自相关 运算 过程 中采用极性 相关法 和峰值估算 法 以提 在
高运算效率．
关 键词 ：自相关 ； 音周期 ； 性相关 基 极 中图分 类号 ：Ｔ ９ ２ ３ Ｍ １．４ 文献标志码 ：Ａ
Ａｕｔ ｃ ｒ ｅａ ｉ ｎ Ｐｉｃ ｒｏ ｔ ｃｉ ｎ ｏ ｉｅ Ｓ ｇ ａ ｏ ｏ ｒ ｌ ｔｏ ｔｈ Ｐｅ ｉ ｄ Ｄｅ ｅ ｔｏ ｆＶｏ ｃ ｉｎ ｌ
语音信号基 音周期检 测 的一般 处理过 程…
收 稿 日期 ：２ １ ０ ０ ０— ７—１ ２
通讯作者简 介 ：赵 萍（ ９ ９一） 女 ， 士 ， 师 ， １７ ， 硕 讲 吉林敦 化人 ． 主要 研 究方 向为 信号 与信息 处 理． － ａ ：ｑｕ２ ３ Ｅ ｍ ｉ ｑ ｙ０ ０ ＠ ｌ
第２ ７卷第 ３期 ２１ ０ １年 ６月
上
海
电 力 学 院 学
报
Ｖｏ ． ７． Ｎｏ ３ ２ １ ．
Ｊｕ ａ ｏ Ｓ ａ ｇ ａ ｏｒ ｌ ｆ ｈｎ ｌｉ ｎ １ Ｕｎｖｒｉ ｏ Ｅｅｔｃ ｏ ｒ ຫໍສະໝຸດ ｅｓｙ ｆ ｌｃ ｉ Ｐ ｗｅ ｔ ｒ
Ｊｎ ２ １ ｕｅ ０１
语音信号的自相关基音周期检测
语音信号的自相关基音周期检测语音信号是一种人类最基本的交流方式，它包含人声的频率、强度和时间三个方面的信息。
语音信号的自相关和基音周期是语音信号分析中的重要技术，对于语音信号分析、识别、合成等应用有着积极的作用。
一、语音信号的自相关语音信号的自相关是指语音信号的样本与样本之间的相关性。
在语音信号中，相邻的样本之间一般都具有相关性，该相关性可以通过计算信号的自相关函数获得。
自相关函数描述了语音信号在不同延迟情况下的相似程度，也就是说，自相关函数可以反映语音信号的周期特征和基音周期。
二、语音信号的基音周期检测语音信号的基音周期是指人语中相邻两个基音周期之间的时间长度。
基音周期检测是一项关键的语音信号分析技术，在音素识别、语音合成、语音编码等领域应用广泛。
基音周期检测方法主要有三种：自相关法、FFT法和LP法。
自相关法是指通过计算信号与自身在不同延迟下的相似度，判断语音信号的基音周期。
具体来说，自相关法首先将语音信号进行预加重和分帧处理，然后计算每一帧的自相关函数，最后采用模板匹配的方法找到最强的周期峰值，从而得到基音周期。
FFT法则是将分帧后的语音信号进行傅里叶变换，提取频谱信息，并通过在频域滑动一个窗口，检测周期性的能量最大值，确定基音周期数。
LP法是通过线性预测，将语音信号分解成具有不同频率的谐波分量，然后利用实验数据验证模型，得到基音周期。
总的来说，不同的基音周期检测方法有其各自的优缺点。
自相关法较为简单但在噪声环境下准确度不高，FFT法可以检测到非周期性的基音，但精度不如自相关法，LP法精度较高但计算复杂度较大。
作为一门复杂的反演问题，语音信号的自相关和基音周期检测一直是语音处理研究中的重要问题，目前的研究主要集中在解决语音信号分析和识别中的实际问题和应用，为提高语音合成、语音编码等方面的应用水平提供技术支持。
Matlab在基音周期检测教学中的应用研究
Matlab在基音周期检测教学中的应用研究
基音周期检测是语音信号处理中的一个重要任务，它用于测量和分析语音信号的基频，也被称为基音周期或基频周期。
在不同的语音处理应用中都需要对语音信号进行基音周期
检测，例如语音合成、声音识别以及语音增强等。
首先是基于自相关函数的方法。
自相关函数是指语音信号自身与自身经过不同时间延
迟的信号之间的相关性。
通过计算自相关函数的峰值位置和间隔，可以获得语音信号的基
音周期。
在Matlab中，可以使用xcorr函数计算自相关函数，然后通过寻找峰值来确定基音周期。
另外还有一些更高级的基音周期检测算法，例如利用神经网络、隐马尔可夫模型等方
法进行基音周期检测。
在Matlab中，可以使用Neural Network Toolbox和Hidden Markov Model Toolbox等工具箱来实现这些算法。
除了基音周期检测算法，Matlab还提供了许多其他的语音信号处理工具，例如滤波器设计、频谱分析、语音合成等，这些功能都可以与基音周期检测结合使用，以实现更复杂
的语音处理任务。
Matlab在基音周期检测教学中的应用研究非常广泛。
通过使用Matlab提供的工具箱
和函数，可以方便地实现不同的基音周期检测算法，并应用于语音处理的各个领域。
语音信号处理 实验报告用修正的短时自相关检测语音的基音周期
语音信号处理课程实验报告专业班级通信学号姓名指导教师实验名称 用修正的短时自相关检测语音的基音周期 同组人 专业班级通信 学号 姓名 成绩 一、实验目的 1.熟悉前一个实验程序以及中心削波的意义 2.用Matlab 实现用修正的短时自相关检测语音的基音周期。
3.分析修正的短时自相关在基音周期检测中的应用。
4.能够对程序进行重新编制。
二、实验原理 如果x(n)是一个周期为P 的信号，则其自相关函数也是周期为P 的信号，且在信号周期的整数倍处，自相关函数取最大值。
语音的浊音信号具有准周期性，其自相关函数在基音周期的整数倍处取最大值。
计算两相邻最大峰值间的距离，就可以估计出基音周期。
观察浊音信号的自相关函数图，其中真正反映基音周期的只是其中少数几个峰，而其余大多数峰都是由于声道的共振特性引起的。
因此为了突出反映基音周期的信息，同时压缩其他无关信息，减小运算量，有必要对语音信号进行适当预处理后再进行自相关计算以获得基音周期。
第一种方法是先对语音信号进行低通滤波，再进行自相关计算。
因为语音信号包含十分丰富的谐波分量，基音频率的范围分布在50~500Hz 左右，即使女高音升c 调最高也不会超过1kHz ，所以采用1kHz 的低通滤波器先对语音信号进行滤波，保留基音频率；再用2kHz 采样频率进行采样；最后用2~20ms 的滞后时间计算短时自相关，帧长取10~20ms ，即可估计出基音周期。
第二种方法是先对语音信号进行中心削波处理，再进行自相关计算。
本实验采用第二种方法。
且中心削波函数如式（3-1）所示： ……………………………………装………………………………………订…………………………………………线………………………………………()()0()()L L L L L L x x x x f x x x x x x x x ->⎧⎪=-≤≤⎨⎪+<⎩一般削波电平L x 取本帧语音最大幅度的60%~70%。
《语音信号处理》实验2-基音周期估计
华南理工大学《语音信号处理》实验报告实验名称：基音周期估计姓名：学号：班级：10级电信5班日期：2013年5 月15日1.实验目的本次试验的目的是通过matlab编程，验证课本中基音周期估计的方法，本实验采用的方法是自相关法。
2. 实验原理1、基音周期基音是发浊音时声带震动所引起的周期性，而基音周期是指声带震动频率的倒数。
基音周期是语音信号的重要的参数之一，它描述语音激励源的一个重要特征，基音周期信息在多个领域有着广泛的应用，如语音识别、说话人识别、语音分析与综合以及低码率语音编码，发音系统疾病诊断、听觉残障者的语音指导等。
因为汉语是一种有调语言，基音的变化模式称为声调，它携带着非常重要的具有辨意作用的信息，有区别意义的功能，所以，基音的提取和估计对汉语更是一个十分重要的问题。
由于人的声道的易变性及其声道持征的因人而异，而基音周期的范围又很宽，而同—个人在不同情态下发音的基音周期也不同，加之基音周期还受到单词发音音调的影响，因而基音周期的精确检测实际上是一件比较困难的事情。
基音提取的主要困难反映在：①声门激励信号并不是一个完全周期的序列，在语音的头、尾部并不具有声带振动那样的周期性，有些清音和浊音的过渡帧是很难准确地判断是周期性还是非周期性的。
②声道共振峰有时会严重影响激励信号的谐波结构，所以，从语音信号中直接取出仅和声带振动有关的激励信号的信息并不容易。
③语音信号本身是准周期性的(即音调是有变化的)，而且其波形的峰值点或过零点受共振峰的结构、噪声等的影响。
④基音周期变化范围大，从老年男性的50Hz到儿童和女性的450Hz，接近三个倍频程，给基音检测带来了一定的困难。
由于这些困难，所以迄今为止尚未找到一个完善的方法可以对于各类人群(包括男、女、儿童及不向语种)、各类应用领域和各种环境条件情况下都能获得满意的检测结果。
尽管基音检测有许多困难，但因为它的重要性，基音的检测提取一直是一个研究的课题，为此提出了各种各样的基音检测算法，如自相关函数(ACF)法、峰值提取算法(PPA)、平均幅度差函数(AMDF)法、并行处理技术、倒谱法、SIFT、谱图法、小波法等等。
09语言信号处理第九章语音检测分析
�
下图为语音信号对数频谱示意图. 下图为语音信号对数频谱示意图.
频谱包络的慢变分量 虚线所示) (虚线所示) 基音谐波峰值 的快变分量 实线所示) (实线所示)
幅度
频率Hz
通过滤波或再取一次傅立叶反变换,即可将 慢变分量与快变分量分离开. 下图为倒谱的示意图.
幅度
靠近原点的低 倒频部分是频 谐波峰值的变 谱包络的变换换,表示基音. 表示基音.
C[x]
-CL
+CL
x
CL 由语音信号的峰值幅度来确定
(2)中心削波后的语音再计算自相关函数,求第 )中心削波后的语音再计算自相关函数, 一最大峰值点
中心削波示意图
输入语音
中心削 波后的 语音
计算自相关函数的运算量是很大的,可用一些 减少短时自相关运算的有效方法. 如可对中心削波函数进行修正. 如可对中心削波函数进行修正.
第九章 语音检测分析
基音检测
自相关法 倒谱法 简化逆滤波法
共振峰估计
倒谱法 LPC法 LPC法
9.1 基音周期估计
基音是指发浊音时声带振动所引起的周期 性
基音周期是语音信号最重要的参数之一,也称 是语音产生模型中激励源的一个重要参数 对汉语:是一种有调语言,基音的变化模式称为声调, 对汉语:是一种有调语言,基音的变化模式称为声调, 它携带着非常重要的具有辨意作用的信息. 它携带着非常重要的具有辨意作用的信息.
∑| x
w
(m + k ) x w (m) |
也呈现周期性. 对于周期性的 x(n) ,γ n (k ) 也呈现周期性. 与自相关函数相反,在周期的各整数倍点上 γ n (k ) 与自相关函数相反, 具有的是谷值,而不是峰值. 具有的是谷值,而不是峰值.
Matlab在基音周期检测教学中的应用研究
Matlab在基音周期检测教学中的应用研究摘要：目前，基音周期检测是语音信号处理中的一个重要研究方向。
在基音周期检测教学中，基音周期的准确检测和测量对于语音分析和合成至关重要。
本文研究了Matlab在基音周期检测教学中的应用，通过分析语音信号的基频轮廓，探讨了基音周期检测算法的原理和方法，以及基于Matlab的应用实例。
实验结果表明，Matlab在基音周期检测的教学中具有较好的效果和实用价值，有助于学生深入理解和掌握基音周期检测的原理和方法。
关键词：Matlab；基音周期检测；教学；语音信号分析一、引言基音周期检测是语音信号处理中的一个重要研究方向，其主要任务是准确检测和测量语音信号中每个基频周期的起始位置和长度。
基音周期检测的准确性对于语音分析和合成至关重要，因为基音周期特征可以用于声码器的参数估计和语音合成。
目前，基音周期检测技术已经成熟，存在的问题主要是在复杂语音环境中基音周期检测的稳定性和准确性。
在基音周期检测的教学中，如何向学生直观地解释基音周期检测算法的原理和方法，有助于学生理解和掌握基音周期检测的技术。
Matlab是一种强大的数值计算和工程仿真软件，对于语音信号处理有着广泛的应用。
在基音周期检测教学中，Matlab可以方便地实现基音周期检测算法，分析和重构语音信号的基频轮廓，提取基音周期特征。
基音周期检测算法的原理和方法主要包括自相关法、互相关法和基频追踪法。
通过Matlab，可以直观地展示基音周期检测算法的计算流程和实现步骤。
通过Matlab的自相关函数xcorr，可以计算语音信号的自相关函数，进而实现自相关法的基音周期检测算法。
Matlab还提供了一系列的工具箱和函数，如数字信号处理工具箱和声音处理工具箱，可以方便地应用于基音周期检测算法的实现和分析。
基于Matlab的应用实例可以帮助学生深入理解和掌握基音周期检测算法。
学生可以通过Matlab实现基音周期检测算法，并对不同语音信号进行测试和比较。
语音信号基音周期检测的matlab程序
function nmax=find_maxn(r)%寻找峰值最大的n值及基音周期%r,自相关序列%maxn,为峰值最大的nzer=find(r==0); %找第一个零点如果存在jiaocha=0; %找第一近零点ii=1;while (jiaocha<=0)if(r(ii)>0 && r(ii+1)<0 && (ii+1)<length(r))jiaocha=ii;endii=ii+1;if ii==length(r) %没有找到符合要求的点jiaocha=1;endendif length(zer)>0 %检查是否存在零点if zer(1)<jiaocha %存在，则和jiaocha比较大小，用于祛除前点的对基音周期的查找带来的影响jiaocha=zer(1);endendr(1:jiaocha)=0; %祛除影响maxn=max(r); %找最大值temp=find(r==maxn);%返回第一个最大值nmax=temp(1);function jiyinzhouqi(filename,shift)%短时自相关分析%filename语音文件*.wav%zhouqi基音周期shift=10;[signal,fs]=wavread('f:/mywork/1.wav');shift=round(fs*shift); %帧移n1=fix(fs*0.97)+1; %分析起点970ms,帧长30msn2=fix(fs*1)+1;ii=1;for ii=1:(length(signal)-n1)/shift %分析次数if n2<length(signal) %防止溢出data=signal(n1:n2);N=n2-n1+1;R=zeros(1,N); %基音周期（n）多次分析数组for k=1:N-1 %求自相关序列for jj=1:N-kR(k)=R(k)+data(jj)*data(jj+k);endendvalue(ii)=find_maxn(R); %调用基音周期（n）分析函数n1=n1+shift; %移动帧n2=n2+shift;endend%figure(3)%plot(R);%axis([0,1000 -300 300])figure(1)stem(value);axis([0 length(value) 0 1000])len =length(value); %基音周期（n）多次分析数组长度aver=mean(value);index=find(abs((value-aver))>aver/5);value(index)=0; %去除大野点的影响len=len-length(index);for jj=1:3:len/3 %中值平滑，滑动窗口宽度3，精度为中值1/4（剔除野点）average=(value(jj)+value(jj+1)+value(jj+2))/3;for kk=1:3if abs((value(jj-1+kk))-average)>average/4value(jj-1+kk)=0; %将野点置零，同时数组长度减一len=len-1;endendendfigure(2)stem(value);axis([0 length(value) 0 max(value)])Tp=sum(value)/len/fs %求基音周期（Tp）。
基于自相关函数的藏语语音基音周期检测
ｙ）Ｃ（】｛ ０ Ｊ Ｃ （＝ Ｉ ） 【 ， ｘ ： （ ≤
（ ） ￡ ＜－ ＋ｃ ， ） － （ Ｃ￡
（） ５
２ 自相 关 函数 基 音 周 期检 测 的基 本 原 理 ， 对 于离 散 信 号 ＾ 自相 关 函数 定 义为
式 中 为 削 波 电平 。ｘｎ是 原 始 语音 信 号 ，（ 是 削波 后 的语 音 （ ） ｙｎ ） 信 号 。 见 经 中心 削 波 处 理 后 ， 可 由声 道谐 振 特 性 产 生 的 大部 分 阻 ） ） ） （） １ 尼振 荡 已被 削 去 ． 留下 来 的是 包 含 基音 周期 的信 息 。 保 进 一 步分 析 可 知 ． 自相 关 函数 提 取 基 音 周 期 时 ， 心 的是 用 关 Ｒ 表征 了将 信 号 延迟 后 与 它 本 身 的相 似 性 。 １ ｋ是 自相 关 的 延 迟 时 间 。可 以 证 明 ， 自相关 函数 具 有 如下 性 质 。 时 间 ．也 就 是 自相 关 函数 峰 值 出 现 的 位 置 ，而 峰 值 本身 无 关 紧 （） １自相 关 函数 具 有 周 期保 持 性 。 果 ｘｎ 周期 信 号 ， 其 要 。这 样 就 可 以采 用 三 电 平 削 波 法 ， 输入 输 出关 系 为 如 （） 是 则 其 自相 关 数也 具 有 周 期 性 ， 两 者 的 周期 相 同 。 且 自相 关 函 数 提 供 了 ｆ，Ⅲ Ｃ ． ￡ （＞ 估 计 信 号 周期 值 的一 种 方 法 。 ｙ ） Ｃ【 】 ０Ｉｎ Ｃ （ ＝ ’（） ；｛ ，（ Ｊ ￡ ｘ） （） ６ Ｉ１（ ＜ Ｃ 一， ” 一 ￡ ） （１ 号 延迟 后 其 自相 关 函数 保 持 不 变 。 即信 号 ｘｎ的 自相 ２信 （１ 关 函 数 与它 的延 迟 信 号 信 号 的 自相 关 函数 完 全相 同 。 经 三 电平 削 波后 ， 自相关 函数 的计 算 极 为 简单 。用 ｙ ） 削 波 （ 表示 ｎ （ 白噪 声 的 自相 关 函数 具 有 原 点 集 中 性 。 自相 关 函数 用 于 器 的输 出 ． ３ ） 则短 时 自相 关 函数 为 基 音 周 期检 测 时 具 有 一 定 的抗 噪性 尺（） 】＝∑ 【（＋ ）’ ） 玎 ｍ＋ ）’ ＋｝ ，、 ｜ ｝ 玎 ｍｗ ｍ】 （＋ ｗ ｍ ｜】 １ （ 【 （ 】 ） 由语 音 信 号产 生 的数 字 模 型 可 知 ．浊 音 是 由声 门处 形 成 的 脉 冲 串激 励 声 道 产 生 的 ．浊音 等 于 脉 冲 串与 声 道 单 位 冲激 响 应 如果 短 时 窗 是矩 形 窗 ． 上式 变 为 则 Ⅳ 一１ 一ｋ 的卷积。 由于 激 励 源 是 准 周期 性 的 ． 因而 语 音 信 号 也 呈 现 准 周 期 尺 （） ∑ ｙｎ ｍ 玎 ｍ＋ ） ｜ ＝ 】 ｝ （ ＋ ）（ ＋ ｋ ，、 Ｑ 性 ， 重 复 频 率取 决 于 激 励 源 的频 率 （ 其 基音 频 率 ） 但 是 由 于受 到 。 声 道 阻 尼振 荡 的 影 响 。 音 信 号 的 波形 结 构 表 现 得极 为 复 杂 ． 浊 为 式 中 ｙｎ ｍｙｎ ｍ＋ ） 取 值 只 有 一 ， ｌ三 种 情 况 ， 而 不 用 作 （＋ ）（＋ ｋ的 ｌＯ， 因 基 音 周 期 的 自动 检 测 带 来 了极 大 的 困难 。 此 ． 以利 用 自相 关 乘 法 运算 ． 为 可 这就 可 以大 大 节省 计 算 时 间 。 函 数 的 周期 保 持 性 ， 间接 地 来 估 计浊 音 信 号 的 基 音 周 期 。 周 ４ 藏语 语 音 的 基 音周 期 检 测 来 在 ． 期 的 整数 倍 处 。 的 自相 关 函数 可 以达 到极 大 值 ． 而从 自相 关 它 从 图１ 是从藏语语音目 龟 目 １中 气 截取的一段浊 函 数 的第 一 个 峰 值 的 位 置 可 以 间 接地 估 计 出基 音 周 期 ．这 就 是 音 信 号 。采 样 频 率 ｆ ｌｋ 。 ｓ Ｏ Ｈｚ ＝ 用 自相 关 函数 进 行 基 音 周 期 检测 的基 本 原 理 语 音 信 号是 非 平 稳 的 随 机 过 程 ， 特性 是 随 时 间 变 化 的 ． 其 但
基音检测技术及其在语音信号处理中的应用研究
基音检测技术及其在语音信号处理中的应用研究摘要：本文先对自相关函数法、平均幅度差函数法、倒谱法、小波变换法这四种经典的基音检测方法的原理进行分析；然后用MATLAB软件完成了基于短时自相关函数法、倒谱法的基音检测。
通过实验可知，这两种方法针对实验室环境下的一帧语音信号比较准确，但在噪声环境下不够理想。
为了克服这两种方法的不足，本文对它们均做了改进，对于短时自相关函数法增加了预处理，使得基音轨迹曲线更加明显。
在对基音检测的倒谱法进行分析时，提出了一种功率谱二次处理的二次谱基音检测方法，该方法克服了倒谱法基音检测的抗噪能力低的弱点，基音轨迹曲线估计的准确性也得到了改善。
关键词：基音检测；短时能量；自相关函数；倒谱函数；预处理；基音轨迹Abstract：Four typical pitch detection methods are analyzed firstly, i.e, autocorrelation function, average magnitude difference function, cepstrum, and wavelet transform. The pitch detection based on short time energy autocorrelation function and cepstrum functions are implemented with MATLAB. The experimental results show that both methods work well for the speech in the lab. But for the noised signal, the detection performance is not so good. In order to overcome these shortcomings，some improvements are presented in this article .For the short time autocorrelation, there introduce the pro-processing operations ，as a result ，the accuracy of the estimated pitch contour is improved. The cepstrum method of speech pitch detection is analyzed carefully，and also propose a secondary processing power spectrum of secondary spectrum pitch detection methods. It used the power spectrum reprocessing results of speech to extract the pitch contour. The presented method not only overcame the shortcomings of cepstrum method, but also improved the accuracy of the estimated pitch contour.Keywords：pitch detection；short time energy；autocorrelation function；cepstrum function；pro-processing；pitch contour前言基音是指发浊音时声带振动所引起的周期性,而声带振动频率的倒数就是基音周期。
基于自相关的语音基音周期检测方法研究
为信号 ｘｎ和 Ｙ（） （） ｎ的互相关函数 。如果 Ｙ（）ｘ（）则互相关 函 ｎ ＝ ｎ， 数 变 成 了 自相 关 函数 ：
Ｒ（ Ⅱ ｍ）
Ｘ（） ｎ ｍ） ｎｘ（＋
式 中 ， 为 信 号 的延 迟 点 数 。 ｍ 自相 关 函 数反 映 了信 号 ｘｆ 和 自身 ｎ １ 作 了一 段 延 迟之 后 的 ｘ（ ｍ） 似 程 度 。 ｎ 的相 ＋ 自相 关 函 数 具 有 以下 的 性 质 ： 果 序 列 ｘｎ 具 有 周 期 Ｎ ， 如 （） 。 则 其 自相关 函数 也 是 同 周期 的周 期 函 数 。 ： ｘ 】 （ 即 若 （ ＝ｘｎ＋Ｎ１ ｎ ｎ ， 则 Ｒ ＝Ｒｍ ＋ 。 ㈤ （ 清 音 信 号 没 有 周期 性 ，它 的 自相 关 函 数 也 没 有 周 期 ， ㈤ Ｒ 会 随着 １ 的增 大 迅 速 衰 减 。浊 音 信 号 具 有 准 周 期 性 ． 的 自相 １ １ 它 关 函数 Ｒ 具 有 与 ｘｎ 同 的周 期 。浊 音 信 号 的 自相 关 函数 在 ㈤ （１ 相 基 音 周 期 的 整数 倍 位 置 上 出现 峰 值 ．而 清 音 信 号 的 自相 关 函数 没 有 明显 的 峰值 。 自相 关法 基 音 检 测 就 是 利用 Ｒ ｍ 的这 一 性 质 ｆ） 对 语 音 信 号 进行 基 音 检 测 的 。 于 浊音 信 号 。 要检 测 到 Ｎ 的位 对 只 置 ， 可 以估 计 语音 信 号 的基 音周 期 值 。 就
福
建
电
脑
２ ０ 年 第 １ 期 ０８ １
基 于 自相 关 的语 音 基 音 周 期检 测 方 法研 究