【短时平均过零率】语音信号短时幅值和短时过零率【含Matlab源码 4942期】

⛄一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【短时平均过零率】基于matlab语音信号短时幅值和短时过零率【含Matlab源码 4942期】
点击上面蓝色字体，直接付费下载，即可。

获取代码方式2：
付费专栏Matlab语音处理（初级版）

备注：
点击上面蓝色字体付费专栏Matlab语音处理（初级版），扫描上面二维码，付费29.9元订阅海神之光博客付费专栏Matlab语音处理（初级版），凭支付凭证，私信博主，可免费获得1份本博客上传CSDN资源代码（有效期为订阅日起，三天内有效）；
点击CSDN资源下载链接：1份本博客上传CSDN资源代码

⛄二、部分源代码

% 语音信号时域波形、短时幅值和短时过零率%%%%%%%%%%%%%%%%
%%语音信号变量y_sound; 短时幅值m1_y;短时能量m2_y;过零率zero_y%%%%

clear;
clc;
close all;
% [y_sound,Fs]=audioread( ‘C:\Users\lenovo\Desktop\大三春季学期课程与项目\语音信号处理-黄青华\实验\实验一\my_voice\t.wav’);%打开声音信号
[FileName,PathName]=uigetfile('.wav’); % uigetfile动态选择文件
dir=[PathName,FileName]; % 实质就是文件的名字
[y_sound,Fs]=audioread(dir); % 读取特定路径下的语音文件
y_sound=y_sound(:,1); % 选择语音的单通道进行画图，如果不是那么画的是双通道的
y_length=length(y_sound); % 语音信号数字化之后的长度
T=0:(1/Fs):(y_length-1)/Fs; % 声音信号的时长
Ny=10Fs/1000; %短时帧长，10ms
if Ny/2>fix(Ny/2) % 将Ny变成偶数
Ny=Ny-1;
else
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if y_length>Ny
frame_all=floor((y_length-(Ny/2))/(Ny/2)); %frame_moving equal to a half of frame_length, get number of all frame,to use “enframe” too.
y2_sound=y_sound(1:(frame_all+1)*(Ny/2));
y3=reshape(y2_sound,(Ny/2),[]);
y4=[y3(:,1:frame_all);y3(:,2:(frame_all+1))];
else
frame_all=1;
y4=y_sound;
end
%%%%%%%加方窗
for i=1:(frame_all)
m1_y(i)=sum(abs(y4(:,i))); %计算短时幅值
m2_y(i)=sum(y4(:,i).*y4(:,i)); %计算短时能量

zero_y(i)=sum(abs(sign(y4(2:Ny,i))-sign(y4(1:(Ny-1),i))))/2;   %计算短时过零率

end
%%%%%%%加哈明窗
win_h=hamming(Ny); %win_y=y_sound.*win_h;
for i=1:(frame_all)
mh1_y(i)=sum(abs(y4(:,i).*win_h)); %计算短时幅值

mh2_y(i)=sum((y4(:,i).win_h).(y4(:,i).*win_h)); %计算短时能量
end

⛄三、运行结果

⛄四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1]韩纪庆,张磊,郑铁然.语音信号处理（第3版）[M].清华大学出版社，2019.
[2]柳若边.深度学习:语音识别技术实践[M].清华大学出版社，2019.
[3]宋云飞,姜占才,魏中华.基于MATLAB GUI的语音处理界面设计[J].科技信息. 2013,(02)

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除