Matlab提取pcm文件的MFCC参数输出为excel表格

参考网上最常见版本的MATLAB提取MFCC参数代码做了简单的修改,修改的内容有两个:一个是读取文件改成了pcm文件,一个是输出参数保存为excel表格,可以成功运行,希望能对其他人有所帮助。参考地址:https://blog.csdn.net/jojozhangju/article/details/18678861%在提取MFCC参数之前需要加载并使用VOICEBOX工具包fs=16000;...
摘要由CSDN通过智能技术生成

参考网上最常见版本的MATLAB提取MFCC参数代码做了简单的修改,修改的内容有两个:一个是读取文件改成了pcm文件,一个是输出参数保存为excel表格,可以成功运行,希望能对其他人有所帮助。

参考地址:https://blog.csdn.net/jojozhangju/article/details/18678861


%在提取MFCC参数之前需要加载并使用VOICEBOX工具包
fs=16000;
fileId = fopen('F:\sample.pcm','r');
x = fread(fileId,inf,'int16');
%sound(x,fs);
bank=melbankm(24,256,fs,0,0.4,'t');%Mel滤波器的阶数为24,fft变换的长度为256,采样频率为16000Hz
%归一化mel滤波器组系数
bank=full(bank);
bank=bank/max(bank(:));
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参数
for i=1:size(xx,1)
y=xx(i,:);
s=y'.*hamming(256);
	t=abs(fft(s));%fft快速傅立叶变换
t=t.^2;
	c1=dctcoef*log(b
  • 0
    点赞
  • 6
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
MFCC(Mel-Frequency Cepstral Coefficient)是一种常用于语音信号处理的特征参数,可以用于语音识别、语音合成、说话人识别等应用中。MATLAB提取MFCC特征参数的步骤如下: 1. 读取音频数据并进行预处理:读取音频文件,进行预加重、分帧、加窗等预处理操作。 2. 计算功率谱密度:对每一帧信号进行傅里叶变换,计算功率谱密度。 3. 应用梅尔滤波器组:将功率谱密度通过一组梅尔滤波器组,将其转换为梅尔频率谱。 4. 取对数:对梅尔频率谱取对数。 5. 进行离散余弦变换:对取对数后的梅尔频率谱进行离散余弦变换,得到MFCC系数。 MATLAB中可以使用mfcc函数来提取MFCC特征参数,示例代码如下: ```matlab % 读取音频数据 [x, Fs] = audioread('example.wav'); % 预处理 x = preemphasis(x, Fs); % 预加重 frameLen = round(25 * Fs / 1000); % 帧长 overlap = round(15 * Fs / 1000); % 帧重叠 frames = enframe(x, frameLen, overlap); % 分帧 window = hamming(frameLen); % 窗函数 frames = frames .* repmat(window, 1, size(frames, 2)); % 加窗 % 计算MFCC特征参数 coeffs = mfcc(frames, Fs, 'NumCoeffs', 13); % 取13个MFCC系数 ``` 其中,preemphasis函数用于进行预加重,enframe函数用于将信号分帧,hamming函数用于生成汉明窗函数。mfcc函数的参数说明如下: - frames:分帧后的信号矩阵,每一列为一帧信号。 - Fs:采样率。 - 'NumCoeffs':MFCC系数个数,默认为 20。 - 'NoDCT':是否进行离散余弦变换,默认为 false,即进行离散余弦变换。 - 'LogEnergy':是否计算对数能量,默认为 true。 - 'MelSlope':梅尔滤波器组斜率,默认为 1。 - 'MelSpacing':梅尔滤波器组间隔,默认为 'Log',即对数间隔。 - 'Window':分帧所用的窗函数,默认为 'Hamming'。 - 'OverlapLength':帧重叠长度,默认为 round(2.5 * frameLen / 10)。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值