基于MATLAB算法数字变声器

基于MATLAB算法的数字变声器【摘要】变声器的原理是通过改变输入声音频率，进行改变声音的音色、音调，使输出声音在感官上与原声音不同。变声器是借助对声音音色和音调的双重复合改变，实现输出声音的改变。通过自己发声，共振峰频率的改变是基于重采样实现的。同时用LPC倒谱分析法分析共振峰的范围，通过MATLAB编写程序并构置GUI界面。实验中通过MATLAB软件对采集来的语言信号进行频谱分析，讨论不同人之间的噪音源参数和共振峰参数，通过对不同人语言信号频率和幅度的改变来实现不同人之间的语音转换。

【关键词】短时自相关法；LPC倒谱；语音信号；matlab GUI；

1引言

随着生活水平的提高，科技的不断进步，很多人为了娱乐，从而希望改变自己的声音；还有如今的许多的访问节目为了保护被访问者，都对声音进行了相应的处理。本设计通过编写matlab程序，借助对声音音色和音调的双重复合改变，实现输出声音的改变。共振峰频率的改变是基于重采样实现的，从重采样原理知道，这也同时印发了基频的变化，为保证基频变化和共振峰频率变化的独立、互不相关，在基频移动时必须考虑抵消重采样带来的偏移，理论上只要基频检测足够精确，确实可以保证基频改变和共振峰频率改变间的互不相关。保证变声效果的自然度主要是没有采用基因检测将基因移动和共振峰变化彻底隔离的缘故。

本次课程设计就是运用我们所学到的理论知识，用MATLAB软件来实现对语音信号的变声处理，理论联系实际，从而更好地掌握以及运用所学习的知识。

2数字变声器的原理与算法

2.1基本原理

语音科学家将人类发声过程视作一个由声门源输送的气流经以声道、口、鼻腔组成的滤波器调制而成的。人类语音可分为有声语音和无声语音，前者是由声带振动激励的脉冲信号经声腔调制变成不同的音，它是人类语言中元音的基础，声带振动的频率称为基频。无声语音则是声带保持开启状态，禁止振动引发的。一般来说，由声门振动决定的基频跟说话人的性别特征有关，如下表，而无声语音则没有体现这个特征。说话人的个性化音色和语音的另外一个声学参数——共振峰频率的分布有关。儿童由于声道短，其共振峰频率高于成年人，成年女性的声道一般短于成年男性，所以女性的共振峰频率一般高于男性。

在进行性别变声时，主要考虑基音周期、基频和共振峰频率的变化。其中男生、女生和童声的基频、共振峰的关系如图1所示；基音周期改变时，基频、共振峰同时变化，若伸展既有男变女、女变童，反之亦可。本实验是基于打开一种声音进行相关参数提取，修改

接近于女声、男声或童声，实现声音的变换。

2.2实现过程及算法

采用线性预测参数合成法。线性预测参数合成法利用LPC语音分析方法，通过分析自然语音样本，计算出LPC系数，就可以建立信号产生模型，从而合成出语音。大致流程如图2所示。

2.2.1基于短时自相关法的基音周期估值

对语言信号进行低通滤波，然后进行自相关计算。在低通滤波时，采用巴特沃斯滤波器。

根据人的说话特征设定相应指标参数，对本段语音设计算出巴特沃斯模拟滤波器的阶数

N为5，3dB截止频率Ωc,Ωc=(100.1a p−1)−1

2n,算出

Ω

c

为0.175，归一化低通原型系统函数

为：

G a（p）=

1

p5+b4p4+b3p3+b2p2+b1p+b0

其中b0=1.0000,b1=3.2361,b2=5.2361,b3=5.2361,b4=3.2361将p=s/Ωc带人G a（p）中，得到低通滤波器，

H s=

Ωc5

s5+b4Ωc2s4+b2Ωc3s2+b1Ωc4s+b0Ωc5

根据设定的滤波器编写matlab程序，当信号经过低通滤波器后，对原始信号滤波产生结果如下图2所示，低通滤波后，保留基音频率，然后再用8kHz采样频率进行采样，采样序列为x(n)，然后进行下一步的自相关计算。

2.2.2语音信号的短时自相关函数图3 原始信号低通滤波

图2 制作流程

定义语音信号自相关函数如下：

R n (k )=

∑

x n (n )x n (n +k )n =k −n −1

n =0

其中k 为信号延迟点数: x n (n )为语音信号；N 为语音帧长度。经过低通滤波之后，取160个样点数，帧长取10ms ，对每帧语音求短时自相关，取得自相关最大点数，自相关函数在基音周期处表现为峰值，这些峰值点之间的间隔的平均值就是基音周期。

2.2.3 LPC 倒谱法提取共振峰

通过线性预测分析得到合成滤波器的系统函数为

H (Z )=1

1−∑a i z −i

p

i =1 其冲击响应为h(n),a i 为预测系数。然后求h(n)的倒谱h ˆ

(n),首先根据同态分析方法，

有H

ˆ(z)=logH(z),因为H （z ）是最小相位的，即在单位圆内是解析的，所以H ˆ(z)可以展开成级数形式，即(z)=∑h ˆ(n )z −n ∞n =1

,也就是说(z)的逆变换(n)是存在的，设(n)=0，将等式两端分别对z −1求导，得到：

∑n h

ˆ(n )∞n =1z

−n +1

=∑u

i z −i +1p

i =1

1−∑a i z

−i

p

i =1

有（1-∑a i z −i p

i =1）∑n (n )∝n =1z

−n +1=∑u i z −i +1∝n =1, 令式左右两边Z 的各次幂前系数分别相等得到：

(1)=1 (n)=a n +∑(1−i n )a i (n −i )n −1i =1 1≤n ≤p (n)=∑(1−i n

)a i h ˆ(n −i )n −1i =1 n>p 在本实验中，取预测阶数为10，语音的样点数为160，按照上式可直接从预测系数a i 求

得倒谱(n)，这个倒谱是根据线性预测模型得到，即称为LPC 倒谱。 2.2.4线性预测语音信号合成

根据线性预测的基本思想，用过去M 个样点值来预测现在或未来的样点值：

y ˆ（n ）=∑a i y (n −i )M i =1

H ˆH ˆh ˆh ˆh ˆh

ˆh ˆh ˆh ˆh

ˆ