view userProgramsTim/speechmanipulation.m @ 38:c2204b18f4a2 tip

End nov big change

function varargout = speechmanipulation(varargin)
% SPEECHMANIPULATION MATLAB code for speechmanipulation.fig
%      SPEECHMANIPULATION, by itself, creates a new SPEECHMANIPULATION or raises the existing
%      singleton*.
%
%      H = SPEECHMANIPULATION returns the handle to a new SPEECHMANIPULATION or the handle to
%      the existing singleton*.
%
%      SPEECHMANIPULATION('CALLBACK',hObject,eventData,handles,...) calls the local
%      function named CALLBACK in SPEECHMANIPULATION.M with the given input arguments.
%
%      SPEECHMANIPULATION('Property','Value',...) creates a new SPEECHMANIPULATION or raises the
%      existing singleton*.  Starting from the left, property value pairs are
%      applied to the GUI before speechmanipulation_OpeningFcn gets called.  An
%      unrecognized property name or invalid value makes property application
%      stop.  All inputs are passed to speechmanipulation_OpeningFcn via varargin.
%
%      *See GUI Options on GUIDE's Tools menu.  Choose "GUI allows only one
%      instance to run (singleton)".
%
% See also: GUIDE, GUIDATA, GUIHANDLES

% Edit the above text to modify the response to help speechmanipulation

% Last Modified by GUIDE v2.5 05-Oct-2011 11:21:58

% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name',       mfilename, ...
                   'gui_Singleton',  gui_Singleton, ...
                   'gui_OpeningFcn', @speechmanipulation_OpeningFcn, ...
                   'gui_OutputFcn',  @speechmanipulation_OutputFcn, ...
                   'gui_LayoutFcn',  [] , ...
                   'gui_Callback',   []);
if nargin && ischar(varargin{1})
    gui_State.gui_Callback = str2func(varargin{1});
end

