comparison multithreshold 1.46/testOME.m @ 22:45f28c49461e master

removing duplicate changes

 
 figure(2), clf, subplot(2,1,1)
 set(2,'position',[5   349   268   327])
 semilogx(HuberFrequencies, 20*log10(HuberDisplacementAt80dBSPL/1e-10),...
     'ko', 'MarkerFaceColor','k', 'Marker','o', 'markerSize',6)
-% Generate test stimulus .................................................................
+hold on
 
-%% independent test using discrete frequencies
+%% Generate test stimulus .................................................................
+
+% independent test using discrete frequencies
 peakResponses=[];
 peakTMpressure=[];
 frequencies=[200 400 HuberFrequencies 10000];
 for toneFrequency=frequencies
     inputSignal=amp*sin(2*pi*toneFrequency*time);

     AN_spikesOrProbability='probability';
     % switch off AR & MOC (Huber's patients were deaf)
     paramChanges{1}='OMEParams.rateToAttenuationFactorProb=0;';
     paramChanges{2}='DRNLParams.rateToAttenuationFactorProb = 0;';
 
-    global OMEoutput  OMEextEarPressure TMoutput ARAttenuation
+    global OMEoutput  OMEextEarPressure TMoutput ARattenuation
     % BF is irrelevant
     MAP1_14(inputSignal, sampleRate, -1, ...
         MAPparamsName, AN_spikesOrProbability, paramChanges);
 
     peakDisplacement=max(OMEoutput(floor(end/2):end));
     peakResponses=[peakResponses peakDisplacement];
 
     peakTMpressure=[peakTMpressure max(OMEextEarPressure)];
-    disp([' greatest AR attenuation:   ' num2str(min(ARAttenuation))])
+    disp([' AR attenuation (dB):   ' num2str(20*log10(min(ARattenuation)))])
 end
 
 %% Report
 disp('frequency displacement(m)')
 % disp(num2str([frequencies' peakResponses']))
 fprintf('%6.0f \t%10.3e\n',[frequencies' peakResponses']')
 
 % stapes peak displacement
 figure(2), subplot(2,1,1), hold on
 semilogx(frequencies, 20*log10(peakResponses/1e-10), 'r', 'linewidth',4)
+set(gca,'xScale','log')
 % ylim([1e-11 1e-8])
 xlim([100 10000]), ylim([0 30])
 grid on
 title(['stapes at ' num2str(leveldBSPL) ' (NB deaf)'])
 ylabel('disp: dB re 1e-10m')