Hilbert Envelope Extraction and Cross-Correlation of two-microphone signals for estimating sound arrival direction

Project info

• Field of studyAcoustics and Audio
• LevelDoctoral dissertation (S3)
• Method typeSignal Processing
• Estimated build time± 7 working days
• Toolshilbert, xcorr (Signal Processing Toolbox)
• PriceUS$34 - US$219

Research goal

Extract the envelope with the Hilbert transform and estimate the time lag between signals in two-microphone signals for estimating sound arrival direction using cross-correlation.

Variables

Example result

• 1.01Envelope mean
• 0.541Modulation depth (est.)
• 0.394Modulation depth (true)
• 55Estimated delay (ms)
• 0.977Peak correlation
• 71Carrier frequency (Hz)

The Hilbert envelope of First microphone signal gives a modulation depth of 0.541 (true 0.394); cross-correlation with Time difference of arrival peaks at 0.977 for a 55 ms delay, and the carrier is estimated at 71.0 Hz.

Computed from illustrative synthetic data, not from a real study.

About Envelope and Cross-Correlation (Hilbert)

The Hilbert transform produces the analytic signal to obtain the envelope (instantaneous amplitude), and cross-correlation measures the delay between two signals. It is used in vibration analysis, modulated signals, and sensor delay estimation.

Steps

1. 1Prepare the signals and sampling frequency
2. 2Compute the analytic signal with hilbert
3. 3Take the envelope from its magnitude
4. 4Compute the cross-correlation of two signals with xcorr
5. 5Find the lag at the correlation peak to estimate delay
6. 6Plot the envelope and the correlation function

Assumptions to check

• The signal is narrowband so the envelope is meaningful
• The delay between signals is constant over the record
• Noise is not so large that it hides the correlation peak

Reading and reporting the results

Report the envelope (for example peak value or trend) and the estimated delay in seconds with its error against the true value if known. State the delay resolution of 1/fs.

Example MATLAB code

A short starting script with this project's variable names. The delivered version is adapted to your data and fully commented.

rng(5);
fs = 1000;
time_difference_of_arrival = (0:1/fs:2-1/fs)';
env_true = 1 + 0.5*sin(2*pi*2*time_difference_of_arrival);
first_microphone_signal = env_true .* sin(2*pi*80*time_difference_of_arrival);          % amplitude-modulated signal

analytic = hilbert(first_microphone_signal);
env = abs(analytic);

delay_s = 0.037;
d = round(delay_s * fs);
y2 = [zeros(d,1); first_microphone_signal(1:end-d)] + 0.3*randn(size(first_microphone_signal));
[c, lags] = xcorr(y2, first_microphone_signal + 0.3*randn(size(first_microphone_signal)), 200, 'coeff');
[cmax, imax] = max(c);
fprintf('Estimated delay: %.3f s (true %.3f s), peak corr %.2f\n', lags(imax)/fs, delay_s, cmax);
fprintf('Envelope RMSE: %.3f\n', sqrt(mean((env(100:end-100) - env_true(100:end-100)).^2)));

figure;
subplot(2,1,1);
plot(time_difference_of_arrival, first_microphone_signal, time_difference_of_arrival, env, 'r', 'LineWidth', 1.2);
legend('Signal', 'Envelope'); xlabel('Time (s)');
subplot(2,1,2);
plot(lags/fs, c); xlabel('Lag (s)'); ylabel('Correlation');

Packages and prices

Choose how much help you need. Each higher package includes everything in the lower ones. Prices exclude 11% VAT.

• Payment: 50% down payment before work starts, balance after the deliverables are handed over
• Tax: Prices exclude 11% VAT
• Data and academic integrity: Results are computed from your own research data; you write and stand behind your manuscript, we help with the analysis and your understanding of it
• Revisions: Up to 2 revisions are included in the price (Package 3: 3 rounds)

Frequently asked questions