短时时域处理（一）

最新推荐文章于 2022-05-13 10:12:25 发布

齐燕博

最新推荐文章于 2022-05-13 10:12:25 发布

阅读量346

点赞数

CC 4.0 BY-SA版权

分类专栏： matlab 语音

本文链接：https://blog.csdn.net/lovekobe1997/article/details/119026166

matlab 同时被 2 个专栏收录

17 篇文章

订阅专栏

语音

13 篇文章

订阅专栏

短时能量

根据公式自己敲的短时能量代码：

clc;
clear;
[x,fs]=audioread('C3_2_y.wav');
inc=100;
win=hamming(200);
fz=enframe(x,win,inc);
fn=size(fz,1);
for i=1:fn
    u=fz(i,:);
    u2=u.*u;
    power(i)=sum(u2);
end

在这里插入图片描述
短时平均幅度:

$[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-KQAharBd-1627008170073)(C:\Users\Qyb\AppData\Roaming\Typora\typora-user-images\image-20210722161120119.png)]$

%短时平均幅度计算函数
function para=STMn(x,win,inc)
    X=enframe(x,win,inc)';     % 分帧
    fn=size(X,2);              % 求出帧数
    for i=1 : fn
        u=X(:,i);              % 取出一帧
        para(i)=sum(abs(u))/100;         % 对一帧累加求和
    end
 end

$[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-Cfs5TB5k-1627008170073)(C:\Users\Qyb\AppData\Roaming\Typora\typora-user-images\image-20210722192815761.png)]$
短时平均过零率:

在这里插入图片描述

clc;
clear;
[x,fs]=audioread('C3_2_y.wav');%读取音频
inc=100;%帧移
win=hamming(200);%汉明窗和窗长
fz=enframe(x,win,inc);%分帧
fn=size(fz,1);%帧数
zc=size(fz,2);%帧长
para=zeros(1,fn);
for i=1:fn
    u=fz(i,:);%取一帧
    for j=1:(zc-1)   %每一帧的数据，用于比较
        if u(j)*u(j+1)<0%如果相邻两个数相乘小于0，则存放过零率的数组在对应帧加一
            para(i)=para(i)+1;
        end
    end
end

在这里插入图片描述