close all; clc;
load('whitenoize.mat'); %generated white noise
load('whitenoiserec.mat'); %recorded white nosie
x = orgwhite(20000:50000); % extract the white noise part only
y = recwhite(20000:50000)'; % extract the corresponding response to white noise part
h2 = ifft(fft(y)./fft(x)); %to get impulse response
plot(abs(h2));
Best Answer