annotate userProgramsRM/testDPOAE.m @ 38:c2204b18f4a2 tip

End nov big change
author Ray Meddis <rmeddis@essex.ac.uk>
date Mon, 28 Nov 2011 13:34:28 +0000
parents
children
rev   line source
rmeddis@38 1 % testDPOAE
rmeddis@38 2
rmeddis@38 3 addpath (['..' filesep 'testPrograms'])
rmeddis@38 4
rmeddis@38 5 leveldB=60;
rmeddis@38 6 f1=3000;
rmeddis@38 7 frequencyDiffs=20:20:1000;
rmeddis@38 8 result=[];
rmeddis@38 9 frequenciesSoFar=[];
rmeddis@38 10 for f2=f1+frequencyDiffs
rmeddis@38 11 [frequencies fft_ampdB]=testDPOAE (leveldB, [f1 f2]);
rmeddis@38 12 dpFreq=2*f1-f2;
rmeddis@38 13 [a idx]=min((frequencies-dpFreq).^2);
rmeddis@38 14 result=[result fft_ampdB(idx)];
rmeddis@38 15 frequenciesSoFar=[frequenciesSoFar dpFreq];
rmeddis@38 16 figure(4), plot(frequenciesSoFar, result)
rmeddis@38 17 title(['F1= ' num2str(f1) ' F2= ' ...
rmeddis@38 18 num2str(f1+ [min(frequencyDiffs) max(frequencyDiffs)])...
rmeddis@38 19 ' leveldB= ' num2str(leveldB)])
rmeddis@38 20 xlabel('DP (2f1- f2) frequency'), ylim([0 100])
rmeddis@38 21 end
rmeddis@38 22
rmeddis@38 23 grid on
rmeddis@38 24
rmeddis@38 25 disp(result)