用于计算MFCC的MATLAB代码

Question

我有一个问题,如果没关系.我最近在寻找计算MFCC的算法.我找到了一个很好的教程,而不是代码,所以我试着自己编写代码.我仍然觉得我错过了一件事.在下面的代码中,我对信号进行FFT,计算归一化功率,使用三角形形状对信号进行滤波,并最终求和每个组的能量以获得MFCC.

function output = mfcc(x,M,fbegin,fs)
    MF = @(f) 2595.*log10(1 + f./700);
    invMF = @(m) 700.*(10.^(m/2595)-1);

    M = M+2; % number of triangular filers
    mm = linspace(MF(fbegin),MF(fs/2),M); % equal space in mel-frequency
    ff = invMF(mm); % convert mel-frequencies into frequency

    X = fft(x);
    N = length(X); % length of a short time window
    N2 = max([floor(N+1)/2 floor(N/2)+1]); %
    P = abs(X(1:N2,:)).^2./N; % NoFr no. of periodograms
    mfccShapes = triangularFilterShape(ff,N,fs); %

    output = log(mfccShapes'*P);
end

function [out,k] = triangularFilterShape(f,N,fs)
    N2 = max([floor(N+1)/2 floor(N/2)+1]);
    M = length(f);
    k = linspace(0,fs/2,N2);
    out = zeros(N2,M-2);
    for m=2:M-1
        I = k >= f(m-1) & k <= f(m);
        J = k >= f(m) & k <= f(m+1);
        out(I,m-1) = (k(I) - f(m-1))./(f(m) - f(m-1));
        out(J,m-1) = (f(m+1) - k(J))./(f(m+1) - f(m));
    end
end

如果我犯了错误,有人可以确认这是对的,也可以指示我>我用简单的纯音测试它,在我看来,它给了我合理的答案.

任何帮助非常感谢:)

PS.我正在研究如何应用矢量化的Cosinus变换.看起来我需要一个MxM的变换系数矩阵,但我找不到任何可以解释如何做的源.

Answer 1

您可以通过将结果与其他类似的实现进行比较来自行测试,您可以在这里 找到完全可配置的matlab工具箱.MFCC甚至是将MFCC反转回时间信号的功能,这对于测试目的非常方便:

melfcc.m - 从声音波形计算PLP和MFCC的主要功能,支持许多选项.

invmelfcc.m - 主要功能,用于从倒谱系数反转回频谱图和(噪声激励)波形,选项与melfcc精确匹配(反转处理).

页面本身有很多关于包的用法的信息.