diff userProgramsTim/speech_stimuli_check.m @ 38:c2204b18f4a2 tip

End nov big change

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/userProgramsTim/speech_stimuli_check.m	Mon Nov 28 13:34:28 2011 +0000
@@ -0,0 +1,168 @@
+%skript to test several stimuli with MAP1_14
+% Tim Juergens, September 2011
+
+%addpath('C:\Documents and Settings\tjurgens\My Documents\Dropbox\matlab');
+addpath('C:\MAP\userPrograms');
+parameterfile='Normal';
+
+%Read the wavfiles
+[da,sfreq] = wavread('..\wavFileStore\new-da-44khz.wav'); %artificial stimulus 
+ba = wavread('..\wavFileStore\new-ba-44khz.wav'); %artificial stimulus
+
+%set level of speech stimuli (RMS 1 equals 94 dB)
+da_69dB = da./sqrt(mean(da.^2)).*10^(-(94-69)/20);
+da_49dB = da./sqrt(mean(da.^2)).*10^(-(94-49)/20);
+ba_69dB = ba./sqrt(mean(ba.^2)).*10^(-(94-69)/20);
+%20*log10(sqrt(mean(ba_69dB.^2))/20e-6) %reference pressure: 20uPa
+
+
+cd ..
+cd MAP
+
+%MAP and store AN output of HSR fibers in variables
+MAP1_14(da_69dB,sfreq,-1,parameterfile,'probability');
+global ANprobRateOutput
+AN_da_69dB = ANprobRateOutput(1:size(ANprobRateOutput,1)/2,:);
+
+ MAP1_14(da_49dB,sfreq,-1,parameterfile,'probability');
+ global ANprobRateOutput
+ AN_da_49dB = ANprobRateOutput(1:size(ANprobRateOutput,1)/2,:);
+% 
+ MAP1_14(ba_69dB,sfreq,-1,parameterfile,'probability');
+ global ANprobRateOutput
+ global savedBFlist
+ AN_ba_69dB = ANprobRateOutput(1:size(ANprobRateOutput,1)/2,:);
+
+%Do the IPIH analysis
+ [iih_da_69dB,IPIhisttime_da_69dB,IPIhistweight_da_69dB] = track_formants_from_IPI_guy(AN_da_69dB,sfreq);
+ %[iih_da_69dB,IPIhisttime_da_69dB,IPIhistweight_da_69dB] = IPIHextract(AN_da_69dB,sfreq);
+ poolIPI_across_channels(IPIhisttime_da_69dB,IPIhistweight_da_69dB);
+ map_iih_onto_log(iih_da_69dB,30,sfreq);
+  [iih_da_49dB,IPIhisttime_da_49dB,IPIhistweight_da_49dB] = track_formants_from_IPI_guy(AN_da_49dB,sfreq);
+%  poolIPI_across_channels(IPIhisttime_da_49dB,IPIhistweight_da_49dB);
+%  
+  [iih_ba_69dB,IPIhisttime_ba_69dB,IPIhistweight_ba_69dB] = track_formants_from_IPI_guy(AN_ba_69dB,sfreq);
+%  [tmp,ctr]=poolIPI_across_channels(IPIhisttime_ba_69dB,IPIhistweight_ba_69dB);
+%  
+ 
+ 
+ 
+ %%%%%% OLLO stimuli %%%%%%%%%%%%%%%%%
+%  OLLOwavfiles = {'dahd','bahb','dehd','behb','died','bieb','dohd','bohb','duhd','buhb','atta','ascha','assa'};
+%  OLLO_level = 69; %level of OLLO files in dB SPL
+%  paramChanges = [];%{'DRNLParams.rateToAttenuationFactorProb = 0;'};
+%  
+% for iCounter = 1:length(OLLOwavfiles)
+%     
+%     %read the stimuli
+%     eval(['[' OLLOwavfiles{iCounter} ',sfreq_OLLO] = wavread([''..\wavFileStore\S02M_L' ...
+%         sprintf('%3.3i',sloga2iloga(OLLOwavfiles{iCounter})) '_V6_M1_N2_CS0.wav'']);'])
+%     
+%     %delete preceding and subsequent silence and resample to 44100 Hz
+%     %sampling frequency
+%     eval([OLLOwavfiles{iCounter} ' = cutsignal(' OLLOwavfiles{iCounter} ',sfreq_OLLO,''d_d'');']);
+%     eval([OLLOwavfiles{iCounter} ' = resample(' OLLOwavfiles{iCounter} ',sfreq,sfreq_OLLO);']);
+%     
+%     %set level
+%     partfilename = [OLLOwavfiles{iCounter} num2str(OLLO_level) 'dB'];
+%     eval([partfilename ' = ' OLLOwavfiles{iCounter} ...
+%         './sqrt(mean(' OLLOwavfiles{iCounter} '.^2)).*10^(-(94-' num2str(OLLO_level) ')/20);']);
+%     
+%     %use MAP
+%     eval(['MAP1_14(' partfilename ',sfreq,-1,parameterfile,''probability'', paramChanges );']);
+%     global ANprobRateOutput
+%     eval(['AN_' partfilename ' = ANprobRateOutput(1:size(ANprobRateOutput,1)/2,:);']);
+%     
+%     %Do the IPIH analysis
+%     eval(['[iih_' partfilename ',IPIhisttime_' partfilename ',IPIhistweight_' ...
+%         partfilename '] = track_formants_from_IPI_guy(AN_' partfilename ',sfreq);']);
+%     eval(['poolIPI_across_channels(IPIhisttime_' partfilename ',IPIhistweight_' partfilename ');']);
+%     title(partfilename);
+%     %set(gca,'Title',partfilename);
+%     xlabel('Interval (ms)')
+%     ylabel('Stimulus time (ms)');
+%     
+%     eval(['map_iih_onto_log(iih_' partfilename ',30,sfreq);']);
+%     title(partfilename);
+% end
+
+
+%%%%%%%%%%%%% OLLO stimuli from different speakers
+%  OLLOwavfiles = {'S01F_L111_V6_M1_N2_CS0.wav','S02M_L111_V6_M1_N2_CS0.wav'};
+%  OLLO_level = 69; %level of OLLO files in dB SPL
+% for iCounter = 1:length(OLLOwavfiles)
+%     
+%     partfilename = [OLLOwavfiles{iCounter}(1:end-4) num2str(OLLO_level) 'dB'];
+%     %read the stimuli
+%     eval(['[' partfilename ',sfreq_OLLO] = wavread([''..\wavFileStore\' OLLOwavfiles{iCounter} ''']);'])
+%     
+%     %delete preceding and subsequent silence and resample to 44100 Hz
+%     %sampling frequency
+%     eval([partfilename ' = cutsignal(' partfilename ',sfreq_OLLO,''d_d'');']);
+%     eval([partfilename ' = resample(' partfilename ',sfreq,sfreq_OLLO);']);
+%     
+%     %set level
+%     
+%     eval([partfilename ' = ' partfilename ...
+%         './sqrt(mean(' partfilename '.^2)).*10^(-(94-' num2str(OLLO_level) ')/20);']);
+%     
+%     %use MAP
+%     eval(['MAP1_14(' partfilename ',sfreq,-1,parameterfile,''probability'');']);
+%     global ANprobRateOutput
+%     eval(['AN_' partfilename ' = ANprobRateOutput(1:size(ANprobRateOutput,1)/2,:);']);
+%     
+%     %Do the IPIH analysis
+%     eval(['[iih_' partfilename ',IPIhisttime_' partfilename ',IPIhistweight_' ...
+%         partfilename '] = track_formants_from_IPI_guy(AN_' partfilename ',sfreq);']);
+%     eval(['poolIPI_across_channels(IPIhisttime_' partfilename ',IPIhistweight_' partfilename ');']);
+%     title(partfilename);
+%     %set(gca,'Title',partfilename);
+%     xlabel('Interval (ms)')
+%     ylabel('Stimulus time (ms)');
+%     
+%     eval(['map_iih_onto_log(iih_' partfilename ',30,sfreq);']);
+%     title(partfilename);
+% end
+
+
+
+%%%%%% da stimuli with different pitches %%%%%%%%%%%%%%%%%
+ 
+  dawavfiles = { '200ms_da_080Hz.wav','200ms_da_100Hz.wav','200ms_da_120Hz.wav','200ms_da_140Hz.wav', ... 
+      '200ms_da_160Hz.wav','200ms_da_180Hz.wav','200ms_da_200Hz.wav','200ms_da_220Hz.wav', ...
+      '200ms_da_240Hz.wav', ...
+      'noise.wav', ...
+      'da_whispered.wav', ...
+      };
+      
+ da_level = 69; %level of OLLO files in dB SPL
+for iCounter = 1:length(dawavfiles)
+    
+    %read the stimuli
+    partfilename = ['da' dawavfiles{iCounter}(1:end-4)];
+    eval([ partfilename ' = wavread(''..\wavFileStore\' dawavfiles{iCounter} ''');'])
+    
+    %set level
+    eval([partfilename ' = ' partfilename ...
+        './sqrt(mean(' partfilename '.^2)).*10^(-(94-' num2str(da_level) ')/20);']);
+    
+    %use MAP
+    eval(['MAP1_14(' partfilename ',sfreq,-1,parameterfile,''probability'');']);
+    global ANprobRateOutput
+    eval(['AN_' partfilename ' = ANprobRateOutput(1:size(ANprobRateOutput,1)/2,:);']);
+    
+    %Do the IPIH analysis
+    eval(['[iih_' partfilename ',IPIhisttime_' partfilename ',IPIhistweight_' ...
+            partfilename '] = IPIHextract(AN_' partfilename ',sfreq);']);%partfilename '] = track_formants_from_IPI_guy(AN_' partfilename ',sfreq);']);%
+    eval(['poolIPI_across_channels(IPIhisttime_' partfilename ',IPIhistweight_' partfilename ');']);
+    title(partfilename);
+    %set(gca,'Title',partfilename);
+    xlabel('Interval (ms)')
+    ylabel('Stimulus time (ms)');
+    
+    eval(['map_iih_onto_log(iih_' partfilename ',30,sfreq);']);
+    title(partfilename);
+end
+% 
+% 
+