rmeddis@38: % testDPOAE
rmeddis@38: 
rmeddis@38: addpath (['..' filesep 'testPrograms'])
rmeddis@38: 
rmeddis@38: leveldB=60;
rmeddis@38: f1=3000;
rmeddis@38: frequencyDiffs=20:20:1000;
rmeddis@38: result=[];
rmeddis@38: frequenciesSoFar=[];
rmeddis@38: for f2=f1+frequencyDiffs
rmeddis@38:     [frequencies fft_ampdB]=testDPOAE (leveldB, [f1 f2]);
rmeddis@38:     dpFreq=2*f1-f2;
rmeddis@38:     [a idx]=min((frequencies-dpFreq).^2);
rmeddis@38:     result=[result fft_ampdB(idx)];
rmeddis@38:     frequenciesSoFar=[frequenciesSoFar dpFreq];
rmeddis@38:     figure(4), plot(frequenciesSoFar, result)
rmeddis@38:     title(['F1= ' num2str(f1) '  F2= ' ...
rmeddis@38:         num2str(f1+ [min(frequencyDiffs) max(frequencyDiffs)])...
rmeddis@38:         '  leveldB= ' num2str(leveldB)])
rmeddis@38:     xlabel('DP (2f1- f2) frequency'), ylim([0 100])
rmeddis@38: end
rmeddis@38: 
rmeddis@38: grid on
rmeddis@38: 
rmeddis@38: disp(result)