查看: 92|回复: 1|关注: 0
打印 上一主题 下一主题

快乐扑克3对子走势图: [未答复] 提取mfcc后幅值非常大

[复制链接]

新手

5 麦片

财富积分


快乐扑克3如何跟对子 www.vqxik.tw 050


4

主题

8

帖子

0

最佳答案
跳转到指定楼层
1#
发表于 2019-11-24 11:58:17 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
我参照网上的程序提取的mfcc绘制图形出来后发现幅值比他的大很多很多,不知道哪里出了问题。第一个图片是我的,第二个是网上别人的。
function Mfcc
clc;close all;
tic
[x,fs] = audioread('E:\\1matlab\\filter.wav');
bank=melbankm(24,256,fs,0,0.4,'t');%Mel滤波器的阶数为24,fft变换的长度为256,采样频率为8000Hz
%归一化mel滤波器组系数
bank=full(bank);
bank=bank/max(bank(:));
dctcoef(12,24) = 0;
for k=1:12 %归一化mel滤波器组系数
    n=0:23;
    dctcoef(k,:)=cos((2*n+1)*k*pi/(2*24));
end
w=1+6*sin(pi*(1:12)./12);%归一化倒谱提升窗口
w=w/max(w);%预加重滤波器
xx=double(x);
xx=filter((1-0.9375),1,xx);%语音信号分帧
xx=enframe(xx,256,80);%对x 256点分为一帧
%计算每帧的MFCC参数
m(size(xx,1),12) = 0;
for i=1:size(xx,1)
    y=xx(i,:);
    s=y'.*hamming(256);
    t=abs(fft(s));%fft快速傅立叶变换
    t=t.^2;
    c1=dctcoef*log(bank*t(1:129));
    c2=c1.*w';
    m(i,:)=c2';
end
%求取一阶差分系数
dtm=zeros(size(m));
for i=3:size(m,1)-2
    dtm(i,:)=-2*m(i-2,:)-m(i-1,:)+m(i+1,:)+2*m(i+2,:);
end
dtm=dtm/3;
%求取二阶差分系数
dtmm=zeros(size(dtm));
for i=3:size(dtm,1)-2
    dtmm(i,:)=-2*dtm(i-2,:)-dtm(i-1,:)+dtm(i+1,:)+2*dtm(i+2,:);
end
dtmm=dtmm/3;
%合并mfcc参数和一阶、二阶差分mfcc参数
ccc=[m dtm dtmm];
%去除首尾两帧,因为这两帧的一阶差分参数为0
ccc=ccc(3:size(m,1)-2,:);
toc
subplot(211)
ccc_1=ccc(:,1);
plot(ccc_1);title('MFCC');ylabel('幅值');
[~,w]=size(ccc);
A=size(ccc);
subplot(212)
plot([1,w],A);
xlabel('维数');
ylabel('幅值');
title('维数与幅值的关系')   end




新手

5 麦片

财富积分


050


4

主题

8

帖子

0

最佳答案
2#
 楼主| 发表于 2019-11-24 13:35:48 | 只看该作者
我的这段音频是有10分钟左右,我换了一段几秒的音频出来的mfcc幅值就比较小了,这是为什么呢
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /3 下一条

快速回复 快乐扑克3如何跟对子 返回列表