Mercurial > hg > map
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) |