一、元音频谱图

1、理论知识

不同元音的频谱不同。其峰值位置称为共振峰。 不同元音的识别取决于共振峰区域的存在。对于不同频率的相同元音,共振峰非常相似,无论音高如何都能识别它。

2、实验结果

声音采样频率为8000Hz,FFT长度为65536。FFT变换后取0—4000Hz进行频谱观察。

Y=fft(y,nfft);                    %用FFT计算频谱
Y_abs=abs(Y(1:nfft2));

五元音频谱图如下所示:

纵轴幅度化为dB后的频谱图:

二、前三共振峰测试结果

2.1 利用倒谱法计算出共振峰

流程如下:

具体步骤:

1、对语音信号进行预加重,再进行加窗和分帧。预加重的滤波器常设为:

H(z)=1−az−1H(z)=1−az−1

H(z)=1-az^{-1}
帧长wlen=256,帧移inc=256,这里将信号分帧后取其中的六帧进行平均后分析。窗函数为汉明窗,长度为帧长。

xx=enframe(x,256,256);                          %对x 256点分为一帧
%a1多帧求平均
x1=xx(1,:);x2=xx(3,:);x3=xx(5,:);x4=xx(7,:);x5=xx(9,:);x6=xx(11,:);
x2=(x1+x2+x3+x4+x5+x6)/6;
u=filter([1 -.99],1,x2);                         % 预加重
u2=u'.*hamming(wlen);                            % 信号加窗函数

2、做傅里叶变换

U=fft(u2); 

3、取幅值后再取对数log(||)log(||)log(| |)

U_abs=log(abs(U(1:wlen2)));                      % 对数幅值谱(实数部分)

4、进行逆傅里叶变换,得到倒谱序列

Cepst=ifft(U_abs);                               % 求取倒谱

5、在倒频率轴上设置一个低通的窗函数,一般设置为矩形窗

式中n0n0n_0是窗函数的宽度。与倒频率的分辨率有关,即和采样频率及FFT的长度有关。

6、傅里叶变换后取实部部分得到频谱包络线。

7、在包络线上寻找出极大值,就可以获得相应的共振峰参数。

% 寻找峰值
j=1;
for i=2:wlen2-1if spect(i)>spect(i-1)&&spect(i)>spect(i+1)Loc(j)=i;Val(j)=spect(i);j=j+1;end
end
2.2 测试结果

低通窗函数n0=7n0=7n_0=7。测试a1画出波形如图所示:

基于matlab的元音共振峰的估算相关推荐

  1. matlab 共振峰检测,基于matlab的语音共振峰的估计.doc

    基于matlab的语音共振峰的估计.doc 题目基于MATLAB的语音共振峰的估计英文题目MATLABBASEDESTIMATESOFANTS院系电子工程学院专业通信工程姓名年级二零零六级指导教师武良 ...

  2. 共振峰检测matlab百度文库,毕业设计论文-基于MATLAB的语音共振峰的估计.doc

    毕业设计论文-基于MATLAB的语音共振峰的估计 题 目: 基于MATLAB的语音共振峰的估计 英文题目: MATLAB-BASED ESTIMATES OF FORMANTS 院 系: 电子工程学院 ...

  3. 共振峰检测matlab,基于 LPC 系数的共振峰估计

    基于 LPC 系数的共振峰估计 此示例说明如何使用线性预测编码 (LPC) 来估计元音共振峰频率.共振峰频率可以通过求预测多项式的根来计算. 此示例使用语音样本 mtlb.mat,它包含在 Signa ...

  4. 【语音合成】基于matlab线性预测共振峰检测和基音参数语音合成【含Matlab源码 562期】

    ⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[语音合成]基于matlab线性预测共振峰检测和基音参数语音合成[含Matlab源码 562期] 点击上面蓝色字体,直接付费下载,即可. ...

  5. 【特征提取】基于matlab共振峰估计【含Matlab源码 550期】

    ⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[特征提取]基于matlab共振峰估计[含Matlab源码 550期] 点击上面蓝色字体,直接付费下载,即可. 获取代码方式2: 付费专栏 ...

  6. 【毕业设计系列】033:基于MATLAB语言的实时变声器系统【含Matlab源码】

    Date:2022.5.5 文章目录 前言 1.算法原理 2.实现效果 2.1.男声变女声 2.2.女声变男声 2.3.男声变大叔声 3.部分matlab代码实现 前言 在大学毕业设计中,要求实现一个 ...

  7. 基于MATLAB的说话人识别系统

    一.课题介绍 说话是人类相互沟通交流最方便.最快捷的一种方式,世界上每一个说话人都拥有自己特定的语音,正如每个人的指纹一样,都是绝无仅有的.说话人识别应用广泛,现已应用到通信.消费电子产品等各个领域. ...

  8. 【毕业设计/Matlab系列】基于MATLAB语言的实时变声器系统

    Date:2022.5.5 文章目录 前言 1.算法原理 2.实现效果 2.1.男声变女声 2.2.女声变男声 2.3.男声变大叔声 3.部分matlab代码实现 前言 在大学毕业设计中,要求实现一个 ...

  9. 《数字语音处理》- 1. 基于MATLAB的语音信号时域特征分析

    声明 本文仅在CSDN发布,其他均为盗版.请支持正版! 正版链接: https://blog.csdn.net/meenr/article/details/117629753 基于MATLAB的语音信 ...

最新文章

  1. plt转pdf软件_无须转Word也能轻松编辑PDF的软件,它来了!
  2. FD.io/VPP — 流量追踪
  3. 锁Lock 那点事儿
  4. java类里的包路径报错_[ERROR] 致命错误: 在类路径或引导类路径中找不到程序包 java.lang...
  5. 设计模式之--单例模式
  6. NHibernate直接执行SQL进行插入
  7. 动手拆解众筹明星产品安全智能锁KeyWe后,发现它不安全且无解
  8. python 线程-threding示例使用
  9. SQL注入漏洞-POST注入
  10. 聚类-KMeans算法(图解算法原理)
  11. 计算快递费系统(java版)
  12. 白化(Whitening) PCA白化 ZCA白化
  13. python 股票估值_隐藏价值的角落:限售股AAP估值及PYTHON实现方法(上)
  14. 彩虹网盘外链程序源码V5.1|网盘外链源码
  15. 普渡大学计算机科学本科,美国普渡大学计算机科学CS本科申请条件及案例
  16. java中的gc是什么意思的缩写_gc是什么意思(gc的缩写是什么意思)
  17. Google账号登录后直接跳转百度首页,登陆不上
  18. 修改Egg.js端口号
  19. 这是广告clearview ai用来将您的脸卖给警察的广告
  20. 节约粮食海报|节约粮食就是对袁隆平院士最好的怀念!

热门文章

  1. java,阳历转阴历
  2. JavaScript-原生Array对象方法详解(一)
  3. android连接sqlite进行简单的增删改查和事务管理
  4. Android 使用Toolbar+DrawerLayout快速实现仿“知乎APP”侧滑导航效果
  5. Linux操作系统上lsof命令详解
  6. 将客户端计算机加入域前建立计算机账户的必要性
  7. android 新浪财经截屏分享,[图]Android 12截图曝光:启用全新UI 增强隐私保护功能...
  8. Xilinx FPGA部分重配置操作流程
  9. pytorch手写数字识别验证四流网络
  10. lc412. Fizz Buzz