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