rmeddis@28: function testRP(BFs,MAPparamsName,paramChanges) rmeddis@28: % testIHC used for IHC I/O function rmeddis@28: rmeddis@28: global experiment method inputStimulusParams rmeddis@28: global stimulusParameters IHC_VResp_VivoParams IHC_cilia_RPParams rmeddis@28: savePath=path; rmeddis@28: addpath (['..' filesep 'utilities'],['..' filesep 'MAP']) rmeddis@28: dbstop if error rmeddis@28: rmeddis@28: figure(4), clf, rmeddis@28: set (gcf, 'name', ['IHC']) rmeddis@28: set(gcf,'position',[613 354 360 322]) rmeddis@28: drawColors='rgbkmcy'; rmeddis@28: drawnow rmeddis@28: rmeddis@28: if nargin<3, paramChanges=[]; end rmeddis@28: rmeddis@28: levels=-20:10:100; rmeddis@28: nLevels=length(levels); rmeddis@28: toneDuration=.05; rmeddis@28: silenceDuration=.01; rmeddis@28: sampleRate=50000; rmeddis@28: dt=1/sampleRate; rmeddis@28: rmeddis@28: allIHC_RP_peak=[]; rmeddis@28: allIHC_RP_dc=[]; rmeddis@28: rmeddis@28: changes{1}=[]; rmeddis@28: % changes{1}='IHC_cilia_RPParams.Et= 0.100;'; rmeddis@28: % changes{2}='IHC_cilia_RPParams.Et= 0.070;'; rmeddis@28: rmeddis@28: for changeNo=1:length(changes) rmeddis@28: if isempty(changes{1}) rmeddis@28: paramChanges=[]; rmeddis@28: else rmeddis@28: paramChanges{1}=changes{changeNo}; rmeddis@28: end rmeddis@28: for BFno=1:length(BFs) rmeddis@28: BF=BFs(BFno); rmeddis@28: targetFrequency=BF; rmeddis@28: % OR rmeddis@28: %Patuzzi and Sellick test (see ELP & AEM, 2006) rmeddis@28: % targetFrequency=100; rmeddis@28: rmeddis@28: IHC_RP_peak=zeros(nLevels,1); rmeddis@28: IHC_RP_min=zeros(nLevels,1); rmeddis@28: IHC_RP_dc=zeros(nLevels,1); rmeddis@28: rmeddis@28: time=dt:dt:toneDuration; rmeddis@28: rmeddis@28: rampDuration=0.004; rmeddis@28: rampTime=dt:dt:rampDuration; rmeddis@28: ramp=[0.5*(1+cos(2*pi*rampTime/(2*rampDuration)+pi)) ... rmeddis@28: ones(1,length(time)-length(rampTime))]; rmeddis@28: ramp=ramp.*fliplr(ramp); rmeddis@28: rmeddis@28: silence=zeros(1,round(silenceDuration/dt)); rmeddis@28: rmeddis@28: toneStartptr=length(silence)+1; rmeddis@28: toneMidptr=toneStartptr+round(toneDuration/(2*dt)) -1; rmeddis@28: toneEndptr=toneStartptr+round(toneDuration/dt) -1; rmeddis@28: rmeddis@28: levelNo=0; rmeddis@28: for leveldB=levels rmeddis@28: levelNo=levelNo+1; rmeddis@28: % replicate at all levels rmeddis@28: amp=28e-6*10^(leveldB/20); rmeddis@28: rmeddis@28: %% create signal (leveldB/ frequency) rmeddis@28: inputSignal=amp*sin(2*pi*targetFrequency'*time); rmeddis@28: inputSignal= ramp.*inputSignal; rmeddis@28: inputSignal=[silence inputSignal silence]; rmeddis@28: inputStimulusParams.sampleRate=1/dt; rmeddis@28: % global IHC_ciliaParams rmeddis@28: rmeddis@28: %% disable efferent for fast processing rmeddis@28: method.DRNLSave=1; rmeddis@28: method.IHC_cilia_RPSave=1; rmeddis@28: method.IHCpreSynapseSave=1; rmeddis@28: method.IHC_cilia_RPSave=1; rmeddis@28: method.segmentDuration=-1; rmeddis@28: moduleSequence=1:4; rmeddis@28: rmeddis@28: %% run the model rmeddis@28: global DRNLoutput IHC_cilia_output IHCrestingCiliaCond IHCrestingV... rmeddis@28: IHCoutput rmeddis@28: AN_spikesOrProbability='probability'; rmeddis@28: rmeddis@28: MAP1_14(inputSignal, sampleRate, BF, ... rmeddis@28: MAPparamsName, AN_spikesOrProbability, paramChanges); rmeddis@28: rmeddis@28: % DRNL rmeddis@28: DRNLoutput=DRNLoutput; rmeddis@28: DRNL_peak(levelNo,1)=max(DRNLoutput(toneMidptr:toneEndptr)); rmeddis@28: DRNL_min(levelNo,1)=min(DRNLoutput(toneMidptr:toneEndptr)); rmeddis@28: DRNL_dc(levelNo,1)=mean(DRNLoutput(toneMidptr:toneEndptr)); rmeddis@28: rmeddis@28: % cilia rmeddis@28: IHC_ciliaData=IHC_cilia_output; rmeddis@28: IHC_ciliaData=IHC_ciliaData; rmeddis@28: IHC_cilia_peak(levelNo,1)=... rmeddis@28: max(IHC_ciliaData(toneMidptr:toneEndptr)); rmeddis@28: IHC_cilia_min(levelNo,1)=... rmeddis@28: min(IHC_ciliaData(toneMidptr:toneEndptr)); rmeddis@28: IHC_cilia_dc(levelNo,1)=... rmeddis@28: mean(IHC_ciliaData(toneMidptr:toneEndptr)); rmeddis@28: rmeddis@28: % RP rmeddis@28: IHC_RPData=IHCoutput; rmeddis@28: IHC_RPData=IHC_RPData; rmeddis@28: IHC_RP_peak(levelNo,1)=... rmeddis@28: max(IHC_RPData(toneMidptr:toneEndptr)); rmeddis@28: IHC_RP_min(levelNo,1)=... rmeddis@28: min(IHC_RPData(toneMidptr:toneEndptr)); rmeddis@28: IHC_RP_dc(levelNo,1)=... rmeddis@28: mean(IHC_RPData(toneMidptr:toneEndptr)); rmeddis@28: rmeddis@28: end % level rmeddis@28: rmeddis@28: rmeddis@28: disp(['parameter file was: ' MAPparamsName]) rmeddis@28: fprintf('\n') rmeddis@28: rmeddis@28: %% plot DRNL rmeddis@28: subplot(2,2,1) rmeddis@28: % referenceDisp= 9e-9*1000/targetFrequency; rmeddis@28: % plot(levels,20*log10(DRNL_peak/referenceDisp), drawColors(BFno), ... rmeddis@28: % 'linewidth',2), hold on rmeddis@28: referenceDisp=10e-9; rmeddis@28: plot(levels,20*log10(DRNL_peak/referenceDisp), drawColors(BFno), ... rmeddis@28: 'linewidth',2), hold on rmeddis@28: title([' DRNL peak: ' num2str(BFs) ' Hz']) rmeddis@28: ylabel ('log10DRNL(m)'), xlabel('dB SPL') rmeddis@28: xlim([min(levels) max(levels)]), ylim([-10 50]) rmeddis@28: grid on rmeddis@28: rmeddis@28: %% plot cilia displacement rmeddis@28: figure(4) rmeddis@28: subplot(2,2,2) rmeddis@28: restingIHC_cilia=IHCrestingCiliaCond; rmeddis@28: plot(levels, IHC_cilia_peak,'k', 'linewidth',2), hold on rmeddis@28: plot(levels, IHC_cilia_min,'r', 'linewidth',2) rmeddis@28: hold on, rmeddis@28: plot([min(levels) max(levels)], ... rmeddis@28: [restingIHC_cilia restingIHC_cilia], 'g') rmeddis@28: title(' IHC apical cond.') rmeddis@28: ylabel ('IHCcilia(conductance)'), xlabel('dB SPL') rmeddis@28: xlim([min(levels) max(levels)]) rmeddis@28: grid on rmeddis@28: rmeddis@28: %% plot receptor potentials rmeddis@28: figure(4) rmeddis@28: subplot(2,2,3) rmeddis@28: % RP I/O function min and max rmeddis@28: restingRP=IHC_RP_peak(1); rmeddis@28: toPlot= [fliplr(IHC_RP_min(:,1)') IHC_RP_peak(:,1)']; rmeddis@28: microPa= 28e-6*10.^(levels/20); rmeddis@28: microPa=[-fliplr(microPa) microPa]; rmeddis@28: plot(microPa,toPlot, drawColors(BFno), 'linewidth',2) rmeddis@28: % ylim([0 300]) rmeddis@28: rmeddis@28: %% Dallos and Harris data rmeddis@28: dallosx=[-0.9 -0.1 -0.001 0.001 0.01 0.9]; rmeddis@28: dallosy=[-8 -7.8 -6.5 11 16.5 22]/1000 + restingRP; rmeddis@28: hold on, plot(dallosx,dallosy, 'o') rmeddis@28: plot([-1 1], [restingRP restingRP], 'r') rmeddis@28: title(' Dallos(86) data at 800 Hz') rmeddis@28: ylabel ('receptor potential(V)'), xlabel('Pa') rmeddis@28: ylim([-0.08 -0.02]), xlim([-1 1]) rmeddis@28: grid on rmeddis@28: rmeddis@28: %% RP I/O function min and max rmeddis@28: figure(4) rmeddis@28: subplot(2,2,4) rmeddis@28: restingRP=IHC_RP_peak(1); rmeddis@28: peakRP=max(IHC_RP_peak); rmeddis@28: plot(levels, IHC_RP_peak,drawColors(BFno), 'linewidth',2) rmeddis@28: hold on rmeddis@28: plot(levels, IHC_RP_dc, [drawColors(BFno) ':'], 'linewidth',2) rmeddis@28: hold on, rmeddis@28: plot([min(levels) max(levels)], [restingRP restingRP], 'r') rmeddis@28: xlim([min(levels) max(levels)]) rmeddis@28: grid on rmeddis@28: title(['Et= ' num2str(IHC_cilia_RPParams.Et) ': RP (AC- / DC:']) rmeddis@28: ylabel ('RP(V)'), xlabel('dB SPL') rmeddis@28: ylim([-0.08 -0.02]) rmeddis@28: allIHC_RP_peak=[allIHC_RP_peak IHC_RP_peak]; rmeddis@28: allIHC_RP_dc=[allIHC_RP_dc IHC_RP_dc]; rmeddis@28: rmeddis@28: fprintf('level\t peak\t DC\n') rmeddis@28: UTIL_printTabTable([levels' IHC_RP_peak IHC_RP_dc]) rmeddis@28: % disp(['restingIHC_cilia= ' num2str(restingIHC_cilia)]) rmeddis@28: fprintf('peakRP= \t%6.3f', peakRP) rmeddis@28: fprintf('\nrestingRP= \t%6.3f', restingRP) rmeddis@28: fprintf('\ndifference= \t%6.3f\n', (peakRP-restingRP)) rmeddis@28: drawnow rmeddis@28: end rmeddis@28: end rmeddis@28: % UTIL_showStruct(IHC_VResp_VivoParams, 'IHC_VResp_VivoParams') rmeddis@28: UTIL_showStruct(IHC_cilia_RPParams, 'IHC_cilia_RPParams') rmeddis@28: fprintf('level\t peak\n') rmeddis@28: UTIL_printTabTable([levels' allIHC_RP_peak]) rmeddis@28: fprintf('level\t DC\n') rmeddis@28: UTIL_printTabTable([levels' allIHC_RP_dc]) rmeddis@28: rmeddis@28: path(savePath); rmeddis@28: disp(paramChanges)