Mercurial > hg > map
diff 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 |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/userProgramsRM/testDPOAE.m Mon Nov 28 13:34:28 2011 +0000 @@ -0,0 +1,25 @@ +% testDPOAE + +addpath (['..' filesep 'testPrograms']) + +leveldB=60; +f1=3000; +frequencyDiffs=20:20:1000; +result=[]; +frequenciesSoFar=[]; +for f2=f1+frequencyDiffs + [frequencies fft_ampdB]=testDPOAE (leveldB, [f1 f2]); + dpFreq=2*f1-f2; + [a idx]=min((frequencies-dpFreq).^2); + result=[result fft_ampdB(idx)]; + frequenciesSoFar=[frequenciesSoFar dpFreq]; + figure(4), plot(frequenciesSoFar, result) + title(['F1= ' num2str(f1) ' F2= ' ... + num2str(f1+ [min(frequencyDiffs) max(frequencyDiffs)])... + ' leveldB= ' num2str(leveldB)]) + xlabel('DP (2f1- f2) frequency'), ylim([0 100]) +end + +grid on + +disp(result)