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