if nargout
    [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
    gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT


% --- Executes just before speechmanipulation is made visible.
function speechmanipulation_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject    handle to figure
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
% varargin   command line arguments to speechmanipulation (see VARARGIN)
global wavfilescell;
global rateaxeshandle; %handle to axis6 is needed, if radiobuttons12 and 13 change


% Choose default command line output for speechmanipulation
handles.output = hObject;

% Update handles structure
guidata(hObject, handles);

% UIWAIT makes speechmanipulation wait for user response (see UIRESUME)
% uiwait(handles.figure1);
wavfiles = dir(['..' filesep 'wavFileStore' filesep '*.wav']);
for iCounter = 1:length(wavfiles)
    wavfilescell{iCounter} = wavfiles(iCounter).name;
end
set(handles.popupmenu1,'String',wavfilescell);
set(handles.popupmenu2,'String',wavfilescell);
set(handles.uipanel8,'SelectionChangeFcn',@radiobuttonselected);
rateaxeshandle = handles.axes6;

% --- Outputs from this function are returned to the command line.
function varargout = speechmanipulation_OutputFcn(hObject, eventdata, handles) 
% varargout  cell array for returning output args (see VARARGOUT);
% hObject    handle to figure
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Get default command line output from handles structure
varargout{1} = handles.output;

function edit1_Callback(hObject, eventdata, handles)
function edit2_Callback(hObject, eventdata, handles)
function edit3_Callback(hObject, eventdata, handles)
function edit4_Callback(hObject, eventdata, handles)
function edit5_Callback(hObject, eventdata, handles)
function edit6_Callback(hObject, eventdata, handles)
function edit7_Callback(hObject, eventdata, handles)
function edit8_Callback(hObject, eventdata, handles)
function edit9_Callback(hObject, eventdata, handles)
function edit10_Callback(hObject, eventdata, handles)
function edit11_Callback(hObject, eventdata, handles)
function edit12_Callback(hObject, eventdata, handles)


% --- Executes on selection change in popupmenu1.
function popupmenu1_Callback(hObject, eventdata, handles)
% hObject    handle to popupmenu1 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Hints: contents = cellstr(get(hObject,'String')) returns popupmenu1 contents as cell array
%        contents{get(hObject,'Value')} returns selected item from popupmenu1
%% CHOOSE AND LOAD WAVEFORM
global wavfilescell
global originalsignal
global actualsignal

%get the number from the popupmenu
selectedfileno = get(handles.popupmenu1,'Value');

%load it
[waveform,sfreq] = wavread(['..' filesep 'wavFileStore' filesep wavfilescell{selectedfileno}]);

%if it is from OLLO then cut it
if regexp(wavfilescell{selectedfileno}, 'S[0-9]*[A-Z]_L[0-9]*_V[0-9]_M[0-9]_N[0-9]_CS0.wav') %pattern for OLLO-files
    waveform = cutsignal(waveform,sfreq,'a_a');
end

%store it in global variables
actualsignal.waveform = waveform;
actualsignal.sfreq = sfreq;
originalsignal.waveform = waveform;
originalsignal.sfreq = sfreq;

%plot it
multipleplot(actualsignal,handles);

% --- Executes during object creation, after setting all properties.
function popupmenu1_CreateFcn(hObject, eventdata, handles)
% hObject    handle to popupmenu1 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: popupmenu controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end


% --- Executes on button press in pushbutton1.
function pushbutton1_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton1 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
%% LOAD WAVEFORM
global wavfilescell
global originalsignal
global actualsignal

%get the number from the popupmenu
selectedfileno = get(handles.popupmenu1,'Value');

%load it
[waveform,sfreq] = wavread(['..' filesep 'wavFileStore' filesep wavfilescell{selectedfileno}]);

%if it is from OLLO then cut it
if regexp(wavfilescell{selectedfileno}, 'S[0-9]*[A-Z]_L[0-9]*_V[0-9]_M[0-9]_N[0-9]_CS0.wav') %pattern for OLLO-files
    waveform = cutsignal(waveform,sfreq,'a_a');
end

%store it in global variables
actualsignal.waveform = waveform;
actualsignal.sfreq = sfreq;
originalsignal.waveform = waveform;
originalsignal.sfreq = sfreq;

%plot it
multipleplot(actualsignal,handles)


% --- Executes during object creation, after setting all properties.
function edit1_CreateFcn(hObject, eventdata, handles)
% hObject    handle to edit1 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end


% --- Executes on button press in pushbutton2.
function pushbutton2_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton2 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
%% PLAY THE ACTUAL (MANIPULATED) SIGNAL
global actualsignal
global toplay

toplay = audioplayer(actualsignal.waveform,actualsignal.sfreq);
play(toplay);


% --- Executes on button press in pushbutton3.
function pushbutton3_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton3 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

%% MANIPULATE SIGNAL
global wavfilescell
global originalsignal
global actualsignal

actualsignal = originalsignal; %reset

%% speed factor change
speedfactor = str2num(get(handles.edit1,'String'));
newfakesfreq=actualsignal.sfreq*speedfactor;
if speedfactor == 1
else
    actualsignal.waveform = resample(actualsignal.waveform,actualsignal.sfreq,newfakesfreq);
end

%% low and high pass filter (butterworth)
% get the desired filter from radio buttons
radiofilterhandles=get(handles.uipanel2,'Children');
for iCounter = 1:length(radiofilterhandles)
    if get(radiofilterhandles(iCounter),'Value') == 1
        filtertype = get(radiofilterhandles(iCounter),'String');
    end
end
% read parameters
tmp = get(handles.edit3,'String');
cutofffreq=str2num(tmp{1});
tmp = get(handles.edit4,'String');
filterorder=str2num(tmp{1});
%do the filter calculations
switch filtertype
    case 'none'
    case 'High pass'
        if isempty(cutofffreq) || isempty(filterorder)
            error('Please specify cutoff frequency and/or filter order')
        else
            [a,b] = butter(filterorder,cutofffreq*2/actualsignal.sfreq,'high');
            actualsignal.waveform=filter(a,b,actualsignal.waveform);
        end
    case 'Low pass'
        if isempty(cutofffreq)|| isempty(filterorder)
            error('Please specify cutoff frequency')
        else
            [a,b] = butter(filterorder,cutofffreq*2/actualsignal.sfreq,'low');
            actualsignal.waveform=filter(a,b,actualsignal.waveform);
        end
end

%additive noise
% get the desired filter from radio buttons
radionoisehandles=get(handles.uipanel4,'Children');
for iCounter = 1:length(radionoisehandles)
    if get(radionoisehandles(iCounter),'Value') == 1
        noisetype = get(radionoisehandles(iCounter),'String');
    end
end      
% read parameters
tmp = get(handles.edit5,'String');
snr=str2num(tmp{1});
%do the noise calculations
switch noisetype
    case 'none'
    case 'White noise'
        if isempty(snr) 
            error('Please specify SNR')
        else
            whitenoise = rand(length(actualsignal.waveform),1)-0.5;
            %set level
            levelsignal = 20*log10(sqrt(mean(actualsignal.waveform.^2)));
            levelnoise = levelsignal-snr;
            whitenoise = whitenoise./sqrt(mean(whitenoise.^2)).*10^(levelnoise/20);
            %20*log10(sqrt(mean(whitenoise.^2)))
            actualsignal.waveform=actualsignal.waveform+whitenoise;
        end
    case 'From file'
        if isempty(snr) 
            error('Please specify SNR')
        else
            selectedfileno = get(handles.popupmenu2,'Value');

            [noisewaveform,noisesfreq] = wavread(['..' filesep 'wavFileStore' filesep wavfilescell{selectedfileno}]);
            if noisesfreq == actualsignal.sfreq
            else
                noisewaveform = resample(noisewaveform,actualsignal.sfreq,noisesfreq);
                noisesfreq = actualsignal.sfreq;
            end
            if length(noisewaveform) < length(actualsignal.waveform)
                warning('Noise waveform too short. Noise is looped without fading at the endings!');
                noisewaveform = repmat(noisewaveform, ceil(length(actualsignal.waveform)/length(noisewaveform)),1);
            end
            noisewaveform=noisewaveform(1:length(actualsignal.waveform));
            %set level
            levelsignal = 20*log10(sqrt(mean(actualsignal.waveform.^2)));
            levelnoise = levelsignal-snr;
            noisewaveform = noisewaveform./sqrt(mean(noisewaveform.^2)).*10^(levelnoise/20);
            %20*log10(sqrt(mean(noisewaveform.^2)))
            %finally add the two signals
            actualsignal.waveform = actualsignal.waveform + noisewaveform;
        end
end

%plot the resulting manipulated waveform
multipleplot(actualsignal,handles);

% --- Executes during object creation, after setting all properties.
function edit3_CreateFcn(hObject, eventdata, handles)
% hObject    handle to edit3 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end


% --- Executes during object creation, after setting all properties.
function edit4_CreateFcn(hObject, eventdata, handles)
% hObject    handle to edit4 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end


% --- Executes during object creation, after setting all properties.
function edit5_CreateFcn(hObject, eventdata, handles)
% hObject    handle to edit5 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end


% --- Executes during object creation, after setting all properties.
function popupmenu2_CreateFcn(hObject, eventdata, handles)
% hObject    handle to popupmenu2 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: popupmenu controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end


% --- Executes on button press in pushbutton4.
function pushbutton4_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton4 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
%% PLAY THE ORIGINAL WAVEFORM
global originalsignal
global toplay
toplay = audioplayer(originalsignal.waveform,originalsignal.sfreq);
play(toplay)


% --- Executes on button press in pushbutton5.
function pushbutton5_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton5 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
global toplay
stop(toplay);


function edit6_Callback(hObject, eventdata, handles)
% hObject    handle to edit6 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Hints: get(hObject,'String') returns contents of edit6 as text
%        str2double(get(hObject,'String')) returns contents of edit6 as a double


% --- Executes during object creation, after setting all properties.
function edit6_CreateFcn(hObject, eventdata, handles)
% hObject    handle to edit6 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end

% --- Executes during object creation, after setting all properties.
function edit7_CreateFcn(hObject, eventdata, handles)
% hObject    handle to edit7 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end


% --- Executes during object creation, after setting all properties.
function edit8_CreateFcn(hObject, eventdata, handles)
% hObject    handle to edit8 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end


% --- Executes on button press in pushbutton6.
function pushbutton6_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton6 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

%% CALCULATE MAP 1.14 auditory nerve probability output
addpath(['..' filesep 'MAP']);
addpath(['..' filesep 'utilities']);
global actualsignal
global AN_HSRoutput

tmp = get(handles.edit6,'String');
level4MAP=str2num(tmp);
if isempty(level4MAP)
    error('Please specify level');
end

tmp = get(handles.edit7,'String');
parameterfile=tmp;
if isempty(parameterfile)
    error('Please specify parameter file');
end

% identify model parameter changes if any
paramChanges=get(handles.edit8,'string');
if ~strcmp(paramChanges, ';'), paramChanges=[paramChanges ';']; end
eval(paramChanges);

%set level
actualsignal.waveform = actualsignal.waveform./sqrt(mean(actualsignal.waveform.^2)).*10^(-(94-level4MAP)/20);
%20*log10(sqrt(mean(actualsignal.waveform.^2))/20e-6) %reference pressure: 20uPa

MAP1_14(actualsignal.waveform,actualsignal.sfreq,-1,parameterfile, ...
    'probability',paramChanges);
global ANprobRateOutput savedBFlist

%take only the HSR fibers
AN_HSRoutput = ANprobRateOutput(size(ANprobRateOutput)/2+1:end,:);

%plot the auditory nerve firing probability as line plot
start_time = size(AN_HSRoutput,2)/2/actualsignal.sfreq-0.025; %start time to plot in s
end_time = start_time+0.05; %plot 50ms
plotIFRAN(AN_HSRoutput,start_time,end_time,actualsignal.sfreq,savedBFlist,handles.axes4);
set(handles.edit11,'String',num2str(1000*size(AN_HSRoutput,2)/2/actualsignal.sfreq));

%plot the fourierhistogram as image plot
formantpattern = fourierautocorrelationhistogram_direct(AN_HSRoutput,actualsignal.sfreq,handles.axes5);
caxis([0 2000]);
colorbar;

%plot the rate output
plotrateOutput(AN_HSRoutput,actualsignal.sfreq,handles.axes6,[33.7 300]);
set(handles.radiobutton12,'Value',1);

%calculate the IPIH
ipih=track_formants_from_IPI_guy(AN_HSRoutput, actualsignal.sfreq);
%the following code assumes that the bin width is 1/actualsignal.sfreq and 
%the temporal step size is 3ms

ipihfreqaxis=1./[1/actualsignal.sfreq:1/actualsignal.sfreq:size(ipih,1)/actualsignal.sfreq];
ipihtimeaxis=[0:3:size(ipih,2)*3];
set(gcf,'CurrentAxes',handles.axes7);
YTickIdx = 1:floor(numel(ipihfreqaxis)/6):numel(ipihfreqaxis);
XTickIdx = 1:floor(numel(ipihtimeaxis)/6):numel(ipihtimeaxis);
imagesc(ipih);
set(gca, 'YTick', YTickIdx);
set(gca, 'YTickLabel', num2str(  ipihfreqaxis(YTickIdx)', '%0.0f' ));
ylabel('best frequency (Hz)')
set(gca, 'XTick', XTickIdx);
set(gca, 'XTickLabel', XTickIdx.*3);
xlabel('Time (ms)');
caxis([0 8e4]); %set color 
colorbar;

%calculate and plot summarized autocorrelation, code from Ray
method.dt=1/actualsignal.sfreq;
method.segmentNo=1;
method.nonlinCF=savedBFlist;

minPitch=	80; maxPitch=	4000; numPitches=100;    % specify lags
pitches=10.^ linspace(log10(minPitch), log10(maxPitch),numPitches);
pitches=fliplr(pitches);
filteredSACFParams.lags=1./pitches;     % autocorrelation lags vector
filteredSACFParams.acfTau=	.003;       % time constant of running ACF
filteredSACFParams.lambda=	0.12;       % slower filter to smooth ACF
filteredSACFParams.lambda=	0.01;       % slower filter to smooth ACF

filteredSACFParams.plotACFs=0;          % special plot (see code)
filteredSACFParams.plotFilteredSACF=0;  % 0 plots unfiltered ACFs
filteredSACFParams.plotMoviePauses=.5;%.3          % special plot (see code)

filteredSACFParams.usePressnitzer=0; % attenuates ACF at  long lags
filteredSACFParams.lagsProcedure=  'useAllLags';
filteredSACFParams.criterionForOmittingLags=3;
filteredSACFParams.plotACFsInterval=50;%200;

% compute ACF
%switch saveAN_spikesOrProbability
%    case 'probability'
        inputToACF=ANprobRateOutput.^0.5;
%    otherwise
%        inputToACF=ANoutput;
%end

disp ('computing ACF...')
t=method.dt*(1:length(actualsignal.waveform));
[P, BFlist, sacf, boundaryValue] = ...
    filteredSACF(inputToACF, method, filteredSACFParams);
P = real(P); %dont know why sometimes P can be (very slightly) complex
disp(' ACF done.')

% SACF
set(gcf,'CurrentAxes',handles.axes8);
imagesc(P)
ylabel('periodicities (Hz)')
xlabel('time (s)')
%title(['running smoothed (root) SACF. ' saveAN_spikesOrProbability ' input'])
pt=[1 get(gca,'ytick')]; % force top xtick to show
set(gca,'ytick',pt)
set(gca,'ytickLabel', round(pitches(pt)))
tt=get(gca,'xtick');
tt=tt(tt<length(t));
set(gca,'xtickLabel', round(100*t(tt))/100)
colorbar;

% --- Executes on button press in pushbutton7.
function pushbutton7_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton7 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

%% SHOW THE  (PREVIOUSLY STORED) AUDITORY PROFILE CORRESPONDING TO THE
%% PARAMETER FILE
addpath(['..' filesep 'multithreshold 1.46' filesep]);
% read the parameter file
tmp = get(handles.edit7,'String');
parameterfile=tmp;
if isempty(parameterfile)
    error('Please specify parameter file');
end

%what kind of profiles are in the profiles directory?
files_in_profiles = dir(['..' filesep 'profiles' filesep '*.m']);%if exist('',file)
done = 0;

%try to map the specified parameter file to these found in the directory
%and
for iCounter = 1:length(files_in_profiles)
    if strfind(files_in_profiles(iCounter).name,parameterfile)
        %ask the user if this is the correct file - because you'll never
        %know ;)
        stringtoask = ['Is ' files_in_profiles(iCounter).name ' the correct profile file?'];
        Answer = questdlg(stringtoask,'Possible match found','Yes','No','Yes');
        switch Answer,
            case 'Yes'
                plotprofile(files_in_profiles(iCounter).name(1:end-2),'profile_CMA_L');
                done = 1;
            case 'No'
        end %switch
    end %if
end
if ~done
    msgbox('No match for this parameter file found.');
end

function multipleplot(actualsignal,handles)
%% FUNCTION TO DO 3 PLOTS
time_axis = [0:1/actualsignal.sfreq:(length(actualsignal.waveform)-1)/actualsignal.sfreq];

%waveform plot (temporal)
plot(handles.axes1,time_axis,actualsignal.waveform);
set(get(handles.axes1,'XLabel'),'String','Time (s)');
set(get(handles.axes1,'YLabel'),'String','Amplitude');
set(handles.axes1,'XLim',[0 length(actualsignal.waveform)/actualsignal.sfreq]);
highestabsampl = max(abs(actualsignal.waveform))*1.05;
set(handles.axes1,'YLim',[-highestabsampl highestabsampl]);

%average spectrum plot
complspectrum = fft(actualsignal.waveform)/(length(actualsignal.waveform));
frequency = [0:1/time_axis(end):1/(2*(time_axis(2)-time_axis(1)))];
plot(handles.axes2,frequency,20*log10(sqrt(2)*abs(complspectrum(1:round(length(complspectrum)/2)))));
set(get(handles.axes2,'XLabel'),'String','frequency (Hz)');
set(get(handles.axes2,'YLabel'),'String','fourier amplitude (dB)');
set(handles.axes2,'XLim',[100 10000]);
set(handles.axes2,'XScale','log');

%spectrogram plot (10 ms temporal resolution)
[s,f,t] = spectrogram(actualsignal.waveform,hann(round(0.01*actualsignal.sfreq)),[],[],actualsignal.sfreq); %10ms short term windows
set(gcf,'CurrentAxes',handles.axes3);
if max(f)>10000
    [tmp,idx] = min(abs(f-10000));
    imagesc(t,f(1:idx),20*log10(abs(s(1:idx,:))));
else
    imagesc(t,f,20*log10(abs(s)));
end
axis xy;
set(get(handles.axes3,'YLabel'),'String','frequency (Hz)');
set(get(handles.axes3,'XLabel'),'String','Time (s)');


function plotrateOutput(ratepattern,sfreq,axeshandle,colorrange)
%% calculate rate representation and plot, code from Nick
global savedBFlist
%calculate rate representation
ANsmooth = [];%Cannot pre-allocate a size as it is unknown until the enframing
hopSize = 10; %ms
winSize = 25; %ms
winSizeSamples = round(winSize*sfreq/1000);
hann = hanning(winSizeSamples);
hopSizeSamples = round(hopSize*sfreq/1000);
for chan = 1:size(ratepattern,1)
    f = enframe(ratepattern(chan,:), hann, hopSizeSamples);
    ANsmooth(chan,:) = mean(f,2)';
end

%plot rate representation
time_axis_rate=[0:hopSize/1000:size(ratepattern,2)/sfreq];
set(gcf,'CurrentAxes',axeshandle);
YTickIdx = 1:floor(numel(savedBFlist)/6):numel(savedBFlist);
XTickIdx = 1:floor(numel(time_axis_rate)/6):numel(time_axis_rate);
imagesc(ANsmooth);
axis xy;
set(gca, 'YTick', YTickIdx);
set(gca, 'YTickLabel', num2str(  savedBFlist(YTickIdx)', '%0.0f' ));
ylabel('best frequency (Hz)')
set(gca, 'XTick', XTickIdx);
set(gca, 'XTickLabel', XTickIdx.*10);
xlabel('Time (ms)');
caxis(colorrange); %set color from average spontaneous rate to a maximum of 600
colorbar;


% --- Executes if the radiobuttons hsr and lsr are changed
function radiobuttonselected(source, eventdata)

global rateaxeshandle
global ANprobRateOutput
global actualsignal

selected = get(get(source,'SelectedObject'),'String');

if isempty(ANprobRateOutput)
    msgbox('Please calculate AN pattern first!');
elseif strcmp(selected,'HSR')
    AN_HSRoutput = ANprobRateOutput(size(ANprobRateOutput)/2+1:end,:);
    colorrange = [33.7 300];
    plotrateOutput(AN_HSRoutput,actualsignal.sfreq,rateaxeshandle,colorrange);
else
    AN_LSRoutput = ANprobRateOutput(1:size(ANprobRateOutput)/2,:);
    colorrange = [10.7 100];
    plotrateOutput(AN_LSRoutput,actualsignal.sfreq,rateaxeshandle,colorrange);
end


% --- Executes during object creation, after setting all properties.
function edit11_CreateFcn(hObject, eventdata, handles)
% hObject    handle to edit11 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end


% --- Executes on button press in pushbutton8.
function pushbutton8_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton8 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
%plot the auditory nerve firing probability as line plot after grabbing the
%time value from edit11
global AN_HSRoutput
global actualsignal
global savedBFlist

middletime = str2num(get(handles.edit11,'String'))/1000;

start_time = min([size(AN_HSRoutput,2)/actualsignal.sfreq-0.01 max([0 middletime-0.025])]); %start time to plot in s
end_time = max([start_time+0.01 min([start_time+0.05 size(AN_HSRoutput,2)/actualsignal.sfreq-0.01])]); %plot max 50ms
plotIFRAN(AN_HSRoutput,start_time,end_time,actualsignal.sfreq,savedBFlist,handles.axes4);