3,423
社区成员
发帖
与我相关
我的任务
分享
clear all;
[data,fs]=audioread('A1.wav');
xx=double(data);
xx=filter([1 -0.9375],1,xx); % 语音信号预加重
y=enframe(xx,256); %xx=enframe(xx,256,80);
guss=size(y,1); %size(xx,1)返回xx的维数
c=eye(guss);
% 计算每帧的MFCC参数
for i=1:guss
h = y(i,:);
s = h' .* hamming(256);%乘窗
t = rfft(s);
p=melbankm(2,510,44100);
%whos p;
%whos s;
for j=1:2
z=log(p(j,:)*(abs(s).^2));
c(i,:)=dct(z);
end
end
fid=fopen('Mel.txt','w+');
fprintf(fid,'%6.2f',c);
fclose(fid);