diff testPrograms/test2toneMovie.m @ 38:c2204b18f4a2 tip

End nov big change

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/testPrograms/test2toneMovie.m	Mon Nov 28 13:34:28 2011 +0000
@@ -0,0 +1,151 @@
+% function test2toneMovie
+%
+% Demonstration of two-tone suppression in the AN using a
+%   F1 tone suppressed by a 1.5*F1 kHz tone
+%
+
+global  ANprobRateOutput DRNLoutput
+dbstop if error
+restorePath=path;
+addpath (['..' filesep 'MAP'],    ['..' filesep 'wavFileStore'], ...
+    ['..' filesep 'utilities'])
+figure(5), clf
+
+disp('')
+primaryToneFrequency=8000;
+
+probedBs=10 :10: 50;
+nProbeLevels=length(probedBs);
+
+% disp(['F1 F1 level= ' num2str([primaryToneFrequency probedB])])
+
+suppressorLevels=0:10:80;
+
+lowestBF=250; 	highestBF= 8000; 	numChannels=21;
+BFlist=round(logspace(log10(lowestBF), log10(highestBF), numChannels));
+% fix so that probe is middle range and is at a BF
+BFlist=BFlist*primaryToneFrequency/2000;
+suppressorFrequencies=BFlist;
+BFchannel=find(round(BFlist)==round(primaryToneFrequency));
+
+sampleRate= 100000; % Hz (higher sample rate needed for BF>8000 Hz)
+dt=1/sampleRate; % seconds
+
+duration=.050;		      % seconds
+    tone2Duration=duration/2; % s
+startSilenceDuration=0.010;
+startSilence=...
+    zeros(1,startSilenceDuration*sampleRate);
+suppressorStartsPTR=...
+    round((startSilenceDuration+duration/2)*sampleRate);
+
+probedBCount=0;
+for probedB=probedBs
+    probedBCount=probedBCount+1;
+    
+peakResponse= zeros(length(suppressorLevels),length(suppressorFrequencies));
+suppFreqCount=0;
+for suppressorFrequency=suppressorFrequencies
+    suppFreqCount=suppFreqCount+1;
+    
+    suppressorLevelCount=0;
+for suppressorDB=suppressorLevels
+suppressorLevelCount=suppressorLevelCount+1;
+
+    % primary + suppressor
+    primaryLevelDB=probedB;
+    suppressorLevelDB=suppressorDB;
+
+
+% primary BF tone
+    time1=dt: dt: duration;
+    amp=10^(primaryLevelDB/20)*28e-6;
+    inputSignal=amp*sin(2*pi*primaryToneFrequency*time1);
+    rampDuration=.005; rampTime=dt:dt:rampDuration;
+    ramp=[0.5*(1+cos(2*pi*rampTime/(2*rampDuration)+pi)) ones(1,length(time1)-length(rampTime))];
+    inputSignal=inputSignal.*ramp;
+    inputSignal=inputSignal.*fliplr(ramp);
+
+    % suppressor
+    time2= dt: dt: tone2Duration;
+    % B: tone on
+    amp=10^(suppressorLevelDB/20)*28e-6;
+    inputSignal2=amp*sin(2*pi*suppressorFrequency*time2);
+    rampDuration=.005; rampTime=dt:dt:rampDuration;
+    ramp=[0.5*(1+cos(2*pi*rampTime/(2*rampDuration)+pi)) ones(1,length(time2)-length(rampTime))];
+    inputSignal2=inputSignal2.*ramp;
+    inputSignal2=inputSignal2.*fliplr(ramp);
+    % A: initial silence (delay to suppressor)
+    silence=zeros(1,length(time2));
+    inputSignal2=[silence inputSignal2];
+
+    % add tone and suppressor components
+    inputSignal=inputSignal+inputSignal2;
+
+    inputSignal=...
+        [startSilence inputSignal ];
+
+    % run MAP
+    MAPparamsName='Normal';
+    AN_spikesOrProbability='probability';
+    paramChanges={'IHCpreSynapseParams.tauCa=80e-6;'};
+    MAP1_14(inputSignal, sampleRate, BFlist, ...
+        MAPparamsName, AN_spikesOrProbability, paramChanges);
+
+    response=ANprobRateOutput;
+peakResponse(suppressorLevelCount,suppFreqCount)=...
+    mean(response(BFchannel,suppressorStartsPTR:end));
+disp(['F2 level= ', num2str([suppressorFrequency suppressorDB ...
+    peakResponse(suppressorLevelCount,suppFreqCount)])])
+
+
+%% make movie
+    figure(6), subplot(4,1,2)
+    axis tight
+    set(gca,'nextplot','replacechildren');
+    plot(dt:dt:dt*length(inputSignal), inputSignal, 'k')
+    ylim([-.5 .5])
+
+    figure(6), subplot(2,1,2)
+    PSTHbinWidth=0.010;
+    PSTH= UTIL_PSTHmakerb(...
+        response, dt, PSTHbinWidth);
+    [nY nX]=size(PSTH);
+    time=PSTHbinWidth*(0:nX-1);
+    surf(time, BFlist, PSTH)
+    zlim([0 500]),
+    caxis([0 500])
+    shading interp
+    colormap(jet)
+    set(gca, 'yScale','log')
+    xlim([0 max(time)+dt])
+    ylim([0 max(BFlist)])
+    view([0 90]) % view([-8 68])
+    title([' primary + suppressor '])
+
+end
+end
+
+%% plot matrix
+
+peakResponse=peakResponse/peakResponse(1,1);
+
+figure(5), subplot(2,nProbeLevels, probedBCount)
+surf(suppressorFrequencies,suppressorLevels,peakResponse)
+shading interp
+view([0 90])
+set(gca,'Xscale','log')
+xlim([min(suppressorFrequencies) max(suppressorFrequencies)])
+set(gca,'Xtick', [1000  4000],'xticklabel',{'1000', '4000'})
+title([num2str(primaryToneFrequency) ' Hz'])
+
+subplot(2,nProbeLevels,nProbeLevels+probedBCount)
+contourf(suppressorFrequencies,suppressorLevels,peakResponse, ...
+    [.1:.1:.9 0.99] )
+set(gca,'Xscale','log')
+set(gca,'Xtick', [1000  4000],'xticklabel',{'1000', '4000'})
+set(gcf, 'name',['primaryToneFrequency= ' num2str(primaryToneFrequency)])
+title([num2str(probedB) ' dB'])
+end
+
+path=restorePath;