comparison userProgramsRM/test_MAP1_14RAMworks.m @ 38:c2204b18f4a2 tip

End nov big change

+function test_MAP1_14RAM
+
+global dt  DRNLoutput 
+global  DRNLParams TMoutput
+global peakOutputDisp alignedOutput toneFrequencyList levels leveldBSPL logOutput
+
+
+dbstop if error
+restorePath=path;
+addpath (['..' filesep 'MAP'],    ['..' filesep 'wavFileStore'], ...
+    ['..' filesep 'utilities'])
+figure(2), clf
+
+%% # BFlist is the BF of the filter to be assessed
+BFlist=8000;
+numChannels=1;
+
+% define probe frequencies
+numFs=6;
+lowestF=BFlist/5; 	highestF= BFlist*1.4;
+lowestF=4000; 	highestF= 12600;
+toneFrequencyList=round(logspace(log10(lowestF), log10(highestF), numFs));
+
+%%  # parameter file name
+MAPparamsName='Normal';
+
+%% # probability representation
+AN_spikesOrProbability='probability';
+
+%% # tone duration
+sampleRate= 100000;
+duration=0.0200;                 % seconds
+rampDuration=.005;              % raised cosine ramp (seconds)
+beginSilence=0.050;
+endSilence=0.050;
+
+%% # levels
+% probe details
+levels=0:10:80;
+% levels=80;
+
+%% # change model parameters
+% Parameter changes can be used to change one or more model parameters
+%  *after* the MAPparams file has been read
+
+% adjust linear path gain (previously g=100)
+% switch off all efferent effects
+paramChanges={...
+    'DRNLParams.rateToAttenuationFactorProb = 0.00; ',...
+    'OMEParams.rateToAttenuationFactorProb=0.0;', ...
+    'DRNLParams.ctBMdB = -10;'...
+    'DRNLParams.g=200;'...
+    'DRNLParams.linCFs=6500;'...
+    'DRNLParams.linBWs=2000;'...
+    };
+%     'DRNLParams.a=0;'...
+
+
+%% delare 'showMap' options to control graphical output
+showMapOptions.printModelParameters=0;   % prints all parameters
+showMapOptions.showModelOutput=1;       % plot of all stages
+showMapOptions.printFiringRates=0;      % prints stage activity levels
+showMapOptions.showACF=0;               % shows SACF (probability only)
+showMapOptions.showEfferent=0;          % tracks of AR and MOC
+showMapOptions.surfProbability=0;       % 2D plot of HSR response
+showMapOptions.surfSpikes=0;            % 2D plot of spikes histogram
+showMapOptions.ICrates=0;               % IC rates by CNtauGk
+showMapOptions.fileName=[];
+
+
+%% now vary level and frequency
+peakOutputDisp=zeros(length(levels),length(toneFrequencyList));
+peakStapesDisp=zeros(length(levels),length(toneFrequencyList));
+levelNo=0;
+for leveldBSPL=levels
+    levelNo=levelNo+1;
+disp(['level: ' num2str(leveldBSPL)])
+    freqNo=0;
+    for toneFrequency=toneFrequencyList
+        freqNo=freqNo+1;
+
+        %% Generate stimuli
+        dt=1/sampleRate;
+        time=dt: dt: duration;
+        inputSignal=sum(sin(2*pi*toneFrequency'*time), 1);
+        amp=10^(leveldBSPL/20)*28e-6;   % converts to Pascals (peak)
+        inputSignal=amp*inputSignal;
+        % apply ramps
+        if rampDuration>0.5*duration, rampDuration=duration/2; end
+        rampTime=dt:dt:rampDuration;
+        ramp=[0.5*(1+cos(2*pi*rampTime/(2*rampDuration)+pi)) ...
+            ones(1,length(time)-length(rampTime))];
+        inputSignal=inputSignal.*ramp;
+        ramp=fliplr(ramp);
+        inputSignal=inputSignal.*ramp;
+        % add silence
+        intialSilence= zeros(1,round(beginSilence/dt));
+        finalSilence= zeros(1,round(endSilence/dt));
+        inputSignal= [intialSilence inputSignal finalSilence];
+
+        %% run the model
+%         fprintf('\n')
+%         disp([num2str(numChannels) ' channel model: ' AN_spikesOrProbability])
+%         disp('Computing ...')
+
+        MAP1_14(inputSignal, sampleRate, BFlist, ...
+            MAPparamsName, AN_spikesOrProbability, paramChanges);
+
+        peakOutputDisp(levelNo,freqNo)=max(DRNLoutput);
+        peakStapesDisp(levelNo,freqNo)=max(TMoutput);
+
+        %% the model run is complete. Now display the results
+%         UTIL_showMAP(showMapOptions, paramChanges)
+    end % probe frequencies
+figure(2), loglog(toneFrequencyList, peakOutputDisp), hold on
+xlabel('frequency')
+ylabel('peak DRNL displacement (m)')
+    
+end  % levels
+%% alignmed plot
+
+
+x=peakOutputDisp./peakStapesDisp;
+logOutput=20*log10(x.*repmat((2*pi*toneFrequencyList),length(levels),1));
+alignedOutput=logOutput-repmat(logOutput(:,1), 1, length(toneFrequencyList));
+figure(3), clf, semilogx(toneFrequencyList, alignedOutput), hold on
+ylim([-20 80])
+xlim([2000 20000])
+% legend(num2str(levels'))
+xlabel('frequency')
+ylabel('peak DRNL velocity (dB)')
+title(['level ' num2str(leveldBSPL) ' dB SPL'])
+
+disp(['level=' num2str(leveldBSPL)])
+disp(['toneFrequency=' num2str(toneFrequency)])
+
+disp(['attenuation factor =' ...
+    num2str(DRNLParams.rateToAttenuationFactor, '%5.3f') ])
+disp(['attenuation factor (probability)=' ...
+    num2str(DRNLParams.rateToAttenuationFactorProb, '%5.3f') ])
+disp(AN_spikesOrProbability)
+disp(paramChanges)
+
+%     UTIL_printTabTable([toneFrequencyList' peakOutputDisp'])
+
+path(restorePath)
+