0
$\begingroup$

I have been experimenting with cross-correlation function to verify the presence of speech in a recorded file wrt a source file. I tried the following in Matlab:

source = '/home/skrowten-hermit/Programs/male_8k.wav'
silenced = '/home/skrowten-hermit/Programs/male_8k_silence.wav'
halved = '/home/skrowten-hermit/Programs/male_8k_half.wav'
attenuated = '/home/skrowten-hermit/Programs/male_8k_attenuated.wav'
[y1, fs1] = audioread(source)
[y2, fs2] = audioread(halved)
[y3, fs3] = audioread(silenced)
[y4, fs4] = audioread(attenuated)
plot(y1)
plot(y2)
plot(y3)
plot(y4)


%% autocorrelation function wrt source

% calculate autocorrelation
[Rx1, lags1] = xcorr(y1, 'coeff')
tau1 = lags1/fs1

% plot the signal autocorrelation function
figure(6)
plot(tau1, Rx1, 'r')


%% crosscorrelation function wrt source

% calculate correlation and time axis
[Rx2, lags2] = xcorr(y1, y2, 'none')
tau2 = lags2/fs2
[Rx3, lags3] = xcorr(y1, y3, 'none')
tau3 = lags3/fs3
[Rx4, lags4] = xcorr(y1, y4, 'none')
tau4 = lags4/fs4

% plot the signal correlation function
figure(6)
plot(tau2, Rx2, 'r')
figure(6)
plot(tau3, Rx3, 'r')
figure(6)
plot(tau4, Rx4, 'r')

[pRx1, idx1] = max(Rx1)
[pRx2, idx2] = max(Rx2)
[pRx3, idx3] = max(Rx3)
[pRx4, idx4] = max(Rx4)

fprintf('Peak value of Rx1 is %f at %f', pRx1, lags(idx1))
fprintf('Peak value of Rx2 is %f at %f', pRx2, lags(idx2))
fprintf('Peak value of Rx3 is %f at %f', pRx3, lags(idx3))
fprintf('Peak value of Rx4 is %f at %f', pRx4, lags(idx4))

The following are the waveforms of my input files as generated by Matlab (male_8k.wav is my source from which others are generated or recorded):

male_8k.wav male_8k_half.wav male_8k_silence.wav male_8k_attenuated.wav

The auto-correlation generates the following output:

ACF of male_8k.wav

The cross-correlation of the source with the other targets generates the following outputs:

CCF of male_8k.wav with male_8k_half.wav CCF of male_8k.wav with male_8k_silence.wav CCF of male_8k.wav with male_8k_attenuated.wav

The output of peak value and the indices (location) of the peaks are as follows:

Peak value of Rx1 is 1.000000 at 0.000000
Peak value of Rx2 is 10.634055 at 0.000000
Peak value of Rx3 is 0.000905 at -21325.000000
Peak value of Rx4 is 48.637631 at 7516.000000

Since I intend to use the above with recorded files over a network by playing the source file (male_8k.wav) at the transmitter (Tx) end and record at a reciever (Rx) in order to verify if there is some speech detected at Rx and calculate the delay (in ms), I would like to quantify them as success or failure for verification and convert the indices (i.e., the time sample) into a value in ms. I understand that the result (i.e., the peak value) could never be 1 as in ACF, but is it possible to fix a threshold for peak and convert sample number index in such a way that:

  1. I could distinguish between silence and some speech data (attenuated is fine - just need to check if data samples exist at Rx).
  2. I could determine there is a delay of d ms at Rx.

The output values of peaks reading 10.634055 for half the speech data samples and 48.637631 for attenuated speech data samples left me a bit confused. How can I do this effectively/efficiently?

Improve this question
$\endgroup$
4
  • 1
    $\begingroup$ Consider using the Wiener-Hopf equations to measure the delay which will also reveal if you are dealing with echo conditions. I detail that approach at this post: dsp.stackexchange.com/questions/63141/… $\endgroup$
    – Dan Boschen
    Sep 18, 2020 at 15:12
  • $\begingroup$ @DanBoschen This is definitely interesting, but the computed equalizer that specifies the relationship between the source and the sink is to be used to convert the source sequence to a desired sequence right? How can it measure similarity or do a simple verification? Please help me understand. $\endgroup$
    – skrowten_hermit
    Sep 21, 2020 at 8:56
  • 1
    $\begingroup$ I show that specifically here where the delay is computed from the determined channel coefficients (solving for the channel instead of the equalizer by swapping input and output): dsp.stackexchange.com/questions/63141/… $\endgroup$
    – Dan Boschen
    Sep 21, 2020 at 14:03
  • $\begingroup$ Okay. Got it. My tx is going to be y1 and rx would be y2, y3 or y4. How do I choose ntaps, given that I have let's say a sample size of 72000 in tx and rx? And once I get the coefficients, how do I get the delay? $\endgroup$
    – skrowten_hermit
    Oct 6, 2020 at 9:36

0

Reset to default

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.