comparison testPrograms/testRP.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 3ea506487b3b
children
comparison
equal deleted inserted replaced
37:771a643d5c29 38:c2204b18f4a2
1 function testRP(BFs,MAPparamsName,paramChanges) 1 function testRP(BFs,MAPparamsName,paramChanges)
2 % testIHC used for IHC I/O function 2 % testIHC evaluates IHC I/O function
3 % multiple BFs can be used but only one is easier to interpret. 3 % multiple BFs can be used but only one is easier to interpret.
4 % e.g. testRP(1000,'Normal',{}); 4 % e.g. testRP(1000,'Normal',{});
5 5
6 global experiment method inputStimulusParams 6 global experiment method inputStimulusParams
7 global stimulusParameters IHC_VResp_VivoParams IHC_cilia_RPParams 7 global stimulusParameters IHC_VResp_VivoParams IHC_cilia_RPParams
15 drawColors='rgbkmcy'; 15 drawColors='rgbkmcy';
16 drawnow 16 drawnow
17 17
18 if nargin<3 18 if nargin<3
19 paramChanges=[]; 19 paramChanges=[];
20 end
21 if nargin<2
22 MAPparamsName='Normal';
23 end
24 if nargin<3
25 BFs=800;
20 end 26 end
21 27
22 levels=-20:10:100; 28 levels=-20:10:100;
23 nLevels=length(levels); 29 nLevels=length(levels);
24 toneDuration=.05; 30 toneDuration=.05;
131 %% plot cilia displacement 137 %% plot cilia displacement
132 figure(4) 138 figure(4)
133 subplot(2,2,2) 139 subplot(2,2,2)
134 restingIHC_cilia=IHCrestingCiliaCond; 140 restingIHC_cilia=IHCrestingCiliaCond;
135 plot(levels, IHC_cilia_peak,'k', 'linewidth',2), hold on 141 plot(levels, IHC_cilia_peak,'k', 'linewidth',2), hold on
136 plot(levels, IHC_cilia_min,'r', 'linewidth',2) 142 % plot(levels, IHC_cilia_min,'r', 'linewidth',2)
137 hold on, 143 % hold on,
138 plot([min(levels) max(levels)], ... 144 % plot([min(levels) max(levels)], ...
139 [restingIHC_cilia restingIHC_cilia], 'g') 145 % [restingIHC_cilia restingIHC_cilia], 'g')
140 title(' IHC apical cond.') 146 title(' IHC apical cond.')
141 ylabel ('IHCcilia(conductance)'), xlabel('dB SPL') 147 ylabel ('IHCcilia(conductance)'), xlabel('dB SPL')
142 xlim([min(levels) max(levels)]) 148 xlim([min(levels) max(levels)])
143 grid on 149 grid on
144 150