unlock/analysis/andres/Code/getUnlockSpectrum.m File Reference

Functions
Calculating the spectrogram if	isfield (SpectInfo,'nEpochs')&&SpectInfo.doEpochs%For epochs
Calculating the spectrogram if if isfield(unlockInfo.epoch,'nOdds')%eegData(	eegData (:,:,:, 2)
	save (saveFile,'dataSpec','SpecInfo','-v7.3')
Variables
	function [spectData, SpectInfo]
Getting the size of	spectrogram [firstSpec,~] = mtspectrumc(squeeze(eegData(:,:,1,1)),SpectInfo.params)
Initial vble	eegDataSpec = nan(size(firstSpec,1),nChs,2)
Getting the spectrogram tic for	iTrial
end toc Super structure spectData	data = mtspectrumc(squeeze(eegData(:,:,iTrial,1)),SpectInfo.params) mtspectrumc(squeeze(eegData(:,:,iTrial,2)),SpectInfo.params) eegDataSpec
spectData	fspec = fspec
SpectInfo	specData = {'Freq.','chs','epoch','stimFreq'}
SpectInfo	size = size(spectData.data)
SpectInfo	FreqRes = spectData.fspec(3) - spectData.fspec(2)
SpectInfo	hz2spec = length(spectData.fspec)/(SpectInfo.freqResolut)
warning('CHECK IF freqREsolut IS BETTER THAN Fs/2 TO GET PROPER FREQ.limits when plotting!')%%%Saving matrix and info structure%if exist(fullfile(BCIinfo.dirs.saveFilename(1 SpecInfo	saveFile = saveFile

---

Function Documentation

Calculating the spectrogram if if isfield (unlockInfo.epoch,'nOdds') % eegData( eegData	(	:	,
		:	,
		:	,
		2
	)

Calculating the spectrogram if isfield	(	SpectInfo	,
		'nEpochs'
	)

save	(	saveFile	,
		'dataSpec'	,
		'SpecInfo'	,
		'-v7.3'
	)

---

Variable Documentation

else for all the data[spectData.data, spectData.fspec] = mtspectrumc(squeeze(eegData(:,:,iTrial,1)),SpectInfo.params) mtspectrumc(squeeze(eegData(:,:,iTrial,2)),SpectInfo.params) eegDataSpec

Initial vble eegDataSpec = nan(size(firstSpec,1),nChs,2)

SpectInfo FreqRes = spectData.fspec(3) - spectData.fspec(2)

end Spec analysis transformation factors SpectInfo fspec = fspec

function[spectData, SpectInfo]

Initial value:

 getUnlockSpectrum(eegData,SpectInfo)
%
% Calculates the spectrogram using chronux's multitaper functions 
%
% INPUT
% eegData:      matrix. This is the EEG raw data. Has to be in the form [samples x channels]
%
% OUTPUT
% 
%

nChs = SpectInfo.nChs

SpectInfo hz2spec = length(spectData.fspec)/(SpectInfo.freqResolut)

Getting the spectrogram tic for iTrial

Initial value:

 1:SpectInfo.nOdd       % number of epochs
        fprintf('Processing epoch %i\n',iTrial)

warning ('CHECK IF freqREsolut IS BETTER THAN Fs/2 TO GET PROPER FREQ. limits when plotting!') % %% Saving matrix and info structure % if exist(fullfile(BCIinfo.dirs.saveFilename(1 SpecInfo saveFile = saveFile

SpectInfo size = size(spectData.data)

SpectInfo specData = {'Freq.','chs','epoch','stimFreq'}

Getting the size of spectrogram[firstSpec,~] = mtspectrumc(squeeze(eegData(:,:,1,1)),SpectInfo.params)