Matlab中如何实现语音信号短时平均过零率特征提取的源码？

本文共计527个文字，预计阅读时间需要3分钟。

1. 简介

2.部分代码

- 清除所有变量：`clear all;` - 清屏：`clc;` - 关闭所有图形窗口：`close all;` - 设置路径：`filedir=[];` - 设置文件名：`filename='s.wav';` - 构成完整路径和文件名：`fle=[filedir filename];` - 读取数据文件：`[xx, Fs]=audioread(fle);` - 去除趋势：`x=detrend(xx);` - 消除去除趋势

1 简介

2 部分代码

clear all; clc; close all;
filedir=[]; % 设置路径
filename='s.wav'; % 设置文件名
fle=[filedir filename]; % 构成完整的路径和文件名
[xx,Fs]=audioread(fle); % 读入数据文件
x=detrend(xx); % 消除直流分量
wlen=200; inc=80; % 设置帧长、帧移
win=hanning(wlen); % 窗函数
N=length(x); % 求数据长度
X=enframe(x,win,inc)'; % 分帧
fn=size(X,2); % 获取帧数
zcr1=zeros(1,fn); % 初始化
for i=1:fn
z=X(:,i); % 取得一帧数据
for j=1: (wlen-1); % 在一帧内寻找过零点
if z(j)* z(j+1)< 0 % 判断是否为过零点
zcr1(i)=zcr1(i)+1; % 是过零点，记录1次
end
end
end
time=(0:N-1)/Fs; % 计算时间坐标
frameTime=frame2time(fn,wlen,inc,Fs); % 求出每帧对应的时间
% 作图
subplot 211; plot(time,x,'k'); grid;
title('语音波形');
ylabel('幅值'); xlabel(['时间/s' 10 '(a)']);
subplot 212; plot(frameTime,zcr1,'k'); grid;
title('短时平均过零率');
ylabel('幅值'); xlabel(['时间/s' 10 '(b)']);

3 仿真结果

4 参考文献

[1]卓嘎, 边巴旺堆, and 姜军. "基于短时平均能量和短时过零率的藏语语音端点检测研究." 电脑知识与技术：学术版 11(2014):4.

博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真，相关matlab代码问题可私信交流。

部分理论引用网络文献，若有侵权联系博主删除。