Комплексный кепстр — основная оценка частоты

В этом примере показано, как оценить основную частоту докладчика с помощью комплексного кепстра. Пример также оценивает основную частоту с помощью метода пересечения нулем и сравнивает результаты.

Загрузите речевой сигнал. Запись имеет женщину, говорящую "MATLAB". Частота дискретизации составляет 7 418 Гц. Следующий код загружает речевую форму волны, mtlb, и частота дискретизации, Fs, в рабочую область MATLAB®.

load mtlb

Используйте спектрограмму, чтобы идентифицировать речевой сегмент для анализа.

segmentlen = 100;
noverlap = 90;
NFFT = 128;

spectrogram(mtlb,segmentlen,noverlap,NFFT,Fs,'yaxis')

Извлеките сегмент от 0,1 до 0,25 секунд для анализа. Извлеченный сегмент соответствует примерно первому гласному,/æ/, в "MATLAB".

dt = 1/Fs;
I0 = round(0.1/dt);
Iend = round(0.25/dt);
x = mtlb(I0:Iend);

Получите комплексный кепстр.

c = cceps(x);

Выберите область значений времени между 2 и 10 мс, соответствуя частотному диапазону приблизительно 100 - 500 Гц. Идентифицируйте самый высокий пик кепстра в выбранной области значений. Найдите частоту, соответствующую пику. Используйте пик в качестве оценки основной частоты.

t = 0:dt:length(x)*dt-dt;

trng = t(t>=2e-3 & t<=10e-3);
crng = c(t>=2e-3 & t<=10e-3);

[~,I] = max(crng);

fprintf('Complex cepstrum F0 estimate is %3.2f Hz.\n',1/trng(I))

Complex cepstrum F0 estimate is 239.29 Hz.

Постройте кепстр в выбранное время, располагаются и накладывают пик.

plot(trng*1e3,crng)
xlabel('ms')

hold on
plot(trng(I)*1e3,crng(I),'o')
hold off

Используйте zerocrossrate функция на отфильтрованной lowpass и исправленной форме гласного, чтобы оценить основную частоту.

[b0,a0] = butter(2,325/(Fs/2));
xin = abs(x);
xin = filter(b0,a0,xin);
xin = xin-mean(xin);
zc = zerocrossrate(xin);
F0 = 0.5*Fs*zc;
fprintf('Zero-crossing F0 estimate is %3.2f Hz.\n',F0)

Zero-crossing F0 estimate is 234.94 Hz.

Смотрите также

cceps | icceps | rceps