comparison 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
comparison
equal deleted inserted replaced
37:771a643d5c29 38:c2204b18f4a2
1 % testDPOAE
2
3 addpath (['..' filesep 'testPrograms'])
4
5 leveldB=60;
6 f1=3000;
7 frequencyDiffs=20:20:1000;
8 result=[];
9 frequenciesSoFar=[];
10 for f2=f1+frequencyDiffs
11 [frequencies fft_ampdB]=testDPOAE (leveldB, [f1 f2]);
12 dpFreq=2*f1-f2;
13 [a idx]=min((frequencies-dpFreq).^2);
14 result=[result fft_ampdB(idx)];
15 frequenciesSoFar=[frequenciesSoFar dpFreq];
16 figure(4), plot(frequenciesSoFar, result)
17 title(['F1= ' num2str(f1) ' F2= ' ...
18 num2str(f1+ [min(frequencyDiffs) max(frequencyDiffs)])...
19 ' leveldB= ' num2str(leveldB)])
20 xlabel('DP (2f1- f2) frequency'), ylim([0 100])
21 end
22
23 grid on
24
25 disp(result)