【信号处理教程】基于倒谱图判断浊音的基音周期附MATLAB代码

简介: 【信号处理教程】基于倒谱图判断浊音的基音周期附MATLAB代码

✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,matlab项目合作可私信。

🍎个人主页:Matlab科研工作室

🍊个人信条:格物致知。

更多Matlab仿真内容点击👇

智能优化算法       神经网络预测       雷达通信      无线传感器        电力系统

信号处理              图像处理               路径规划       元胞自动机        无人机

⛄ 内容介绍

自相关基音检测算法是语音信号处理的关键技术,算法的效率直接影响语音信号实时处理的质量。在对自相关基音检测算法基本原理进行分析的基础上,设计了Matlab算法实现方案,通过对一段具体语音时域信号采样值进行滤波、分帧、求短时自相关函数,得到了浊音语音的基音周期。

⛄ 代码

clcclear allclose allwaveFile = 'qinghua.wav ';[y, fs] =audioread(waveFile);time=(1:length(y))/fs;frameSize=floor(40*fs/1000);     %帧长startIndex=round(15000);         %起始序号endIndex=startIndex+frameSize-1; %结束序号frame = y(startIndex:endIndex);  %取出该帧frameSize=length(frame);frame2=frame.*hamming(length(frame));  % 加hamming窗rwy = rceps(frame2);                   % 求倒谱ylen=length(rwy);cepstrum=rwy(1:ylen/2); %基音检测LF=floor(fs/500);     %设置基音搜索的范围HF=floor(fs/70);cn=cepstrum(LF:HF);[mx_cep ind]=max(cn); %设置门限,找到峰值位置if mx_cep > 0.08 & ind >LF    a= fs/(LF+ind);else    a=0;endfigure(1);plot(time, y); title(waveFile); axis tightylim=get(gca, 'ylim');line([time(startIndex), time(startIndex)], ylim, 'color', 'r');line([time(endIndex), time(endIndex)], ylim, 'color', 'r');title('语音波形');figure(2);subplot(2,1,1);plot(frame);title('取出帧的波形');subplot(2,1,2);plot(cepstrum);title('倒谱图');[x,sr]=audioread('qinghua.wav');  meen=mean(x);x= x - meen;updRate=floor(20*sr/1000);          %每20ms更新fRate=floor(40*sr/1000);            %40ms一帧n_samples=length(x);nFrames=floor(n_samples/updRate)-1; %帧数k=1;pitch=zeros(1,nFrames);f0=zeros(1,nFrames);LF=floor(sr/500);HF=floor(sr/70);m=1;avgF0=0;for t=1:nFrames        yin=x(k:k+fRate-1);        cn1=rceps(yin);        cn=cn1(LF:HF);        [mx_cep ind]=max(cn);        if mx_cep > 0.08 & ind >LF           a= sr/(LF+ind);        else            a=0;        end         f0(t)=a;        if t>2 & nFrames>3   %中值滤波对基音轨迹图进行平滑           z=f0(t-2:t);           md=median(z);           pitch(t-2)=md;           if md > 0             avgF0=avgF0+md;             m=m+1;           end        else            if nFrames<=3            pitch(t)=a;            avgF0=avgF0+a;            m=m+1;           end         end     k=k+updRate;endfigure(3)subplot(211);plot((1:length(x))/sr, x);ylabel('幅度');xlabel('时间');subplot(212);xt=1:nFrames;xt=20*xt;plot(xt,pitch)xlim([0,3]);axis([xt(1) xt(nFrames) 0 max(pitch)+50]);ylabel('基音频率/HZ');xlabel('时间');

⛄ 运行结果

⛄ 参考文献

[1] 郭武, 王仁华, 戴礼荣. 基于基音周期与清浊音信息的梅尔倒谱参数[J]. 数据采集与处理, 2007, 22(2):5.

[2] 宋建华, 朱勇, 徐军,等. 基于Matlab的一种基音周期检测算法[J]. 信息技术, 2009(1):5.

[3] 段继鹏, 李春泉, 熊殷. 基于MATLAB的语音信号基音周期检测的实现[J]. 中国科技博览, 2008(21):2.

⛳️ 代码获取关注我

❤️部分理论引用网络文献,若有侵权联系博主删除
❤️ 关注我领取海量matlab电子书和数学建模资料


相关文章
|
3天前
|
数据安全/隐私保护
地震波功率谱密度函数、功率谱密度曲线,反应谱转功率谱,matlab代码
地震波格式转换、时程转换、峰值调整、规范反应谱、计算反应谱、计算持时、生成人工波、时频域转换、数据滤波、基线校正、Arias截波、傅里叶变换、耐震时程曲线、脉冲波合成与提取、三联反应谱、地震动参数、延性反应谱、地震波缩尺、功率谱密度
|
3天前
|
数据安全/隐私保护
耐震时程曲线,matlab代码,自定义反应谱与地震波,优化源代码,地震波耐震时程曲线
地震波格式转换、时程转换、峰值调整、规范反应谱、计算反应谱、计算持时、生成人工波、时频域转换、数据滤波、基线校正、Arias截波、傅里叶变换、耐震时程曲线、脉冲波合成与提取、三联反应谱、地震动参数、延性反应谱、地震波缩尺、功率谱密度
|
3天前
|
数据安全/隐私保护
matlab批量计算地震加速度反应谱、速度谱、位移谱,伪速度谱、伪加速度谱;反应谱对比图
地震波格式转换、时程转换、峰值调整、规范反应谱、计算反应谱、计算持时、生成人工波、时频域转换、数据滤波、基线校正、Arias截波、傅里叶变换、耐震时程曲线、脉冲波合成与提取、三联反应谱、地震动参数、延性反应谱、地震波缩尺、功率谱密度
基于混合整数规划的微网储能电池容量规划(matlab代码)
基于混合整数规划的微网储能电池容量规划(matlab代码)
|
3天前
|
算法 调度
面向配电网韧性提升的移动储能预布局与动态调度策略(matlab代码)
面向配电网韧性提升的移动储能预布局与动态调度策略(matlab代码)
|
4天前
|
算法 调度
含多微网租赁共享储能的配电网博弈优化调度(含matlab代码)
含多微网租赁共享储能的配电网博弈优化调度(含matlab代码)
|
4天前
|
运维 算法
基于改进遗传算法的配电网故障定位(matlab代码)
基于改进遗传算法的配电网故障定位(matlab代码)
|
3天前
|
Serverless
基于Logistic函数的负荷需求响应(matlab代码)
基于Logistic函数的负荷需求响应(matlab代码)
|
4天前
|
供应链 算法
基于分布式优化的多产消者非合作博弈能量共享(Matlab代码)
基于分布式优化的多产消者非合作博弈能量共享(Matlab代码)
|
4天前
|
算法 调度
基于多目标粒子群算法冷热电联供综合能源系统运行优化(matlab代码)
基于多目标粒子群算法冷热电联供综合能源系统运行优化(matlab代码)

热门文章

最新文章