%**************************************************************************
%逐浈帧求取语音信号的共振峰
%所得s为一个矩阵,存储了每帧语音信号所对应的共振峰值
%**************************************************************************
function [s,s_bw,s_amp]=pick_peak(sound,samplerate)
%sex:1--男声,2--女声
% sex=1;
%对语音信号进行预处理,得到分帧后的信号
Sw=pre(sound,samplerate);
[frame_number,point]=size(Sw);
%逐帧修正,EST的初始值为每帧初始值的平均
for t=1:1:frame_number
[ppp,bw,amp]=find_peak(Sw(t,:),samplerate,1);
pp(t,:)=[ppp,zeros(1,4-length(ppp))];
end
for t=1:1:4
a=find(pp(:,t)==0);
if isempty(a)
l=0;
else
l=length(a);
end
EST(t)=sum(pp(:,t)/(frame_number-l));
end
% if sex==1
% EST(1)=320;
% EST(2)=1440;
% EST(3)=2760;
% EST(4)=3200;
% elseif sex==2
% EST(1)=480;
% EST(2)=1760;
% EST(3)=3200;
% EST(4)=3520;
% end
% EST=mean(p,1);
%s用来记录当前帧的共振峰估计值
s=zeros(frame_number,4);
s_bw=zeros(frame_number,4);
s_amp=zeros(frame_number,4);
%开始逐帧处理
for t=1:1:frame_number
r=1;
while(s(t,1)==0 || s(t,2)==0 || s(t,3)==0)
[p,bw,amp]=find_peak(Sw(t,:),samplerate,r);
if r>0.88
r=r-0.02;
else
break;
end
peak_number=length(p);
%p_mark用来标注被分配的峰值,已经分配的为0,未被分配的为1;
p_mark=ones(1,peak_number);
%e用来存放峰值与预测值的差
e=zeros(1,peak_number);
%将提取出的峰值对应的填到每个共振峰的位置上,要使得对应的共振峰值与EST值最为接近
for x=1:1:4
for y=1:1:peak_number
e(y)=abs(p(y)-EST(x));
end
[a,b]=min(e);
s(t,x)=p(b);
s_bw(t,x)=bw(b);
s_amp(t,x)=amp(b);
p_mark(b)=0;
end
%如果两个共振峰填入的峰值相同,则将该值保留在与EST值差距最小的位置,而将另外一个共振峰值除去
for x=1:1:3
for y=(x+1):1:4
if s(t,x)==s(t,y)
e(x)=abs(s(t,x)-EST(x));
e(y)=abs(s(t,y)-EST(y));
if e(x)>e(y)
s(t,x)=0;
else
s(t,y)=0;
end
end
end
end
%处理没有分配的峰值
for k=1:1:peak_number
if p_mark(k)~=0
if s(t,k)==0
s(t,k)=p(k);
p_mark(k)=0;
elseif ((amp(k)>=s_amp(t,k)*0.5) && k~=4)
if (s(t,k+1)==0)
s(t,k+1)=s(t,k);
s_bw(t,k+1)=s_bw(t,k);
s_amp(t,k+1)=s_amp(t,k);
s(t,k)=p(k);
s_bw(t,k)=bw(k);
s_amp(t,k)=amp(k);
p_mark(k)=0;
end
elseif (amp(k)>=s_amp(t,k)*0.5) && (k~=1)
if (s(t,k-1)==0)
s(t,k-1)=s(t,k);
s_bw(t,k-1)=s_bw(t,k);
s_amp(t,k-1)=s_amp(t,k);
s(t,k)=p(k);
s_bw(t,k)=bw(k);
s_amp(t,k)=amp(k);
p_mark(k)=0;
end
end
end
end
for k=1:1:4
if (s(t,k)~=0)
EST(k)=s(t,k);
end
end
end
end
IT狂飙
- 粉丝: 4842
- 资源: 2650
最新资源
- 基于BS模式的冷链物流系统pf-springboot毕业项目,适合计算机毕-设、实训项目、大作业学习.zip
- 基于simulink实现的110kV继电保护设计实现
- 这个程序是一个简单的git代理切换器,它的主要功能是帮助用户在需要代理的网络环境和不需要代理的网络环境之间快速切换
- java javacv jt1078 rtp流转flv
- iClient for Classic加载wmts山东天地图示例demo
- lvgl-v7版本代码示例
- 基于模型预测控制的无人驾驶车辆轨迹跟踪问题研究,matlab 附有MATLAB程序与详细的建模过程,研究车辆转向的同学可以作为参考
- 基于SpringBoot+vue的名城小区物业管理系统(含报告).zip
- Xmanager Power Suite 8.0.0005
- 阅后即焚平台系统PHP源码
- C#读写西门子PLC.OPC.数据库.Socket 1、PLC数据通信读写; 2、联合OPC; 3、联合Socket; 4、联合数据库;
- 德普微一级代理 DP023N10TGN TOLL DPMOS N-MOSFET 100V 300A 1.75mΩ
- 在线教育网校系统社区版+带安装说明
- 基于Spring+vue的智慧生活商城系统.zip
- 高可用k8s集群离线部署(五)
- 基于Java的考试系统代码pf-springboot毕业项目,适合计算机毕-设、实训项目、大作业学习.rar
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