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)